A classic approach to computational fluid dynamics is to perform simulations with a fixed set of variables in order to account for parameters and boundary conditions. However, experiments and real-life performance are subject to variability in their conditions. In recent years, the interest of performing simulations under uncertainty is increasing, but this is not yet a common rule, and simulations with lack of information are still taking place. This procedure could be missing details such as whether sources of uncertainty affect dramatic parts in the simulation of the flow. One of the reasons of avoiding to quantify uncertainties is that they usually require to run an unaffordable number of CFD simulations to develop the study. To face this problem, Non-Intrusive Uncertainty Quantification (UQ) has been applied to 3D Reynolds-Averaged Navier-Stokes simulations of an under-expanded jet from an aircraft exhaust with the Spalart-Allmaras turbulent model, in order to assess the impact of inaccuracies and quality in the simulation. To save a large number of computations, sparse grids are used to compute the integrals and built surrogates for UQ. Results show that some regions of the jet plume can be more sensitive than others to variance in both physical and turbulence model parameters. The Spalart-Allmaras turbulent model is demonstrated to have an accurate performance with respect to other turbulent models in RANS, LES and experimental data, and the contribution of a large variance in its parameter is analysed. This investigation explicitly outlines, exhibits and proves the details of the relationship between diverse sources of input uncertainty, the sensitivity of different quantities of interest to said uncertainties and the spatial distribution arising due to their propagation in the simulation of the high-speed jet flow. This analysis represents first numerical study that provides evidence for this heuristic observation.