Simulate Illumina next-gen reads to determine the settings that will allow detection of rare variants (freq 0.001).