Off topic:Tool for computing Tumor Mutational Burden?

Question

[Pin.Bioinf]

Hello,

I know that Tumor Mutational Burden is the 'number of variants per megabase' but, how can I compute this? are there any bcftools or samtools commands or R packages that can obtain this information from the vcf files? I have no idea about this, sorry.

Thank you !