Question: how to visualize rna-seq coverage track in igv?
gravatar for Zee_S
13 months ago by
Zee_S60 wrote:

Hello biostars community,

how does one visualize an rna-seq coverage track alongside chip-seq bigwig tracks in igv?

the problem i have is that, when i zoom into the rna-seq coverage track (becaus eit is a bam file), i lose the resolution on my bigwig tracks (too zommed in) and i can no longer visualize everything in parallel.

the other question is that, coverage track is not normalized, how do i normalize it so that two rna-seq coverage tracks can be compared to one another for assessing the difference in the level of expression?

many thanks!

ADD COMMENTlink modified 13 months ago by jared.andrews078.0k • written 13 months ago by Zee_S60

You can display bam files in larger zoom out than the default. Somewhere in the settings.

ADD REPLYlink written 13 months ago by Asaf8.5k

what do you mean by larger zoom out? can you please elaborate?

ADD REPLYlink written 13 months ago by Zee_S60

You can redefine the point it writes "Zoom in to see alignment", make it wider (will slow down obviously), can't find where at the moment. As for you second question you can right click the coverage track and unmark "Autoscale" and click on "Set Data Range" to set the desired range.

ADD REPLYlink written 13 months ago by Asaf8.5k
gravatar for jared.andrews07
13 months ago by
Memphis, TN
jared.andrews078.0k wrote:

As for your second question, the short answer is that you can't.

Using coverage tracks (even if normalized via RPKM/FPKM/what have you) to try to determine differential gene expression has all sorts of issues. They should only be used as a visual aid to show expression changes between samples that were identified with robust, statistically-backed methods using gene counts like edgeR, DESeq2, limma, etc. If you just want them to serve as that sort of visual aid, then deepTools is an excellent option for creating depth/control-normalized tracks.

ADD COMMENTlink written 13 months ago by jared.andrews078.0k

this is very useful feedback thank you very much.

ADD REPLYlink written 13 months ago by Zee_S60

Hello! I would like to come back to this post to ask: when we visualise RNA-seq tracks, do you include both unique and multi mappers? or why is it that some people only use unique mappers? which is the correct way?

Thank you!

ADD REPLYlink modified 8 months ago • written 8 months ago by Zee_S60
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: 1768 users visited in the last hour