[BioC] how to load RangedData for all chromosomes in one tab in UCSC browser using browserView from rtracklayer
Yue Li
gorillayue at gmail.com
Thu May 31 19:00:25 CEST 2012
Dear List,
I wonder how to load RangedData for all chromosomes in a SINGLE tab in UCSC browser using browserView from rtracklayer.
Currently, if I have a GRanges object for "mm9" build:
> alignGR
GRanges with 238161 ranges and 0 elementMetadata cols:
seqnames ranges strand
<Rle> <IRanges> <Rle>
SRR039212.1000031 chr19 [ 8790316, 8790351] -
SRR039212.1000085 chr5 [106579844, 106579879] +
SRR039212.1000087 chr8 [109778747, 109778782] +
SRR039212.1000088 chr8 [ 93777537, 93777572] +
SRR039212.1000132 chr3 [128910749, 128910784] +
SRR039212.1000149 chr8 [127433402, 127433437] +
SRR039212.1000170 chr15 [ 93546853, 93546888] +
SRR039212.1000174 chr18 [ 32056273, 32056308] -
SRR039212.1000177 chr7 [ 90292453, 90292474] -
... ... ... ...
SRR039212.999792 chr2 [162907151, 162907186] -
SRR039212.999805 chr6 [ 44021338, 44021373] +
SRR039212.999810 chr4 [121106682, 121106717] -
SRR039212.999844 chr19 [ 60848841, 60848876] +
SRR039212.999848 chr5 [117644397, 117644432] -
SRR039212.999854 chr6 [132445007, 132445042] -
SRR039212.999855 chr7 [108392362, 108392397] -
SRR039212.999892 chr9 [ 20946884, 20946919] +
SRR039212.999901 chr2 [168845152, 168845187] +
---
seqlengths:
chr1 chr10 chr11 chr12 chr13 chr14 ... chr8 chr9 chrM chrX chrY
197195432 129993255 121843856 121257530 120284312 125194864 ... 131738871 124076172 16299 166650296 15902555
and I do this:
session <- browserSession()
track(session, "read alignments") <- RangedData(alignGR)
# launch browser view
browserView(session, alignGR)
This will open 22 tabs in my browser corresponding to 22 chromosomes in the GRanges object. I wonder if it would be possible to have just one single tab open for all 22 chromosomes to view on UCSC browser.
Thanks much,
Yue
More information about the Bioconductor
mailing list