% \bibitem{Sos2016}Sos BC, Fung HL, Gao DR, et al. (2016) Characterization of chromatin accessibility with a transposome hypersensitive sites sequencing (THS-seq) assay. \textit{Genome Biol} \textbf{17}:20. % doi: 10.1186/s13059-016-0882-7.
% \bibitem{Ou2018}Ou J, Liu H, Yu J, et al. (2018) ATACseqQC: a Bioconductor package for post-alignment quality assessment of ATAC-seq data.\textit{BMC Genomics} \textbf{19}(1):169. % doi: 10.1186/s12864-018-4559-3.
% \bibitem{Raney2014}Raney BJ, Dreszer TR, Barber GP, et al. (2014) Track data hubs enable visualization of user-defined genome-wide annotations on the UCSC Genome Browser. \textit{Bioinformatics} \textbf{30}, 1003--1005.
% \bibitem{Koehler2011}Koehler R, Issac H, Cloonan N, Grimmond SM. (2011) The uniqueome: a mappability resource for short-tag sequencing. \textit{Bioinformatics} \textbf{27}, 272--274.
% \bibitem{LeeSchatz2012}Lee H, Schatz MC. (2012) Genomic dark matter: the reliability of short read mapping illustrated by the genome mappability score. \textit{Bioinformatics} \textbf{28}, 2097--2105.
% \bibitem{Derrien2012}Derrien T, Estell\'e J, Marco Sola S, et al. (2012) Fast computation and applications of genome mappability. \textit{PLoS One} \textbf{7}, e30377.

\bibitem{Wu1980}Wu C. (1980) The 5$^{\prime}$ ends of \textit{Drosophila} heat shock genes in chromatin are hypersensitive to DNase I. \textit{Nature} \textbf{286}(5776):854--860.

\bibitem{Keene1981}Keene MA, Corces V, Lowenhaupt K, et al. (1981) DNase I hypersensitive sites in Drosophila chromatin occur at the 5$^{\prime}$ ends of regions of transcription. \textit{Proc Natl Acad Sci U S A} \textbf{78}, 143--146.

\bibitem{McGhe1981}McGhee JD, Wood WI, Dolan M, et al. (1981) A 200 base pair region at the 5$^{\prime}$ end of the chicken adult $\beta$-globin gene is accessible to nuclease digestion. \textit{Cell} \textbf{27}, 45--55.

\bibitem{Dorschner2004}Dorschner MO, Hawrylycz M, Humbert R, et al. (2004) High-throughput localization of functional elements by quantitative chromatin profiling. \textit{Nat Methods} \textbf{1}, 219--225.

\bibitem{Sabo2004}Sabo PJ, Humbert R, Hawrylycz M, et al. (2004) Genome-wide identification of DNaseI hypersensitive sites using active chromatin sequence libraries. \textit{Proc Natl Acad Sci U S A} \textbf{101}, 4537-4542.

\bibitem{Sabo2006}Sabo PJ, Kuehn MS, Thurman R, et al. (2006) Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays. \textit{Nat Methods} \textbf{3}, 511--518.

\bibitem{Crawford2006}Crawford GE, Holt IE, Whittle J, et al. (2006) Genome-wide mapping of DNase hypersensitive sites using massively parallel signature sequencing (MPSS). \textit{Genome Res} \textbf{16}, 123--131. 

\bibitem{Boyle2008}Boyle AP, Davis S, Shulha HP, et al. (2008) High-resolution mapping and characterization of open chromatin across the genome. \textit{Cell} \textbf{132}(2):311--322. % doi: 10.1016/j.cell.2007.12.014.

\bibitem{Thurman2012}Thurman RE, Rynes E, Humbert R, et al. (2012) The accessible chromatin landscape of the human genome. \textit{Nature} \textbf{489}(7414):75-82.

\bibitem{Kelly2012}Kelly TK, Liu Y, Lay FD, et al. (2012) Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules. \textit{Genome Res} \textbf{22}(12):2497--2506. % doi: 10.1101/gr.143008.112. % NOME-seq

\bibitem{Krebs2017}Krebs AR, Imanci D, Hoerner L, Gaidatzis D, et al. (2017) Genome-wide Single-Molecule Footprinting Reveals High RNA Polymerase II Turnover at Paused Promoters. \textit{Mol Cell} \textbf{67}(3):411--422.e4. % doi: 10.1016/j.molcel.2017.06.027. 

