[Biopython-dev] Notification: incoming/77

Jeffrey Chang jchang at smi.stanford.edu
Tue Aug 6 15:03:02 EDT 2002


Hi Courtney,

Are you using the version of NCBIWWW from the CVS, or from release
1.00a4?  I believe these issues have been solved in the CVS version.
Please let me know if that doesn't work.

Jeff




On Thu, Aug 01, 2002 at 05:13:06PM -0400, biopython-bugs at bioperl.org wrote:
> JitterBug notification
> 
> new message incoming/77
> 
> Message summary for PR#77
> 	From: zakano at hotmail.com
> 	Subject: Blast bug
> 	Date: Thu, 1 Aug 2002 17:13:06 -0400
> 	0 replies 	0 followups
> 
> ====> ORIGINAL MESSAGE FOLLOWS <====
> 
> >From zakano at hotmail.com Thu Aug  1 17:13:06 2002
> Received: from localhost (localhost [127.0.0.1])
> 	by pw600a.bioperl.org (8.12.2/8.12.2) with ESMTP id g71LD6JD009191
> 	for <biopython-bugs at pw600a.bioperl.org>; Thu, 1 Aug 2002 17:13:06 -0400
> Date: Thu, 1 Aug 2002 17:13:06 -0400
> Message-Id: <200208012113.g71LD6JD009191 at pw600a.bioperl.org>
> From: zakano at hotmail.com
> To: biopython-bugs at bioperl.org
> Subject: Blast bug
> 
> Full_Name: Courtney Harper
> Module: NCBIWWW.py
> Version: ?
> OS: windows
> Submission from: leopard.cgl.ucsf.edu (128.218.21.28)
> 
> 
> THERE ARE TWO PROBLEMS CAUSED BY A CHANGE IN BLAST OUTPUT SINCE THE RELEASE OF
> THE BIOPYTHON VERSION I DOWNLOADED LAST WEEK.
> 
> PROBLEM #1, IN:
>     def _scan_database_info(self, uhandle, consumer):
> 
> THE FOLLOWING IS EXPECTED:
> <b>Database:</b> All GenBank+EMBL+DDBJ+PDB sequences 
>            1,357,022 sequences; 6,227,440,795 total letters
> 
> BLAST OUTPUT HAS BEEN CHANGED TO THE FOLLOWING:
> <b>Database:</b> All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS,
> GSS, or phase 0, 1 or 2 HTGS sequences) 
>            1,357,022 sequences; 6,227,440,795 total letters
> 
> TO FIX THIS, I CHANGED THIS CODE:
>         read_and_call(uhandle, consumer.database_info, contains='Database')
>         read_and_call(uhandle, consumer.database_info, contains='sequences')
>         read_and_call(uhandle, consumer.noevent, blank=1)
> 
> TO THIS CODE:
>         read_and_call(uhandle, consumer.database_info, contains='Database')
>         read_and_call(uhandle, consumer.database_info, contains='sequences')
>         read_and_call(uhandle, consumer.database_info, contains='total')
>         read_and_call(uhandle, consumer.noevent, blank=1)
> 
> 
> PROBLEM #2, IN:
>     def _scan_database_report(self, uhandle, consumer):
> 
> THE FOLLOWING IS EXPECTED:
> <PRE>
>    Database: Non-redundant SwissProt sequences
>      Posted date:  Dec 18, 1999  8:26 PM
> 
> BLAST OUTPUT HAS BEEN CHANGED TO THE FOLLOWING:
> <PRE>
>   Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, GSS,
>   or phase 0, 1 or 2 HTGS sequences)
>     Posted date:  Aug 1, 2002 12:00 AM
> 
> TO FIX THIS, I CHANGED THIS CODE:
>         read_and_call(uhandle, consumer.noevent, start='<PRE>')
>         read_and_call(uhandle, consumer.noevent, start='  Database')
>         read_and_call(uhandle, consumer.posted_date, start='    Posted')
> 
> TO THIS CODE:
>         read_and_call(uhandle, consumer.noevent, start='<PRE>')
>         read_and_call(uhandle, consumer.noevent, start='  Database')
>         read_and_call(uhandle, consumer.database, start='  or')
>         read_and_call(uhandle, consumer.posted_date, start='    Posted')
> 
> 
> 
> 
> 
> _______________________________________________
> Biopython-dev mailing list
> Biopython-dev at biopython.org
> http://biopython.org/mailman/listinfo/biopython-dev



More information about the Biopython-dev mailing list