Hi everybody,
I've a little question. In R, I've two sequences ( DNA sequences ) and I want to align them and known the mismatch positions.
For example :
seq1 <- "ATGC"
seq2 <- "ATGG"
mismatch.pos <- aCoolFunction(seq1,seq2)
mismatch.pos
--> 4
Anyone knows such thing ?
if it doesn't exist in R, is it easily possible in bioperl ?
I've a little question. In R, I've two sequences ( DNA sequences ) and I want to align them and known the mismatch positions.
For example :
seq1 <- "ATGC"
seq2 <- "ATGG"
mismatch.pos <- aCoolFunction(seq1,seq2)
mismatch.pos
--> 4
Anyone knows such thing ?
if it doesn't exist in R, is it easily possible in bioperl ?
Comment