Variational inference for Markov jump processes