Question: Align protein sequence to Pfam file using HMMER3.0 (hmmalign)
0
gravatar for Xiaokang ZH
10 months ago by
Xiaokang ZH40
Xiaokang ZH40 wrote:

I'm using HMMER 3.0 (hmmalign to be specific) to align the whole proteome of Human to Pfam file. I got the proteome from https://www.ensembl.org/info/data/ftp/index.html, got the Pfam file from: ftp://ftp.ebi.ac.uk/pub/databases/Pfam/releases/Pfam31.0/, and downloaded the file Pfam-A.hmm.gz

With the command

hmmalign -o human.out Pfam-A.hmm Homo_sapiens.GRCh38.pep.all.fa

I got the error:

Error: HMM file Pfam-A.hmm does not contain just one HMM
pfam protein hmmer • 595 views
ADD COMMENTlink modified 8 months ago by Biostar ♦♦ 20 • written 10 months ago by Xiaokang ZH40

You might want to try hmmsearch or hmmscan (for the latter, you need to hmmpress your Pfam-A first).

ADD REPLYlink written 10 months ago by cschu1811.4k

Thank you. But why hmmscan, instead of hmmsearch? Since I found an article comparing them here. It seems that they have the same output, but hmmsearch is faster?

ADD REPLYlink written 10 months ago by Xiaokang ZH40

I have always asked this myself. I know, the output is not identical, at least not when you use the tblout/domtblout options. hmmscan is domain/hmm-oriented, whereas hmmsearch is sequence oriented. In the end, the speed issue as mentioned in your link (nice read by the way!) would convince me to use hmmsearch and then reformat the output.

ADD REPLYlink written 10 months ago by cschu1811.4k

Thank you very much!

ADD REPLYlink written 10 months ago by Xiaokang ZH40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1490 users visited in the last hour