Optimization of ribosomal RNA profile alignments.
BACKGROUND
Large alignments of ribosomal RNA sequences are maintained at various sites. New sequences are added to these alignments using a combination of manual and automatic methods. We examine the use of profile alignment methods for rRNA alignment and try to optimize the choice of parameters and sequence weights.


RESULTS
Using a large alignment of eukaryotic SSU rRNA sequences as a test case, we empirically compared the performance of various sequence weighting schemes over a range of gap penalties. We developed a new weighting scheme which gives most weight to the sequences in the profile that are most similar to the new sequence. We show that it gives the most accurate alignments when combined with a more traditional sequence weighting scheme.


BACKGROUND
The source code of all software is freely available by anonymous ftp from chah.ucc.ie in the directory /home/ftp/pub/emmet,in the compressed file PRNAA.tar:


BACKGROUND
emmet@chah.ucc.ie, des@chah.ucc.ie
O'Brien
E A
EA
Department of Biochemistry, University College, Cork, Ireland and EMBL-European Bioinformatics Institute, Hinxton, Cambridge CB10 1RQ, UK. emmet@chah.ucc.ie
Notredame
C
C
Higgins
D G
DG
D000818:Animals; D006801:Humans; D017433:Protein Structure, Secondary; D012335:RNA, Ribosomal / Q000737:chemistry*; D016415:Sequence Alignment; D012984:Software
D012335:RNA, Ribosomal