Bioinformatics Advance Access published online on December 10, 2004
Bioinformatics, doi:10.1093/bioinformatics/bti198
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Centre for Bioinformatics and Biological Computing, School of Information Technology, Murdoch University, Murdoch, WA, Australia, 6150
* To whom correspondence should be addressed.
Motivation: Although pairwise sequence alignment is essential in comparative genomic sequence analysis, it has proven difficult to precisely determine the gap penalties for a given pair of sequences. A common practice is to employ default penalty values. However, there are a number of problems associated with using gap penalties. First, alignment results can vary depending on the gap penalties, making it difficult to explore appropriate parameters. Second, the statistical significance of an alignment score is typically based on a theoretical model of non-gapped alignments, which may be misleading. Third, and finally, there is no way to control the number of gaps for a given pair of sequences, even if the number of gaps is known in advance. Results: In this paper, we develop and evaluate the performance of an alignment technique that allows the researcher to assign a priori set of the number of allowable gaps, rather than using gap penalties. We compare this approach to the Smith and Waterman and Needleman and Wunsch techniques on a set of structurally aligned protein sequences. We demonstrate that this approach outperforms the other techniques, especially for short sequence (56-133 residues) with low similarity (<25%). Further, by employing a statistical measure, we show that it can be used to assess the quality of the alignment in relation to the true alignment with the associated optimal number of gaps. Availability: The implementation of the described methods SANK_AL is available at http://cbbc.murdoch.edu.au/.
Received July 11, 2004
Revised November 17, 2004
Accepted November 29, 2004
Article
Statistical evaluation and comparison of a pairwise alignment algorithm that a priori assigns the number of gaps rather than employing gap penalties
Matthew Bellgard, E-mail: m.bellgard{at}cbbc.murdoch.edu.au
![]()
Abstract ![]()
CiteULike
Connotea
Del.icio.us What's this?