Entering edit mode
10.7 years ago
Pappu
★
2.1k
I am wondering if it is possible to output full length sequences of hmmsearch matches. Otherwise I have to extract the sequence ids and download the full length sequences from uniprot. "-A" option only outputs the alignment of the matching parts of sequences.
Your question is not clear. Which sequences do you want? Are you interested in the sequences that you use to query the model file? You should already have those. Note that a HMMER search does not match any individual sequence. Rather, it matches a profile which is constructed from a multiple alignment of (typically) several sequences.
I want to output the full length sequences of the hits.
I think perhaps you should investigate exactly what a Hidden Markov Model is.