Find Segregating Sites in DNA Sequences
This function gives the indices of segregating (polymorphic) sites in a sample of DNA sequences.
seg.sites(x, strict = FALSE, trailingGapsAsN = TRUE)
| x | a matrix or a list which contains the DNA sequences. | 
| strict | a logical value; if  | 
| trailingGapsAsN | a logical value; if  | 
If the sequences are in a list, they must all be of the same length.
If strict = FALSE (the default), the following rule is used to
determine if a site is polymorphic or not in the presence of ambiguous
bases: ‘A’ and ‘R’ are not interpreted as different, ‘A’ and ‘Y’ are
interpreted as different, and ‘N’ and any other base (ambiguous or
not) are interpreted as not different. If strict = TRUE, all
letters are considered different.
Alignment gaps are considered different from all letters except for
the leading and trailing gaps if trailingGapsAsN = TRUE (which
is the default).
A numeric (integer) vector giving the indices of the segregating sites.
Emmanuel Paradis
base.freq, theta.s, nuc.div (last two in pegas)
data(woodmouse) y <- seg.sites(woodmouse) y length(y)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.