Commit 11260ff8 by Riccardo Vicedomini

fixed a bug in parsing hmmsearch output when a '== domain' block does not contain a hit/alignment

parent d851d812
...@@ -41,7 +41,7 @@ def parseHMMER(fh,hmmDict): ...@@ -41,7 +41,7 @@ def parseHMMER(fh,hmmDict):
# process domain table # process domain table
line = fh.readline() line = fh.readline()
hitTable = [] hitTable = []
while line and not line.lstrip().startswith("== domain "): while line and not line.lstrip().startswith("== domain ") and not line.startswith(">> ") and not line.startswith("Query: "):
tokens = line.split() tokens = line.split()
## domain table row: ## domain table row:
## domNum [!?] score bias c-Evalue i-Evalue hmmfrom hmm-to .. alifrom ali-to .. envfrom env-to .. acc ## domNum [!?] score bias c-Evalue i-Evalue hmmfrom hmm-to .. alifrom ali-to .. envfrom env-to .. acc
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment