[Bioc-devel] rtracklayer::import.bed(genome = Seqinfo)

Bhagwat, Aditya Ad|ty@@Bh@gw@t @end|ng |rom mp|-bn@mpg@de
Tue Oct 22 18:50:38 CEST 2019

Dear Michael,

Sorry for my incomplete email - the send button got hit too fast. Better this time.

rtracklayer::import.bed mentions the argument "genome" to be either a genome identifier (like 'mm10') or a Seqinfo object.

I notice that the second option does not work on my BED file (in attach).

# This works
rtracklayer::import.bed('SRF.bed', genome = 'mm10') # this works

# But this doesn't
seqinfo1<- GenomeInfoDb::seqinfo(BSgenome.Mmusculus.UCSC.mm10::BSgenome.Mmusculus.UCSC.mm10)
rtracklayer::import.bed('SRF.bed', genome = seqinfo1)

So I am requesting feedback.
I thought to use this channel


More information about the Bioc-devel mailing list