HDplot_simulation

R code to simulate data for HDplot. This code also runs HDplot internally and produces plots showing expected distributions based on simulation parameters. Simulation parameters include the number of singleton, duplicate, and diverged duplicate loci, the total population size, the sampled population size, average read depth per locus, statistical distribution of reads per locus to sample from and distribution of reads per allele to sample from.