| seqEqual {alakazam} | R Documentation |
seqEqual checks if two DNA sequences are identical.
seqEqual(seq1, seq2, ignore = as.character(c("N", "-", ".", "?")))
seq1 |
character string containing a DNA sequence. |
seq2 |
character string containing a DNA sequence. |
ignore |
vector of characters to ignore when testing for equality. Default is to ignore c("N",".","-","?") |
Returns TRUE if sequences are equal and FALSE if they are not.
Sequences of unequal length will always return FALSE regardless of
their character values.
Used by pairwiseEqual within collapseDuplicates. See seqDist for calculation Hamming distances between sequences.
# Ignore gaps
seqEqual("ATG-C", "AT--C")
seqEqual("ATGGC", "ATGGN")
seqEqual("AT--T", "ATGGC")
# Ignore only Ns
seqEqual("ATG-C", "AT--C", ignore="N")
seqEqual("ATGGC", "ATGGN", ignore="N")
seqEqual("AT--T", "ATGGC", ignore="N")