下图蓝色块是我想研究的结构域,按老师讲的脚本应该$a[9]==3,可是为什么提取不到这部分的结构域序列呢?请问哪里是不是错了呢?我该怎样调整呢?之前做的这个文件geneid2mRNAid.txt也是空的。

attachments-2019-11-Ike0GcBS5dc7d2f55b007.pngattachments-2019-11-mTaR3Kgg5dc7d32165e05.pngattachments-2019-11-8axwusoe5dc7d34c129c8.png

请先 登录 后评论

1 个回答

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

GFF文件里面的ID看看,和拟南芥蛋白质序列里面的ID是否一样,

怀疑你的ID里面的gene:  等没有删除;

请先 登录 后评论