Abstract We assess whether tabular foundation models can be used as off-the-shelf probabilistic photometric-redshift estimators for quasars in the 12-band S-PLUS DR6 survey, where color–redshift degeneracies produce multimodal posteriors, and spectroscopic training sets are shifted relative to the photometric population. TabPFN 2.5, RealTabPFN 2.5, and TabICL are benchmarked against eight task-specific baselines, including linear conditional Gaussians, FlexZBoost, mixture-density networks, normalizing flows, random forests, and gradient-boosted trees, with training sets from 500 to 121,626 quasars, using both density and point-prediction metrics, together with importance-weighted scores that approximate deployment on the photometric target sample. TabPFN 2.5 is best or statistically tied for best on all metrics except the unweighted CDE loss, on which the normalizing flow is statistically tied and attains the lowest mean value; its largest gains occur for small training sets and in difficult regimes (very bright and faint sources, high redshift), while retaining near-nominal calibration under a covariate shift. Its main practical cost is inference: with frozen weights, large support and target catalogs require substantial GPU/accelerator memory, and full-catalog deployment may need support-set subsampling or distillation. SHAP attributions identify the Wide-field Infrared Survey Explorer W1/W2 as the strongest individual predictors, with UV and optical bands offering nonnegligible refinements. We conclude that TabPFN 2.5 is a strong default for probabilistic quasar photo- z estimation, particularly when training data are limited or when calibration under a covariate shift is critical.