The comprehension and structural modeling of masonry constructions is fundamental to safeguard the integrity of built cultural assets and intervene through adequate actions, especially in earthquake-prone regions. Despite the availability of several modeling strategies and modern computing power, modeling masonry remains a great challenge because of still demanding computational efforts, constraints in performing destructive or semi-destructive in-situ tests, and material uncertainties. This paper investigates the shear behavior of masonry walls by applying a plane-stress FE continuum model with the Modified Masonry-like Material (MMLM). Epistemic uncertainty affecting input parameters of the MMLM is considered in a probabilistic framework. After appointing a suitable probability density function to input quantities according to prior engineering knowledge, uncertainties are propagated to outputs relying on gPCE-based surrogate models to considerably speed up the forward problem-solving. The sensitivity of the response to input parameters is evaluated through the computation of Sobol' indices pointing out the parameters more worthy to be further investigated, when dealing with the seismic assessment of masonry buildings. Finally, masonry mechanical properties are calibrated in a probabilistic setting with the Bayesian approach to the inverse problem based on the available measurements obtained from the experimental load-displacement curves provided by shear compression in-situ tests.