hmmsearch to show all sequences
1
1
Entering edit mode
5.6 years ago
srx3001 ▴ 10

Hello all,

I am using hmmsearch to go over a small database of peptides (basically just a multifasta file) to search with certain hmmprofile's. Now that is all working great, but since I am visualizing the data using plots - like an hmm score for each hit compared to function of the peptide I need something else.

I am curious if its possible to show ALL sequences from the multifasta file instead of just the ones it "hits". If I set the E-value to +100 (from the default inclusion of 0.01) it still only shows only a fraction of the sequences. It would be fine even if the score isn't informative , I would just like to have an overview since this was asked of me in the lab.

Basically my plot would have a couple of peptide functions plotted on Y-axis, as integers (cathegories) and HMMscore would be on X-axis. So far I was only able to get the hmmscore for 120/400 of my peptides when I run the profile for a "general sequence" of the family.

Appreciate all the help

hmmer hmmsearch multifasta • 1.9k views
ADD COMMENT
0
Entering edit mode

I ran into the same problem too. Does anyone have a solution to this?

ADD REPLY
0
Entering edit mode
2.9 years ago
jscience • 0

I think I figured it out. It can be done by turning off the nobias filter (with the --nobias option).

ADD COMMENT

Login before adding your answer.

Traffic: 1874 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