Given BAM file and its BAI file, how can I calculate the TMB for the patient? This bam file is the result of target sequencing, so I would like to calculate the TMB over the genomic region of interest. Is there any package or software that can do this?