ape (version 4.0)

seg.sites: Find Segregating Sites in DNA Sequences

Description

This function gives the indices of segregating (polymorphic) sites in a sample of DNA sequences.

Usage

seg.sites(x)

Arguments

x

a matrix or a list which contains the DNA sequences.

Value

A numeric (integer) vector giving the indices of the segregating sites.

Details

If the sequences are in a list, all the sequences must be of the same length. Ambiguous nucleotides are ignored.

See Also

base.freq, theta.s, nuc.div

Examples

data(woodmouse)
y <- seg.sites(woodmouse)
y
length(y)