def consensus_quality(qual_list):
    """
    Create a consensus quality using the average of the qualities
    """
    qual_avg = [int(sum(qual_score) / len(qual_list)) for qual_score in zip(*qual_list)]
    return [x if x < 41 else 41 for x in qual_avg]