def burden_by_gene(args):
    query = ("SELECT gene from variants WHERE "
             "is_coding=1 and (impact_severity = 'HIGH' or "
             "polyphen_pred = 'probably_damaging')")
    _summarize_by_gene_and_sample(args, query)
calculates per sample the total genetic burden for each gene
dataset/ETHPy150Open arq5x/gemini/gemini/tool_burden_tests.py/burden_by_gene