I'm a data scientist trying to perform predictive modeling on clinical trial patients based on TMB. However, the data provided was raw fast-q files. We have succeeded in generating BAM alignment files. How do I measure non-synonymous mutations per Mb in these files? Is there any R package I can use for the same?