Directories and files
Expand all   Collapse all
FilesPath
Alignment
Replicate 1
bam ./align/rep1/SRR1370876.bam
filtered bam ./align/rep1/SRR1370876.nodup.bam
tagalign ./align/rep1/SRR1370876.nodup.tagAlign.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
tagalign ./align/pseudo_reps/rep1/pr1/SRR1370876.nodup.pr1.tagAlign.gz
Pseudo-replicate 2
tagalign ./align/pseudo_reps/rep1/pr2/SRR1370876.nodup.pr2.tagAlign.gz
Peaks
IDR
Pseudo-replicates
Replicate 1
IDR peak ./peak/idr/pseudo_reps/rep1/rep1-pr.IDR0.05.filt.narrowPeak.gz
Unthresholded IDR peak ./peak/idr/pseudo_reps/rep1/rep1-pr.unthresholded-peaks.txt.gz
SPP
Replicate 1
peak ./peak/spp/rep1/SRR1370876.nodup.tagAlign_x_ControlAllNew.sub50MilPlus.tagAlign.regionPeak.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
peak ./peak/spp/pseudo_reps/rep1/pr1/SRR1370876.nodup.pr1.tagAlign_x_ControlAllNew.sub50MilPlus.tagAlign.regionPeak.gz
Pseudo-replicate 2
peak ./peak/spp/pseudo_reps/rep1/pr2/SRR1370876.nodup.pr2.tagAlign_x_ControlAllNew.sub50MilPlus.tagAlign.regionPeak.gz
MACS2
Naive overlap
peak ./peak/spp/overlap/SRR1370876.nodup.tagAlign_x_ControlAllNew.sub50MilPlus.tagAlign.naive_overlap.regionPeak.gz
QC and logs
Replicate 1
dup. QC ./qc/rep1/SRR1370876.dup.qc
map QC ./qc/rep1/SRR1370876.nodup.flagstat.qc
PBC QC ./qc/rep1/SRR1370876.nodup.pbc.qc
cross-corr. QC ./qc/rep1/SRR1370876.nodup.15M.cc.qc
IDR QC ./qc/IDR_final.qc

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

Workflow diagram (rep : Replicate, ctl : Control, pr : Pseudo-replicate, ppr : Pooled pseudo-replicate)


Flagstat QC (all)
 Reads (QC-passed)Reads (QC-failed)Dupes (QC-passed)Dupes (QC-failed)Mapped Reads% Mapped
rep1
173475990001706172798.35

Dup QC (all)
 Unpaired ReadsPaired ReadsUnmapped ReadsUnpaired DupesPaired DupesPaired Opt. Dupes% Dupes
rep1
14353009002219191000.154615

Flagstat QC (all, filtered)
 Reads (QC-passed)Reads (QC-failed)Dupes (QC-passed)Dupes (QC-failed)Mapped Reads% Mapped
rep1
1213381800012133818100.00

PBC QC (all)
 Total Read PairsDistinct Read PairsOne Read PairTwo Read PairsNRF = Distinct/TotalPBC1 = OnePair/DistinctPBC2 = OnePair/TwoPair
rep1
14339806122291341039392015955450.8528100.8499316.514338

NRF (non redundant fraction)
PBC1 (PCR Bottleneck coefficient 1)
PBC2 (PCR Bottleneck coefficient 2)
PBC1 is the primary measure. Provisionally

  • 0-0.5 is severe bottlenecking
  • 0.5-0.8 is moderate bottlenecking
  • 0.8-0.9 is mild bottlenecking
  • 0.9-1.0 is no bottlenecking


Cross-correlation QC (all)
 numReadsestFragLencorr_estFragLenPhantomPeakcorr_phantomPeakargmin_corrmin_corrphantomPeakCoefrelPhantomPeakCoefQualityTag
rep1
121258801350.167513264057672500.147894915000.14522311.1534898.342542

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
Important columns highlighted, but all/whole file can be stored for display


rep1
 

IDR QC (idr)
IDR_final.qc
NtN1Npconservative_setoptimal_setrescue_ratioself_consistency_ratioreproducibility
011010N/AN/ANaN1.01

rep1 pseudo-reps