R/singleBAMreads.R

.singleBAMreads <- function(bam, se, len, n) {
  w <- vector(mode = "numeric", length = 0)
  for (i in seq_len(len)) {
    se <- seq(i, n, len)
    w <- cbind(w, bam[se])
  }
  return(w)
}
skyhorsetomoon/Trumpet documentation built on May 30, 2019, 3:04 a.m.