Directories and files
Expand all   Collapse all
FilesPath
Alignment
Replicate 1
trimmed fastq1 ./align/rep1/ATAC13-trimmed-pair1.trim.fastq.gz
trimmed fastq2 ./align/rep1/ATAC13-trimmed-pair2.trim.fastq.gz
bam ./align/rep1/ATAC13-trimmed-pair1.trim.PE2SE.bam
filtered bam ./align/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.bam
sorted bam ./align/rep1/ATAC13-trimmed-pair1.trim.PE2SE.srt.bam
bedpe ./align/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.bedpe.gz
tagalign ./align/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tagAlign.gz
Replicate 2
trimmed fastq1 ./align/rep2/ATAC14-trimmed-pair1.trim.fastq.gz
trimmed fastq2 ./align/rep2/ATAC14-trimmed-pair2.trim.fastq.gz
bam ./align/rep2/ATAC14-trimmed-pair1.trim.PE2SE.bam
filtered bam ./align/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.bam
sorted bam ./align/rep2/ATAC14-trimmed-pair1.trim.PE2SE.srt.bam
bedpe ./align/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.bedpe.gz
tagalign ./align/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tagAlign.gz
Signal tracks
MACS2
Replicate 1
p-value ./signal/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.fc.signal.bigwig
fold enrichment ./signal/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval.signal.bigwig
Replicate 2
p-value ./signal/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.fc.signal.bigwig
fold enrichment ./signal/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval.signal.bigwig
Pooled replicate
p-value ./signal/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.fc.signal.bigwig
fold enrichment ./signal/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval.signal.bigwig
Peaks
IDR
True replicates
Rep. 1 vs Rep. 2
IDR peak ./peak/idr/true_reps/rep1-rep2/rep1-rep2.IDR0.1.filt.narrowPeak.gz
Unthresholded IDR peak ./peak/idr/true_reps/rep1-rep2/rep1-rep2.unthresholded-peaks.txt.gz
Optimal set
IDR peak ./peak/idr/optimal_set/rep1-rep2.IDR0.1.filt.narrowPeak.gz
Conservative set
IDR peak ./peak/idr/conservative_set/rep1-rep2.IDR0.1.filt.narrowPeak.gz
MACS2
Replicate 1
gapped peak (p-val thld: 0.01) ./peak/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.gappedPeak.gz
gapped peak ./peak/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.gappedPeak.gz
peak ./peak/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.narrowPeak.gz
peak (p-val thld: 0.01) ./peak/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.narrowPeak.gz
filtered peak ./peak/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.500000.narrowPeak.gz
filtered gapped peak ./peak/macs2/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.500000.gappedPeak.gz
Replicate 2
gapped peak (p-val thld: 0.01) ./peak/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.gappedPeak.gz
gapped peak ./peak/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.gappedPeak.gz
peak ./peak/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.narrowPeak.gz
peak (p-val thld: 0.01) ./peak/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.narrowPeak.gz
filtered peak ./peak/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.500000.narrowPeak.gz
filtered gapped peak ./peak/macs2/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.17.tn5.pf.pval0.1.500000.gappedPeak.gz
Pooled replicate
gapped peak (p-val thld: 0.01) ./peak/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.gappedPeak.gz
gapped peak ./peak/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval0.1.gappedPeak.gz
peak ./peak/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval0.1.narrowPeak.gz
peak (p-val thld: 0.01) ./peak/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.narrowPeak.gz
filtered peak ./peak/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval0.1.500000.narrowPeak.gz
filtered gapped peak ./peak/macs2/pooled_rep/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval0.1.500000.gappedPeak.gz
Naive overlap
peak ./peak/macs2/overlap/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval0.1.500000.naive_overlap.narrowPeak.gz
gapped peak ./peak/macs2/overlap/ATAC13-trimmed-pair1.trim.PE2SE.nodup.17.tn5_pooled.pf.pval0.1.500000.naive_overlap.gappedPeak.gz
QC and logs
Replicate 1
bowtie2 QC ./qc/rep1/ATAC13-trimmed-pair1.trim.PE2SE.align.log
dup. QC ./qc/rep1/ATAC13-trimmed-pair1.trim.PE2SE.dup.qc
map QC ./qc/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.flagstat.qc
PBC QC ./qc/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.pbc.qc
cross-corr. QC ./qc/rep1/ATAC13-trimmed-pair1.trim.PE2SE.nodup.25.cc.qc
ATAQC ./qc/rep1/ATAC13-trimmed-pair1.PE2SE_qc.html
Replicate 2
bowtie2 QC ./qc/rep2/ATAC14-trimmed-pair1.trim.PE2SE.align.log
dup. QC ./qc/rep2/ATAC14-trimmed-pair1.trim.PE2SE.dup.qc
map QC ./qc/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.flagstat.qc
PBC QC ./qc/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.pbc.qc
cross-corr. QC ./qc/rep2/ATAC14-trimmed-pair1.trim.PE2SE.nodup.25.cc.qc
ATAQC ./qc/rep2/ATAC14-trimmed-pair1.PE2SE_qc.html
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)


Bowtie2 QC (all)
 Total Reads% Aligned
rep1
4851428593.06
rep2
4567127992.03

Dup QC (all)
 Unpaired ReadsPaired ReadsUnmapped ReadsUnpaired DupesPaired DupesPaired Opt. Dupes% Dupes
rep1
031490431001272075670546170.403956
rep2
029029646001157012264714780.398562

Flagstat QC (all, filtered)
 Reads (QC-passed)Reads (QC-failed)Dupes (QC-passed)Dupes (QC-failed)Mapped Reads% Mapped
rep1
3753935000037539350100.00
rep2
3491904800034919048100.00

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

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
1800088300.222849954004333650.227639915000.18521981.2031650.88708290
rep2
1676125100.208967806844454650.213672715000.17668451.1827180.87280150

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
 
rep2
 

IDR QC (idr)
IDR_final.qc
NtNpconservative_setoptimal_setrescue_ratioself_consistency_ratioreproducibility
1163120rep1-rep2rep1-rep2InfinityNaN0

true reps (rep1-rep2)