\bibitem{Shipony2018}Shipony Z, Marinov GK, Swaffer MP, et al. (2018) Long-range single-molecule mapping of chromatin accessibility in eukaryotes. \textit{bioRxiv} 504662.

\bibitem{Daley2013}Daley T, Smith AD (2013) Predicting the molecular complexity of sequencing libraries. \textit{Nat Methods} \textbf{10}(4):325--327. % doi: 10.1038/nmeth.2375. 

\bibitem{Wang2019}Wang Y, Wang A, Liu Z, et al. (2019) Single-molecule long-read sequencing reveals the chromatin basis of gene expression. \textit{Genome Res} \textbf{29}(8):1329--1342. % doi: 10.1101/gr.251116.119.  % MeSMLR-seq

\bibitem{Aughey2018}Aughey GN, Estacio Gomez A, Thomson J, et al. (2018) CATaDa reveals global remodelling of chromatin accessibility during stem cell differentiation in vivo. \textit{Elife} \textbf{7}. pii: e32341. % doi: 10.7554/eLife.32341.

\bibitem{Chereji2019}Chereji RV, Eriksson PR, Ocampo J, Clark DJ. (2019) DNA accessibility is not the primary determinant of chromatin-mediated gene regulation \textit{bioRxiv} 639971 % doi: https://doi.org/10.1101/639971

\bibitem{Ponnaluri2017}Ponnaluri VKC, Zhang G, Est\'eve PO, et al. (2017) NicE-seq: high resolution open chromatin profiling. \textit{Genome Biol} \textbf{18}(1):122. % doi: 10.1186/s13059-017-1247-6.

\bibitem{Umeyama2017}Umeyama T, Ito T. (2017) DMS-Seq for In Vivo Genome-wide Mapping of Protein-DNA Interactions and Nucleosome Centers. \textit{Cell Rep} \textbf{21}(1):289--300. % doi: 10.1016/j.celrep.2017.09.035.

\bibitem{Timms2019}Timms RT, Tchasovnikarova IA, Lehner PJ. (2019) Differential viral accessibility (DIVA) identifies alterations in chromatin architecture through large-scale mapping of lentiviral integration sites. \textit{Nat Protoc} \textbf{14}(1):153--170. % doi: 10.1038/s41596-018-0087-5.

\bibitem{Buenrostro2013}Buenrostro JD, Giresi PG, Zaba LC, et al. (2013) Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. \textit{Nat Methods} \textbf{10}, 1213--1218. % doi: 10.1038/nmeth.2688. 

\bibitem{Buenrostro2015}Buenrostro JD, Wu B, Litzenburger UM, et al. (2015) Single-cell chromatin accessibility reveals principles of regulatory variation. \textit{Nature} \textbf{523}(7561):486--490. % doi: 10.1038/nature14590. Epub 2015 Jun 17.

\bibitem{Cusanovich2015}Cusanovich DA, Daza R, Adey A, et al. (2015) Multiplex single cell profiling of chromatin accessibility by combinatorial cellular indexing. \textit{Science} \textbf{348}(6237):910--914. % doi: 10.1126/science.aab1601. Epub 2015 May 7.

\bibitem{ENCODE2012}ENCODE Project Consortium. (2012) An integrated encyclopedia of DNA elements in the human genome. \textit{Nature} \textbf{489}, 57--74.

\bibitem{Amemiya2019}Amemiya HM, Kundaje A, Boyle AP. (2019) The ENCODE Blacklist: Identification of Problematic Regions of the Genome. \textit{Sci Rep} \textbf{9}(1):9354. % doi: 10.1038/s41598-019-45839-z.

