[BioPython] Problem with BLAST parser in NCBIStandalone
Jeffrey Chang
jchang@smi.stanford.edu
Wed, 8 May 2002 22:30:51 -0700
It looks like the format has changed a little bit. Try this patch and
let me know if this doesn't fix the problem.
Jeff
*** NCBIStandalone.py.orig 2002-03-18 13:46:08.056022000 -0800
--- NCBIStandalone.py 2002-05-08 22:25:44.664509000 -0700
***************
*** 199,205 ****
# indicates that no descriptions follow, and we should go straight
# to the alignments.
if not attempt_read_and_call(
! uhandle, consumer.description_header, contains='Score E'):
# Either case 2 or 3. Look for "No hits found".
attempt_read_and_call(uhandle, consumer.no_hits,
contains='No hits found')
--- 199,206 ----
# indicates that no descriptions follow, and we should go straight
# to the alignments.
if not attempt_read_and_call(
! uhandle, consumer.description_header,
! has_re=re.compile(r'Score {4,5}E')):
# Either case 2 or 3. Look for "No hits found".
attempt_read_and_call(uhandle, consumer.no_hits,
contains='No hits found')
On Thu, May 09, 2002 at 12:21:18PM +1000, Chris Francis wrote:
>
>
> Running a blast locally (BLASTP 2.2.3):
>
> blastall -p blastp -d nr.seq -i TyrR.fasta -o TyrR2.out
>
> The resulting file format was unable to be parsed by
> NCBIStandalone.BlastParser, giving the previously posted error.
>
> I have attached the output, as requested.
>
> Thank-you for looking into this.
>
> C.Francis
>
>
>
>
>
>
> _________________________________________________________________
> Send and receive Hotmail on your mobile device: http://mobile.msn.com