Rapid and accurate diagnosis of acute myeloid leukemia (AML) remains a significant challenge, particularly in the context of myelodysplastic syndrome (MDS) or MDS/myeloproliferative neoplasm with NPM1 mutations. This study introduces an innovative approach using holotomography (HT), a 3D label-free quantitative phase imaging technique, to detect NPM1 mutations. We analyzed a dataset of 2073 HT myeloblast images from 48 individuals, including both NPM1 wild-type and mutated samples, to distinguish subcellular morphological changes associated with NPM1 mutations. Employing a convolutional neural network, we analyzed 3D cell morphology, focusing on refractive index distributions. The machine learning model showed high accuracy, with an area under the receiver operating characteristic curve of 0.9375 and a validation accuracy of 76.0%. Our findings reveal distinct morphological differences between the NPM1 wild-type and mutation at the subcellular level. This study demonstrates the potential of HT combined with deep learning for early, efficient, and cost-effective diagnosis of AML, offering a promising alternative to traditional stepwise genetic testing methods and providing additional assistance in morphological myeloblast discrimination. This approach may revolutionize the diagnostic process in leukemia, facilitating early detection and potentially reducing the reliance on extensive genetic testing.