Hello, I am using macs2 for peak calling , on ChIPseq data using the following command. I have done the normalization using the drosophila genome on my data.
macs2 callpeak -t /Normalized_samples_ip/_IP_unique_rmdup_sampled.bam \
-c _001_unique_rmdup_sampled.bam -f BAM --gsize 3.0e9 --broad --broad-cutoff 0.1 --nomodel --extsize 125 --bdg -n peak_call --outdir outputfolder
but the result is generating 0 byte gapped and broad peaks file , is it ok ? or what am i doing wrong
any suggestion?? these are the files generated in my result folder
1.7G Oct 14 08:52 peak_call_control_lambda.bdg 0 Oct 14 08:52 peak_call_peaks.broadPeak 0 Oct 14 08:52 peak_call_peaks.gappedPeak 2.1K Oct 14 08:52 peak_call_peaks.xls 26M Oct 14 08:52 peak_call_treat_pileup.bdg 1.7G Oct 13 15:24 wildcard_control_lambda.bdg 0 Oct 13 15:24 wildcard_peaks.broadPeak 0 Oct 13 15:24 wildcard_peaks.gappedPeak 1.8K Oct 13 15:24 wildcard_peaks.xls 59M Oct 13 15:24 wildcard_treat_pileup.bdg
