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
Aditya
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel
