么有加,你核对一下数据的,在最初的GFF里面搜索一下有没有这个ID的基因;
#设置基因组文件名字变量
genome=Gura_genome.fa
gff=Glycyrrhiza_uralensis_Fisch.gene.gff
species=Glycyrrhiza_uralensis_Fisch
#保留蛋白编码基因
agat_sp_filter_feature_by_attribute_value.pl --gff $gff --attribute biotype --value protein_coding -t '!' -o $species.protein_coding.gff3
#保留最长转录本
agat_sp_keep_longest_isoform.pl --gff Glycyrrhiza_uralensis_Fisch.gene.gff -o $species.longest_isoform.gff3
跑完以上脚本后新的gff文件会把ID换成nbisL1-mrna-,这样好像和我的ID 对不上了