[Biopython] Parsing text output of hmmsearch from HMMER 3.1b1
D. Bradley
db534 at cam.ac.uk
Fri Jun 12 15:52:46 UTC 2015
Hello,
I am trying to parse a text output file from HMMER vesion 3.1b1; the
file was produced using the following command:
* hmmsearch Pkinase.hmm H_sapiens.fasta.gz > PKhits.out
In BioPython, I have tried to begin parsing this file by issuing the
following commands:
* from Bio import SearchIO
* hmmer_qresult = SearchIO.read('PKhits.out', 'hmmer3-text')
However, an error is produced and I get the following output:
"qid = regx.group(1).strip()"
"AttributeError: 'NoneType' object has no attribute 'group'"
I do not know what the problem is. I am a beginner bioinformatician and
this is actually the first time I've tried to parse anything, so maybe
the error is simply in my understanding. However, from a quick web
search, somebody else seems to have had a similar problem
(http://permalink.gmane.org/gmane.comp.python.bio.general/8027).
Any help would be greatly appreciated; I can send my output file
straight away upon request.
Best regards,
David Bradley (PhD student)
More information about the Biopython
mailing list