CmhaDSO MACS2

使用例

ピークコール
macs2 callpeak -c ctrl.bam -t target.bam -g mm -p 1e-5 -f BAM -B -n prefix --outdir outDir

  # -c: control (if needed)
  # -t: target
  # -p: threshold (p < 10^-5)
  # -q: threshold (q-val) FDR-BH
  # -g: species, mm (=mouse), hs (=human)
  # -n: Exp_name
  # -B: bigBed file output
  # --outdir
OUTPUT
*_model.r estimated DNA flagment length
*_peaks.narrowPeak peak region
*_peaks.xls peak region
*_summits.bed summit of peak
*_treat_pileup.bdg read coverage of ChIP data
*_lambda.bdg local bias of read coverage of control data

参考文献