Identify variable sites among sequences.
variableSites(x, bases = c("a", "c", "g", "t", "-"), simplify = TRUE)
a '>gtypes object with sequences, a DNAbin object, or a list of sequences.
DNAbin
character vector of 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 list with:
a DNAbin object composed of variable sites.
a matrix of base pair frequencies by site.
fixedSites
# NOT RUN { data(dolph.haps) variableSites(dolph.haps) # }
Run the code above in your browser using DataCamp Workspace