def nonsynonymous_by_gene(args):
    query = ("SELECT variant_id, gene from variants WHERE "
             "codon_change != 'None'")
    _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/nonsynonymous_by_gene