We present numerical results on the static quark--anti-quark grand canonical potential in full QCD at non-vanishing temperature ($T$) and quark chemical potential ($mu$). Non-zero $mu$-s are reached by means of multi-parameter reweighting. The dynamical staggered simulations were carried out for $n_f=2+1$ flavors with physical quark masses on $4times 12^3$ lattices.