Sites of protein binding to DNA are inferred from footprints or spans of protection against a probing reagent. In most protocols, sites of accessibility to a probe are detected by mapping breaks in DNA strands. As discussed in this unit, such methods obscure molecular heterogeneity by averaging cuts at a given site over all DNA strands in a sample population. The DNA methyltransferase accessibility protocol for individual templates (MAPit), an alternative method described in this unit, localizes protein-DNA interactions by probing with cytosine-modifying DNA methyltransferases followed by bisulfite sequencing. Sequencing individual DNA products after amplification of bisulfite-converted sequences permits assignment of the methylation status of every enzyme target site along a single DNA strand. Use of the GC-methylating enzyme M.CviPI allows simultaneous mapping of chromatin accessibility and endogenous CpG methylation. MAPit is therefore the only footprinting method that can detect subpopulations of molecules with distinct patterns of protein binding or chromatin architecture and correlate them directly with the occurrence of endogenous methylation. Additional advantages of MAPit methylation footprinting as well as considerations for experimental design and potential sources of error are discussed.