[Bioperl-l] [Bioperl]problem with E-value
lidaof
lidaof at gmail.com
Thu Jan 18 05:05:48 UTC 2007
Hi,
the blast report i provided is the last 5000 line of my blastx result
i use "tail -5000 >" to generated that
because the origin result is so huge
and i will update my Bioperl using CVS to check if it works
Cheers
Li
On 1/18/07, Sendu Bala <bix at sendu.me.uk> wrote:
>
> 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'.
>
--
Li
More information about the Bioperl-l
mailing list