Effect of vacancies on the structure of solid molecular parahydrogen studied with variational Monte Carlo simulations