Numerical computations of ITER equilibria in the hybrid scenario using a three-dimensional (3D) magnetohydrodynamic equilibrium code with nested magnetic flux surfaces demonstrate the formation of internal 3D helical cores similar to saturated ideal internal kink modes under reversed magnetic shear conditions when the minimum value of the inverse rotational transform is in the neighbourhood of unity.