Congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency (21OHD) is autosomal recessive disorder of cortisol biosynthesis. Genetic defects in CYP21A2 cause 21OHD. The aim of this study was to determine spectrum of mutations in CYP21A2 in a large cohort and analyze the genotype-phenotype correlation to assess predictive characteristics of genotype. We investigated a total of 113 patients with 21OHD. Next-generation sequencing and Multiplex ligation-dependent probe amplification of the CYP21A2 gene were performed in patients and their parents. The genotypes were categorized into Groups 0, A, B, and C according to the residual 21-hydroxylase activities. In this study, the group A was divided into two subgroups as A1 and A2. Three novel variants were found. The genotype-phenotype correlation of the mutation classification was 91.5%. Positive predictivity of subgroups A1 was higher than groups A and subgroups A2. Our study reports genotype-phenotype correlations in the largest 21OHD cohort in Turkey. This correlation sustained when we analyzed our data in combination with metadata from other published studies. This study confirms that CYP21A2 genotyping with next-generation sequencing and MLPA can accurately and reliably confirm the diagnosis of 21OHD. We propose a new classification by dividing group A into two new subgroups to better predict the phenotype. In light of this very high genotype-phenotype correlation, with their ever-increasing availability, declining cost, and turnaround time, we propose that molecular genetic studies can be more economical and practical alternative to the current initial diagnostic laboratory studies based on assays of intermediary steroid metabolites.