how to call weak peak using macs2
0
0
Entering edit mode
18 months ago

We are interesting in ChIP-seq small peak , we have avidence that the the small or weak peaks are functional sites , how to call weak peak using macs2? I try to use the parameter q 0.05, but It can't call the weaker peaks.

MACS2 ChIP-Seq • 424 views
ADD COMMENT
0
Entering edit mode

There is no gold standard answer for this. You will need a good input control, probably chromatin input is the best sequenced to similar depths as the IP samples. You can use dedicated replicate-aware peak callers such as PePr or use other filtering strategies that increase confidence in your calls such as IDR or custom strategies like "peak has to be in at least 2/3 samples at FDR below a certain threshold". Still, confidence in small peaks requires excellent data quality. If you have a poor antibody then you will have a hard time calling weak binding events. Are you after a transcription factor? How is data quality?

ADD REPLY

Login before adding your answer.

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