Question: Reads detected/"spill over" on introns after using bamCoverage function from deepTools
gravatar for kdcastillo
14 days ago by
kdcastillo0 wrote:

Hello! We recently used the recent version of bamCoverage from the deepTools suite for our ribo-seq and RNA-seq bam files to generate bigwig files. The conversions ran smoothly, however when we loaded the bigwig files into IGV, we noticed some reads spilling over to some exon-intron junctions when we zoomed into each track. We believe these are not genome misannotations because the non-normalized/non-bigwig bam files don't show these junction reads. Any advice is appreciated. Thank you!

Best, Kat

ADD COMMENTlink modified 2 days ago by Biostar ♦♦ 20 • written 14 days ago by kdcastillo0

Hello, it would help if you add code and screenshots to understand the problem. My first guess when it comes to local anomalies is the chosen bin size and extension parameter (-bs -e).

ADD REPLYlink written 14 days ago by ATpoint44k

enter image description hereHello! I agree that it might be due to bin size, I tried the default and bin size of 10, and I am still getting the same issue. But I attached the script (for bin size 10) and the example image of the problem I'm seeing.

bamCoverage -b filename.bam -o --normalizeUsing RPKM -bs 10

Thanks for the response!


ADD REPLYlink modified 12 days ago by i.sudbery10k • written 12 days ago by kdcastillo0
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: 1677 users visited in the last hour