Changes between Version 22 and Version 23 of SOPs/chip_seq_peaks
- Timestamp:
- 01/29/16 15:05:16 (9 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/chip_seq_peaks
v22 v23 39 39 * MACS2 is appropriate for both proteins like transcription factors that may have narrow peaks, as well as histone modifications that may affect broader regions. For broader peaks we recommend using --nomodel,--broad, --nolambda (if there's no control), and using the fragment size calculated on the strand cross correlation analysis. We recommend using macs2 rather than macs14 for broad peaks. Identifying reproducible peaks (across replicates) with IDR requires MACS2 and works best with a relaxed p-value threshold. 40 40 {{{ 41 bsub macs2 callpeak -t IP_reads.mapped_only.bam -c Control_reads.mapped_only.bam -f BAM -g mm -n Name --nomodel -B --extsize "size calculated on t e strand croscorrelation analysis"41 bsub macs2 callpeak -t IP_reads.mapped_only.bam -c Control_reads.mapped_only.bam -f BAM -g mm -n Name --nomodel -B --extsize "size calculated on the strand croscorrelation analysis" 42 42 bsub macs2 callpeak -t IP_reads.mapped_only.bam -c Control_reads.mapped_only.bam -f BAM -g mm -n Name -B -p 1e-3 43 43 }}}