1)使用HMM模型搜索序列数据库(以青蟹蛋白库为例，简写为qingxie.pep),同源参考序列(query.fas)

hmmbuild： 用多重比对序列构建HMM模型；
hmmsearch： 使用HMM模型搜索序列库；

步骤1：
pfam下载多重比对文件的种子序列(PF02898_seed_NOS.txt)
如果没有就用多条参考序列进行多重比对，然后用hmmbuild构建模型（NOS基因为例）
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/hmmbuild PF02898_seed_NOS.hmm PF02898_seed_NOS.txt
步骤2:用hmmsearch去搜索青蟹蛋白库
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/hmmsearch PF02898_seed_NOS.hmm qingxie.pep >NOS.out

2) phmmer和jackhmmer的用法
phmmer:与Blastp类似，使用蛋白质序列(一条或者多条都行),对数据库进行搜索
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/phmmer query.fas qingxie.pep >qingxie.phmmer.out
jackhmmer：与psiBlast类似，蛋白质序列迭代搜索蛋白质序列库
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/jackhmmer query.fas qingxie.pep >qingxie.jackhmmer.out

3) nhmmer（类似于blastn，输入的query可以是fasta格式的DNA或者RNA序列,但是只能是一条。）
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/nhmmer Ssa_sdy.ref.cds(核酸参考序列) Omy_omyt1(核酸库) >nhmmer.out

nhmmer如果是多条请先做多重序列比对（hmmalign），然后使用建立模型(hmmbuild)，然后生成的hmm作为输入）
nhmmer MADE1.hmm Ssa_sdy.ref.cds(核酸参考序列)  Omy_omyt1(核酸库) > nhmmer.out

3)使用蛋白质序列搜索HMM数据库
Pfam-A.hmm数据库路径：/PUBLIC/software/DENOVO/bio/database/Pfam/Pfam-A.hmm

步骤1 ：首先要对HMM数据库进行格式化（包括压缩以及创建索引），很快会完成。
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/hmmpress Pfam-A.hmm  
步骤2：用hmmscan将蛋白序列去搜索HMM数据库
/PUBLIC/software/DENOVO/bio/software/hmmer-3.1b2-linux-intel-x86_64/binaries/hmmscan  Pfam-A.hmm query.fas >hmmcan.out
 
 接下来可以注释出结构域，并提提出结构域序列，建树。

hmm 软件的使用的更多相关文章

  1. 转：从头开始编写基于隐含马尔可夫模型HMM的中文分词器

    http://blog.csdn.net/guixunlong/article/details/8925990 从头开始编写基于隐含马尔可夫模型HMM的中文分词器之一 - 资源篇 首先感谢52nlp的 ...

  2. 在weka中添加libSVM或者HMM等新算法

    转:http://kasy-13.blog.163.com/blog/static/8214691420143226365887/ Weka的全名是怀卡托智能分析环境(Waikato Environm ...

  3. 共有11款Python 中文分词库开源软件

    件过滤: 排序: 收录时间 | 浏览数 Python 中文分词库 Yaha "哑哈"中文分词,更快或更准确,由你来定义.通过简单定制,让分词模块更适用于你的需求. "Ya ...

  4. DBImport V3.7版本发布及软件稳定性（自动退出问题）解决过程分享

    DBImport V3.7介绍: 1:先上图,再介绍亮点功能: 主要的升级功能为: 1:增加(Truncate Table)清表再插入功能: 清掉再插,可以保证两个库的数据一致,自己很喜欢这个功能. ...

  5. nw.js桌面软件开发系列 第0.1节 HTML5和桌面软件开发的碰撞

    第0.1节 HTML5和桌面软件开发的碰撞 当我们谈论桌面软件开发技术的时候,你会想到什么?如果不对技术本身进行更为深入的探讨,在我的世界里,有这么多技术概念可以被罗列出来(请原谅我本质上是一个Win ...

  6. 标准产品+定制开发:专注打造企业OA、智慧政务云平台——山东森普软件，交付率最高的技术型软件公司

    一.公司简介山东森普信息技术有限公司(以下简称森普软件)是一家专门致力于移动互联网产品.企业管理软件定制开发的技术型企业.公司总部设在全国五大软件园之一的济南齐鲁软件园.森普SimPro是由Simpl ...

  7. 结巴分词3--基于汉字成词能力的HMM模型识别未登录词

    作者:zhbzz2007 出处:http://www.cnblogs.com/zhbzz2007 欢迎转载,也请保留这段声明.谢谢! 1 算法简介 在 结巴分词2--基于前缀词典及动态规划实现分词 博 ...

  8. JAVA程序员常用软件整理下载

    ********为了大家学习方便,特意整理软件下载如下:*************Java类软件:-------------------------------JDK7.0:http://pan.ba ...

  9. spring boot 实战：我们的第一款开源软件

    在信息爆炸时代,如何避免持续性信息过剩,使自己变得专注而不是被纷繁的信息所累?每天会看到各种各样的新闻,各种新潮的技术层出不穷,如何筛选出自己所关心的? 各位看官会想,我们是来看开源软件的,你给我扯什 ...

随机推荐

  1. dede模块管理一片空白或没有列表内容的解决办法

    为什么dede后台模块管理,打开之后一片空白,又或者没有列表,插件与其他模块的使用也是正常的. 这主要是因为我们在安装模块,然后又卸载模块,卸载的时候选择了删除安装文件,就会出这个问题. 这里面分为两 ...

  2. springMVC学习(11)-json数据交互和RESTful支持

    一.json数据交互: json数据格式在接口调用中.html页面中较常用,json格式比较简单,解析还比较方便. 比如:webservice接口,传输json数据. springMVC进行json交 ...

  3. Frequently Asked Questions

    转自:http://www.tornadoweb.org/en/stable/faq.html Frequently Asked Questions Why isn’t this example wi ...

  4. Log4net详细说明（全）

    转自:http://www.cnblogs.com/zhangchenliang/p/4546352.htmlhttp://www.cnblogs.com/zhangchenliang/p/45463 ...

  5. 常用Java程序片段

    1.改变数组的大小 package com.js.ai.modules.jsa.test; public class Testxf { private static Object resizeArra ...

  6. 小朋友学Python（2）

    本节学习python的第一个程序:输出“Hello World!”.咱们用两种方法实现 方法(一) 进入python环境,直接使用print方法   hello world.png 方法(二) 先编程 ...

  7. HTML|CSS之布局相关总结

    知识内容: 1.浮动相关 2.display属性 3.居中显示 4.盒模型和box-sizing 5.position属性 6.响应式设计 7.flex布局 8.其他 参考:http://zh.lea ...

  8. Table-Driven Design 表驱动设计

    注:本文所有代码来自 http://www.codeproject.com/Articles/42732/Table-driven-Approach 在许多程序中,经常需要处理那些拥有种种色色不同特性 ...

  9. 5.mybatis实战教程(mybatis in action)之五:与spring3集成(附源码)

    转自:https://blog.csdn.net/nnn9223643/article/details/41962097 在 这一系列文章中,前面讲到纯粹用mybatis 连接数据库, 然后 进行增删 ...

  10. Delphi 解析HTML

    uses mshtml; IHTMLEleMent.ID; IHTMLEleMent.tagName; IHTMLEleMent.title;elmt._className;elmt.getAttri ...