哥,您好,如何检查GFF文件里面的mRNA ID是否对应一致,或者如何改脚本,get_gene_weizhi.pl

attachments-2019-11-TqlsuRFD5dc4cfefe3da0.jpgattachments-2019-11-XQW6O0Bf5dc4cff4d3f9e.jpg

请先 登录 后评论

1 个回答

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

脚本认为ID=与;之间的为gene或mRNA的ID,建议你把GFF文件里面的gene:  CDS:  transcript: 删除; 然后再看看fasta里面的ID与GFF里面的ID是否一致。 这样保持所有地方的ID一致,就不会因为找不到对应的ID而出错。


如果上面按照我说的都改了,脚本这里27行就可以删除gene:    获取ID信息了:

attachments-2019-11-VVHgFlw95dc4d15549883.png

请先 登录 后评论