A, B. Distribution of the max CADD38 PHRED (model for GRCh38) scores for (A) cysteine (n=1,401) and (B) lysine (n=4,363) CpDAAs of low, medium, and high intrinsic reactivities, defined by isoTOP-ABPP ratios, low (R10:1>5), medium (2<R10:1<5), high (R<2) Kruskal-Wallis nonparametric test to examine reactivity group difference, Wilcox test used for pairwise comparisons (BH-adjusted p-values, *p. adj = 0.04, **p. adj = 0.0037, ***p. adj = 0.00013). The boxplot boxes represent the lower and upper quartiles, with the central band as median. Notches show the confidence interval based on median ± 1.58*IQR/sqrt(n), and whiskers mark observations that satisfy quartiles ± 1.5*IQR. Median CADD38 max codon scores with bootstrapped 95% confidence intervals for reactive groups are: low CpD Cys 27.3 [26.9, 28.0], medium CpD Cys 28.55 [27.80, 29.05], high CpD Cys 31 [28.8, 32.0], low CpD Lys 29.5 [29.3, 29.6], medium CpD Lys 29.25 [28.85, 29.50], high CpD Lys 29.05 [28.50, 29.55]. 