这篇文章给大家分享的是有关hmmsearch如何寻找相似序列的内容。小编觉得挺实用的,因此分享给大家做个参考,一起跟随小编过来看看吧。
hmmsearch寻找相似序列
HMMER是基于隐马尔可夫模型,用于生物序列分析工作的一个非常强大的软件包,它的一般用途是识别同源蛋白或核苷酸序列和进行序列比对。与BLAST、FASTA等序列比对和数据库搜索工具相比,HMMER更准确。
再有HMM模型的情况下寻找相似序列可用hmmsearch,用法如下:
hmmsearch --domtblout Ensembl.txt PF01535.hmm Zea_mays.AGPv3.24.pep.all.fa
输入文件PF01535.hmm为HMM模型。 比较序列文件可以是FASTA格式。
--domtblout:domtblout格式输出
结果按照E-values值从小到大排序,形式与blast类似。
其中target name是每个目标序列的名称;
query name是查询序列的名称;
score是比对得分,分值越高说明越相似;
E-value目标序列的期望值(统计意义);
最重要的是E-value值,值越小,越可信,相当于一个统计量。
感谢各位的阅读!关于“hmmsearch如何寻找相似序列”这篇文章就分享到这里了,希望以上内容可以对大家有一定的帮助,让大家可以学到更多知识,如果觉得文章不错,可以把它分享出去让更多的人看到吧!
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。