def N_to_Nbp(N, a, b, l):
    return N * (l / b) * (a / l) ** 2  # number of base pairs (bp)