Directories and files
Expand all   Collapse all
FilesPath
Alignment
Replicate 1
trimmed fastq1 ./align/rep1/Donor-1-day-224-ATAC_read1.trim.fastq.gz
trimmed fastq2 ./align/rep1/Donor-1-day-224-ATAC_read2.trim.fastq.gz
bam ./align/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.bam
filtered bam ./align/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.bam
bedpe ./align/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.bedpe.gz
tagalign ./align/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tagAlign.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
tagalign ./align/pseudo_reps/rep1/pr1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.pr1.tagAlign.gz
Pseudo-replicate 2
tagalign ./align/pseudo_reps/rep1/pr2/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.pr2.tagAlign.gz
Signal tracks
MACS2
Replicate 1
p-value ./signal/macs2/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tn5.pf.fc.signal.bigwig
fold enrichment ./signal/macs2/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tn5.pf.pval.signal.bigwig
Peaks
MACS2
Replicate 1
gapped peak (p-val thld: 0.01) ./peak/macs2/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tn5.pf.gappedPeak.gz
gapped peak ./peak/macs2/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tn5.pf.pval0.1.gappedPeak.gz
peak ./peak/macs2/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tn5.pf.pval0.1.narrowPeak.gz
peak (p-val thld: 0.01) ./peak/macs2/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.tn5.pf.narrowPeak.gz
Pseudo-replicates
Replicate 1
Pseudo-replicate 1
gapped peak ./peak/macs2/pseudo_reps/rep1/pr1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.pr1.tn5.pf.pval0.1.gappedPeak.gz
peak ./peak/macs2/pseudo_reps/rep1/pr1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.pr1.tn5.pf.pval0.1.narrowPeak.gz
Pseudo-replicate 2
gapped peak ./peak/macs2/pseudo_reps/rep1/pr2/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.pr2.tn5.pf.pval0.1.gappedPeak.gz
peak ./peak/macs2/pseudo_reps/rep1/pr2/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.0.pr2.tn5.pf.pval0.1.narrowPeak.gz
QC and logs
Replicate 1
bowtie2 QC ./qc/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.align.log
dup. QC ./qc/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.dup.qc
map QC ./qc/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.flagstat.qc
PBC QC ./qc/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.pbc.qc
cross-corr. QC ./qc/rep1/Donor-1-day-224-ATAC_read1.trim.PE2SE.nodup.25.cc.qc
ATAQC ./qc/rep1/Donor-1-day-224-ATAC_read1.PE2SE_qc.html

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
5810675297.83

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

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

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

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
2274953100.305241919699402550.325388315000.18340721.6642850.85810550

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