Check DNA Alignments
This function performs a series of diagnostics on a DNA alignement.
checkAlignment(x, check.gaps = TRUE, plot = TRUE, what = 1:4)
x |
an object of class |
check.gaps |
a logical value specifying whether to check the distribution of alignment gaps. |
plot |
a logical value specifying whether to do the plots. |
what |
an integer value giving the plot to be done. By default, four plots are done on the same figure. |
This function prints on the console a series of diagnostics on the set a aligned DNA sequences. If alignment gaps are present, their width distribution is analysed, as well as the width of contiguous base segments. The pattern of nucleotide diversity on each site is also analysed, and a relevant table is printed.
If plot = TRUE
, four plots are done: an image of the
alignement, the distribution of gap widths (if present), the Shannon
index of nucleotide diversity along the sequence, and the number of
observed bases along the sequence.
If the sequences contain many gaps, it might be better to set
check.gaps = FALSE
to skip the analysis of contiguous
segments.
NULL
Emmanuel Paradis
data(woodmouse) checkAlignment(woodmouse) layout(1)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.