Question: IGV tracks fail to render even after zooming
gravatar for mark.rose
10 months ago by
United States
mark.rose20 wrote:

Hi All

I have a tool that performs alignment against a reference sequence and then creates an IGV session file based on that alignment and other features of the reference sequence.

genome="dataset_18581.dat" locus="dataset_18581.dat" version="3"
<Resource name="726173-1-1 align" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//726173-1-1-1510602244-bowtie2.bam.sort.rmdup.bam.sort.bam"/>
<Resource name="726173-1-1 variants" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//726173-1-1-1510602244-bowtie2.bam.sort.rmdup.bam.sort.bam.vcf"/>
<Resource name="726173-1-1 zero coverage gaps" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//726173-1-1-1510602244-zero_coverage_gaps.bedgraph"/>
<Resource name="726173-1-1 below threshold coverage gaps" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//726173-1-1-1510602244-below_threshold_coverage_gaps.bedgraph"/>
<Resource name="726173-1-1 zero anchored coverage gaps" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//anchored-726173-1-1-1510602244-zero_coverage_gaps.bedgraph"/>
<Resource name="726173-1-1 below threshold anchored coverage gaps" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//anchored-726173-1-1-1510602244-below_threshold_coverage_gaps.bedgraph"/>
<Resource name="726173-1-1 soft clipping" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//COT102Frag1F4R1.txt-1510602244-soft_clipping.bedgraph"/>
<Resource name="726173-1-1 non-ref frequency" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//COT102Frag1F4R1.txt-1510602244-non-ref_freq.bedgraph"/>
<Resource name="726173-1-1 strand bias" path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//COT102Frag1F4R1.txt-1510602244-strand_bias.bedgraph"/>
<Resource name="726173-1-1 homopolymers " path="http://demeter//project/LEE2705/TEST/TK0000000/171013_M01610_0074_000000000-D30NY/726173-1-1/1510602244//COT102Frag1F4R1.txt-1510602244-homopolymers.bedgraph"/>

I have used this several times and it has worked fine. However, I have a dataset now where several of these tracks are not rendered though they are uploaded (because they are present in the track window) and they are read by IGV (since the displayed data range is correct for each track that contains data (not all do)). For this dataset the reference is small (1426bp) and the reads are 250bp PE. The alignment (bam) and coverage tracks seem to render properly but tracks showing strand bias and homopolymers (which contain data) are empty. I can zoom in as far as possible at locations where features should appear but there is nothing. Likewise no choice of "windowing function" makes them appear either.

Any suggestions


igv track • 462 views
ADD COMMENTlink modified 10 months ago by Jim Robinson220 • written 10 months ago by mark.rose20

Have you indexed the bedgraph files?

ADD REPLYlink written 10 months ago by Devon Ryan84k

I presume you mean with tabix. No, but I have never done that in the past and never, in the many times I've done this, had any issue. Is there a reason why indexing would be needed in some cases but not others? This is a smaller reference than I have used in the past.

ADD REPLYlink written 10 months ago by mark.rose20

I'm not sure if tabix is needed or if IGV uses it's own thing, but either way. Often IGV runs into memory issues if you don't index things like this. If that's the cause then it's an easy enough fix.

ADD REPLYlink written 10 months ago by Devon Ryan84k

Well I looked more online and found some people using tabix for indexing bed files for IGV and so that is what I did (after commenting the header)

sort -k1,1 -k2,2n test.strand_bias.bedgraph | bgzip > test.strand_bias.bedgraph.gz; tabix -s 1 -b 2 -e 3 test.strand_bias.bedgraph.gz

I then substituted the in path to this gzipped file in my session file and reloaded. The only difference is that the line which separates positive and negatives values (100 to -100) tuned blue but still no features.

ADD REPLYlink written 10 months ago by mark.rose20

Do you see any errors printed to the terminal? If you're not starting IGV from a terminal, do so so you can see these errors. At the end of the day that's likely the only way to diagnose this.

ADD REPLYlink written 10 months ago by Devon Ryan84k

I am launching it on a windows machine with a igv.bat file. When I do so a terminal window as well as the viewer appears. There is no error message in the terminal just the usual startup info ![][1]

Also I tried upgrading to version 2.4.4. No difference.

I then tried dropping the number of reads to align to my 1.4kb reference to 100 (250bp PE). Java has 2G of allotted memory to run IGV. I would think this would be way more than enough if there is some memory limitation. While I do now get some tracks with data showing up, they are ones reporting on aspects of coverage, which of course now is low with only 100 reads aligned. The homopolymer and strand bias track remain as before, blank.

ADD REPLYlink modified 10 months ago • written 10 months ago by mark.rose20
gravatar for Jim Robinson
10 months ago by
Jim Robinson220
United States
Jim Robinson220 wrote:

Hi, is your chromosome name really COT102Frag1F4R1.txt?

Its possible there are no records in the bedgraph file that match the sequence names in your reference fasta. That's just a guess though.

ADD COMMENTlink written 10 months ago by Jim Robinson220
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 876 users visited in the last hour