RandomSNPs[trees, theta] generates SNPs by simulating mutations on local trees with scaled mutation rate theta. Returns a list with element consisting of SNP location and allele vector at the location. In the allele vector, 0 deotes ancestral allele and 1 denotes mutant allele. |