Online scalable Gaussian processes with conformal prediction for guaranteed coverage