[Bioperl-l] frac_* methods

aaron.j.mackey at gsk.com aaron.j.mackey at gsk.com
Wed Oct 3 01:40:14 UTC 2007


Let me second Jason's comment that while BLAST is a great search program, 
it is not a very good alignment algorithm.  In this day and age with so 
many good pairwise alignment algorithms out there (customized for the 
context in which the alignment is performed), BLAST-based alignments 
should frankly be ignored.  See: exonerate, pairagon, etc.

Oh, and since ssearch35 (the Smith-Waterman algorithm that comes with the 
FASTA package) is now vector-parallelized on most i386 architectures, it 
is only about 10 times slower than BLAST for complete database searches 
(with superior sensitivity/specificity); add PVM or MPI-based CPU 
parallelization on top of that, and there's almost no reason to even run 
BLAST anymore ...

-Aaron

bioperl-l-bounces at lists.open-bio.org wrote on 10/02/2007 06:22:59 PM:

> I think my answer before was something to the tune of:
> 
> Use an alignment algorithm that finds a single best alignment like 
> FASTA or Smith-Waterman (SW) if what you want is a single number that 
> represents the alignment.  BLAST is great for fast searching but 
> FASTA or SW/SSEARCH are going to be better at creating an alignment. 
> Consider the -postsw option in WUBLAST as well as it will realign the 
> HSPs with SW.
> 
> I personally never use the frac alignment summary stats for the Hit 
> object for this reason unless I know I am going to have a single HSP.
> 
> -jason
> 
> On Oct 2, 2007, at 2:41 PM, Thiago Venancio wrote:
> 
> > Hi all,
> >
> > This topic was discussed before, but I would like to put it on the 
> > list
> > again, maybe someone has an update.
> >
> > The methods frac_identical, frac_conserved, frac_aligned_query and
> > frac_aligned_hit can also be used in the hit context, after HSP 
> > tilling. In
> > my point of view, it is better to use it just in HSPs individually, 
> > because
> > there are some rare/strange kinds of alignments. However, we 
> > frequently need
> > to get one measure of the whole alignment.
> >
> > Any of the BioPerl masters has an update on this topic ? What is 
> > the best
> > current usage ?
> >
> > Best.
> >
> > Thiago
> >
> > -- 
> > "Innovation distinguishes between a leader and a follower."
> >     Steve Jobs
> >
> > ========================
> > Thiago Motta Venancio, MSc
> > PhD student in Bioinformatics
> > University of Sao Paulo
> > ========================
> > _______________________________________________
> > Bioperl-l mailing list
> > Bioperl-l at lists.open-bio.org
> > http://lists.open-bio.org/mailman/listinfo/bioperl-l
> 
> --
> Jason Stajich
> jason at bioperl.org
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
> 





More information about the Bioperl-l mailing list