The effect of differences in nominally equal sectors of a bladed disk (mistuning) is a well-known problem for designers since the forced response may show localized amplification of the blade response with respect to a cyclically symmetric (tuned) configuration. In order to perform a large number of simulations in a reasonable amount of time to characterize the highest blade response, corresponding to the worst mistuning pattern, reduction techniques have been developed, where the mistuning is introduced directly in a compact, reduced order model (ROM) obtained from very large finite element (FE) models. Typically, mistuning is introduced in the ROM in terms of natural frequency perturbations of the blade; nevertheless, a better insight is specifically required in order to correlate mistuning to a specific source (geometrical, material, contact mistuning). In this paper, a reduction technique is presented to take into account mistuning due to the contact uncertainties at the blade root joint, which can be caused by design tolerances, manufacturing process, assembly procedures, wear, etc. The technique takes its basis from the Craig–Bampton component mode synthesis (CB-CMS) applied to the uncoupled blade and disk sector, which is typically included in most of the FE software for easy implementation in standard industrial practice. The full set of master degrees-of-freedom (DOFs) at the random contacts are purposely reduced using an optimal local modal basis based on the Gram–Schmidt interface (GSI) technique developed by the authors. Experimental evidence of actual uncertain contact obtained during joint preloading is used to find an appropriate base to represent typical contact patterns.