This is a test version of Biostars. For the public version, visit https://www.biostars.org.
ChIPseq gapped and braod peaks file are empty

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

chipseq broadpeak gapped peak

0 answers

No answers yet.

Log in to answer this question.