Title: ZNF684

Report generated at Tue May 9 22:42:18 PDT 2017

Command line arguments

-title ZNF684 -nth 5 -species hg38 -url_base http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/HughesZnfData/ZNF684/out -tag /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz -ctl_tag /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -final_stage idr

Pipeline version

Latest git commit SHA1: 22fb22592f96625fe282ca1ab9e189f4dac78157 (Tue May 9 20:51:05 2017)


Directories and files
Expand all   Collapse all
FilesPath
Alignment
Replicate 1
Tag-align ./align/rep1/SRR1370904_1.nodup.tagAlign.gz
Control 1
Tag-align /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
Tag-align ./align/pseudo_reps/rep1/pr1/SRR1370904_1.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
Tag-align ./align/pseudo_reps/rep1/pr2/SRR1370904_1.nodup.pr2.tagAlign.gz
Signal tracks
MACS2
Replicate 1
P-value ./signal/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.bw
Fold enrichment ./signal/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.bw
Peaks
SPP
Replicate 1
Region peak ./peak/spp/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.regionPeak.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
Region peak ./peak/spp/pseudo_reps/rep1/pr1/SRR1370904_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.regionPeak.gz
Pseudo-replicate 2
Region peak ./peak/spp/pseudo_reps/rep1/pr2/SRR1370904_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.regionPeak.gz
Naive overlap
Peak ./peak/spp/overlap/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.regionPeak.gz
IDR
Pseudo-replicates
Replicate 1
IDR peak ./peak/spp/idr/pseudo_reps/rep1/ZNF684_rep1-pr.IDR0.05.filt.narrowPeak.gz
Optimal set
IDR peak ./peak/spp/idr/optimal_set/ZNF684_rep1-pr.IDR0.05.filt.narrowPeak.gz
Conservative set
IDR peak ./peak/spp/idr/conservative_set/ZNF684_rep1-pr.IDR0.05.filt.narrowPeak.gz
MACS2
Replicate 1
Narrow peak ./peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.narrowPeak.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
Narrow peak ./peak/macs2/pseudo_reps/rep1/pr1/SRR1370904_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.narrowPeak.gz
Pseudo-replicate 2
Narrow peak ./peak/macs2/pseudo_reps/rep1/pr2/SRR1370904_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.narrowPeak.gz
QC and logs
Replicate 1
Cross-corr. log ./qc/rep1/SRR1370904_1.nodup.15M.cc.qc
Cross-corr. plot ./qc/rep1/SRR1370904_1.nodup.15M.cc.plot.pdf
IDR
IDR QC log ./qc/ZNF684_IDR_final.qc
Pseudo-replicates
Replicate 1
FRiP ./peak/spp/idr/pseudo_reps/rep1/ZNF684_rep1-pr.IDR0.05.filt.narrowPeak.FRiP.qc
IDR plot ./peak/spp/idr/pseudo_reps/rep1/ZNF684_rep1-pr.unthresholded-peaks.txt.png

Visualization
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)
Reads15000000
Est. Fragment Len.170
Corr. Est. Fragment Len.0.182626302070605
Phantom Peak50
Corr. Phantom Peak0.1748288
Argmin. Corr.1500
Min. Corr.0.1727319
NSC1.057282
RSC4.718593

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

  • 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
 
Number of peaks
rep1300000
rep1-pr1284205
rep1-pr2284262
overlap16340

  • 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 peaks)
Fraction of Reads in Peak0.000755717

  • 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)
Nt0
N1266
N20
Np0
N optimal266
N conservative266
Optimal Setself_pseudo_rep_rep1
Conservative Setself_pseudo_rep_rep1
Rescue Ratio0.0
Self Consistency Ratio1.0
Reproducibility Testpass

  • N1: Replicate 1 self-consistent IDR 0.05 peaks (comparing two pseudoreplicates generated by subsampling Rep1 reads)
  • N2: Replicate 2 self-consistent IDR 0.05 peaks (comparing two pseudoreplicates generated by subsampling Rep2 reads)
  • Nt: True Replicate consistent IDR 0.05 peaks (comparing true replicates Rep1 vs Rep2 )
  • Np: Pooled-pseudoreplicate consistent IDR 0.05 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'


rep1 pseudo-reps