#!/bin/bash -Ceuo pipefail
calculate_tmb.py \
--in-bed xgen-exome-hyb-panel-v2-targets-hg38_short.mrg_chr21.bed \
--in-vcf HCC1395_tumor_vs_HCC1395_BL.tnseq.filtered_snpEff.ann.vcf.gz \
--threshold 10.0 \
--tumor-sample-name 'HCC1395_tumor' \
--normal-sample-name 'HCC1395_BL' \
--ignore-high-ploidy \
--filter-for-pass \
--filter-for-biallelic \
--min-tumor-vaf 0.05 \
--out-summary-file HCC1395_tumor_vs_HCC1395_BL.tnseq.filtered_snpEff.ann.tmb.txt \
--out-tmb-summary-csv HCC1395_tumor_vs_HCC1395_BL.tnseq.filtered_snpEff.ann.tmb.csv \
--out-variant-list-file HCC1395_tumor_vs_HCC1395_BL.tnseq.filtered_snpEff.ann.vars.txt