Contrast-enhanced ultrasound (CEUS) plays a pivotal role in the diagnosis of primary breast cancer and in axillary lymph node (ALN) metastasis. However, the imaging features that are clinically crucial for lymph node metastasis have not been fully elucidated. Hence, we developed a bimodal model to predict ALN metastasis in patients with early breast cancer by integrating CEUS images with the annotated imaging features. The model adopted a light-gradient boosting machine to produce feature importance, enabling the extraction of clinically crucial imaging features. In this retrospective study, the diagnostic performance of the model was investigated using 788 CEUS images of ALNs obtained from 788 patients who underwent breast surgery between 2013 and 2021, with the ground truth defined by the pathological diagnosis. The results indicated that the test cohort had an area under the receiver operating characteristic curve (AUC) value of 0.93 (95% confidence interval: 0.88, 0.98). The model had an accuracy of 0.93, which was higher than the radiologist's diagnosis (accuracy of 0.85). The most important imaging features were heterogeneous enhancement, diffuse cortical thickening, and eccentric cortical thickening. Our model has an excellent diagnostic performance, and the extracted imaging features could be crucial for confirming ALN metastasis in clinical settings.