[Bioperl-l] [BULK] can't parse blast file anymore

Christopher Fields cjfields at uiuc.edu
Mon May 8 02:04:13 UTC 2006


These are debugging lines (not errors); you still have the -verbose flag set.  

Did you follow Jason's advice?  I believe he's right on the money about the issue 
at hand...

Chris

---- Original message ----
>Date: Sun, 07 May 2006 19:41:14 -0600
>From: Hubert Prielinger <hubert.prielinger at gmx.at>  
>Subject: Re: [Bioperl-l] [BULK]  can't parse blast file anymore  
>To: Torsten Seemann <torsten.seemann at infotech.monash.edu.au>, bioperl-
l at bioperl.org, Chris Fields <cjfields at uiuc.edu>, Jason Stajich 
<jason.stajich at duke.edu>
>
>hi,
>I have corrected that and now I finally I got a few error messages:
>
>blast.pm: unrecognized line Reference: Altschul, Stephen F., Thomas L. 
>Madden, Alejandro A. Schäffer,
>blast.pm: unrecognized line Jinghui Zhang, Zheng Zhang, Webb Miller, and 
>David J. Lipman
>blast.pm: unrecognized line (1997), "Gapped BLAST and PSI-BLAST: a new 
>generation of
>blast.pm: unrecognized line protein database search programs", Nucleic 
>Acids Res. 25:3389-3402.
>blast.pm: unrecognized line RID: 
1137529800-24476-151611170370.BLASTQ1
>
>after that line it stops without terminating....
>
>
>Torsten Seemann wrote:
>> Hubert Prielinger wrote:
>>> ok, thanks
>>> I have submitted the bug
>>> bug #1994
>>
>> This is a line from the script you sent to Bugzilla:
>>
>> my $search = new Bio::SearchIO (
>> -verbose => 1,-format => 'blast', -file => $file)
>> or die "could not open blast report" if not defined my $search;
>>
>> Althoygh syntactically correct, I don't think it is doing what you want.
>> Please change it to this:
>>
>> my $search = new Bio::SearchIO(-format => 'blast', -file => $file) or 
>> die "could not open blast report";
>>
>> or alternatively, this:
>>
>> my $search = new Bio::SearchIO(-format => 'blast', -file => $file);
>> if (not defined $search) {
>>   die "could not open blast report";
>> }
>>
>> and let us know what happens.
>>
>> all the example output you have supplied still suggests that 
>> Bio::SearchIO can not load or parse your blast report.
>>
>




More information about the Bioperl-l mailing list