\bibitem{Langmead2009}Langmead B, Trapnell C, Pop M, et al. (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. \textit{Genome Biol} \textbf{10}, R25.

\bibitem{Langmead2012}Langmead B, Salzberg SL. (2012) Fast gapped-read alignment with Bowtie 2. \textit{Nat Methods} \textbf{9}, 357--359.

\bibitem{Li2009a}Li H, Handsaker B, Wysoker A, et al. (2009) The Sequence Alignment/Map format and SAMtools. \textit{Bioinformatics} \textbf{25}, 2078--2079.

\bibitem{Feng2012}Feng J, Liu T, Qin B, et al. (2012) Identifying ChIP-seq enrichment using MACS. \textit{Nat Protoc} \textbf{7}, 1728--1740. 

\bibitem{LiIDR2011}Li Q, Brown J, Huang H, et al. (2011) Measuring reproducibility of high-throughput experiments. \textit{Ann Appl Stat} \textbf{5}, 1752--1779.

\bibitem{Kuhn2013}Kuhn RM, Haussler D, Kent WJ (2013) The UCSC genome browser and associated tools. \textit{Brief Bioinform} \textbf{14}, 144--161.

\bibitem{Kent2010}Kent WJ, Zweig AS, Barber G, et al. (2010) BigWig and BigBed: enabling browsing of large distributed datasets. \textit{Bioinformatics} \textbf{26}, 2204--2207.

\bibitem{Love2014}Love MI, Huber W, Anders S. (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. \textit{Genome Biol} \textbf{15}(12):550.

\bibitem{Schep2017}Schep AN, Wu B, Buenrostro JD, Greenleaf WJ. (2017) chromVAR: inferring transcription-factor-associated accessibility from single-cell epigenomic data. \textit{Nat Methods} \textbf{14}, 975--978. % doi: 10.1038/nmeth.4401.

\bibitem{Schep2015}Schep AN, Buenrostro JD, Denny SK, et al. (2015) Structured nucleosome fingerprints enable high-resolution mapping of chromatin architecture within regulatory regions. \textit{Genome Res} \textbf{25}, 1757--1770. % doi: 10.1101/gr.192294.115.

\bibitem{Ramirez2016}Ram\'irez F, Ryan DP, Gr\"uning B, et al (2016) deepTools2: a next generation web server for deep-sequencing data analysis. \textit{Nucleic Acids Res} \textbf{44}(W1):W160--165. % doi: 10.1093/nar/gkw257. Epub 2016 Apr 13.

\bibitem{Quinlan2010}Quinlan AR, Hall IM (2010) BEDTools: a flexible suite of utilities for comparing genomic features. \textit{Bioinformatics} \textbf{26}(6):841--842. % doi: 10.1093/bioinformatics/btq033. Epub 2010 Jan 28.

\bibitem{Liao2014}Liao Y, Smyth GK, Shi W. (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. \textit{Bioinformatics} \textbf{30}(7):923--930. % doi: 10.1093/bioinformatics/btt656. 

\bibitem{Bolger2014}Bolger AM, Lohse M, Usadel B. (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. \textit{Bioinformatics} \textbf{30}(15):2114--2120. % doi: 10.1093/bioinformatics/btu170. 

\bibitem{Martin2011}Martin M (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. \textit{EMBnet.journal} \textbf{17}(1):10--12.

\bibitem{Li2009b}Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. \textit{Bioinformatics} \textbf{25}, 1754--1760.

\bibitem{Corces2017}Corces MR, Trevino AE, Hamilton EG, et al. (2017) An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. \textit{Nat Methods} \textbf{14}, 959--962. % doi: 10.1038/nmeth.4396. 

\bibitem{HazkaniCovo2010} Hazkani-Covo E, Zeller RM, Martin W. (2010) Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes. \textit{PLoS Genet} \textbf{6}(2):e1000834.

\bibitem{Marinov2014}Marinov GK, Wang YE, Chan D, Wold BJ. (2014) Evidence for site-specific occupancy of the mitochondrial genome by nuclear transcription factors. \textit{PLoS One} \textbf{9}(1):e84713. % doi: 10.1371/journal.pone.0084713. 

\bibitem{Smith2015}Smith DR, Keeling PJ. (2015) Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes. \textit{Proc Natl Acad Sci U S A} \textbf{112}(33):10177--10184. % doi: 10.1073/pnas.1422049112. 

\bibitem{Landt2012}Landt SG, Marinov GK, Kundaje A, et al. (2012) ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. \textit{Genome Res} \textbf{22}(9):1813--1831. % doi: 10.1101/gr.136184.11

\bibitem{Marinov2014b}Marinov GK, Kundaje A, Park PJ, Wold BJ. (2014) Large-scale quality analysis of published ChIP-seq data. \textit{G3 (Bethesda)} \textbf{4}(2):209--223. % doi: 10.1534/g3.113.008680.

\bibitem{Tarbell2019}Tarbell ED, Liu T. (2019) HMMRATAC: a Hidden Markov ModeleR for ATAC-seq. \textit{Nucleic Acids Res} pii: gkz533. % doi: 10.1093/nar/gkz533. [Epub ahead of print]

\bibitem{McCarthy2012}McCarthy DJ, Chen Y, Smyth GK. (2012) Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. \textit{Nucleic Acids Res} \textbf{40}(10):4288--4297. % doi: 10.1093/nar/gks042. 

\bibitem{Ritchie2015}Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. 2015. limma powers differential expression analyses for RNA-sequencing and microarray studies. \textit{Nucleic Acids Res} \textbf{43}(7):e47. % doi: 10.1093/nar/gkv007. 

\bibitem{tSNE}van der Maaten LJP, Hinton GE. (2008) Visualizing High-Dimensional Data Using t-SNE. \textit{J Mach Learn Res} \textbf{9}:2579--2605.

\bibitem{UMAP}Becht E, McInnes L, Healy J, et al. (2018) Dimensionality reduction for visualizing single-cell data using UMAP. \textit{Nat Biotechnol} \textbf{37}:38--44.

\bibitem{Li2019}Li Z, Schulz MH, Look T, et al. (2019) Identification of transcription factor binding sites using ATAC-seq. \textit{Genome Biol} \textbf{20}(1):45. % doi: 10.1186/s13059-019-1642-2.

\bibitem{Hesselberth2009}Hesselberth JR, Chen X, Zhang Z, et al. (2009) Global mapping of protein-DNA interactions in vivo by digital genomic footprinting. \textit{Nat Methods} \textbf{6}(4):283--289. % doi: 10.1038/nmeth.1313. 

\bibitem{Neph2012a}Neph S, Stergachis AB, Reynolds A, et al. (2012) Circuitry and dynamics of human transcription factor regulatory networks. \textit{Cell} \textbf{150}, 1274--1286. % doi: 10.1016/j.cell.2012.04.040. 

\bibitem{Neph2012b}Neph S, Vierstra J, Stergachis AB, et al. (2012) An expansive human regulatory lexicon encoded in transcription factor footprints. \textit{Nature} \textbf{489}, 83--90. % doi: 10.1038/nature11212.

\bibitem{Stergachis2013}Stergachis AB, Neph S, Reynolds A, et al. (2013) Developmental fate and cellular maturity encoded in human regulatory DNA landscapes. \textit{Cell} \textbf{154}, 888--903. % doi: 10.1016/j.cell.2013.07.020.

\bibitem{Pique-Regi2011}Pique-Regi R, Degner JF, Pai AA, et al. (2011) Accurate inference of transcription factor binding from DNA sequence and chromatin accessibility data. \textit{Genome Res} \textbf{21}(3):447--455. % doi: 10.1101/gr.112623.110. 

\bibitem{Cuellar-Partida2012}Cuellar-Partida G, Buske FA, McLeay RC, et al. (2012) Epigenetic priors for identifying active transcription factor binding sites. \textit{Bioinformatics} \textbf{28}(1):56--62. % doi: 10.1093/bioinformatics/btr614. 

\bibitem{Piper2013}Piper J, Elze MC, Cauchy P, et al. (2013) Wellington: a novel method for the accurate identification of digital genomic footprints from DNase-seq data. \textit{Nucleic Acids Res} \textbf{41}(21):e201. % doi: 10.1093/nar/gkt850. 

\bibitem{Sherwood2014}Sherwood RI, Hashimoto T, O'Donnell CW, et al. (2014) Discovery of directional and nondirectional pioneer transcription factors by modeling DNase profile magnitude and shape. \textit{Nat Biotechnol} \textbf{32}(2):171--178. % doi: 10.1038/nbt.2798. 

\bibitem{He2014}He HH, Meyer CA, Hu SS, et al. (2014) Refined DNase-seq protocol and data analysis reveals intrinsic bias in transcription factor footprint identification. \textit{Nat Methods} \textbf{11}, 73--78. % doi: 10.1038/nmeth.2762. 

\bibitem{Sung2014}Sung MH, Guertin MJ, Baek S, Hager GL. 2014. DNase footprint signatures are dictated by factor dynamics and DNA sequence. \textit{Mol Cell} \textbf{56}(2):275--285. % doi: 10.1016/j.molcel.2014.08.016. 

\bibitem{Gusmao2014}Gusmao EG, Dieterich C, Zenke M, Costa IG. (2014) Detection of active transcription factor binding sites with the combination of DNase hypersensitivity and histone modifications. \textit{Bioinformatics} \textbf{30}(22):3143--3151. % doi: 10.1093/bioinformatics/btu519. 

\bibitem{Raj2015}Raj A, Shim H, Gilad Y, et al. (2015). msCentipede: Modeling Heterogeneity across Genomic Sites and Replicates Improves Accuracy in the Inference of Transcription Factor Binding. \textit{PLoS One} \textbf{10}(9):e0138030. % doi: 10.1371/journal.pone.0138030. eCollection 2015.

\bibitem{Yardimci2015}Yardimci GG, Frank CL, Crawford GE, Ohler U. (2015) Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection. \textit{Nucleic Acids Res} \textbf{42}(19):11865--11878. % doi: 10.1093/nar/gku810. 

\bibitem{Gusmao2016}Gusmao EG, Allhoff M, Zenke M, Costa IG. (2016) Analysis of computational footprinting methods for DNase sequencing experiments. \textit{Nat Methods} \textbf{13}(4):303--309. % doi: 10.1038/nmeth.3772. 

\bibitem{Quach20117}Quach B, Furey TS. (2017) DeFCoM: analysis and modeling of transcription factor binding sites using a motif-centric genomic footprinter. \textit{Bioinformatics} \textbf{33}(7):956--963. % doi: 10.1093/bioinformatics/btw740.

\bibitem{Baek2017}Baek S, Goldstein I, Hager GL. (2017) Bivariate Genomic Footprinting Detects Changes in Transcription Factor Activity.  \textit{Cell Rep} \textbf{19}(8):1710--1722. % doi: 10.1016/j.celrep.2017.05.003.

\bibitem{Karabacak2019}Karabacak Calviello A, Hirsekorn A, Wurmus R, et al. (2019) Reproducible inference of transcription factor footprints in ATAC-seq and DNase-seq datasets using protocol-specific bias modeling. \textit{Genome Biol} \textbf{20}(1):42. % doi: 10.1186/s13059-019-1654-y.

\bibitem{Sung2016}Sung MH, Baek S, Hager GL. (2016) Genome-wide footprinting: ready for prime time? \textit{Nat Methods} \textbf{13}(3):222--228. % doi: 10.1038/nmeth.3766.

\bibitem{Vierstra2016}Vierstra J, Stamatoyannopoulos JA. (2016) Genomic footprinting. \textit{Nat Methods} \textbf{13}(3):213--221. %  doi: 10.1038/nmeth.3768. 

\bibitem{Mortazavi2006}Mortazavi A, Leeper Thompson EC, Garcia ST, et al. (2006) Comparative genomics modeling of the NRSF/REST repressor network: from single conserved sites to genome-wide repertoire. \textit{Genome Res} \textbf{16}(10):1208--1221.

\bibitem{Mortazavi2007}Johnson DS, Mortazavi A, Myers RM, Wold B. (2007). Genome-wide mapping of in vivo protein-DNA interactions. \textit{Science} \textbf{316}(5830):1497--1502.

\bibitem{Weirauch2014}Weirauch MT, Yang A, Albu M, et al. (2014) Determination and inference of eukaryotic transcription factor sequence specificity. \textit{Cell} \textbf{158}, 1431--1443. % doi: 10.1016/j.cell.2014.08.009. 

\bibitem{Grant2011}Grant CE, Bailey TL, Noble WS. (2011) FIMO: scanning for occurrences of a given motif. \textit{Bioinformatics} \textbf{27}, 1017--1018. % doi: 10.1093/bioinformatics/btr064.

\bibitem{Henikoff2011}Henikoff JG, Belsky JA, Krassovsky K, et al. (2011) Epigenome characterization at single base-pair resolution. \textit{Proc Natl Acad Sci U S A} \textbf{108}, 18318--18323. % doi: 10.1073/pnas.1110731108. 

\bibitem{Fu2008}Fu Y, Sinha M, Peterson CL, Weng Z. (2008) The insulator binding protein CTCF positions 20 nucleosomes around its binding sites across the human genome. \textit{PLoS Genet} \textbf{4}, e1000138. % doi: 10.1371/journal.pgen.1000138.