ATAqC

Sample Information

Sample Day6-D1-FrozenKera-13B-TGGATCTG-GCGATCTA
Genome hg19
Paired/Single-ended Paired-ended
Read length 76

Summary

Read count from sequencer 23,402,650
Read count successfully aligned 22,535,641
Read count after filtering for mapping quality 18,154,333
Read count after removing duplicate reads 15,742,985
Read count after removing mitochondrial reads (final read count) 12,489,724
Note that all these read counts are determined using 'samtools view' - as such,
these are all reads found in the file, whether one end of a pair or a single 
end read. In other words, if your file is paired end, then you should divide 
these counts by two. Each step follows the previous step; for example, the 
duplicate reads were removed after reads were removed for low mapping quality.
This bar chart also shows the filtering process and where the reads were lost
over the process. Note that each step is sequential - as such, there may
have been more mitochondrial reads which were already filtered because of
high duplication or low mapping quality. Note that all these read counts are 
determined using 'samtools view' - as such, these are all reads found in 
the file, whether one end of a pair or a single end read. In other words, 
if your file is paired end, then you should divide these counts by two.

Alignment statistics

Bowtie alignment log

11701325 reads; of these:
  11701325 (100.00%) were paired; of these:
    615078 (5.26%) aligned concordantly 0 times
    7848259 (67.07%) aligned concordantly exactly 1 time
    3237988 (27.67%) aligned concordantly >1 times
    ----
    615078 pairs aligned concordantly 0 times; of these:
      38538 (6.27%) aligned discordantly 1 time
    ----
    576540 pairs aligned 0 times concordantly or discordantly; of these:
      1153080 mates make up the pairs; of these:
        867009 (75.19%) aligned 0 times
        153040 (13.27%) aligned exactly 1 time
        133031 (11.54%) aligned >1 times
96.30% overall alignment rate

  

Samtools flagstat

23402650 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 duplicates
22535641 + 0 mapped (96.30%:-nan%)
23402650 + 0 paired in sequencing
11701325 + 0 read1
11701325 + 0 read2
22172494 + 0 properly paired (94.74%:-nan%)
22329170 + 0 with itself and mate mapped
206471 + 0 singletons (0.88%:-nan%)
71386 + 0 with mate mapped to a different chr
33405 + 0 with mate mapped to a different chr (mapQ>=5)

  

Filtering statistics

Mapping quality > q30 (out of total) 18,154,333 0.776
Duplicates (after filtering) 2,411,348.0 0.267
Mitochondrial reads (out of total) 727,334 0.055
Final reads (after all filters) 12,489,724 0.534
Mapping quality refers to the quality of the read being aligned to that 
particular location in the genome. A standard quality score is > 30. 
Duplications are often due to PCR duplication rather than two unique reads
mapping to the same location. High duplication is an indication of poor 
libraries. Mitochondrial reads are often high in chromatin accessibility
assays because the mitochondrial genome is very open. A high mitochondrial
fraction is an indication of poor libraries. Based on prior experience, a
final read fraction above 0.70 is a good library.
  

Library complexity statistics

ENCODE library complexity metrics

Metric Result
NRF 0.70558843551 out of range [0.8, inf]
PBC1 0.783076191797 out of range [0.8, inf]
PBC2 4.97267911619 - OK
The non-redundant fraction (NRF) is the fraction of non-redundant mapped reads 
in a dataset; it is the ratio between the number of positions in the genome 
that uniquely mapped reads map to and the total number of uniquely mappable 
reads. The NRF should be > 0.8. The PBC1 is the ratio of genomic locations
with EXACTLY one read pair over the genomic locations with AT LEAST one read 
pair. PBC1 is the primary measure, and the PBC1 should be close to 1. 
Provisionally 0-0.5 is severe bottlenecking, 0.5-0.8 is moderate bottlenecking, 
0.8-0.9 is mild bottlenecking, and 0.9-1.0 is no bottlenecking. The PBC2 is 
the ratio of genomic locations with EXACTLY one read pair over the genomic 
locations with EXACTLY two read pairs. The PBC2 should be significantly 
greater than 1.

Picard EstimateLibraryComplexity

20,322,065

Yield prediction

Preseq performs a yield prediction by subsampling the reads, calculating the
number of distinct reads, and then extrapolating out to see where the 
expected number of distinct reads no longer increases. The confidence interval
gives a gauge as to the validity of the yield predictions.

Fragment length statistics

Metric Result
Fraction of reads in NFR 0.652103851231 - OK
NFR / mono-nuc reads 3.16790402172 - OK
Presence of NFR peak OK
Presence of Mono-Nuc peak OK
Presence of Di-Nuc peak OK
Open chromatin assays show distinct fragment length enrichments, as the cut 
sites are only in open chromatin and not in nucleosomes. As such, peaks 
representing different n-nucleosomal (ex mono-nucleosomal, di-nucleosomal) 
fragment lengths will arise. Good libraries will show these peaks in a 
fragment length distribution and will show specific peak ratios.

Sequence quality metrics

GC bias

Open chromatin assays are known to have significant GC bias. Please take this
into consideration as necessary.

Annotation-based quality metrics

Enrichment plots (TSS)

Open chromatin assays should show enrichment in open chromatin sites, such as
TSS's. An average TSS enrichment is above 6-7. A strong TSS enrichment is 
above 10.
  

Annotated genomic region enrichments

Fraction of reads in universal DHS regions 2,896,770 0.232
Fraction of reads in blacklist regions 42,537 0.003
Fraction of reads in promoter regions 563,151 0.045
Fraction of reads in enhancer regions 3,654,471 0.293
Fraction of reads in called peak regions 2,277,653 0.182
Signal to noise can be assessed by considering whether reads are falling into
known open regions (such as DHS regions) or not. A high fraction of reads 
should fall into the universal (across cell type) DHS set. A small fraction
should fall into the blacklist regions. A high set (though not all) should
fall into the promoter regions. A high set (though not all) should fall into 
the enhancer regions. The promoter regions should not take up all reads, as
it is known that there is a bias for promoters in open chromatin assays.

Comparison to Roadmap DNase

This bar chart shows the correlation between the Roadmap DNase samples to
your sample, when the signal in the universal DNase peak region sets are 
compared. The closer the sample is in signal distribution in the regions
to your sample, the higher the correlation.