This sequence alignment falls on the empirical threshold which distinguishes instances of significant similarity from uncertain alignments [41].
To account for slight positional variations in ncRNA predictions, and minor inaccuracies, gaps and small indels in the genomic alignments, we only infer orthology among similar ncRNA sequences across the genome alignment where the alignment falls within a range of ±80 nucleotides for snoRNAs and ±40 nucleotides for (the generally shorter) miRNAs across the entire alignment.
One approach could be to implement model partitioning by inferring among-site variation from the alignment, for example, using a mixture model in the context of a Bayesian framework, such as that developed by Lartillot and Philippe (2004) in PhyloBayes, where each site in the alignment falls into one of several classes characterized by its own set of frequencies (CAT model).
Despite the high level of detection, more than 80% of the alignments fell outside the coding sequences.
For one of the eight BES pairs, 20 hits were detected with the V. vinifera genome and all of these intertwined alignments fell into the 'no-gapped' category for chromosome 2 of V. vinifera.
A majority of the genes in the 4-way alignments fell into one of the three categories, defined based on the microarray experiment: of the 990 genes, 99 (10%), 94 (9%) and 628 (63%) fell into the sexual, vegetative and constitutive gene categories, respectively.
With this additional condition, the probability of scoring pairs that contain spurious alignments falls to less than 0.01%.
Only the Tetrahymena genome produced a match, which upon various CLUSTALW alignments falls outside of the Group B superclade from which the Myc and Max clades group.
Ten of these residue positions show statistically significant evidence for positive selection, and two of them (sites with coordinates 49 and 50 in the protein alignment) fall within the known Tp1 epitope (Table 2, Fig. 1b).
To ensure that the endpoints of the optimal alignment fall on the N- and C-terminal of some blocks, H i, j may be re-initialized to S(a i, b j ) (replacing initialization to 0) if a i or b j is the first column in its respective block and traceback through H i, j is required to terminate at that position.
