Modified gravity (MG) theories aim to reproduce the observed acceleration of the Universe by reducing the dark sector while simultaneously recovering General Relativity (GR) within dense environments. Void studies appear to be a suitable scenario to search for imprints of alternative gravity models on cosmological scales. Voids cover an interesting range of density scales where screening mechanisms fade out, which reaches from a density contrast δ ≈ −1 close to their centers to δ ≈ 0 close to their boundaries. We present an analysis of the level of distinction between GR and two modified gravity theories, the Hu-Sawicki f (R) and the symmetron theory. This study relies on the abundance, linear bias, and density profile of voids detected in n-body cosmological simulations. We define voids as connected regions made up of the union of spheres with a mean density given by ρ v = 0.2 ρ m , but disconnected from any other voids. We find that the height of void walls is considerably affected by the gravitational theory, such that it increases for stronger gravity modifications. Finally, we show that at the level of dark matter n-body simulations, our constraints allow us to distinguish between GR and MG models with | f R0 | > 10 −6 and z S S B > 1. Differences of best-fit values for MG parameters that are derived independently from multiple void probes may indicate an incorrect MG model. This serves as an important consistency check.