基因家族2.0数据准备,跑完代码后生成新的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 对不上了attachments-2023-12-523IEqC8657ed0119ae72.png

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

么有加,你核对一下数据的,在最初的GFF里面搜索一下有没有这个ID的基因;

请先 登录 后评论
  • 1 关注
  • 0 收藏,773 浏览
  • gz_li 提出于 2023-12-17 18:41

相似问题