R/sq_comp.R

In g-r-eg/fold2: Hairpin bisulphite sequence analyser.

#' Get the complementary DNA sequence.
#'
#' @param sq sequence
#' @export
sq_comp <- function(sq) {
  ret <- character(length(sq))
  ret[sq == "A"] <- "T"
  ret[sq == "T"] <- "A"
  ret[sq == "G"] <- "C"
  ret[sq == "C"] <- "G"
  ret[sq == "N"] <- "N"
  return(ret)
}