Title: humanMonocyteDHS

Report generated at Tue Jul 3 15:42:43 EDT 2018

Command line arguments

-title humanMonocyteDHS -nth 1 -species hg38 -dnase_seq -se -enable_idr -tag1 /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteCrawfordFemaleUnknown/align/rep1/ENCFF000TBL_trimmed.nodup.tagAlign.gz -tag2 /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteCrawfordFemaleUnknown/align/rep2/ENCFF000TBP_trimmed.nodup.tagAlign.gz -tag3 /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteStamFemale34/align/rep1/ENCFF231SIM.nodup.tagAlign.gz -tag4 /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteStamFemaleUnknown/align/rep1/MonocyteStamFemaleUnknown.nodup.tagAlign.gz -tag5 /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteStamMale37/align/rep1/ENCFF690LNW.nodup.tagAlign.gz -out_dir /projects/pfenninggroup/jemmie/ml_human_monocyte_brain/humanMonocyteDHS

Pipeline version

Latest git commit SHA1: 56c48f22e03b2cc4c37b2e29db46ad8658e97c28 (Fri Apr 27 18:28:27 2018)


Directories and files
Expand all   Collapse all
FilesPath
Alignment
Replicate 1
Tag-align /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteCrawfordFemaleUnknown/align/rep1/ENCFF000TBL_trimmed.nodup.tagAlign.gz
Replicate 2
Tag-align /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteCrawfordFemaleUnknown/align/rep2/ENCFF000TBP_trimmed.nodup.tagAlign.gz
Replicate 3
Tag-align /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteStamFemale34/align/rep1/ENCFF231SIM.nodup.tagAlign.gz
Replicate 4
Tag-align /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteStamFemaleUnknown/align/rep1/MonocyteStamFemaleUnknown.nodup.tagAlign.gz
Replicate 5
Tag-align /projects/pfenninggroup/machineLearningForComputationalBiology/regElEvoGrant/HumanDNase/MonocyteStamMale37/align/rep1/ENCFF690LNW.nodup.tagAlign.gz
Pooled replicate
Tag-align ./align/pooled_rep/ENCFF000TBL_trimmed.nodup_pooled.tagAlign.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
Tag-align ./align/pseudo_reps/rep1/pr1/ENCFF000TBL_trimmed.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
Tag-align ./align/pseudo_reps/rep1/pr2/ENCFF000TBL_trimmed.nodup.pr2.tagAlign.gz
Replicate 2
Pseudo-replicate 1
Tag-align ./align/pseudo_reps/rep2/pr1/ENCFF000TBP_trimmed.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
Tag-align ./align/pseudo_reps/rep2/pr2/ENCFF000TBP_trimmed.nodup.pr2.tagAlign.gz
Replicate 3
Pseudo-replicate 1
Tag-align ./align/pseudo_reps/rep3/pr1/ENCFF231SIM.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
Tag-align ./align/pseudo_reps/rep3/pr2/ENCFF231SIM.nodup.pr2.tagAlign.gz
Replicate 4
Pseudo-replicate 1
Tag-align ./align/pseudo_reps/rep4/pr1/MonocyteStamFemaleUnknown.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
Tag-align ./align/pseudo_reps/rep4/pr2/MonocyteStamFemaleUnknown.nodup.pr2.tagAlign.gz
Replicate 5
Pseudo-replicate 1
Tag-align ./align/pseudo_reps/rep5/pr1/ENCFF690LNW.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
Tag-align ./align/pseudo_reps/rep5/pr2/ENCFF690LNW.nodup.pr2.tagAlign.gz
Pooled pseudo-replicates
Pooled pseudo-replicate 1
Tag-align ./align/pooled_pseudo_reps/ppr1/ENCFF000TBL_trimmed.nodup.pr1_pooled.tagAlign.gz
Pooled pseudo-replicate 2
Tag-align ./align/pooled_pseudo_reps/ppr2/ENCFF000TBL_trimmed.nodup.pr2_pooled.tagAlign.gz
Signal tracks
MACS2
Replicate 1
P-value ./signal/macs2/rep1/ENCFF000TBL_trimmed.nodup.pf.pval.signal.bigwig
Fold enrichment ./signal/macs2/rep1/ENCFF000TBL_trimmed.nodup.pf.fc.signal.bigwig
Replicate 2
P-value ./signal/macs2/rep2/ENCFF000TBP_trimmed.nodup.pf.pval.signal.bigwig
Fold enrichment ./signal/macs2/rep2/ENCFF000TBP_trimmed.nodup.pf.fc.signal.bigwig
Replicate 3
P-value ./signal/macs2/rep3/ENCFF231SIM.nodup.pf.pval.signal.bigwig
Fold enrichment ./signal/macs2/rep3/ENCFF231SIM.nodup.pf.fc.signal.bigwig
Replicate 4
P-value ./signal/macs2/rep4/MonocyteStamFemaleUnknown.nodup.pf.pval.signal.bigwig
Fold enrichment ./signal/macs2/rep4/MonocyteStamFemaleUnknown.nodup.pf.fc.signal.bigwig
Replicate 5
P-value ./signal/macs2/rep5/ENCFF690LNW.nodup.pf.pval.signal.bigwig
Fold enrichment ./signal/macs2/rep5/ENCFF690LNW.nodup.pf.fc.signal.bigwig
Pooled replicate
P-value ./signal/macs2/pooled_rep/ENCFF000TBL_trimmed.nodup_pooled.pf.pval.signal.bigwig
Fold enrichment ./signal/macs2/pooled_rep/ENCFF000TBL_trimmed.nodup_pooled.pf.fc.signal.bigwig
Peaks
MACS2
Replicate 1
Narrow peak (p-val thresh=0.01) ./peak/macs2/rep1/ENCFF000TBL_trimmed.nodup.pf.pval0.01.300K.narrowPeak.gz
Replicate 2
Narrow peak (p-val thresh=0.01) ./peak/macs2/rep2/ENCFF000TBP_trimmed.nodup.pf.pval0.01.300K.narrowPeak.gz
Replicate 3
Narrow peak (p-val thresh=0.01) ./peak/macs2/rep3/ENCFF231SIM.nodup.pf.pval0.01.300K.narrowPeak.gz
Replicate 4
Narrow peak (p-val thresh=0.01) ./peak/macs2/rep4/MonocyteStamFemaleUnknown.nodup.pf.pval0.01.300K.narrowPeak.gz
Replicate 5
Narrow peak (p-val thresh=0.01) ./peak/macs2/rep5/ENCFF690LNW.nodup.pf.pval0.01.300K.narrowPeak.gz
Pooled replicate
Narrow peak (p-val thresh=.01) ./peak/macs2/pooled_rep/ENCFF000TBL_trimmed.nodup_pooled.pf.pval0.01.300K.narrowPeak.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
Narrow peak ./peak/macs2/pseudo_reps/rep1/pr1/ENCFF000TBL_trimmed.nodup.pr1.pf.pval0.01.300K.narrowPeak.gz
Pseudo-replicate 2
Narrow peak ./peak/macs2/pseudo_reps/rep1/pr2/ENCFF000TBL_trimmed.nodup.pr2.pf.pval0.01.300K.narrowPeak.gz
Replicate 2
Pseudo-replicate 1
Narrow peak ./peak/macs2/pseudo_reps/rep2/pr1/ENCFF000TBP_trimmed.nodup.pr1.pf.pval0.01.300K.narrowPeak.gz
Pseudo-replicate 2
Narrow peak ./peak/macs2/pseudo_reps/rep2/pr2/ENCFF000TBP_trimmed.nodup.pr2.pf.pval0.01.300K.narrowPeak.gz
Replicate 3
Pseudo-replicate 1
Narrow peak ./peak/macs2/pseudo_reps/rep3/pr1/ENCFF231SIM.nodup.pr1.pf.pval0.01.300K.narrowPeak.gz
Pseudo-replicate 2
Narrow peak ./peak/macs2/pseudo_reps/rep3/pr2/ENCFF231SIM.nodup.pr2.pf.pval0.01.300K.narrowPeak.gz
Replicate 4
Pseudo-replicate 1
Narrow peak ./peak/macs2/pseudo_reps/rep4/pr1/MonocyteStamFemaleUnknown.nodup.pr1.pf.pval0.01.300K.narrowPeak.gz
Pseudo-replicate 2
Narrow peak ./peak/macs2/pseudo_reps/rep4/pr2/MonocyteStamFemaleUnknown.nodup.pr2.pf.pval0.01.300K.narrowPeak.gz
Replicate 5
Pseudo-replicate 1
Narrow peak ./peak/macs2/pseudo_reps/rep5/pr1/ENCFF690LNW.nodup.pr1.pf.pval0.01.300K.narrowPeak.gz
Pseudo-replicate 2
Narrow peak ./peak/macs2/pseudo_reps/rep5/pr2/ENCFF690LNW.nodup.pr2.pf.pval0.01.300K.narrowPeak.gz
Pooled pseudo-replicates
Pooled pseudo-replicate 1
Narrow peak ./peak/macs2/pooled_pseudo_reps/ppr1/ENCFF000TBL_trimmed.nodup.pr1_pooled.pf.pval0.01.300K.narrowPeak.gz
Pooled pseudo-replicate 2
Narrow peak ./peak/macs2/pooled_pseudo_reps/ppr2/ENCFF000TBL_trimmed.nodup.pr2_pooled.pf.pval0.01.300K.narrowPeak.gz
Naive overlap
True replicates
Rep. 1 vs. Rep. 2
Narrow peak ./peak/macs2/overlap/true_reps/rep1-rep2/humanMonocyteDHS_rep1-rep2.naive_overlap.filt.narrowPeak.gz
Rep. 1 vs. Rep. 3
Narrow peak ./peak/macs2/overlap/true_reps/rep1-rep3/humanMonocyteDHS_rep1-rep3.naive_overlap.filt.narrowPeak.gz
Rep. 1 vs. Rep. 4
Narrow peak ./peak/macs2/overlap/true_reps/rep1-rep4/humanMonocyteDHS_rep1-rep4.naive_overlap.filt.narrowPeak.gz
Rep. 1 vs. Rep. 5
Narrow peak ./peak/macs2/overlap/true_reps/rep1-rep5/humanMonocyteDHS_rep1-rep5.naive_overlap.filt.narrowPeak.gz
Rep. 2 vs. Rep. 3
Narrow peak ./peak/macs2/overlap/true_reps/rep2-rep3/humanMonocyteDHS_rep2-rep3.naive_overlap.filt.narrowPeak.gz
Rep. 2 vs. Rep. 4
Narrow peak ./peak/macs2/overlap/true_reps/rep2-rep4/humanMonocyteDHS_rep2-rep4.naive_overlap.filt.narrowPeak.gz
Rep. 2 vs. Rep. 5
Narrow peak ./peak/macs2/overlap/true_reps/rep2-rep5/humanMonocyteDHS_rep2-rep5.naive_overlap.filt.narrowPeak.gz
Rep. 3 vs. Rep. 4
Narrow peak ./peak/macs2/overlap/true_reps/rep3-rep4/humanMonocyteDHS_rep3-rep4.naive_overlap.filt.narrowPeak.gz
Rep. 3 vs. Rep. 5
Narrow peak ./peak/macs2/overlap/true_reps/rep3-rep5/humanMonocyteDHS_rep3-rep5.naive_overlap.filt.narrowPeak.gz
Rep. 4 vs. Rep. 5
Narrow peak ./peak/macs2/overlap/true_reps/rep4-rep5/humanMonocyteDHS_rep4-rep5.naive_overlap.filt.narrowPeak.gz
Pseudo-replicates
Replicate 1
Narrow peak ./peak/macs2/overlap/pseudo_reps/rep1/humanMonocyteDHS_rep1-pr.naive_overlap.filt.narrowPeak.gz
Replicate 2
Narrow peak ./peak/macs2/overlap/pseudo_reps/rep2/humanMonocyteDHS_rep2-pr.naive_overlap.filt.narrowPeak.gz
Replicate 3
Narrow peak ./peak/macs2/overlap/pseudo_reps/rep3/humanMonocyteDHS_rep3-pr.naive_overlap.filt.narrowPeak.gz
Replicate 4
Narrow peak ./peak/macs2/overlap/pseudo_reps/rep4/humanMonocyteDHS_rep4-pr.naive_overlap.filt.narrowPeak.gz
Replicate 5
Narrow peak ./peak/macs2/overlap/pseudo_reps/rep5/humanMonocyteDHS_rep5-pr.naive_overlap.filt.narrowPeak.gz
Pooled pseudo-replicates
Narrow peak ./peak/macs2/overlap/pooled_pseudo_reps/humanMonocyteDHS_ppr.naive_overlap.filt.narrowPeak.gz
Optimal set
Overlapping peak ./peak/macs2/overlap/optimal_set/humanMonocyteDHS_ppr.naive_overlap.filt.narrowPeak.gz
Conservative set
Overlapping peak ./peak/macs2/overlap/conservative_set/humanMonocyteDHS_rep1-rep2.naive_overlap.filt.narrowPeak.gz
IDR
True replicates
Rep. 1 vs. Rep. 2
IDR peak ./peak/macs2/idr/true_reps/rep1-rep2/humanMonocyteDHS_rep1-rep2.IDR0.1.filt.narrowPeak.gz
Rep. 1 vs. Rep. 3
IDR peak ./peak/macs2/idr/true_reps/rep1-rep3/humanMonocyteDHS_rep1-rep3.IDR0.1.filt.narrowPeak.gz
Rep. 1 vs. Rep. 4
IDR peak ./peak/macs2/idr/true_reps/rep1-rep4/humanMonocyteDHS_rep1-rep4.IDR0.1.filt.narrowPeak.gz
Rep. 1 vs. Rep. 5
IDR peak ./peak/macs2/idr/true_reps/rep1-rep5/humanMonocyteDHS_rep1-rep5.IDR0.1.filt.narrowPeak.gz
Rep. 2 vs. Rep. 3
IDR peak ./peak/macs2/idr/true_reps/rep2-rep3/humanMonocyteDHS_rep2-rep3.IDR0.1.filt.narrowPeak.gz
Rep. 2 vs. Rep. 4
IDR peak ./peak/macs2/idr/true_reps/rep2-rep4/humanMonocyteDHS_rep2-rep4.IDR0.1.filt.narrowPeak.gz
Rep. 2 vs. Rep. 5
IDR peak ./peak/macs2/idr/true_reps/rep2-rep5/humanMonocyteDHS_rep2-rep5.IDR0.1.filt.narrowPeak.gz
Rep. 3 vs. Rep. 4
IDR peak ./peak/macs2/idr/true_reps/rep3-rep4/humanMonocyteDHS_rep3-rep4.IDR0.1.filt.narrowPeak.gz
Rep. 3 vs. Rep. 5
IDR peak ./peak/macs2/idr/true_reps/rep3-rep5/humanMonocyteDHS_rep3-rep5.IDR0.1.filt.narrowPeak.gz
Rep. 4 vs. Rep. 5
IDR peak ./peak/macs2/idr/true_reps/rep4-rep5/humanMonocyteDHS_rep4-rep5.IDR0.1.filt.narrowPeak.gz
Pseudo-replicates
Replicate 1
IDR peak ./peak/macs2/idr/pseudo_reps/rep1/humanMonocyteDHS_rep1-pr.IDR0.1.filt.narrowPeak.gz
Replicate 2
IDR peak ./peak/macs2/idr/pseudo_reps/rep2/humanMonocyteDHS_rep2-pr.IDR0.1.filt.narrowPeak.gz
Replicate 3
IDR peak ./peak/macs2/idr/pseudo_reps/rep3/humanMonocyteDHS_rep3-pr.IDR0.1.filt.narrowPeak.gz
Replicate 4
IDR peak ./peak/macs2/idr/pseudo_reps/rep4/humanMonocyteDHS_rep4-pr.IDR0.1.filt.narrowPeak.gz
Replicate 5
IDR peak ./peak/macs2/idr/pseudo_reps/rep5/humanMonocyteDHS_rep5-pr.IDR0.1.filt.narrowPeak.gz
Pooled pseudo-replicates
IDR peak ./peak/macs2/idr/pooled_pseudo_reps/humanMonocyteDHS_ppr.IDR0.1.filt.narrowPeak.gz
Optimal set
IDR peak ./peak/macs2/idr/optimal_set/humanMonocyteDHS_ppr.IDR0.1.filt.narrowPeak.gz
Conservative set
IDR peak ./peak/macs2/idr/conservative_set/humanMonocyteDHS_rep1-rep4.IDR0.1.filt.narrowPeak.gz
QC and logs
Replicate 1
Cross-corr. plot ./qc/rep1/ENCFF000TBL_trimmed.nodup.no_chrM.25M.cc.plot.pdf
Replicate 2
Cross-corr. plot ./qc/rep2/ENCFF000TBP_trimmed.nodup.no_chrM.25M.cc.plot.pdf
Replicate 3
Cross-corr. plot ./qc/rep3/ENCFF231SIM.nodup.no_chrM.25M.cc.plot.pdf
Replicate 4
Cross-corr. plot ./qc/rep4/MonocyteStamFemaleUnknown.nodup.no_chrM.25M.cc.plot.pdf
Replicate 5
Cross-corr. plot ./qc/rep5/ENCFF690LNW.nodup.no_chrM.25M.cc.plot.pdf
Naive overlap
Naive overlap QC log ./qc/humanMonocyteDHS_peak_overlap_final.qc
True replicates
Rep. 1 vs. Rep. 2
FRiP ./peak/macs2/overlap/true_reps/rep1-rep2/humanMonocyteDHS_rep1-rep2.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 1 vs. Rep. 3
FRiP ./peak/macs2/overlap/true_reps/rep1-rep3/humanMonocyteDHS_rep1-rep3.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 1 vs. Rep. 4
FRiP ./peak/macs2/overlap/true_reps/rep1-rep4/humanMonocyteDHS_rep1-rep4.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 1 vs. Rep. 5
FRiP ./peak/macs2/overlap/true_reps/rep1-rep5/humanMonocyteDHS_rep1-rep5.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 2 vs. Rep. 3
FRiP ./peak/macs2/overlap/true_reps/rep2-rep3/humanMonocyteDHS_rep2-rep3.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 2 vs. Rep. 4
FRiP ./peak/macs2/overlap/true_reps/rep2-rep4/humanMonocyteDHS_rep2-rep4.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 2 vs. Rep. 5
FRiP ./peak/macs2/overlap/true_reps/rep2-rep5/humanMonocyteDHS_rep2-rep5.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 3 vs. Rep. 4
FRiP ./peak/macs2/overlap/true_reps/rep3-rep4/humanMonocyteDHS_rep3-rep4.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 3 vs. Rep. 5
FRiP ./peak/macs2/overlap/true_reps/rep3-rep5/humanMonocyteDHS_rep3-rep5.naive_overlap.filt.narrowPeak.FRiP.qc
Rep. 4 vs. Rep. 5
FRiP ./peak/macs2/overlap/true_reps/rep4-rep5/humanMonocyteDHS_rep4-rep5.naive_overlap.filt.narrowPeak.FRiP.qc
Pseudo-replicates
Replicate 1
FRiP ./peak/macs2/overlap/pseudo_reps/rep1/humanMonocyteDHS_rep1-pr.naive_overlap.filt.narrowPeak.FRiP.qc
Replicate 2
FRiP ./peak/macs2/overlap/pseudo_reps/rep2/humanMonocyteDHS_rep2-pr.naive_overlap.filt.narrowPeak.FRiP.qc
Replicate 3
FRiP ./peak/macs2/overlap/pseudo_reps/rep3/humanMonocyteDHS_rep3-pr.naive_overlap.filt.narrowPeak.FRiP.qc
Replicate 4
FRiP ./peak/macs2/overlap/pseudo_reps/rep4/humanMonocyteDHS_rep4-pr.naive_overlap.filt.narrowPeak.FRiP.qc
Replicate 5
FRiP ./peak/macs2/overlap/pseudo_reps/rep5/humanMonocyteDHS_rep5-pr.naive_overlap.filt.narrowPeak.FRiP.qc
Pooled pseudo-replicates
FRiP ./peak/macs2/overlap/pooled_pseudo_reps/humanMonocyteDHS_ppr.naive_overlap.filt.narrowPeak.FRiP.qc
IDR
IDR QC log ./qc/humanMonocyteDHS_IDR_final.qc
True replicates
Rep. 1 vs. Rep. 2
FRiP ./peak/macs2/idr/true_reps/rep1-rep2/humanMonocyteDHS_rep1-rep2.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep1-rep2/humanMonocyteDHS_rep1-rep2.unthresholded-peaks.txt.png
Rep. 1 vs. Rep. 3
FRiP ./peak/macs2/idr/true_reps/rep1-rep3/humanMonocyteDHS_rep1-rep3.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep1-rep3/humanMonocyteDHS_rep1-rep3.unthresholded-peaks.txt.png
Rep. 1 vs. Rep. 4
FRiP ./peak/macs2/idr/true_reps/rep1-rep4/humanMonocyteDHS_rep1-rep4.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep1-rep4/humanMonocyteDHS_rep1-rep4.unthresholded-peaks.txt.png
Rep. 1 vs. Rep. 5
FRiP ./peak/macs2/idr/true_reps/rep1-rep5/humanMonocyteDHS_rep1-rep5.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep1-rep5/humanMonocyteDHS_rep1-rep5.unthresholded-peaks.txt.png
Rep. 2 vs. Rep. 3
FRiP ./peak/macs2/idr/true_reps/rep2-rep3/humanMonocyteDHS_rep2-rep3.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep2-rep3/humanMonocyteDHS_rep2-rep3.unthresholded-peaks.txt.png
Rep. 2 vs. Rep. 4
FRiP ./peak/macs2/idr/true_reps/rep2-rep4/humanMonocyteDHS_rep2-rep4.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep2-rep4/humanMonocyteDHS_rep2-rep4.unthresholded-peaks.txt.png
Rep. 2 vs. Rep. 5
FRiP ./peak/macs2/idr/true_reps/rep2-rep5/humanMonocyteDHS_rep2-rep5.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep2-rep5/humanMonocyteDHS_rep2-rep5.unthresholded-peaks.txt.png
Rep. 3 vs. Rep. 4
FRiP ./peak/macs2/idr/true_reps/rep3-rep4/humanMonocyteDHS_rep3-rep4.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep3-rep4/humanMonocyteDHS_rep3-rep4.unthresholded-peaks.txt.png
Rep. 3 vs. Rep. 5
FRiP ./peak/macs2/idr/true_reps/rep3-rep5/humanMonocyteDHS_rep3-rep5.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep3-rep5/humanMonocyteDHS_rep3-rep5.unthresholded-peaks.txt.png
Rep. 4 vs. Rep. 5
FRiP ./peak/macs2/idr/true_reps/rep4-rep5/humanMonocyteDHS_rep4-rep5.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/true_reps/rep4-rep5/humanMonocyteDHS_rep4-rep5.unthresholded-peaks.txt.png
Pseudo-replicates
Replicate 1
FRiP ./peak/macs2/idr/pseudo_reps/rep1/humanMonocyteDHS_rep1-pr.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/pseudo_reps/rep1/humanMonocyteDHS_rep1-pr.unthresholded-peaks.txt.png
Replicate 2
FRiP ./peak/macs2/idr/pseudo_reps/rep2/humanMonocyteDHS_rep2-pr.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/pseudo_reps/rep2/humanMonocyteDHS_rep2-pr.unthresholded-peaks.txt.png
Replicate 3
FRiP ./peak/macs2/idr/pseudo_reps/rep3/humanMonocyteDHS_rep3-pr.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/pseudo_reps/rep3/humanMonocyteDHS_rep3-pr.unthresholded-peaks.txt.png
Replicate 4
FRiP ./peak/macs2/idr/pseudo_reps/rep4/humanMonocyteDHS_rep4-pr.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/pseudo_reps/rep4/humanMonocyteDHS_rep4-pr.unthresholded-peaks.txt.png
Replicate 5
FRiP ./peak/macs2/idr/pseudo_reps/rep5/humanMonocyteDHS_rep5-pr.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/pseudo_reps/rep5/humanMonocyteDHS_rep5-pr.unthresholded-peaks.txt.png
Pooled pseudo-replicates
FRiP ./peak/macs2/idr/pooled_pseudo_reps/humanMonocyteDHS_ppr.IDR0.1.filt.narrowPeak.FRiP.qc
IDR plot ./peak/macs2/idr/pooled_pseudo_reps/humanMonocyteDHS_ppr.unthresholded-peaks.txt.png

