tests/testthat/test-bamread.R

In zslastman/RiboStan: Estimates Isoform Specific RPF Footprint Densities

test_that("reading bam file works",
{
  data(chr22_anno)
  testbam <- system.file("extdata", "nchr22.bam", package = "Ribostan", mustWork = TRUE)
  rpfs <- get_readgr(testbam, chr22_anno)
  expect_equal(length(rpfs), 63346)
  expect_true(all(seqnames(rpfs) %in% names(chr22_anno$exonsgrl)))
  expect_true(all(colnames(mcols(rpfs)) == "readlen"))
  expect_true(all(strand(rpfs) == "+"))
}
)