MergeIdenticalHaplotypes: Merge Alleles with Identical DNA Sequences
Description
If any alleles within a locus have identical alleleNucleotides values
(including those identical based on IUPAC ambiguity codes),
this function merges those alleles, summing their read depths. This function is
primarily intended to be used internally in cases where tags vary in length
within a locus, resulting in truncated alleleNucleotides.
Usage
MergeIdenticalHaplotypes(object, ...)
Value
A RADdata object identical to object, but with alleles merged.
data(exampleRAD)
# change a haplotype for this exampleexampleRAD$alleleNucleotides[5] <- "GY"
nAlleles(exampleRAD)
exampleRAD <- MergeIdenticalHaplotypes(exampleRAD)
nAlleles(exampleRAD)