Visualization (add -url_base [URL_ROOT_DIR_FOR_OUT_DIR] to the command line.)
Visualize  JSON (datahub)



Workflow diagram (g./b. peak : gapped/broad peak from MACS2, n. peak : narrow peak from MACS2, r. peak : relaxed peak from SPP)


Enrichment QC (strand cross-correlation measures)
 rep1
rep2
rep3
rep4
rep5
Reads2500000025000000191638172500000020027378
Est. Fragment Len.00000
Corr. Est. Fragment Len.0.2865472930698020.3432796196442540.3389168554288890.3504333551294850.338124673851978
Phantom Peak2020404040
Corr. Phantom Peak0.28380290.3369520.36229150.36655890.3650881
Argmin. Corr.15001500150015001500
Min. Corr.0.25551730.27139050.19050850.24107560.1948309
NSC1.121441.2648921.7790121.4536251.735477
RSC1.0970261.0965140.86392950.87149290.8416312

NOTE: Reads from replicates are subsampled to a max of 25M

  • Normalized strand cross-correlation coefficient (NSC) = col9 in outFile
  • Relative strand cross-correlation coefficient (RSC) = col10 in outFile
  • Estimated fragment length = col3 in outFile, take the top value


rep1
 
rep2
 
rep3
 
rep4
 
