difference between treat_pileup and bdgcmp fold enrichment tracks macs2
0
1
Entering edit mode
5 weeks ago
buffealo ▴ 40

Hello,

I created bigwig file from a treat_pileup.bdg file generated by macs2 and also used treat_pileup.bdg and control_lambda.bdg with macs2 bdgcmp.

Here is my codes;

peak calling

macs2 callpeak -t sample.bam -c sample_input.bam -g hs -f BAM -q 0.001 --bdg --outdir /folder -n sample

bedGraphToBigWig sample_treat_pileup.bdg hg19.len sample_treat_pileup.bw

fold enrichment tracks

macs2 bdgcmp -t sample_treat_pileup.bdg -c sample_control_lambda.bdg -o sample_foldEnrichment.bdg -m FE

sort -k 1,1 -k2,2n sample_foldEnrichment.bdg > sample_foldEnrichment_sorted.bdg

bedGraphToBigWig sample_foldEnrichment_sorted.bdg hg19.len sample_foldEnrichment_sorted.bw

When I see the tracks on IGV,

tracks

{UP: macs2 treat_pileup BOTTOM: Fold enrichment} the tracks look very different. Peak locations which are to focus, are different, right? Also highest peaks are changing between two signals.. I read that fold enrichment (using the control_lambda) dedcut the noise. So, Can I assume that the true signal can be seen in fold enrichment tracks?

Thank you.

chipseq chipsignals chipseqtracks chipseqanalysis IGVtracks • 96 views
ADD COMMENT

Login before adding your answer.

Traffic: 1712 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