def standardize_snps(G):
mean = G.mean(<NUM_LIT:0>)<EOL>std = G.std(<NUM_LIT:0>)<EOL>return (G - mean) / std<EOL>
r"""
    Standardize variantes.

    Parameters
    ----------
    G : (`n_inds`, `n_snps`) array
        Genetic data

    Returns
    -------
    G_out : standardized array