Prediction interval of generalized linear model

I am trying to do predictions using a generalized linear model, and my model has a gamma distribution for the response variable and a reciprocal link function. Are there ways to obtain the prediction interval or prediction distribution by simulation? I couldn't find any software with such a built-in function.