rep5
 
Number of peaks
rep1300000
rep1-pr1256104
rep1-pr2256125
rep2281071
rep2-pr1204884
rep2-pr2205027
rep3122352
rep3-pr1101875
rep3-pr2102385
rep4156862
rep4-pr1143268
rep4-pr2144083
rep5135362
rep5-pr1114524
rep5-pr2114419
pooled300000
ppr1300000
ppr2300000

  • ppr1: Raw peaks called on the first pooled pseudoreplicates
  • ppr2: Raw peaks called on the second pooled pseudoreplicates
  • repi: Raw peaks called on true replicate i
  • repi-pr1 : Raw peaks called on the first pseudoreplicate from replicate i
  • repi-pr2 : Raw peaks called on the second pseudoreplicates from replicate i
  • overlap : Overlapped peaks (filtered if blacklist exists)


Enrichment QC (Fraction of reads in overlapping peaks)
 ppr
rep1-pr
rep1-rep2
rep1-rep3
rep1-rep4
rep1-rep5
rep2-pr
rep2-rep3
rep2-rep4
rep2-rep5
rep3-pr
rep3-rep4
rep3-rep5
rep4-pr
rep4-rep5
rep5-pr
Fraction of Reads in Peak0.23270.1245690.2067890.1748910.1891130.174070.2200550.1659060.1780240.1638690.2550030.1794130.1741290.2688490.1808920.264335

  • ppr: Overlapping peaks comparing pooled pseudo replicates
  • rep1-pr: Overlapping peaks comparing pseudoreplicates from replicate 1
  • rep2-pr: Overlapping peaks comparing pseudoreplicates from replicate 2
  • repi-repj: Overlapping peaks comparing true replicates (rep i vs. rep j)


