In this work, we revise the conventional description of J/Psi(1S), Y(1S), Psi(2S) and Y(2S) elastic photo- and electroproduction off a nucleon target within the color dipole picture and carefully study various sources of theoretical uncertainties in calculations of the corresponding electroproduction cross sections. For this purpose, we test the corresponding predictions using a bulk of available dipole cross section parametrisations obtained from deep inelastic scattering data at HERA. Specifically, we provide the detailed analysis of the energy and hard-scale dependencies of quarkonia yields employing the comprehensive treatment of the quarkonia wave functions in the Schroedinger equation based approach for a set of available c-bar{c} and b-bar{b} interquark interaction potentials. Besides, we quantify the effect of Melosh spin rotation, the Q^2-dependence of the diffractive slope and an uncertainty due to charm and bottom quark mass variations.