PUT

ATAC-seq data from putamen

Pipeline version: v1.4.1

Report generated at 2019-06-14 07:55:02

Paired-end: [True, True, True, True]

Pipeline type: ATAC-Seq

Genome: hg38.tsv

Peak caller: MACS2

Alignment

Flagstat (raw BAM)

Marking duplicates (filtered BAM)

Library complexity (filtered non-mito BAM)

Flagstat (filtered/deduped BAM)

Filtered and duplicates removed

Peak calling

IDR (Irreproducible Discovery Rate) plots

Reproducibility QC and peak detection statistics

Enrichment

Strand cross-correlation measures

Performed on subsampled reads (25M)

Fraction of reads in overlapping peaks

Fraction of reads in IDR peaks

ATAQC

Summary table

Replicate 1

Sample Information

Summary

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

27247094 reads; of these:
  27247094 (100.00%) were paired; of these:
    41339 (0.15%) aligned concordantly 0 times
    20628768 (75.71%) aligned concordantly exactly 1 time
    6576987 (24.14%) aligned concordantly >1 times
    ----
    41339 pairs aligned concordantly 0 times; of these:
      6053 (14.64%) aligned discordantly 1 time
    ----
    35286 pairs aligned 0 times concordantly or discordantly; of these:
      70572 mates make up the pairs; of these:
        50402 (71.42%) aligned 0 times
        3820 (5.41%) aligned exactly 1 time
        16350 (23.17%) aligned >1 times
99.91% overall alignment rate

  

Samtools flagstat

95444864 + 0 in total (QC-passed reads + QC-failed reads)
40950676 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
95394462 + 0 mapped (99.95%:-nan%)
54494188 + 0 paired in sequencing
27247094 + 0 read1
27247094 + 0 read2
54411510 + 0 properly paired (99.85%:-nan%)
54439110 + 0 with itself and mate mapped
4676 + 0 singletons (0.01%:-nan%)
5318 + 0 with mate mapped to a different chr
1996 + 0 with mate mapped to a different chr (mapQ>=5)

  

Note that the flagstat command counts alignments, not reads. please 
use the read counts table to get accurate counts of reads at each
stage of the pipeline.

Filtering statistics

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

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

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

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.

Peak statistics

Naive overlap peak file statistics

IDR peak file statistics

For a good ATAC-seq experiment in human, you expect to get 100k-200k peaks
for a specific cell type.

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 in human (hg19) is above 6. A strong TSS enrichment is
above 10. For other references please see https://www.encodeproject.org/atac-seq/
  

Annotated genomic region enrichments

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.

Replicate 2

Sample Information

Summary

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

52712090 reads; of these:
  52712090 (100.00%) were paired; of these:
    129825 (0.25%) aligned concordantly 0 times
    39176936 (74.32%) aligned concordantly exactly 1 time
    13405329 (25.43%) aligned concordantly >1 times
    ----
    129825 pairs aligned concordantly 0 times; of these:
      8677 (6.68%) aligned discordantly 1 time
    ----
    121148 pairs aligned 0 times concordantly or discordantly; of these:
      242296 mates make up the pairs; of these:
        169197 (69.83%) aligned 0 times
        26920 (11.11%) aligned exactly 1 time
        46179 (19.06%) aligned >1 times
99.84% overall alignment rate

  

Samtools flagstat

190505688 + 0 in total (QC-passed reads + QC-failed reads)
85081508 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
190336491 + 0 mapped (99.91%:-nan%)
105424180 + 0 paired in sequencing
52712090 + 0 read1
52712090 + 0 read2
105164530 + 0 properly paired (99.75%:-nan%)
105210878 + 0 with itself and mate mapped
44105 + 0 singletons (0.04%:-nan%)
8712 + 0 with mate mapped to a different chr
2912 + 0 with mate mapped to a different chr (mapQ>=5)

  

Note that the flagstat command counts alignments, not reads. please 
use the read counts table to get accurate counts of reads at each
stage of the pipeline.

