Question: macs2: how to retrieve normalized counts
gravatar for olechnwin
2.5 years ago by
olechnwin0 wrote:

Anyone knows how to get the normalized counts for the ChIP-seq peaks using MACS2 ?

I would like to get the count of reads in the peak's region after normalization in both input and treatment (ChIP pull down).

MACS2 generates the treat_pileup.bdg and the control_lambda.bdg. What exactly is this? Can I use these files to extract the information to get the pileup (normalized read counts) in the entire peak region and not just at the summit which was outputted in the peak.xls file.

Thank you in advance for your help.

chip-seq macs2 • 1.4k views
ADD COMMENTlink modified 2.5 years ago by ATpoint14k • written 2.5 years ago by olechnwin0
gravatar for ATpoint
2.5 years ago by
ATpoint14k wrote:

You could take the bed file (narrowPeak or broadPeak) that contains the peak regions and then use BEDtools intersect or BEDtools coverage together with your sam/bam/bed file that contains the sequencing data to get counts per peak. Then simply normalize these counts by, e.g. reads per million or any appropriate scaling factor.

ADD COMMENTlink written 2.5 years ago by ATpoint14k
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: 899 users visited in the last hour