alakazam (version 0.2.2)

testSeqEqual: Test DNA sequences for equality.

Description

testSeqEqual checks if two DNA sequences are identical.

Usage

testSeqEqual(seq1, seq2, ignore = c("N", "-", ".", "?"))

Arguments

seq1

character string containing a DNA sequence.

seq2

character string containing a DNA sequence.

ignore

vector of characters to ignore when testing for equality.

Value

• 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.

See Also

Used by collapseDuplicates.

Examples

# Ignore gaps
testSeqEqual("ATG-C", "AT--C")
testSeqEqual("ATGGC", "ATGGN")
testSeqEqual("AT--T", "ATGGC")

# Ignore only Ns
testSeqEqual("ATG-C", "AT--C", ignore="N")
testSeqEqual("ATGGC", "ATGGN", ignore="N")
testSeqEqual("AT--T", "ATGGC", ignore="N")