Intermediate-risk acute myeloid leukemia (IR-AML), which accounts for a substantial number of AML cases, is highly heterogeneous. Although several mutations have been identified, the heterogeneity of AML is uncertain because novel mutations have yet to be discovered. Here we applied next generation sequencing (NGS) platform to screen mutational hotspots in 410 genes relevant to hematological malignancy. IR-AML samples (N=95) were sequenced by Illumina Hiseq and mutations in 101 genes were identified. Only seven genes (CEBPA, NPM1, DNMT3A, FLT3-ITD, NRAS, IDH2 and WT1) were mutated in more than 10% of patients. Genetic interaction analysis identified several cooperative and exclusive patterns of overlapping mutations. Mutational analysis indicated some correlation between genotype and phenotype. FLT3-ITD mutations were identified as independent factors of poor prognosis, while CEBPA mutations were independent favorable factors. Co-occurrence of FLT3-ITD, NPM1 and DNMT3A mutations was identified with associated with specific clinical AML features and poor outcomes. Furthermore, by integrating multiple mutations in the survival analysis, 95 IR-AML patients could be stratified into three distinct risk groups allowing reductions in IR-AML by one-third. Our study offers deep insights into the molecular pathogenesis and biology of AML and indicated that the prognosis of IR-AML could be further stratified by different mutation combinations which may direct future treatment intervention.