[BioC] how to load RangedData for all chromosomes in one tab in UCSC browser using browserView from rtracklayer

Steve Lianoglou mailinglist.honeypot at gmail.com
Fri Jun 1 16:29:11 CEST 2012


Hi Yue,

On Thu, May 31, 2012 at 7:23 PM, Yue Li <gorillayue at gmail.com> wrote:
> Not sure if you were asking a question or implying that it's not feasible with UCSC browser. Anyhow, I could just upload to UCSC the bed/bam/bedGraph/gff file containing ranges for multiple chromosomes. Then go to the browser-view to visualize any chromosome in a single window with my browser?

I was just trying to get some clarification is all -- it sounded like
you wanted to see all of the different regions at once when you said:

"""I wonder how to load RangedData for all chromosomes in a SINGLE tab
in UCSC browser using browserView from rtracklayer."""

-steve

>
>
> On 2012-05-31, at 2:12 PM, Steve Lianoglou wrote:
>
>> Hi,
>>
>> On Thu, May 31, 2012 at 1:00 PM, Yue Li <gorillayue at gmail.com> wrote:
>>> 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.
>>
>> Is there anyway that you know of to have the genome browser do what
>> you are asking "directly"?
>>
>> I mean, if you were just navigating w/ the browser alone (not using
>> rtracklayer), how would you get it to do what you are asking?
>>
>> --
>> Steve Lianoglou
>> Graduate Student: Computational Systems Biology
>>  | Memorial Sloan-Kettering Cancer Center
>>  | Weill Medical College of Cornell University
>> Contact Info: http://cbio.mskcc.org/~lianos/contact
>



-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact



More information about the Bioconductor mailing list