Parallel kinetic Monte Carlo simulation of Coulomb glasses