Remove the intron peaks from the genome browser
1
0
Entering edit mode
8 months ago
bright602 ▴ 50

Hi, I finish my mapping and upload the bigwig files to IGV. However, I found lots of intron regions also have peaks, could someone tell me how to remove those intron peaks?

Thanks for your help! enter image description here

RNAseq • 464 views
ADD COMMENT
0
Entering edit mode

What kind of data is this? How did you generate your tracks? If it's RNA Seq data, and you used bedtools to generate coverage, you might try the -split option. There's a million ways to generate coverage tracks - guessing what you did, or why you want to hide introns without knowing what kind of data it is, is sort of useless. Please tell us what you have, what you did, how you did it, and why you did it.

ADD REPLY
0
Entering edit mode

Thanks for your reply. Yes, they are RNAseq data and I tried using the following two commands to get the bigwig files. But I don't know why there are intron peaks on those RNAseq tracks.

bedtools genomecov -ibam ${file}.unique.bam -bg | sort -k1,1 -k2,2n > ${file}.bg

${script}/bedGraphToBigWig ${file}.bg ${script}/${genome}.chrom.sizes ${file}.bw

ADD REPLY
3
Entering edit mode
8 months ago
seidel 9.6k

Given that this is RNA Seq data, when you create your coverage track using bedtools, supply the -split option, which will account for the space between exons and not add coverage to display. The figure shows the difference with and without split. Modifying your own code:

bedtools genomecov -ibam input.bam -bg -split | sort -k1,1 -k2,2n > output.bg

enter image description here

ADD COMMENT

Login before adding your answer.

Traffic: 692 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6