Add script to estimate theta at 18 loci for bad mutations samples