Entering edit mode
Lori nicely explained how a BSgenome
can be converted into a GenomicRanges::GRanges
object.
Would it be useful to turn this into an exported S4 method in GenomicRanges
or BSgenome
?
methods::setAs( "BSgenome",
"GRanges",
function(from) as(GenomeInfoDb::seqinfo(from), 'GRanges')
Thank you Michael .