hmmscan using Pfam databases

Question

[dimitrischat]

Hello,

wanted to ask if you could do a specific search using pfam databases, f.e. plants ( taxid: 33090 )?

Thanks

hmmscan --cpu 8 --domtblout pfam.domtblout /path/to/Pfam-A.hmm transdecoder_dir/longest_orfs.pep

edit: or maybe i could edit the database Pfam-A.hmm and keep only taxonomy 33090 (viridiplantae) ?

[Fatima]

Could you elaborate on your question?

How To Ask A Good Question:

[[ Please read before posting a question ::: How To Ask A Good Question ]]

HMMER Tutorial

https://saf.bio.caltech.edu/saf_manuals/hmmer/v2_3_2.pdf

Answer 1

[h.mon]

Go to the Pfam search page, select the taxonomy tab, enter Viridiplantae as query term and toggle the Find domains unique to query term box. It should return Viridiplantae-specific domains, but I am unable to confirm, as I get the following message:

There was a system error on your last request.

This error has been recorded and will be investigated. If you would like to let us know what you were doing when the error occurred, please file a bug report.

You might like to go back to the previous page and see if it works this time...

I don't know if this is a temporary issue, a bug or some error from my part. I tried the above search strategy with Viridiplantae, Embryophyta, Tracheophyta, Spermatophyta, all produced the same error message. You may try contacting Pfam to further investigate the issue. If one does not select the Find domains unique to query term box, the search does indeed returns results - but this is not what you want.

edit: as Pfam is organized in terms of protein domains, it may well be there is no domain exclusive to the above search terms. If this is the case, the Pfam error message could be improved.

[dimitrischat]

Thanks a lot for your help. Unfortunately i get the same error when toggling the "..unique.." option. I contacted them via email, hope they reply.

[dimitrischat]

what do you mean by that "as Pfam is organized in terms of protein domains, it may well be there is no domain exclusive to the above search terms. If this is the case, the Pfam error message could be improved." ?