[BioPython] Problem with blastx output parsing =~

Christof Winter winter at biotec.tu-dresden.de
Mon Jun 4 17:08:34 UTC 2007


Hi Maia:

Could you post the begin and end of your blastx output? I think you can omit the 
query and hits in between...

Cheers,
Christof


Italo Maia wrote:
> Well, i have a perfectly fine blastx output that throws an error when parsed
> by biopython.
> It gives me this output:
> 
> Traceback (most recent call last):
>  File "<stdin>", line 1, in <module>
>  File "/var/lib/python-support/python2.5/Bio/Blast/NCBIStandalone.py", line
> 624, in parse
>    self._scanner.feed(handle, self._consumer)
>  File "/var/lib/python-support/python2.5/Bio/Blast/NCBIStandalone.py", line
> 99, in feed
>    self._scan_parameters(uhandle, consumer)
>  File "/var/lib/python-support/python2.5/Bio/Blast/NCBIStandalone.py", line
> 570, in _scan_parameters
>    has_re=re.compile(r"[Ll]ength of \s*[Dd]atabase"))
>  File "/var/lib/python-support/python2.5/Bio/ParserSupport.py", line 300,
> in read_and_call
>    raise SyntaxError, errmsg
> SyntaxError: Line does not match regex '[Ll]ength of \s*[Dd]atabase':
> Number of HSP's gapped: 136690
> 
> What could i do??? I'm using ubuntu feisty here.
> 
> 



More information about the Biopython mailing list