The dynamics of charge carriers in doped graphene, i.e., graphene with a gap in the energy spectrum depending on the substrate, in the presence of a Coulomb impurity with charge Z is considered within the effective two-dimensional Dirac equation. The wave functions of carriers with conserved angular momentum J = M + 1/2 are determined for a Coulomb potential modified at small distances. This case, just as any two-dimensional physical system, admits both integer and half-integer quantization of the orbital angular momentum in plane, M = 0, ±1, ±2, …. For J = 0, ±1/2, ±1, critical values of the effective charge Zcr(J, n) are calculated for which a level with angular momentum J and radial quantum numbers n = 0 and n = 1 reaches the upper boundary of the valence band. For Z Zcr (J, n = 0), scattering phases are calculated as a function of hole energy for several values of supercriticality, as well as the positions e0 and widths γ of quasistationary states as a function of supercriticality. The values of e0* and width γ* are pointed out for which quasidiscrete levels may show up as Breit–Wigner resonances in the scattering of holes by a supercritical impurity. Since the phases are real, the partial scattering matrix is unitary, so that the radial Dirac equation is consistent even for Z > Zcr. In this single-particle approximation, there is no spontaneous creation of electron–hole pairs, and the impurity charge cannot be screened by this mechanism.