This function compares two haplotypes and returns a summary of the differences.
diffHaplo(h, a = 1, b = 2)
an object of class
two integers (or character strings) giving the indices (or labels) of the two haplotypes to be compared.
This function prints the number of transitions and transversions between both sequences, and returns a data frame with three columns giving the positions of the differences and the nucleotides in each sequence at these positions.
a data frame with three columns named
pos (position of the
differences) and the labels of the two haplotypes compared.
1 2 3 4 5