Neuroblastoma is the most prevalent solid tumor in early childhood, with a 5-year overall survival rate of 40-60% in high-risk cases. Therefore, the identification of novel biomarkers for the diagnosis, prognosis, and therapy of neuroblastoma is crucial for improving the clinical outcomes of these patients. In this study, we conducted the whole-exome sequencing of 48 freshly frozen tumor samples obtained from the Biobank. Somatic variants were identified and selected using a bioinformatics analysis pipeline. The mutational signatures were determined using the Mutalisk online tool. Cancer driver genes and druggable mutations were predicted using the Cancer Genome Interpreter. The most common mutational signature was single base substitution 5. MUC4, MUC16, and FLG were identified as the most frequently mutated genes. Using the Cancer Genome Interpreter, we identified five recurrent cancer driver mutations spanning MUC16, MUC4, ALK, and CTNND1, with the latter being novel and containing a missense mutation, R439C. We also identified 11 putative actionable mutations including NF1 Q1798*, Q2616*, and S636X, ALK F1174L and R1275Q, SETD2 P10L and Q1829E, BRCA1 R612S, NOTCH1 D1670V, ATR S1372L, and FGFR1 N577K. Our findings provide a comprehensive overview of the novel information relevant to the underlying molecular pathogenesis and therapeutic targets of neuroblastoma.