10 保留蛋白编码gff

#保留蛋白编码基因

  agat_sp_filter_feature_by_attribute_value.pl --gff  $gff --attribute biotype --value protein_coding -t '!' -o $species.protein_coding.gff3

请问老师为什么运行这个命令的时候保留的gff文件里是biotype=rRNA/tRNA这些非编码基因(第一个图),而discard的却是protein_coding(第二个图)

attachments-2023-10-uExR5Fcx65309daf4ad58.pngattachments-2023-10-AxnhEUGO65309dcd5fb08.png

请先 登录 后评论
  • 0 关注
  • 0 收藏,1359 浏览
  • MenG 提出于 2023-10-19 11:09

相似问题