def q2q3_mean(coverage):
    """ Calculate the mean of the 2nd and 3rd quartiles """
    if len(coverage) < 4:
        return coverage.mean()
    coverage.sort()
    quartile_size = int(len(coverage) / 4)
    return coverage[quartile_size:-quartile_size].mean()