Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Jane M
    Senior Member
    • Aug 2011
    • 239

    #1

    Visualize unionBedGraphs file in UCSC

    Dear all,

    I have generated bedgraphs that I can visualize in UCSC. For example:
    Code:
    head BIO_31-T1-ARN.bedGraph
    track type=bedGraph name=BIO_31-T1-ARN
    chr5	170000000	170000600	0.00
    chr5	170000600	170000700	2.00
    chr5	170000700	170000850	0.00
    chr5	170000850	170000950	2.00
    chr5	170000950	170001000	0.00
    chr5	170001000	170001100	2.00
    chr5	170001100	170001250	0.00
    chr5	170001250	170001400	2.00
    chr5	170001400	170002200	0.00
    I want to compare bam of samples from 2 conditions by visualizing them on UCSC. In one condition, I have 5 samples and in the other condition, there are 16 samples. I was told that it is possible to visualize only one "summarizing sample" per condition.
    I ran unionBedGraphs from bedtools on the 5 samples of condition 1:
    unionBedGraphs -i BIO-12-T1-ARN/tophat_out/BIO-12-T1-ARN.bedGraph BIO_26-T1-ARN/tophat_out/BIO_26-T1-ARN.bedGraph BIO_28-T1-ARN/tophat_out/BIO_28-T1-ARN.bedGraph BIO_052_T1_ARN/tophat_out/BIO_052_T1_ARN.bedGraph BIO_057_T1_ARN/tophat_out/BIO_057_T1_ARN.bedGraph -header > Cond1.bedGraph
    I got a merged file of 1.4Go:
    Code:
    head Cond1.bedGraph
    chrom	start	end
    chr11	60000000	60008500	0	0	0	0	0.00
    chr11	60008500	60008600	0	0	0	0	1.00
    chr11	60008600	60009850	0	0	0	0	0.00
    chr11	60009850	60010050	0	0	0	0	2.00
    chr11	60010050	60010400	0	0	0	0	0.00
    chr11	60010400	60010450	0	0	0	0	1.00
    chr11	60010450	60016500	0	0	0	0	0.00
    chr11	60016500	60016600	0	0	0	0	1.00
    chr11	60016600	60047950	0	0	0	0	0.00
    When I try to upload this file in UCSC, I get the following error:
    Error File 'Cond1.bedGraph' - Unrecognized format line 1 of file: chrom start end
    I changed the header line with:
    track type=bedGraph name=Cond1
    but I got the following error:
    needMem: trying to allocate 1414098818 bytes (limit: 500000000)
    Then, I took only the 500,000 first lines of Cond1.bedGraph with "track type=bedGraph name=Cond1" as Header, but again, I got an error:
    Error File 'Cond1_500000.bedGraph' - Error line 2 of custom track: Expecting + or - in strand
    Finally I tried another header:
    Code:
    head Cond1_500000_v3.bedGraph
    chrom	start	end	S1	S2	S3	S4	S5
    chr11	60000000	60008500	0	0	0	0	0.00
    chr11	60008500	60008600	0	0	0	0	1.00
    chr11	60008600	60009850	0	0	0	0	0.00
    chr11	60009850	60010050	0	0	0	0	2.00
    chr11	60010050	60010400	0	0	0	0	0.00
    chr11	60010400	60010450	0	0	0	0	1.00
    chr11	60010450	60016500	0	0	0	0	0.00
    chr11	60016500	60016600	0	0	0	0	1.00
    chr11	60016600	60047950	0	0	0	0	0.00
    ending with:
    Error File Cond1_500000_v3.bedGraph' - Unrecognized format line 1 of file: chrom start end S1 S2 S3 S4 S5 (note: chrom names are case sensitive, e.g.: correct: 'chr1', incorrect: 'Chr1', incorrect: '1')
    Now, I wonder if it is possible to visualize this merged file in UCSC.
    Can you please tell me if it is possible to visualize it somehow? Should I modify something in the header?

    Thank you for your help,
    Jane
    Last edited by Jane M; 11-25-2016, 01:45 AM.

Latest Articles

Collapse

  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM
  • SEQadmin2
    New Genomics Technologies Take Aim at Long-Standing Limits
    by SEQadmin2


    Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

    We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
    ...
    09-28-2026, 10:25 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
29 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-15-2026, 10:52 AM
0 responses
39 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-14-2026, 05:39 AM
0 responses
45 views
0 reactions
Last Post SEQadmin2  
Working...