#保留蛋白编码基因
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(第二个图)
刚刚测试过没问题的,你自己搞错了吧;
微信里点“发现”,扫一下
二维码便可将本文分享至朋友圈。
刚刚测试过没问题的,你自己搞错了吧;