tests/testthat/test-createGeneAnnotation.R

window <- GRanges("chr16:11348649-11349648")
gene_anno <- createGeneAnnotation(
    window = window,
    TxDb = TxDb.Hsapiens.UCSC.hg19.knownGene::TxDb.Hsapiens.UCSC.hg19.knownGene
)

test_that("Output of createGeneAnnotation is correct", {
  expect_s4_class(gene_anno, "GRanges")
  expect_equal(unique(strand(gene_anno)), factor(c("+", "-"), levels=c("+", "-", "*")))
  expect_equal(unique(gene_anno$type), c("GENE", "EXON"))
})

Try the UMI4Cats package in your browser

Any scripts or data that you put into this service are public.

UMI4Cats documentation built on Dec. 31, 2020, 2:01 a.m.