Prediction interval of generalized linear models

When doing predictions using a generalized linear model, how to get the prediction interval or prediction distribution?

My model has a gamma distribution and a reciprocal link function in particular.

Can anyone please help me on solving this by either analytical methods or simulation?