Statistical evaluation and comparison of a pairwise alignment algorithm that a priori assigns the number of gaps rather than employing gap penalties.
BACKGROUND
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. 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 with the Smith-Waterman and Needleman-Wunsch techniques on a set of structurally aligned protein sequences. We demonstrate that this approach outperforms the other techniques, especially for short sequences (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.


BACKGROUND
The implementation of the described methods SANK_AL is available at http://cbbc.murdoch.edu.au/


BACKGROUND
matthew@cbbc.murdoch.edu.au.
Nozaki
Yasuyuki
Y
Centre for Bioinformatics and Biological Computing, Murdoch University, Murdoch, WA 6150, Australia.
Bellgard
Matthew
M
D000465:Algorithms; D017124:Conserved Sequence; D003627:Data Interpretation, Statistical; D008956:Models, Chemical; D008957:Models, Genetic; D015233:Models, Statistical; D011506:Proteins / Q000737:chemistry*; D016415:Sequence Alignment / Q000379:methods*; D020539:Sequence Analysis, Protein / Q000379:methods*; D017386:Sequence Homology, Amino Acid; D012984:Software; D012986:Software Validation
D011506:Proteins