  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 $\delta
\approx -1$ close to their centers to $\delta \approx 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 {\it \textup{mean}} density
given by $\overline\rho_v=0.2\,\overline\rho_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_{SSB} > 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.
