[Bioperl-l] RemoteBlast with same e-val (10 vs 1e-10) works
Madeleine Lemieux
mlemieux at bioinfo.ca
Fri Nov 12 00:40:27 EST 2004
Affan,
I ran your code with the expected value set to 10, as it is by default
on the NCBI Blast form, and got exactly the same results as when I
blasted your sequence on the NCBI server.
Madeleine
----- Output from your code calling RemoteBlast with e-val = 10
Query Name: affan
hit name is gb|AE016803.1|
score is 22
hit name is gb|AC124687.4|
score is 20
hit name is gb|AC112701.6|
score is 20
hit name is gb|AC016932.17|
score is 20
hit name is gb|AE008564.1|AE008564
score is 20
hit name is gb|AC073964.3|AC073964
score is 20
----- NCBI Blast output
Score
E
Sequences producing significant alignments: (bits)
Value
gi|27361305|gb|AE016803.1| Vibrio vulnificus CMCP6 chromoso... 44
0.28
gi|27356737|gb|AC124687.4| Mus musculus BAC clone RP24-447D... 40
4.4
gi|22296769|gb|AC112701.6| Mus musculus BAC clone RP23-171A... 40
4.4
gi|15789223|gb|AC016932.17| Homo sapiens 3 BAC RP11-166C10 ... 40
4.4
gi|15459696|gb|AE008564.1|AE008564 Streptococcus pneumoniae... 40
4.4
More information about the Bioperl-l
mailing list