How to run your code on a pair of bam files?

[bagavi]

Hi

I had a pair of normal and tumor bam files, and I have a bed file containing ~90,000 MS loci (as per your nature paper).
What is the simplest way to run your code?

Thanks in advance