The prolonged course of the COVID-19 pandemic necessitates sustained surveillance of emerging variants. This study aimed to develop a multiplex real-time polymerase chain reaction (rt-PCR) suitable for the real-time tracking of Omicron subvariants in clinical and wastewater samples. Plasmids containing variant-specific mutations were used to develop a MeltArray assay. After a comprehensive evaluation of both analytical and clinical performance, the established assay was used to detect Omicron variants in clinical and wastewater samples, and the results were compared with those of next-generation sequencing (NGS) and droplet digital PCR (ddPCR). The MeltArray assay identified 14 variant-specific mutations, enabling the detection of five Omicron sublineages (BA.2*, BA.5.2*, BA.2.75*, BQ.1*, and XBB.1*) and eight subvariants (BF.7, BN.1, BR.2, BQ.1.1, XBB.1.5, XBB.1.16, XBB.1.9, and BA.4.6). The limit of detection (LOD) of the assay was 50 copies/reaction, and no cross-reactivity was observed with 15 other respiratory viruses. Using NGS as the reference method, the clinical evaluation of 232 swab samples exhibited a clinical sensitivity of > 95.12% (95% CI 89.77–97.75%) and a specificity of > 95.21% (95% CI, 91.15–97.46%). When used to evaluate the Omicron outbreak from late 2022 to early 2023, the MeltArray assay performed on 1408 samples revealed that the epidemic was driven by BA.5.2* (883, 62.71%) and BF.7 (525, 37.29%). Additionally, the MeltArray assay demonstrated potential for estimating variant abundance in wastewater samples. The MeltArray assay is a rapid and scalable method for identifying SARS-CoV-2 variants. Integrating this approach with NGS and ddPCR will improve variant surveillance capabilities and ensure preparedness for future variants.