Gaussian process modelling in approximate Bayesian computation to estimate horizontal gene transfer in bacteria