How To Reconstruct A Consensus From A Hmmer Hmm Profile
1
1
Entering edit mode
9.3 years ago

Hi,

I am trying to (and understand how to) extract the consensus of a protein alignment from its HMM profile (produced by HMMer2). I understand that from the "match emission line" I can get the corresponding amino acid (the one with the higher match emission score) when the most probable transition is 'm->m'. Am I right in assuming that a "dot" has to be inserted in the consensus when the most probable transition is 'm->d' or 'd->d'? what happens with 'd->', 'm->i', and 'i->i'? Any help would be appreciated.

Cheers,

M;

hmm hmmer alignment consensus • 2.9k views
ADD COMMENT
4
Entering edit mode
9.3 years ago
robert.d.finn ▴ 130

Try using the program hmmemit from the HMMER suite. I think that you need the -c option to get the consensus... The -h will list all options.

And yes a . Represents an insert respect to the model and a dash a delete state.

Best

Rob

ADD COMMENT
0
Entering edit mode

hmmemit makes it. I can get all the information I need to understand the process from the source code. Thanks!

ADD REPLY

Login before adding your answer.

Traffic: 1131 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6