Explicitly correlated configuration interaction wavefunctions using spherical Gaussians. Formulation and initial application to LiH