Reproducibility QC and Peak Detection Statistics (Overlapping peaks)
Nt178376
N1157157
N2147838
Np223202
N optimal223202
N conservative178376
Optimal Setpooled_pseudo_rep
Conservative Setrep1-rep2
Rescue Ratio1.2513006234022515
Self Consistency Ratio2.209961610394724
Reproducibility Testborderline

  • N1: Replicate 1 self-consistent overlapping peaks (comparing two pseudoreplicates generated by subsampling Rep1 reads)
  • N2: Replicate 2 self-consistent overlapping peaks (comparing two pseudoreplicates generated by subsampling Rep2 reads)
  • Nt: True Replicate consisten overlapping peaks (comparing true replicates Rep1 vs Rep2 )
  • Np: Pooled-pseudoreplicate consistent overlapping peaks (comparing two pseudoreplicates generated by subsampling pooled reads from Rep1 and Rep2 )
  • Self-consistency Ratio: max(N1,N2) / min (N1,N2)
  • Rescue Ratio: max(Np,Nt) / min (Np,Nt)
  • Reproducibility Test: If Self-consistency Ratio >2 AND Rescue Ratio > 2, then 'Fail' else 'Pass'


Enrichment QC (Fraction of reads in IDR peaks)
 ppr
rep1-pr
rep1-rep2
rep1-rep3
rep1-rep4
rep1-rep5
rep2-pr
rep2-rep3
rep2-rep4
rep2-rep5
rep3-pr
rep3-rep4
rep3-rep5
rep4-pr
rep4-rep5
rep5-pr
Fraction of Reads in Peak0.1923150.09500920.1572810.1662650.1674520.1570310.1658150.1465150.1424880.1427820.2248590.1515730.1466060.2436880.1534850.233109

  • ppr: IDR peaks comparing pooled pseudo replicates
  • rep1-pr: IDR peaks comparing pseudoreplicates from replicate 1
  • rep2-pr: IDR peaks comparing pseudoreplicates from replicate 2
  • repi-repj: IDR peaks comparing true replicates (rep i vs. rep j)


