怎么在下一代测序数据bam文件中做indel标记？

我通过bwa得到了sam数据，并经过samtools转换得到bam文件，我现在需要做indel的模拟，
请问我应该怎么做模拟呢？