Hi All,
I have HMM profiles but the raw sequences or multiple sequence alignments are not available. I saw the tool hmmemit, but this randomly generates sequences that match a given profile HMM. I do not want random sequences. Is there a way to obtain sequences or multiple sequence alignment that generated the HMM profile? Additionally, I know skylign generates hmm logos, is there a way to obtain the consensus sequences from the logo? The other option would be to programmatically parse the HMM profile to try to generate the consensus sequence.
Any input would be appreciated!
Thanks in advance.
•
written
by
frcamacho
2 answers
I do not think it is possible to get those sequences as HMM algorithm does not save them, however, you can use hmmemit as you mentioned. It is possible to tweak the hmmemit parameters to match the profile. Please refer to the HMMER documentation.
Unless the original sequences were saved somewhere, you can't get them back from the HMM alone. hmmemit -c will give you a consensus sequence.
Log in to answer this question.
-
Scan multiple sequences on multiple hmm profiles
written by biobiu
-
minimum number of protein sequences for a sequence logo
written by Stefano •
-
Generate HMM profile without an MSA (hmmer)
written by elenatea.russo •
-
Multiple Sequence Alignment
written by Natasha •
-
Running HMMER via Python
written by Ghoti •
-
Produce Multiple Sequence Alignments using Profile HMMs
written by bordieszter •
-
Strategy for RBP sequence logo alignments to other RNA sequences
written by syntonic.c •
-
Is there a way to perform a multiple sequence alignment in HMMER (hmmalign) using multiple profiles?
written by marinaga •
-
Creating Sequence Logos From The Command Line
written by Qdjm
-
Method For Randomly Selecting Subsequences In An Alignment ?
written by Spyros •
hmmemitis your only option as far as I'm aware. It doesn't randomly emit a sequence, it produces the most likely sequence given the hidden states of the Markov Model profile (a kind of consensus sequence).