There's a good answer already but I wanted to mention some other possibilities.
If you're prepared to make a parametric assumption, any parametric model should let you get an approximate quantile.
For example, consider a GLM. At any combination of the predictors, you can calculate the conditional (fitted) distribution.
If you're doing out-of-sample prediction you can get an approximate predicted value -- though incorporating parameter estimation uncertainty is a bit more tricky, in cases where it's not tractable it can be done via approximation, such as treating the difference between the true parameters and the estimated parameters as if the asymptotic case applied (multivariate normal) and then simulating.
There's also bootstrapping approaches to obtaining quantiles for predictions in the cases where you don't want to make a parametric assumption. (The book by Davison & Hinkley Bootstrap methods and their application, has information on how to do bootstrapped prediction intervals. These days this would be regarded as a form of parametric bootstrap, though it's not necessarily to think of it that way)