Identify fixed sites among sequences.
fixedSites(x, bases = c("a", "c", "g", "t", "-"), simplify = TRUE)
a gtypes object with sequences, a list of sequences, or a consensus sequence. Sequences must be aligned.
gtypes
a character vector of valid bases to consider.
if there is a single locus, return result in a simplified form? If FALSE a list will be returned wth one element per locus.
FALSE
a vector of fixed sites. Element names are site positions in the original sequence.
variableSites
# NOT RUN { data(dolph.haps) fixedSites(dolph.haps) # }
Run the code above in your browser using DataCamp Workspace