HMMER是一个软件包,它提供了制作蛋白质和DNA序列域家族概率模型的工具,称为轮廓隐马尔可夫模型、轮廓HMM或仅轮廓,并使用这些轮廓来注释新序列、搜索序列数据库以寻找其他同源物,以及进行深度多重序列比对。HMMER是已知蛋白质和DNA序列域家族的几种综合比对和图谱的基础,包括Pfam数据库。
网页版 https://www.ebi.ac.uk/Tools/hmmer/文章来源:https://www.toymoban.com/news/detail-642737.html
conda install hmmer
# 各种安装方式。http://hmmer.org/documentation.html
# install_path/hmmer-3.3.2/tutorial test文件位置
### 1. hmmsearch:用profile文件搜索序列数据库
## Step 1: build a profile with hmmbuild
# globins4.sto 四个珠蛋白序列的MSA比对结果,Stockholm format
hmmbuild globins4.hmm globins4.sto
# 统计对序列比对(MSA) 的HMM(隐马尔科夫) profile
hmmstat globins4.hmm
## Step 2: search the sequence database with hmmsearch
# 下载数据库文件如 uniprot_sprot.fasta, globins45.fa为演示数据库
hmmsearch globins4.hmm globins45.fa > globins4.out
# 输出globins4.out为sequence top hits list,BLAST-like style
# E-value 期望的假阳性率
### 2 在序列数据库中迭代搜索
jackhmmer HBB_HUMAN uniprot_sprot.fasta
### 3. 用序列文件搜索profile数据库
## Step 1: create a profile database file
hmmbuild globins4.hmm globins4.sto
hmmbuild fn3.hmm fn3.sto
hmmbuild Pkinase.hmm Pkinase.sto
cat globins4.hmm fn3.hmm Pkinase.hmm > minifam
## Step 2: compress and index the flatfile with hmmpress
hmmpress minifam
## Step 3: search the profile database with hmmscan
hmmscan minifam 7LESS_DROME
参考:
http://eddylab.org/software/hmmer/Userguide.pdf文章来源地址https://www.toymoban.com/news/detail-642737.html
到了这里,关于HMMER-序列分析软件介绍的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!