Combining Genetic Algorithms, Gibbs Sampler, and Gradient-based Inversion to Estimate Uncertainty in 2D FWI