Methods and compositions are provided for identifying any of the presence, location and phasing of methylated and/or hydroxymethylated cytosines in nucleic acids including long stretches of DNA. In some embodiments, the method may comprise reacting a first portion (aliquot) of a nucleic acid sample with a dioxygenase and optionally a glucosyltransferase in a reaction mixture containing the nucleic acid followed by a reaction with a cytidine deaminase to detect and optionally map 5mC in a DNA. Optionally, a second portion can be reacted with glucosyltransferase followed by reaction with a cytidine deaminase to detect and optionally map 5hmC in a DNA.