Filtering statistics

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

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

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

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.

Peak statistics

Naive overlap peak file statistics

IDR peak file statistics

For a good ATAC-seq experiment in human, you expect to get 100k-200k peaks
for a specific cell type.

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 in human (hg19) is above 6. A strong TSS enrichment is
above 10. For other references please see https://www.encodeproject.org/atac-seq/
  

Annotated genomic region enrichments

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.

Replicate 3

Sample Information

Summary

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

49027823 reads; of these:
  49027823 (100.00%) were paired; of these:
    48693 (0.10%) aligned concordantly 0 times
    38514404 (78.56%) aligned concordantly exactly 1 time
    10464726 (21.34%) aligned concordantly >1 times
    ----
    48693 pairs aligned concordantly 0 times; of these:
      4789 (9.84%) aligned discordantly 1 time
    ----
    43904 pairs aligned 0 times concordantly or discordantly; of these:
      87808 mates make up the pairs; of these:
        64814 (73.81%) aligned 0 times
        4520 (5.15%) aligned exactly 1 time
        18474 (21.04%) aligned >1 times
99.93% overall alignment rate

  

Samtools flagstat

163936258 + 0 in total (QC-passed reads + QC-failed reads)
65880612 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
163871444 + 0 mapped (99.96%:-nan%)
98055646 + 0 paired in sequencing
49027823 + 0 read1
49027823 + 0 read2
97958260 + 0 properly paired (99.90%:-nan%)
97984506 + 0 with itself and mate mapped
6326 + 0 singletons (0.01%:-nan%)
6528 + 0 with mate mapped to a different chr
2773 + 0 with mate mapped to a different chr (mapQ>=5)

  

Note that the flagstat command counts alignments, not reads. please 
use the read counts table to get accurate counts of reads at each
stage of the pipeline.

Filtering statistics

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

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

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

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.

Peak statistics

Naive overlap peak file statistics

IDR peak file statistics

For a good ATAC-seq experiment in human, you expect to get 100k-200k peaks
for a specific cell type.

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 in human (hg19) is above 6. A strong TSS enrichment is
above 10. For other references please see https://www.encodeproject.org/atac-seq/
  

Annotated genomic region enrichments

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.

Replicate 4

Sample Information

Summary

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

38371500 reads; of these:
  38371500 (100.00%) were paired; of these:
    35866 (0.09%) aligned concordantly 0 times
    27760554 (72.35%) aligned concordantly exactly 1 time
    10575080 (27.56%) aligned concordantly >1 times
    ----
    35866 pairs aligned concordantly 0 times; of these:
      5296 (14.77%) aligned discordantly 1 time
    ----
    30570 pairs aligned 0 times concordantly or discordantly; of these:
      61140 mates make up the pairs; of these:
        34088 (55.75%) aligned 0 times
        4660 (7.62%) aligned exactly 1 time
        22392 (36.62%) aligned >1 times
99.96% overall alignment rate

  

Samtools flagstat

146180152 + 0 in total (QC-passed reads + QC-failed reads)
69437152 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
146146064 + 0 mapped (99.98%:-nan%)
76743000 + 0 paired in sequencing
38371500 + 0 read1
38371500 + 0 read2
76671268 + 0 properly paired (99.91%:-nan%)
76701606 + 0 with itself and mate mapped
7306 + 0 singletons (0.01%:-nan%)
5342 + 0 with mate mapped to a different chr
1606 + 0 with mate mapped to a different chr (mapQ>=5)

  

Note that the flagstat command counts alignments, not reads. please 
use the read counts table to get accurate counts of reads at each
stage of the pipeline.

Filtering statistics

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

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

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

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.

Peak statistics

Naive overlap peak file statistics

IDR peak file statistics

For a good ATAC-seq experiment in human, you expect to get 100k-200k peaks
for a specific cell type.

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 in human (hg19) is above 6. A strong TSS enrichment is
above 10. For other references please see https://www.encodeproject.org/atac-seq/
  

Annotated genomic region enrichments

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.