bigWig tracks inconsistent with MACS2 callpeak

0

Entering edit mode

3.8 years ago

lida ▴ 10

Hi, all

Recently, I have been dealing with some H3K27me3 ChIP-seq data. I used MACS2 to callpeak (with --broad ) and visualized the bw files (converted from bam) in IGV.

For me, this entire window seems possessed H3K27me3 modification, but the MACS2 just take the high signal region as H3K27me3 peak.

What I want to do is to get genes with H3K27me3 modification. If I take MACS2 as reference it may overlook some genes with low level H3K27me3 modification.

My question is how do yours handling this broad peak issues? Do we have better peakcallers for broad peaks like H3K27me3?

Thanks in advance

Lida

ChIP-Seq MACS2 • 1.3k views

ADD COMMENT • link 3.8 years ago by lida ▴ 10

0

Entering edit mode

Can you post your full MACS2 command, and are any of those tracks input or some sort of negative control?

ADD REPLY • link 3.8 years ago by rpolicastro 13k

0

Entering edit mode

Sure, my MACS2 command is macs2 callpeak -t <IP.bam> -c <Input.bam> -f BAM -g dm -p 0.01 -B --outdir macs_res --min-length 3000 --nomodel --extsize 300 --broad

This is input 1-5.scale.bw) and IP 11-15.scale.bw) samples

ADD REPLY • link 3.8 years ago by lida ▴ 10

Login before adding your answer.