We use relativistic multiconfiguration Dirac-Hartree-Fock and configuration interaction calculations to study 5s2nl and 5s5p2 configurations of In I, Sn II, Sb III, and Te IV. Energies, transition amplitudes, Land´e gJ-factors, and hyperfine constants are calculated using a correlation model that accounts for valence and core-valence correlation. Also spin- and orbital polarization effects are accounted for by single excitations from all core-shells to an increasing set of active orbitals. Transformed to the LSJ-coupling scheme, using the new features of the GRASP2K program, the calculated wave functions shed light on the difficulties in labeling some states due to the extensive 5s25d and 5s5p2 configuration interaction. Our results are compared with experimental values and values from relativistic many-body perturbation theory (RMBPT) and all-order single-double (SD) calculations. The theoretical work is complemented with experiment, and new hyperfine interaction constants are derived for several states in In I from high resolution Fourier Transform Spectra.