No peaks detected in chipseq analysis-macs3
0
0
Entering edit mode
3 months ago
bioinfo_ga ▴ 70

Hi .. I am doing peak analysis using macs3 for ralstonia samples with the below command for paired end data:

macs3 callpeak -t t.bed -c c.bed --outdir ../Peak/c --nomodel -q 0.01

i tried inc the q value still do not get any peaks what could be the reason?

Thanks

chipseq • 325 views
ADD COMMENT
1
Entering edit mode

The most obvious reason is (please do not feel offended, it's not meant to be mocking you) that there are no peaks. Have you checked data on the IGV? Make a bigwig track, e.g. using bamCoverage from deeptools, and just "look" at the data.

ADD REPLY
0
Entering edit mode

Side issue: MACS3 can accept BAM as an input: MACS3 docs so there is no need to create often giant BED files from BAMs

ADD REPLY

Login before adding your answer.

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