One of the pervasive features of cancer is the diversity of mutations found in malignant cells within the same tumor; a phenomenon called clonal diversity or intratumor heterogeneity. Clonal diversity allows tumors to adapt to the selective pressure of treatment and likely contributes to the development of treatment resistance and cancer recurrence. Thus, the ability to precisely delineate the clonal substructure of a tumor, including the evolutionary history of its development and the co-occurrence of its mutations, is necessary to understand and overcome treatment resistance. However, DNA sequencing of bulk tumor samples cannot accurately resolve complex clonal architectures. Here, we performed high-throughput single-cell DNA sequencing to quantitatively assess the clonal architecture of acute myeloid leukemia (AML). We sequenced a total of 556,951 cells from 77 patients with AML for 19 genes known to be recurrently mutated in AML. The data revealed clonal relationship among AML driver mutations and identified mutations that often co-occurred (e.g., NPM1/FLT3-ITD, DNMT3A/NPM1, SRSF2/IDH2, and WT1/FLT3-ITD) and those that were mutually exclusive (e.g., NRAS/KRAS, FLT3-D835/ITD, and IDH1/IDH2) at single-cell resolution. Reconstruction of the tumor phylogeny uncovered history of tumor development that is characterized by linear and branching clonal evolution patterns with latter involving functional convergence of separately evolved clones. Analysis of longitudinal samples revealed remodeling of clonal architecture in response to therapeutic pressure that is driven by clonal selection. Furthermore, in this AML cohort, higher clonal diversity (≥4 subclones) was associated with significantly worse overall survival. These data portray clonal relationship, architecture, and evolution of AML driver genes with unprecedented resolution, and illuminate the role of clonal diversity in therapeutic resistance, relapse and clinical outcome in AML.