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. DEDECMS ShowMsg()样式修改 提示信息的修改以及美化

    织梦DedeCMS系统,处处都在用到提示信息,但是这个提示框,前台后台一层不变,太死板了,可能有很多人都有过去修改它的想法,只是苦于不知道去哪里 改.今天我就来说说这个吧,DedeCMS的所有提示信息 ...

  2. 安卓秘钥生成命令以及SHA1值获取办法

    切换到秘钥所在目录,例如:cd C:\Program Files\Java\jdk1.8.0_171\bin 执行命令:keytool -genkey -alias demo.keystore -ke ...

  3. C/C++基础----类

    IO类属于不能被拷贝的类型,因此只能通过引用来传递.同时读取和写入操作都会改变流的内容,所以接收的是普通引用. 类内的友元声明仅仅指定了访问的权限,需要在友元声明之外再专门对函数进行一次声明. 可变数 ...

  4. 【Spring学习笔记-MVC-13】Spring MVC之文件上传

    作者:ssslinppp       1. 摘要 Spring MVC为文件上传提供了最直接的支持,这种支持是通过即插即用的MultipartResolve实现的.Spring使用Jakarta Co ...

  5. C++ 获取类成员函数地址方法 浅析

    C语言中可以用函数地址直接调用函数:   void print ()   {   printf ("function print");   }   typdef void (*fu ...

  6. tsql 循环id读取

    declare @IDList as varchar(max) declare @ID as int declare @i as int set @IDList='' )) + ',' from ta ...

  7. Dubbo与Zookeeper、Spring整合使用

    Dubbo与Zookeeper.Spring整合使用 Dubbo采用全spring配置方式,透明化接入应用,对应用没有任何API侵入,只需用Spring加载Dubbo的配置即可,Dubbo基于Spri ...

  8. php变量详细讲解

    变量是用于存储信息的"容器". 定义一个变量的语法: $变量名 = 值; 使用变量的例子: <?php $x=5; $y=6; $z=$x+$y; echo $z; ?> ...

  9. Linux Performance Analysis and Tools(Linux性能分析和工具)

    首先来看一张图: 上面这张神一样的图出自国外一个Lead Performance Engineer(Brendan Gregg)的一次分享,几乎涵盖了一个系统的方方面面,任何人,如果没有完善的计算系统 ...

  10. python的分支循环

    知识内容: 1.if-else分支结构 2.while循环 3.for循环 4.循环结构综述 5.break和continue语句 一.if-else分支结构 1.单分支选择结构 if 表达式: 语句 ...