构建泛基因列表,提取最长转录本这一步,运行几个小时之后没有结果是怎么回事?

老师,我在构建泛基因集列表时,做提取最长转录本这一步,for i in `cat ID.txt`;do    $script/agat_sp_keep_longest_isoform.pl --gff $datadir/${i}.gff -o gff/${i}.longest_isoform.gff3 ; done     由于我的基因组较多,113个,所以先做了这一步。运行了好几个小时,结果里有所有基因组的agat.log文件,但查看bed,cds,gff目录下只有4个基因组的名称,内容为空,是怎么回事呢?
attachments-2025-02-rovhXBIi67b48a786b148.pngattachments-2025-02-KHSu3tCA67b48a92b4d0a.png

attachments-2025-02-zcDh55Pr67b48ad5241d2.pngattachments-2025-02-dmjWQ97i67b48afdd1ee7.pngattachments-2025-02-ILtsdG6467b48b1c80a35.png

请先 登录 后评论
  • 0 关注
  • 0 收藏,268 浏览
  • Connie 提出于 2025-02-18 21:33

相似问题