Temozolomide kills cancer cells by formingO6-methylguanine (O6-MeG), which leads to apoptosis by mismatch-repair overload, as an important basis for glioblastoma therapy. However,O6-MeG repair byO6-methylguanine-DNA methyltransferase (MGMT) strongly contributes to drug resistance. Genomic profiles ofO6-MeG may be useful for understanding howO6-MeG accumulation can be influenced by repair mechanisms, but there are no methods to map genomic locations ofO6-MeG. Here, we devised an immunoprecipitation– and polymerase-stalling-based strategy to locateO6-MeG across the whole genome at single-nucleotide resolution. We applied this strategy, termedO6-MeG-seq, to analyzeO6-MeG formation and repair with regards to sequence contexts and functional genomic regions in glioblastoma-derived cell lines deficient in MGMT and evaluated the impact of MGMT transfection on these profiles.O6-MeG signatures from the sequencing data were highly similar to mutational signatures extracted from cancer genomes of patients previously treated with temozolomide, consistent withO6-MeG initiating mutations. Furthermore, it appears that MGMT does not preferentially repairO6-MeG with respect to sequence context, chromatin state or gene expression level, however, may protect oncogenes from mutations. Finally, an MGMT-independent strand bias inO6-MeG accumulation in highly expressed genes suggests alternative repair influences. These data provide high resolution insight on howO6-MeG formation and repair is impacted by genome structure and regulation.O6-MeG-seq also provides a strategy for future studies of DNA modification signatures as diagnostic markers for addressing drug resistance and preventing secondary cancers.