We have extended the fully variational molecular orbital (FVMO) method to the full-configuration interaction (CI) wave function (full-CI FVMO). All variational parameters in the full-CI scheme, i.e., exponents and centers in Gaussian-type function (GTF) basis set, and nuclear positions, as well as the CI coefficients, are simultaneously optimized by using their analytical gradients. We have applied the full-CI FVMO method to the ground and electronic excited states of hydrogen molecule. In the ground state, the total energy (-1.174 015 hartree) and the internuclear distance (1.4016 bohr) obtained by the full-CI FVMO calculation with [8s4p2d] GTFs agree very well with the high-level calculation by the 249 term expansion in elliptic coordinates (-1.174 476 hartree and 1.4010 bohr, respectively). The excitation energies to the (1)Sigma(u)(+), (1)Pi(u), (3)Sigma(g)(+), and (3)Pi(u) Rydberg states calculated by the full-CI FVMO method with [8s4p2d] GTFs coincide with the experimental values within 52 cm(-1). The present result can not be obtained with the conventional basis set approach because of the fact that our full-CI FVMO calculation gives an extremely accurate wave function with a relatively small number of basis functions owing to the extension of flexibility in the variational space.