Reproducibility QC and Peak Detection Statistics (Irreproducible Discovery Rate)
Nt91131
N186979
N277554
Np126727
N optimal126727
N conservative91131
Optimal Setpooled_pseudo_rep
Conservative Setrep1-rep4
Rescue Ratio1.3906025392018084
Self Consistency Ratio1.8730134803393772
Reproducibility Testpass

  • N1: Replicate 1 self-consistent IDR 0.1 peaks (comparing two pseudoreplicates generated by subsampling Rep1 reads)
  • N2: Replicate 2 self-consistent IDR 0.1 peaks (comparing two pseudoreplicates generated by subsampling Rep2 reads)
  • Nt: True Replicate consistent IDR 0.1 peaks (comparing true replicates Rep1 vs Rep2 )
  • Np: Pooled-pseudoreplicate consistent IDR 0.1 peaks (comparing two pseudoreplicates generated by subsampling pooled reads from Rep1 and Rep2 )
  • Self-consistency Ratio: max(N1,N2) / min (N1,N2)
  • Rescue Ratio: max(Np,Nt) / min (Np,Nt)
  • Reproducibility Test: If Self-consistency Ratio >2 AND Rescue Ratio > 2, then 'Fail' else 'Pass'


true reps (rep1-rep2)
 
true reps (rep1-rep3)
 
true reps (rep1-rep4)
 
true reps (rep1-rep5)
 
true reps (rep2-rep3)
 
true reps (rep2-rep4)
 
true reps (rep2-rep5)
 
true reps (rep3-rep4)
 
true reps (rep3-rep5)
 
true reps (rep4-rep5)
 
pooled pseudo-reps
 
rep1 pseudo-reps
 
rep2 pseudo-reps
 
rep3 pseudo-reps
 
rep4 pseudo-reps
 
rep5 pseudo-reps