[Bioperl-l] [Bioperl]problem with E-value

Sendu Bala bix at sendu.me.uk
Wed Jan 17 16:54:27 UTC 2007


Chris Fields wrote:
> I pretty sure I know exactly what the problem is and how to fix it (if 
> you haven't done it already).  Looks like the parser trashes the rest of 
> the BLAST results data since it's not catching the next BLAST report 
> header (and not breaking out of the while() loop).  That may be what is 
> triggering the e-value error.  I have been parsing new BLAST reports 
> recently w/o that one popping up, but it may be a difference between the 
> web BLAST report and the executable (wouldn't be the first time that has 
> happened.
> 
> Did you want me to take a look?

Sure. In case you didn't notice, the reason it isn't catching the next 
BLAST report is the lack of alignments. For the few results that do have 
alignments, that's where it 'works'.



More information about the Bioperl-l mailing list