Skip to content

[Bioc-devel] Error on Windows when selecting a single chromosome from a BSgenome

1 message · Hervé Pagès

#
Hi Bernat,

Sorry for the delay in getting back to you. Like Dan, I cannot
reproduce this either. Without more input from you (e.g. sessionInfo()),
there is not much we can do.

Note that an easy way to get the coordinates of the masked regions in
a data.frame is with something like this:

   path <- system.file("extdata", "chr1.masks.rda",
                     package="BSgenome.Hsapiens.UCSC.hg19.masked")
   load(path)
   as.data.frame(chr1.masks[["AGAPS"]])

This is for the assembly gaps on chr1. You can do the same for any
of the 3 other types of masks (intra-contig ambiguities, RepeatMasker,
Tandem Repeats Finder) on any other chromosome/scaffold.

Cheers,
H.
On 03/20/2015 02:50 PM, Dan Tenenbaum wrote: