Large-scale mutagenesis of target DNA sequences allows researchers to comprehensively assess the effects of single-nucleotide changes. Here we demonstrate the construction of a systematic allelic series (SAS) using massively parallel single-nucleotide mutagenesis with reversibly terminated deoxyinosine triphosphates (rtITP). We created a mutational library containing every possible single-nucleotide mutation surrounding the active site of the TEM-1 β-lactamase gene. When combined with high-throughput functional assays, SAS mutational libraries can expedite the functional assessment of genetic variation.