Wow- that is a significant speedup!<br>
<br>
I don't mind the minor score differences either- does anyone else?<br>
<br>
Marty<br><br><div><span class="gmail_quote">On 12/16/05, <b class="gmail_sendername">Mike Muratet</b> <<a href="mailto:muratem@eng.uah.edu">muratem@eng.uah.edu</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Greetings<br><br>I ran a test on the the 2.2.13 BLAST that I thought the list might find<br>interesting. (Thanks to Joe Landman for putting it on the list, I thought<br>I was on that particular NCBI distribution but never saw the message until
<br>Joe put it here.)<br><br>I made a database out of the TIGR gene index for pig, SSGI.021105<br><br>I made a test query out of 200 70mers based on 5 sequences from<br>SSGI.021105.<br><br>I ran on a simple test on a Solaris/Sparc64 using
<br><br>time blastall [-V F] -p blastn -i test.fa -d SSGI.021105 > test.blast<br><br>yielding<br><br>          
2.2.12       2.2.13 no
-V      2.2.13 with -V F<br>user      1m
39.8s     1m
50.3s           0m
9.2s<br>real      1m
46.9s     1m
57.7            0m
9.5s<br>sys      
0m  6.3s     0m
6.3s            0m
0.3s<br><br>A diff on the output between v.12 and v.13(no V -F) produces only<br>differences in the header because the version numbers are different.<br><br>The diff between v.13(no V -F) and v.13 produces small differences in the
<br>scores and e-values. For example:<br><br>< NP1115201 GB|AY277629.1|AAQ24490.1 type XVII collagen [Sus scrofa]<br>139   4e-33<br>< AJ656543<br>34   0.18<br>< BP160342<br>32   0.73<br>< TC237977<br>32   0.73
<br>< TC212301<br>32   0.73<br>< CJ009554<br>30   2.9<br>< AJ657806<br>30   2.9<br>< BI344938<br>30   2.9<br>< AJ651764<br>30   2.9<br>< TC220684 UP|S3A3_HUMAN (Q12874) Splicing factor 3A subunit 3 (Sp...
<br>30   2.9<br>---<br>> NP1115201 GB|AY277629.1|AAQ24490.1 type XVII collagen [Sus scrofa]<br>138   6e-33<br>> AJ656543<br>34   0.19<br>> BP160342<br>32   0.75<br>> TC237977<br>32   0.75<br>> TC212301<br>32   
0.75<br>> CJ009554<br>30   3.0<br>> AJ657806<br>30   3.0<br>> BI344938<br>30   3.0<br>> AJ651764<br>30   3.0<br>> TC220684 UP|S3A3_HUMAN (Q12874) Splicing factor 3A subunit 3 (Sp...<br>30   3.0<br><br>So, it looks to be significantly faster.
<br><br>Cheers<br><br>Mike<br>_______________________________________________<br>Bioinformatics.Org general forum  -  <a href="mailto:BiO_Bulletin_Board@bioinformatics.org">BiO_Bulletin_Board@bioinformatics.org</a><br><a href="https://bioinformatics.org/mailman/listinfo/bio_bulletin_board">
https://bioinformatics.org/mailman/listinfo/bio_bulletin_board</a><br></blockquote></div><br><br clear="all"><br>-- <br>-- <br>Martin Gollery<br>Associate Director<br>Center For Bioinformatics<br>University of Nevada at Reno
<br>Dept. of Biochemistry / MS330<br>775-784-7042<br>-----------