Nothing
# Create example granges
bsgenome <- BSgenome.Mmusculus.UCSC.mm10::BSgenome.Mmusculus.UCSC.mm10
granges <- GenomicRanges::GRanges('chr1',
'100-200',
strand = '-',
seqinfo = seqinfo(bsgenome))
chrlength <- GenomeInfoDb::seqlengths(bsgenome)[['chr1']]
# Test
context('extend')
test_that('extend works', {
expect_equal(
start(extend(granges, -5, 5)),
start(granges) - 5)
expect_equal(
end(extend(granges, -5, 5)),
end(granges) + 5)
})
test_that('extend warns for coordinates < 1', {
expect_warning(extend(granges, -500))
})
test_that('extend warns for coordinates > chrlength', {
expect_warning(extend(granges, 1, 1 + chrlength, plot = FALSE))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.