基因组注释信息GFF文件里面这些信息都有的,自己提取出来就行, 如果需要批量提取建议写perl或者python脚本: 可学习: perl入门到精通、perl语言高级、
回答于 2019-07-24 17:56
目前没有这样的课程,可参考文献:https://www.omicsclass.com/article/331
回答于 2019-07-24 09:57
SMART用法错了吧,参考这里:https://www.omicsclass.com/article/681
回答于 2019-07-23 17:37
可能有的转录本没有基因ID,所以这个脚本报错了,你可以在gff里面搜索对应的mRna行把他删除再运行脚本。
回答于 2019-07-23 10:04
如果全长序列差异太大建议用结构域序列构建进化树,检查很短的6条序列是否具有完整结构域,酌情删除。
回答于 2019-07-23 10:02
高版本的hmmsearch 会报这个错误,建议删除这个参数--cut_tc hmmer 3.1b2 不会报错,
回答于 2019-07-22 10:45