PMC:1635699 / 12044-13846
Annnotations
{"target":"https://pubannotation.org/docs/sourcedb/PMC/sourceid/1635699","sourcedb":"PMC","sourceid":"1635699","source_url":"https://www.ncbi.nlm.nih.gov/pmc/1635699","text":"Each alignment program has to use a substitution matrix for replacement of characters during the alignment process. Traditionally these matrices differentiate between transitions (purine to purine and pyrimidine to pyrimidine substitutions) and transversions (purine to pyrimidine and vice versa), but more complex matrices have been described in the literature. An example for the latter are the RIBOSUM matrices [30] used by RSEARCH to score alignments of single-stranded regions. To address the question whether incorporating RIBOSUM matrices results in a significant performance change, we used the RIBOSUM 85–60 4 × 4 matrix as substitution matrix for CLUSTALW, ALIGN-M and POA, as these programs allow an easy integration of non-default substitution matrices via command line options. Since gap-costs and substitution matrix values are interdependent we adjusted the original RIBOSUM values to the range of the default values. We applied Wilcoxon tests to test whether using the RIBOSUM matrix (instead of the simpler default matrices) yields a statistical significant performance change. Results are summarized in Table 4. POA and ALIGN-M perform significantly better, only CLUSTALW's performance suffers from RIBOSUM utilization. The reason for CLUSTALW's performance loss is not obvious to us; it might be that CLUSTALW's dynamic variation of gap penalties in a position and residue specific manner [27] works optimally only with CLUSTALW's default matrix. Furthermore, the RIBOSUM 4 × 4 matrix is based on nucleotide substitutions in single-stranded regions whereas we used it as a general substitution matrix. Other matrices, based on base-paired as well as loop regions from a high-quality alignment of ribosomal RNA [31], gave, however, no significantly different results (data not shown).","tracks":[{"project":"2_test","denotations":[{"id":"17062125-14499004-1688900","span":{"begin":415,"end":417},"obj":"14499004"},{"id":"17062125-7984417-1688901","span":{"begin":1409,"end":1411},"obj":"7984417"},{"id":"17062125-15860776-1688902","span":{"begin":1730,"end":1732},"obj":"15860776"}],"attributes":[{"subj":"17062125-14499004-1688900","pred":"source","obj":"2_test"},{"subj":"17062125-7984417-1688901","pred":"source","obj":"2_test"},{"subj":"17062125-15860776-1688902","pred":"source","obj":"2_test"}]}],"config":{"attribute types":[{"pred":"source","value type":"selection","values":[{"id":"2_test","color":"#ecbd93","default":true}]}]}}