[Bioperl-l] What does Expect(2) mean in a blast result?
sac at bioperl.org
Sat Nov 10 17:59:28 EST 2007
The Bioperl blast parser should extract that value and you can obtain
it from an HSP object, via the HSPI::n() method, documented here:
Dave's basically correct in his explanation. It's a result of the
application of sum statistics by the blast algorithm. You can read all
about it in Korf et al's BLAST book. Here's the relevant section:
On Nov 10, 2007 3:04 AM, Dave Messina <David.Messina at sbc.su.se> wrote:
> Hi Amir,
> I don't have my BLAST book handy, and my memory is a little fuzzy, but I
> think the Expect(2) you're seeing is the E-value based on both HSPs
> combined. And I think this is why you see the same Expect value for both --
> because it is shared between them (which sounds like what you wanted).
> Again, this is just from memory, but I think this is an option that has to
> be turned on rather than something which Blast decides to do on its own.
> I don't know whether BioPerl reports this or not. Would you mind e-mailing
> me a entire BLAST report as a sample? When I have some time I'd like to play
> around with this a bit.
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
More information about the Bioperl-l