I want to test my SNV calling pipeline using exome sequencing data. So I need simulation data with mutation positions known. I used Wessim to produce simulated reads. But I don't know how to put artificial mutations to the genome positions as I like. Look forward to your suggestion. Thx