The parasite drag calculation in VSPAERO is not suitable for this kind of thing. It is a very simple model based on NACA 0012 data. As you change the airfoil, the drag model will not change.
That said, I am confident that a genetic algorithm will be very inefficient at designing a wing.
Families of airfoils are often parameterized in terms of design lift coefficient. You should first choose a family based on your needs (supercritical, laminar flow, etc). Then, consider the cl distribution required to obtain your target cl*c distribution given your planform c distribution. (If that sentence doesn't make sense, you need to study wing design until it does). Next, choose the airfoil design cl distribution to match the target cl distribution. Finally, adjust the twist distribution to get any final adjustments to your lift distribution. There are techniques to do this that will give an answer much faster than optimization. The result will be a wing that is very difficult to beat -- for almost zero computational cost.
Rob