Biostar Beta. Not for public use.
How to caclulate ChIP-Seq coverage?
1
Entering edit mode
2.4 years ago
Tori • 80
Finland

I have a ChIP-seq sample along with input control. I want to know the percentage of genome being enriched and coverage. I used bamFingerprint from deepTools and CHANCE. They gave conflicting results. How can I do it manuallly in a reliable way?

ChIP-Seq • 2.3k views
ADD COMMENTlink
0
Entering edit mode
13 months ago
Freiburg, Germany

Presumably you have a BED or similar file with the position of the significant peaks. You can sum the widths of the peaks and divide by the total genome size to get the fraction of the genome enriched. For the coverage, samtools depth is a convenient enough method. Just average the per-site depths. You can also supply the BED file of your peaks to get the average coverage inside the peaks, if that's helpful.

ADD COMMENTlink
0
Entering edit mode

I only have SAM file.

ADD REPLYlink
0
Entering edit mode

Then you need to use something to first call enriched regions. This is not something you want to do/write yourself unless you're into method development.

ADD REPLYlink
0
Entering edit mode

If I already have narrowPeak file, which contains the enriched region, how can I convert it to RPM or RPKM or tag density? Thanks.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1