Current implementation is deployed to the DNAnexus platform.

1. Map reads with BWA, mark duplicates Picard, and remove duplicates.
2. Estimate library complexity and calculate calculate NRF (non-redundant fraction), PBC1, PBC2 (PCR bottleneck coefficient).
3. Calculate cross-correlation analysis with spp/phantompeakqualtools.
4. Generate p-value and fold-over-control signal tracks for each replicate and replicates pooled with MACS2.

1. Call peaks with MACS2.
2. Calculate and report overlapping peaks from both replicates.

1. Call peaks with SPP.
2. Threshold peaks with IDR.
3. Report IDR-thresholded peak sets, self-consistency ratio, rescue ratio, reproducibility test.