Nuclear deformation effects are theoretically investigated in terms of deformation corrections of the electronic binding and transition energies, $g$ factor, and hyperfine splitting constant. By solving the Dirac equation twice, with the nuclear potential calculated from Fermi and deformed Fermi nuclear density distributions, we separate the deformation effect in binding energies and wavefunctions. The parameters for both models are determined from experimental data. The considered corrections are of interest for spectral analysis and are numerically calculated for the widest possible range of nuclei, consisting over 1100 different samples. The subtleties between different sources of measured data and the corresponding results are discussed. In addition, the importance of deformation effects for the search of new physics with singly-charged ions is examined.