def get_Qs(alpha, beta, seq_dict, gs):
    """Calculates the observed amino acid frequency with pseudocounts,
    see formula 5 of NucAcRes paper"""
    return {aa: (alpha * seq_dict[aa] + beta * gs[aa]) / (alpha + beta) for aa in gs.keys()}