The study of complex turbulent flows by means of large-eddy simulation approaches has become increasingly popular in many scientific and engineering applications. The underlying filtering operation of the approach enables to significantly reduce the spatial and temporal resolution requirements by means of representing only large-scale motions. However, the small-scale stresses and their effects on the resolved flow field are not negligible, and therefore require additional modeling. As a consequence, the assumptions made in the closure formulations become potential sources of model-form uncertainty that can impact the quantities of interest. The objective of this work, thus, is to perform a model-form sensitivity analysis in large-eddy simulations of an axisymmetric turbulent jet following an eigenspace-based strategy recently proposed. The approach relies on introducing perturbations to the decomposed subgrid-scale stress tensor within a range of physically plausible values. These correspond to discrepancy in magnitude (trace), anisotropy (eigenvalues) and orientation (eigenvectors) of the normalized, small-scale stresses with respect to a given tensor state, such that propagation of their effects can be assessed. The generality of the framework with respect to the six degrees of freedom of the small-scale stress tensor makes it