[Biopython] NCBIWWW.qblast Expect value / E-value

Jasper Koehorst jasperkoehorst at gmail.com
Wed May 20 13:52:22 UTC 2009


Yes we have, at the moment are we trying to identify a lot of small peptides
in several organisms for a project. What we would like to do is to ignore
the e-value for this moment. But a solution has yet to be found.
jasper koehorst


2009/5/20 Peter Cock <p.j.a.cock at googlemail.com>

> On Wed, May 20, 2009 at 1:51 PM, Jasper Koehorst
> <jasperkoehorst at gmail.com> wrote:
> > Im running against problems when i try to identify short peptide
> sequences
> > of ±10 - 20 peptides. When i run this at the NCBI Blast website i get
> > results, these all have
> > high E-Values that i dont really care at the moment. The problem is when
> i
> > do this in biopython as stated below, i will not get any results...
> >
> > I believe this is due to the fact that biopython will not show results
> with
> > a "high" E-Value. Is there a way to change this? So it will allow results
> > with an E-value of ±1500 or more?
>
> Nothing in Biopython limits the expectation values - our qblast function
> defaults to 10, but you can set this to what you like. However, the NCBI
> may be imposing their own limit. Are you sure using anything more than
> 10 is actually meaningful?
>
> > I tried the sript below but that does not quit work...
> >
> > Anybody has an idea?
> >
> >        result_handle =
> > NCBIWWW.qblast("blastp","nr",peptide[1],entrez_query=i, expect=2000,
> > matrix_name='BLOSUM80')
> >        blast_records = NCBIXML.parse(result_handle)
> >        blast_record = blast_records.next()
> >        E_VALUE_TRESH = 10000000
> >        for alignment in blast_record.alignments:
> >            for hsp in alignment.hsps:
> >
>
> I would guess (from previous examples) that this is due to the NCBI
> website and QBLAST API using different default parameters - the
> NCBI likes to change the defaults on the website from time to time,
> and these may differ from what you are getting via their QBLAST
> API. I would start by checking the gap parameters.
>
> See also:
> http://lists.open-bio.org/pipermail/biopython/2008-May/004252.html
> http://lists.open-bio.org/pipermail/biopython/2007-August/003679.html
>
> Peter
>




More information about the Biopython mailing list