A novel variational Bayes multiple locus Z-statistic for genome-wide association studies with Bayesian model averaging