Parsing Jackhmmer Output

Question

[Pawel Szczesny]

What are the options available to parse output (full output with alignments) from jackhmmer? As far as I can test, BioPerl hmmer module isn't able to read it correctly. Are there any ready-to-use libraries (language doesn't matter), or I should write something on my own?

[barkh22g]

Is there any update on this? Has a parsing program been updated or written to handle jackhmmer output?

Answer 1

[Pawel Szczesny]

Well, apparently the correct answer is that no library exists that parses full jackhmmer output. Its authors even write:

The output format is designed to be human-readable, but is often so voluminous that reading it is impractical, and parsing it is a pain. The --tblout and --domtblout options save output in simple tabular formats that are concise and easier to parse.

In case one needs sequences of hits, there are two options. One is to painfully write a parser. The other one is to extract sequence based on coordinates provided by tabular formats.