Bioinformatics Group

My Home Page

Marker Selection Methods

Description of Method (PDF)

Input File Format

Running the Programs

Output

Whole Chromosome Analysis

Download C Source Code

Web Server


Wellcome Trust Centre for Human Genetics

SNP Selection Page

Output Format

Both span_haplotype, greedy_haplotype have the same output format, but big_haplotype has a different output format

span_haplotype, greedy_haplotype

  • The sequence of best marker subsets of size m=1,2,3...max, found by the program are printed, separated by **************************
  • The start of each sectioin looks like this:
    
              **************************
              total_haps = 48 len = 17 done = 1 used = 3 n_haps = 7 entropy 2.911532e+00 5.003353e+00 SNPS used: 10 15 2 
    
    

    The header tells us there are 48 haplotypes of length 17 markers; the best subset containing 3 snps is reported, which divides the data into 7 haplotypes. The entropy is 2.911, compared with 5.00 when all markers are used. The SNPs used are numbered 10, 15, 2

  • The header is followed by the haplotypes sorted into the merged haplotypes:
    
                              -*-------*----*--
                 1         10 11212121112112112 111 1.834862e-02 1.834862e-02 3.243802e+00
                 1          6 11211112212112112 111 9.174312e-03 2.752294e-02 3.243802e+00
    

    The columns are: the merged haplotype id, the original haplotype id, the original haplotype, the merged haplotype, the probability of the haplotype, the cumulated probability for the merged haplotype, the cumulative entropy.

An example of the output for the input data is given here:

	      
        
        
          **************************
          total_haps = 48 len = 17 done = 1 used = 1 n_haps = 1 entropy 9.994535e-01 5.003353e+00 SNPS used: 10 
                          ---------*-------
             1         14 11212121121112111 2 9.174312e-03 9.174312e-03 9.629617e-01
             1         20 11221112212122211 1 9.174312e-03 1.834862e-02 9.629617e-01
             1         26 12111122212122211 1 9.174312e-03 2.752294e-02 9.629617e-01
             1          7 11211122212122211 1 9.174312e-03 3.669725e-02 9.629617e-01
             1          4 11121122212122211 1 3.669725e-02 7.339450e-02 9.629617e-01
             1         40 21121122212122211 1 1.834862e-02 9.174312e-02 9.629617e-01
             1         46 22121122212122211 1 3.669725e-02 1.284404e-01 9.629617e-01
             1         22 11221122212122211 1 9.174312e-03 1.376147e-01 9.629617e-01
             1         44 21221122212122211 1 1.834862e-02 1.559633e-01 9.629617e-01
             1         27 12121112212222211 1 1.834862e-02 1.743119e-01 9.629617e-01
             2         11 11212121121111112 2 9.174312e-03 9.174312e-03 1.402273e+00
             2         48 22221211121112112 2 9.174312e-03 1.834862e-02 1.402273e+00
             2         29 12121121121112112 2 9.174312e-03 2.752294e-02 1.402273e+00
             2         15 11212121121112112 2 2.752294e-02 5.504587e-02 1.402273e+00
             2         18 11212221121112112 2 9.174312e-03 6.422018e-02 1.402273e+00
             2         10 11212121112112112 1 1.834862e-02 8.256881e-02 1.402273e+00
             2          6 11211112212112112 1 9.174312e-03 9.174312e-02 1.402273e+00
             2          3 11121122212112112 1 9.174312e-03 1.009174e-01 1.402273e+00
             2         39 21121122212112112 1 3.669725e-02 1.376147e-01 1.402273e+00
             2         33 12121122212112112 1 0.000000e+00 1.376147e-01 1.402273e+00
             2         34 12121122212112112 1 1.009174e-01 2.385321e-01 1.402273e+00
             2         25 12111121122112112 2 9.174312e-03 2.477064e-01 1.402273e+00
             2         38 21121121122112112 2 4.587156e-02 2.935780e-01 1.402273e+00
             2         31 12121121122112112 2 7.339450e-02 3.669725e-01 1.402273e+00
             2         42 21221121122112112 2 9.174312e-03 3.761468e-01 1.402273e+00
             2          2 11112121122112112 2 9.174312e-03 3.853211e-01 1.402273e+00
             2         17 11212121122112112 2 9.174312e-02 4.770642e-01 1.402273e+00
             2         41 21121122212212112 1 1.834862e-02 4.954128e-01 1.402273e+00
             2         36 12121122212212112 1 4.587156e-02 5.412844e-01 1.402273e+00
             2         37 12221122212212112 1 9.174312e-03 5.504587e-01 1.402273e+00
             2         32 12121121122122112 2 9.174312e-03 5.596330e-01 1.402273e+00
             2          8 11212111121111212 2 1.834862e-02 5.779817e-01 1.402273e+00
             2         47 22221211121111212 2 4.587156e-02 6.238532e-01 1.402273e+00
             2         28 12121121121111212 2 9.174312e-03 6.330275e-01 1.402273e+00
             2         12 11212121121111212 2 9.174312e-03 6.422018e-01 1.402273e+00
             2         45 22121122212111212 1 9.174312e-03 6.513761e-01 1.402273e+00
             2         19 11221112212112212 1 9.174312e-03 6.605505e-01 1.402273e+00
             2         21 11221122212112212 1 2.752294e-02 6.880734e-01 1.402273e+00
             2         43 21221121122112212 2 1.834862e-02 7.064220e-01 1.402273e+00
             2         24 11221221122212212 2 9.174312e-03 7.155963e-01 1.402273e+00
             2         23 11221122212122212 1 1.834862e-02 7.339450e-01 1.402273e+00
             2          9 11212121111111222 1 9.174312e-03 7.431193e-01 1.402273e+00
             2          5 11211111121111222 2 9.174312e-03 7.522936e-01 1.402273e+00
             2          1 11112121121111222 2 9.174312e-03 7.614679e-01 1.402273e+00
             2         13 11212121121111222 2 3.669725e-02 7.981651e-01 1.402273e+00
             2         16 11212121122111222 2 0.000000e+00 7.981651e-01 1.402273e+00
             2         30 12121121121112222 2 2.752294e-02 8.256881e-01 1.402273e+00
             2         35 12121122212112222 1 0.000000e+00 8.256881e-01 1.630436e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 2 n_haps = 3 entropy 1.975919e+00 5.003353e+00 SNPS used: 10 15 
                          ---------*----*--
             1         20 11221112212122211 12 9.174312e-03 9.174312e-03 2.186196e+00
             1         26 12111122212122211 12 9.174312e-03 1.834862e-02 2.186196e+00
             1          7 11211122212122211 12 9.174312e-03 2.752294e-02 2.186196e+00
             1          4 11121122212122211 12 3.669725e-02 6.422018e-02 2.186196e+00
             1         40 21121122212122211 12 1.834862e-02 8.256881e-02 2.186196e+00
             1         46 22121122212122211 12 3.669725e-02 1.192661e-01 2.186196e+00
             1         22 11221122212122211 12 9.174312e-03 1.284404e-01 2.186196e+00
             1         44 21221122212122211 12 1.834862e-02 1.467890e-01 2.186196e+00
             1         27 12121112212222211 12 1.834862e-02 1.651376e-01 2.186196e+00
             2         10 11212121112112112 11 1.834862e-02 1.834862e-02 2.615266e+00
             2          6 11211112212112112 11 9.174312e-03 2.752294e-02 2.615266e+00
             2          3 11121122212112112 11 9.174312e-03 3.669725e-02 2.615266e+00
             2         39 21121122212112112 11 3.669725e-02 7.339450e-02 2.615266e+00
             2         33 12121122212112112 11 0.000000e+00 7.339450e-02 2.615266e+00
             2         34 12121122212112112 11 1.009174e-01 1.743119e-01 2.615266e+00
             2         41 21121122212212112 11 1.834862e-02 1.926606e-01 2.615266e+00
             2         36 12121122212212112 11 4.587156e-02 2.385321e-01 2.615266e+00
             2         37 12221122212212112 11 9.174312e-03 2.477064e-01 2.615266e+00
             2         45 22121122212111212 12 9.174312e-03 2.568807e-01 2.615266e+00
             2         19 11221112212112212 12 9.174312e-03 2.660550e-01 2.615266e+00
             2         21 11221122212112212 12 2.752294e-02 2.935780e-01 2.615266e+00
             2         23 11221122212122212 12 1.834862e-02 3.119266e-01 2.615266e+00
             2          9 11212121111111222 12 9.174312e-03 3.211009e-01 2.615266e+00
             2         35 12121122212112222 12 0.000000e+00 3.211009e-01 2.615266e+00
             3         14 11212121121112111 21 9.174312e-03 9.174312e-03 3.141519e+00
             4         11 11212121121111112 21 9.174312e-03 9.174312e-03 3.203612e+00
             4         48 22221211121112112 21 9.174312e-03 1.834862e-02 3.203612e+00
             4         29 12121121121112112 21 9.174312e-03 2.752294e-02 3.203612e+00
             4         15 11212121121112112 21 2.752294e-02 5.504587e-02 3.203612e+00
             4         18 11212221121112112 21 9.174312e-03 6.422018e-02 3.203612e+00
             4         25 12111121122112112 21 9.174312e-03 7.339450e-02 3.203612e+00
             4         38 21121121122112112 21 4.587156e-02 1.192661e-01 3.203612e+00
             4         31 12121121122112112 21 7.339450e-02 1.926606e-01 3.203612e+00
             4         42 21221121122112112 21 9.174312e-03 2.018349e-01 3.203612e+00
             4          2 11112121122112112 21 9.174312e-03 2.110092e-01 3.203612e+00
             4         17 11212121122112112 21 9.174312e-02 3.027523e-01 3.203612e+00
             4         32 12121121122122112 21 9.174312e-03 3.119266e-01 3.203612e+00
             4          8 11212111121111212 22 1.834862e-02 3.302752e-01 3.203612e+00
             4         47 22221211121111212 22 4.587156e-02 3.761468e-01 3.203612e+00
             4         28 12121121121111212 22 9.174312e-03 3.853211e-01 3.203612e+00
             4         12 11212121121111212 22 9.174312e-03 3.944954e-01 3.203612e+00
             4         43 21221121122112212 22 1.834862e-02 4.128440e-01 3.203612e+00
             4         24 11221221122212212 22 9.174312e-03 4.220183e-01 3.203612e+00
             4          5 11211111121111222 22 9.174312e-03 4.311927e-01 3.203612e+00
             4          1 11112121121111222 22 9.174312e-03 4.403670e-01 3.203612e+00
             4         13 11212121121111222 22 3.669725e-02 4.770642e-01 3.203612e+00
             4         16 11212121122111222 22 0.000000e+00 4.770642e-01 3.203612e+00
             4         30 12121121121112222 22 2.752294e-02 5.045872e-01 3.701551e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 3 n_haps = 7 entropy 2.911532e+00 5.003353e+00 SNPS used: 10 15 2 
                          -*-------*----*--
             1         10 11212121112112112 111 1.834862e-02 1.834862e-02 3.243802e+00
             1          6 11211112212112112 111 9.174312e-03 2.752294e-02 3.243802e+00
             1          3 11121122212112112 111 9.174312e-03 3.669725e-02 3.243802e+00
             1         39 21121122212112112 111 3.669725e-02 7.339450e-02 3.243802e+00
             1         33 12121122212112112 112 0.000000e+00 7.339450e-02 3.243802e+00
             1         34 12121122212112112 112 1.009174e-01 1.743119e-01 3.243802e+00
             1         41 21121122212212112 111 1.834862e-02 1.926606e-01 3.243802e+00
             1         36 12121122212212112 112 4.587156e-02 2.385321e-01 3.243802e+00
             1         37 12221122212212112 112 9.174312e-03 2.477064e-01 3.243802e+00
             2         20 11221112212122211 121 9.174312e-03 9.174312e-03 3.742508e+00
             2         26 12111122212122211 122 9.174312e-03 1.834862e-02 3.742508e+00
             2          7 11211122212122211 121 9.174312e-03 2.752294e-02 3.742508e+00
             2          4 11121122212122211 121 3.669725e-02 6.422018e-02 3.742508e+00
             2         40 21121122212122211 121 1.834862e-02 8.256881e-02 3.742508e+00
             2         46 22121122212122211 122 3.669725e-02 1.192661e-01 3.742508e+00
             2         22 11221122212122211 121 9.174312e-03 1.284404e-01 3.742508e+00
             2         44 21221122212122211 121 1.834862e-02 1.467890e-01 3.742508e+00
             2         27 12121112212222211 122 1.834862e-02 1.651376e-01 3.742508e+00
             3         45 22121122212111212 122 9.174312e-03 9.174312e-03 4.171578e+00
             3         19 11221112212112212 121 9.174312e-03 1.834862e-02 4.171578e+00
             3         21 11221122212112212 121 2.752294e-02 4.587156e-02 4.171578e+00
             3         23 11221122212122212 121 1.834862e-02 6.422018e-02 4.171578e+00
             3          9 11212121111111222 121 9.174312e-03 7.339450e-02 4.171578e+00
             3         35 12121122212112222 122 0.000000e+00 7.339450e-02 4.171578e+00
             4         14 11212121121112111 211 9.174312e-03 9.174312e-03 4.448142e+00
             5         11 11212121121111112 211 9.174312e-03 9.174312e-03 4.510236e+00
             5         48 22221211121112112 212 9.174312e-03 1.834862e-02 4.510236e+00
             5         29 12121121121112112 212 9.174312e-03 2.752294e-02 4.510236e+00
             5         15 11212121121112112 211 2.752294e-02 5.504587e-02 4.510236e+00
             5         18 11212221121112112 211 9.174312e-03 6.422018e-02 4.510236e+00
             5         25 12111121122112112 212 9.174312e-03 7.339450e-02 4.510236e+00
             5         38 21121121122112112 211 4.587156e-02 1.192661e-01 4.510236e+00
             5         31 12121121122112112 212 7.339450e-02 1.926606e-01 4.510236e+00
             5         42 21221121122112112 211 9.174312e-03 2.018349e-01 4.510236e+00
             5          2 11112121122112112 211 9.174312e-03 2.110092e-01 4.510236e+00
             5         17 11212121122112112 211 9.174312e-02 3.027523e-01 4.510236e+00
             5         32 12121121122122112 212 9.174312e-03 3.119266e-01 4.510236e+00
             6          8 11212111121111212 221 1.834862e-02 1.834862e-02 5.034498e+00
             6         47 22221211121111212 222 4.587156e-02 6.422018e-02 5.034498e+00
             6         28 12121121121111212 222 9.174312e-03 7.339450e-02 5.034498e+00
             6         12 11212121121111212 221 9.174312e-03 8.256881e-02 5.034498e+00
             6         43 21221121122112212 221 1.834862e-02 1.009174e-01 5.034498e+00
             6         24 11221221122212212 221 9.174312e-03 1.100917e-01 5.034498e+00
             6          5 11211111121111222 221 9.174312e-03 1.192661e-01 5.034498e+00
             6          1 11112121121111222 221 9.174312e-03 1.284404e-01 5.034498e+00
             6         13 11212121121111222 221 3.669725e-02 1.651376e-01 5.034498e+00
             6         16 11212121122111222 221 0.000000e+00 1.651376e-01 5.034498e+00
             6         30 12121121121112222 222 2.752294e-02 1.926606e-01 5.492233e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 4 n_haps = 14 entropy 3.570538e+00 5.003353e+00 SNPS used: 10 15 2 1 
                          **-------*----*--
             1         10 11212121112112112 1111 1.834862e-02 1.834862e-02 4.566240e+00
             1          6 11211112212112112 1111 9.174312e-03 2.752294e-02 4.566240e+00
             1          3 11121122212112112 1111 9.174312e-03 3.669725e-02 4.566240e+00
             1         39 21121122212112112 1112 3.669725e-02 7.339450e-02 4.566240e+00
             1         41 21121122212212112 1112 1.834862e-02 9.174312e-02 4.566240e+00
             2         33 12121122212112112 1121 0.000000e+00 0.000000e+00 4.882410e+00
             2         34 12121122212112112 1121 1.009174e-01 1.009174e-01 4.882410e+00
             2         36 12121122212212112 1121 4.587156e-02 1.467890e-01 4.882410e+00
             2         37 12221122212212112 1121 9.174312e-03 1.559633e-01 4.882410e+00
             3         20 11221112212122211 1211 9.174312e-03 9.174312e-03 5.300504e+00
             3          7 11211122212122211 1211 9.174312e-03 1.834862e-02 5.300504e+00
             3          4 11121122212122211 1211 3.669725e-02 5.504587e-02 5.300504e+00
             3         40 21121122212122211 1212 1.834862e-02 7.339450e-02 5.300504e+00
             3         22 11221122212122211 1211 9.174312e-03 8.256881e-02 5.300504e+00
             3         44 21221122212122211 1212 1.834862e-02 1.009174e-01 5.300504e+00
             4         19 11221112212112212 1211 9.174312e-03 9.174312e-03 5.634415e+00
             4         21 11221122212112212 1211 2.752294e-02 3.669725e-02 5.634415e+00
             4         23 11221122212122212 1211 1.834862e-02 5.504587e-02 5.634415e+00
             4          9 11212121111111222 1211 9.174312e-03 6.422018e-02 5.634415e+00
             5         26 12111122212122211 1221 9.174312e-03 9.174312e-03 5.888780e+00
             5         46 22121122212122211 1222 3.669725e-02 4.587156e-02 5.888780e+00
             5         27 12121112212222211 1221 1.834862e-02 6.422018e-02 5.888780e+00
             6         45 22121122212111212 1222 9.174312e-03 9.174312e-03 6.143145e+00
             6         35 12121122212112222 1221 0.000000e+00 9.174312e-03 6.143145e+00
             7         14 11212121121112111 2111 9.174312e-03 9.174312e-03 6.205239e+00
             8         11 11212121121111112 2111 9.174312e-03 9.174312e-03 6.267332e+00
             8         15 11212121121112112 2111 2.752294e-02 3.669725e-02 6.267332e+00
             8         18 11212221121112112 2111 9.174312e-03 4.587156e-02 6.267332e+00
             8         38 21121121122112112 2112 4.587156e-02 9.174312e-02 6.267332e+00
             8         42 21221121122112112 2112 9.174312e-03 1.009174e-01 6.267332e+00
             8          2 11112121122112112 2111 9.174312e-03 1.100917e-01 6.267332e+00
             8         17 11212121122112112 2111 9.174312e-02 2.018349e-01 6.267332e+00
             9         48 22221211121112112 2122 9.174312e-03 9.174312e-03 6.733319e+00
             9         29 12121121121112112 2121 9.174312e-03 1.834862e-02 6.733319e+00
             9         25 12111121122112112 2121 9.174312e-03 2.752294e-02 6.733319e+00
             9         31 12121121122112112 2121 7.339450e-02 1.009174e-01 6.733319e+00
             9         32 12121121122122112 2121 9.174312e-03 1.100917e-01 6.733319e+00
            10          8 11212111121111212 2211 1.834862e-02 1.834862e-02 7.083765e+00
            10         12 11212121121111212 2211 9.174312e-03 2.752294e-02 7.083765e+00
            10         43 21221121122112212 2212 1.834862e-02 4.587156e-02 7.083765e+00
            10         24 11221221122212212 2211 9.174312e-03 5.504587e-02 7.083765e+00
            10          5 11211111121111222 2211 9.174312e-03 6.422018e-02 7.083765e+00
            10          1 11112121121111222 2211 9.174312e-03 7.339450e-02 7.083765e+00
            10         13 11212121121111222 2211 3.669725e-02 1.100917e-01 7.083765e+00
            10         16 11212121122111222 2211 0.000000e+00 1.100917e-01 7.083765e+00
            11         47 22221211121111212 2222 4.587156e-02 4.587156e-02 7.434212e+00
            11         28 12121121121111212 2221 9.174312e-03 5.504587e-02 7.434212e+00
            11         30 12121121121112222 2221 2.752294e-02 8.256881e-02 7.731316e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 5 n_haps = 21 entropy 3.927351e+00 5.003353e+00 SNPS used: 10 15 2 1 3 
                          ***------*----*--
             1         10 11212121112112112 11112 1.834862e-02 1.834862e-02 5.456897e+00
             1          6 11211112212112112 11112 9.174312e-03 2.752294e-02 5.456897e+00
             1          3 11121122212112112 11111 9.174312e-03 3.669725e-02 5.456897e+00
             2         39 21121122212112112 11121 3.669725e-02 3.669725e-02 5.631876e+00
             2         41 21121122212212112 11121 1.834862e-02 5.504587e-02 5.631876e+00
             3         33 12121122212112112 11211 0.000000e+00 0.000000e+00 5.862145e+00
             3         34 12121122212112112 11211 1.009174e-01 1.009174e-01 5.862145e+00
             3         36 12121122212212112 11211 4.587156e-02 1.467890e-01 5.862145e+00
             3         37 12221122212212112 11212 9.174312e-03 1.559633e-01 5.862145e+00
             4         20 11221112212122211 12112 9.174312e-03 9.174312e-03 6.280239e+00
             4          7 11211122212122211 12112 9.174312e-03 1.834862e-02 6.280239e+00
             4          4 11121122212122211 12111 3.669725e-02 5.504587e-02 6.280239e+00
             4         22 11221122212122211 12112 9.174312e-03 6.422018e-02 6.280239e+00
             5         19 11221112212112212 12112 9.174312e-03 9.174312e-03 6.534604e+00
             5         21 11221122212112212 12112 2.752294e-02 3.669725e-02 6.534604e+00
             5         23 11221122212122212 12112 1.834862e-02 5.504587e-02 6.534604e+00
             5          9 11212121111111222 12112 9.174312e-03 6.422018e-02 6.534604e+00
             6         40 21121122212122211 12121 1.834862e-02 1.834862e-02 6.788970e+00
             6         44 21221122212122211 12122 1.834862e-02 3.669725e-02 6.788970e+00
             7         26 12111122212122211 12211 9.174312e-03 9.174312e-03 6.963949e+00
             7         27 12121112212222211 12211 1.834862e-02 2.752294e-02 6.963949e+00
             8         35 12121122212112222 12211 0.000000e+00 0.000000e+00 7.106606e+00
             9         46 22121122212122211 12221 3.669725e-02 3.669725e-02 7.106606e+00
            10         45 22121122212111212 12221 9.174312e-03 9.174312e-03 7.281586e+00
            11         14 11212121121112111 21112 9.174312e-03 9.174312e-03 7.343679e+00
            12         11 11212121121111112 21112 9.174312e-03 9.174312e-03 7.405772e+00
            12         15 11212121121112112 21112 2.752294e-02 3.669725e-02 7.405772e+00
            12         18 11212221121112112 21112 9.174312e-03 4.587156e-02 7.405772e+00
            12          2 11112121122112112 21111 9.174312e-03 5.504587e-02 7.405772e+00
            12         17 11212121122112112 21112 9.174312e-02 1.467890e-01 7.405772e+00
            13         38 21121121122112112 21121 4.587156e-02 4.587156e-02 7.812111e+00
            13         42 21221121122112112 21122 9.174312e-03 5.504587e-02 7.812111e+00
            14         29 12121121121112112 21211 9.174312e-03 9.174312e-03 8.042380e+00
            14         25 12111121122112112 21211 9.174312e-03 1.834862e-02 8.042380e+00
            14         31 12121121122112112 21211 7.339450e-02 9.174312e-02 8.042380e+00
            14         32 12121121122122112 21211 9.174312e-03 1.009174e-01 8.042380e+00
            15         48 22221211121112112 21222 9.174312e-03 9.174312e-03 8.376291e+00
            16          8 11212111121111212 22112 1.834862e-02 1.834862e-02 8.438385e+00
            16         12 11212121121111212 22112 9.174312e-03 2.752294e-02 8.438385e+00
            16         24 11221221122212212 22112 9.174312e-03 3.669725e-02 8.438385e+00
            16          5 11211111121111222 22112 9.174312e-03 4.587156e-02 8.438385e+00
            16          1 11112121121111222 22111 9.174312e-03 5.504587e-02 8.438385e+00
            16         13 11212121121111222 22112 3.669725e-02 9.174312e-02 8.438385e+00
            16         16 11212121122111222 22112 0.000000e+00 9.174312e-02 8.438385e+00
            17         43 21221121122112212 22122 1.834862e-02 1.834862e-02 8.754555e+00
            18         28 12121121121111212 22211 9.174312e-03 9.174312e-03 8.860393e+00
            18         30 12121121121112222 22211 2.752294e-02 3.669725e-02 8.860393e+00
            19         47 22221211121111212 22222 4.587156e-02 4.587156e-02 9.239329e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 6 n_haps = 25 entropy 4.196385e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 
                          ***------**---*--
             1          3 11121122212112112 111112 9.174312e-03 9.174312e-03 5.901820e+00
             2         10 11212121112112112 111122 1.834862e-02 1.834862e-02 5.963913e+00
             2          6 11211112212112112 111122 9.174312e-03 2.752294e-02 5.963913e+00
             3         39 21121122212112112 111212 3.669725e-02 3.669725e-02 6.106571e+00
             3         41 21121122212212112 111212 1.834862e-02 5.504587e-02 6.106571e+00
             4         33 12121122212112112 112112 0.000000e+00 0.000000e+00 6.336840e+00
             4         34 12121122212112112 112112 1.009174e-01 1.009174e-01 6.336840e+00
             4         36 12121122212212112 112112 4.587156e-02 1.467890e-01 6.336840e+00
             5         37 12221122212212112 112122 9.174312e-03 9.174312e-03 6.743179e+00
             6          4 11121122212122211 121112 3.669725e-02 3.669725e-02 6.805272e+00
             7         20 11221112212122211 121122 9.174312e-03 9.174312e-03 6.980252e+00
             7          7 11211122212122211 121122 9.174312e-03 1.834862e-02 6.980252e+00
             7         22 11221122212122211 121122 9.174312e-03 2.752294e-02 6.980252e+00
             8         19 11221112212112212 121122 9.174312e-03 9.174312e-03 7.122909e+00
             8         21 11221122212112212 121122 2.752294e-02 3.669725e-02 7.122909e+00
             8         23 11221122212122212 121122 1.834862e-02 5.504587e-02 7.122909e+00
             8          9 11212121111111222 121121 9.174312e-03 6.422018e-02 7.122909e+00
             9         40 21121122212122211 121212 1.834862e-02 1.834862e-02 7.377274e+00
            10         44 21221122212122211 121222 1.834862e-02 1.834862e-02 7.483113e+00
            11         26 12111122212122211 122112 9.174312e-03 9.174312e-03 7.588951e+00
            11         27 12121112212222211 122112 1.834862e-02 2.752294e-02 7.588951e+00
            12         35 12121122212112222 122112 0.000000e+00 0.000000e+00 7.731608e+00
            13         46 22121122212122211 122212 3.669725e-02 3.669725e-02 7.731608e+00
            14         45 22121122212111212 122212 9.174312e-03 9.174312e-03 7.906588e+00
            15          2 11112121122112112 211112 9.174312e-03 9.174312e-03 7.968681e+00
            16         14 11212121121112111 211121 9.174312e-03 9.174312e-03 8.030774e+00
            17         11 11212121121111112 211121 9.174312e-03 9.174312e-03 8.092868e+00
            17         15 11212121121112112 211121 2.752294e-02 3.669725e-02 8.092868e+00
            17         18 11212221121112112 211121 9.174312e-03 4.587156e-02 8.092868e+00
            17         17 11212121122112112 211122 9.174312e-02 1.376147e-01 8.092868e+00
            18         38 21121121122112112 211212 4.587156e-02 4.587156e-02 8.486624e+00
            19         42 21221121122112112 211222 9.174312e-03 9.174312e-03 8.690581e+00
            20         29 12121121121112112 212111 9.174312e-03 9.174312e-03 8.752674e+00
            20         25 12111121122112112 212112 9.174312e-03 1.834862e-02 8.752674e+00
            20         31 12121121122112112 212112 7.339450e-02 9.174312e-02 8.752674e+00
            20         32 12121121122122112 212112 9.174312e-03 1.009174e-01 8.752674e+00
            21         48 22221211121112112 212221 9.174312e-03 9.174312e-03 9.086585e+00
            22          1 11112121121111222 221111 9.174312e-03 9.174312e-03 9.148678e+00
            23          8 11212111121111212 221121 1.834862e-02 1.834862e-02 9.210772e+00
            23         12 11212121121111212 221121 9.174312e-03 2.752294e-02 9.210772e+00
            23         24 11221221122212212 221122 9.174312e-03 3.669725e-02 9.210772e+00
            23          5 11211111121111222 221121 9.174312e-03 4.587156e-02 9.210772e+00
            23         13 11212121121111222 221121 3.669725e-02 8.256881e-02 9.210772e+00
            23         16 11212121122111222 221122 0.000000e+00 8.256881e-02 9.210772e+00
            24         43 21221121122112212 221222 1.834862e-02 1.834862e-02 9.507876e+00
            25         28 12121121121111212 222111 9.174312e-03 9.174312e-03 9.613714e+00
            25         30 12121121121112222 222111 2.752294e-02 3.669725e-02 9.613714e+00
            26         47 22221211121111212 222221 4.587156e-02 4.587156e-02 9.992650e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 7 n_haps = 29 entropy 4.403736e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 
                          ***------***--*--
             1          3 11121122212112112 1111121 9.174312e-03 9.174312e-03 6.262890e+00
             2         10 11212121112112112 1111221 1.834862e-02 1.834862e-02 6.324983e+00
             2          6 11211112212112112 1111221 9.174312e-03 2.752294e-02 6.324983e+00
             3         39 21121122212112112 1112121 3.669725e-02 3.669725e-02 6.467641e+00
             3         41 21121122212212112 1112122 1.834862e-02 5.504587e-02 6.467641e+00
             4         33 12121122212112112 1121121 0.000000e+00 0.000000e+00 6.697910e+00
             4         34 12121122212112112 1121121 1.009174e-01 1.009174e-01 6.697910e+00
             4         36 12121122212212112 1121122 4.587156e-02 1.467890e-01 6.697910e+00
             5         37 12221122212212112 1121222 9.174312e-03 9.174312e-03 7.104249e+00
             6          4 11121122212122211 1211121 3.669725e-02 3.669725e-02 7.166342e+00
             7          9 11212121111111222 1211211 9.174312e-03 9.174312e-03 7.341322e+00
             8         20 11221112212122211 1211221 9.174312e-03 9.174312e-03 7.403415e+00
             8          7 11211122212122211 1211221 9.174312e-03 1.834862e-02 7.403415e+00
             8         22 11221122212122211 1211221 9.174312e-03 2.752294e-02 7.403415e+00
             9         19 11221112212112212 1211221 9.174312e-03 9.174312e-03 7.546073e+00
             9         21 11221122212112212 1211221 2.752294e-02 3.669725e-02 7.546073e+00
             9         23 11221122212122212 1211221 1.834862e-02 5.504587e-02 7.546073e+00
            10         40 21121122212122211 1212121 1.834862e-02 1.834862e-02 7.776342e+00
            11         44 21221122212122211 1212221 1.834862e-02 1.834862e-02 7.882180e+00
            12         26 12111122212122211 1221121 9.174312e-03 9.174312e-03 7.988018e+00
            12         27 12121112212222211 1221122 1.834862e-02 2.752294e-02 7.988018e+00
            13         35 12121122212112222 1221121 0.000000e+00 0.000000e+00 8.130676e+00
            14         46 22121122212122211 1222121 3.669725e-02 3.669725e-02 8.130676e+00
            15         45 22121122212111212 1222121 9.174312e-03 9.174312e-03 8.305655e+00
            16          2 11112121122112112 2111121 9.174312e-03 9.174312e-03 8.367748e+00
            17         14 11212121121112111 2111211 9.174312e-03 9.174312e-03 8.429842e+00
            18         11 11212121121111112 2111211 9.174312e-03 9.174312e-03 8.491935e+00
            18         15 11212121121112112 2111211 2.752294e-02 3.669725e-02 8.491935e+00
            18         18 11212221121112112 2111211 9.174312e-03 4.587156e-02 8.491935e+00
            19         17 11212121122112112 2111221 9.174312e-02 9.174312e-02 8.695892e+00
            20         38 21121121122112112 2112121 4.587156e-02 4.587156e-02 9.012062e+00
            21         42 21221121122112112 2112221 9.174312e-03 9.174312e-03 9.216019e+00
            22         29 12121121121112112 2121111 9.174312e-03 9.174312e-03 9.278112e+00
            23         25 12111121122112112 2121121 9.174312e-03 9.174312e-03 9.340206e+00
            23         31 12121121122112112 2121121 7.339450e-02 8.256881e-02 9.340206e+00
            23         32 12121121122122112 2121121 9.174312e-03 9.174312e-02 9.340206e+00
            24         48 22221211121112112 2122211 9.174312e-03 9.174312e-03 9.656376e+00
            25          1 11112121121111222 2211111 9.174312e-03 9.174312e-03 9.718469e+00
            26          8 11212111121111212 2211211 1.834862e-02 1.834862e-02 9.780563e+00
            26         12 11212121121111212 2211211 9.174312e-03 2.752294e-02 9.780563e+00
            26          5 11211111121111222 2211211 9.174312e-03 3.669725e-02 9.780563e+00
            26         13 11212121121111222 2211211 3.669725e-02 7.339450e-02 9.780563e+00
            27         24 11221221122212212 2211222 9.174312e-03 9.174312e-03 1.005713e+01
            27         16 11212121122111222 2211221 0.000000e+00 9.174312e-03 1.005713e+01
            28         43 21221121122112212 2212221 1.834862e-02 1.834862e-02 1.011922e+01
            29         28 12121121121111212 2221111 9.174312e-03 9.174312e-03 1.022506e+01
            29         30 12121121121112222 2221111 2.752294e-02 3.669725e-02 1.022506e+01
            30         47 22221211121111212 2222211 4.587156e-02 4.587156e-02 1.060399e+01
          **************************
          total_haps = 48 len = 17 done = 1 used = 8 n_haps = 32 entropy 4.562160e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 
                          ***---*--***--*--
             1          3 11121122212112112 11111212 9.174312e-03 9.174312e-03 6.562034e+00
             2         10 11212121112112112 11112212 1.834862e-02 1.834862e-02 6.624128e+00
             2          6 11211112212112112 11112211 9.174312e-03 2.752294e-02 6.624128e+00
             3         39 21121122212112112 11121212 3.669725e-02 3.669725e-02 6.766785e+00
             4         41 21121122212212112 11121222 1.834862e-02 1.834862e-02 6.941765e+00
             5         33 12121122212112112 11211212 0.000000e+00 0.000000e+00 7.047603e+00
             5         34 12121122212112112 11211212 1.009174e-01 1.009174e-01 7.047603e+00
             6         36 12121122212212112 11211222 4.587156e-02 4.587156e-02 7.381514e+00
             7         37 12221122212212112 11212222 9.174312e-03 9.174312e-03 7.585470e+00
             8          4 11121122212122211 12111212 3.669725e-02 3.669725e-02 7.647564e+00
             9          9 11212121111111222 12112112 9.174312e-03 9.174312e-03 7.822543e+00
            10         20 11221112212122211 12112211 9.174312e-03 9.174312e-03 7.884637e+00
            10          7 11211122212122211 12112212 9.174312e-03 1.834862e-02 7.884637e+00
            10         22 11221122212122211 12112212 9.174312e-03 2.752294e-02 7.884637e+00
            11         19 11221112212112212 12112211 9.174312e-03 9.174312e-03 8.027294e+00
            11         21 11221122212112212 12112212 2.752294e-02 3.669725e-02 8.027294e+00
            11         23 11221122212122212 12112212 1.834862e-02 5.504587e-02 8.027294e+00
            12         40 21121122212122211 12121212 1.834862e-02 1.834862e-02 8.257563e+00
            13         44 21221122212122211 12122212 1.834862e-02 1.834862e-02 8.363401e+00
            14         26 12111122212122211 12211212 9.174312e-03 9.174312e-03 8.469240e+00
            15         35 12121122212112222 12211212 0.000000e+00 0.000000e+00 8.531333e+00
            16         27 12121112212222211 12211221 1.834862e-02 1.834862e-02 8.531333e+00
            17         46 22121122212122211 12221212 3.669725e-02 3.669725e-02 8.637171e+00
            18         45 22121122212111212 12221212 9.174312e-03 9.174312e-03 8.812150e+00
            19          2 11112121122112112 21111212 9.174312e-03 9.174312e-03 8.874244e+00
            20         14 11212121121112111 21112112 9.174312e-03 9.174312e-03 8.936337e+00
            21         11 11212121121111112 21112112 9.174312e-03 9.174312e-03 8.998431e+00
            21         15 11212121121112112 21112112 2.752294e-02 3.669725e-02 8.998431e+00
            21         18 11212221121112112 21112112 9.174312e-03 4.587156e-02 8.998431e+00
            22         17 11212121122112112 21112212 9.174312e-02 9.174312e-02 9.202387e+00
            23         38 21121121122112112 21121212 4.587156e-02 4.587156e-02 9.518558e+00
            24         42 21221121122112112 21122212 9.174312e-03 9.174312e-03 9.722514e+00
            25         29 12121121121112112 21211112 9.174312e-03 9.174312e-03 9.784608e+00
            26         25 12111121122112112 21211212 9.174312e-03 9.174312e-03 9.846701e+00
            26         31 12121121122112112 21211212 7.339450e-02 8.256881e-02 9.846701e+00
            26         32 12121121122122112 21211212 9.174312e-03 9.174312e-02 9.846701e+00
            27         48 22221211121112112 21222111 9.174312e-03 9.174312e-03 1.016287e+01
            28          1 11112121121111222 22111112 9.174312e-03 9.174312e-03 1.022497e+01
            29          8 11212111121111212 22112111 1.834862e-02 1.834862e-02 1.028706e+01
            29         12 11212121121111212 22112112 9.174312e-03 2.752294e-02 1.028706e+01
            29          5 11211111121111222 22112111 9.174312e-03 3.669725e-02 1.028706e+01
            29         13 11212121121111222 22112112 3.669725e-02 7.339450e-02 1.028706e+01
            30         16 11212121122111222 22112212 0.000000e+00 0.000000e+00 1.056362e+01
            31         24 11221221122212212 22112222 9.174312e-03 9.174312e-03 1.056362e+01
            32         43 21221121122112212 22122212 1.834862e-02 1.834862e-02 1.062572e+01
            33         28 12121121121111212 22211112 9.174312e-03 9.174312e-03 1.073155e+01
            33         30 12121121121112222 22211112 2.752294e-02 3.669725e-02 1.073155e+01
            34         47 22221211121111212 22222111 4.587156e-02 4.587156e-02 1.111049e+01
          **************************
          total_haps = 48 len = 17 done = 1 used = 9 n_haps = 37 entropy 4.719923e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 
                          ***---*--****-*--
             1          3 11121122212112112 111112121 9.174312e-03 9.174312e-03 6.787643e+00
             2          6 11211112212112112 111122111 9.174312e-03 9.174312e-03 6.849736e+00
             3         10 11212121112112112 111122121 1.834862e-02 1.834862e-02 6.911829e+00
             4         39 21121122212112112 111212121 3.669725e-02 3.669725e-02 7.017668e+00
             5         41 21121122212212112 111212221 1.834862e-02 1.834862e-02 7.192647e+00
             6         33 12121122212112112 112112121 0.000000e+00 0.000000e+00 7.298485e+00
             6         34 12121122212112112 112112121 1.009174e-01 1.009174e-01 7.298485e+00
             7         36 12121122212212112 112112221 4.587156e-02 4.587156e-02 7.632396e+00
             8         37 12221122212212112 112122221 9.174312e-03 9.174312e-03 7.836353e+00
             9          4 11121122212122211 121112122 3.669725e-02 3.669725e-02 7.898446e+00
            10          9 11212121111111222 121121121 9.174312e-03 9.174312e-03 8.073425e+00
            11         20 11221112212122211 121122112 9.174312e-03 9.174312e-03 8.135519e+00
            12         19 11221112212112212 121122111 9.174312e-03 9.174312e-03 8.197612e+00
            13          7 11211122212122211 121122122 9.174312e-03 9.174312e-03 8.259706e+00
            13         22 11221122212122211 121122122 9.174312e-03 1.834862e-02 8.259706e+00
            14         21 11221122212112212 121122121 2.752294e-02 2.752294e-02 8.365544e+00
            14         23 11221122212122212 121122122 1.834862e-02 4.587156e-02 8.365544e+00
            15         40 21121122212122211 121212122 1.834862e-02 1.834862e-02 8.569501e+00
            16         44 21221122212122211 121222122 1.834862e-02 1.834862e-02 8.675339e+00
            17         26 12111122212122211 122112122 9.174312e-03 9.174312e-03 8.781177e+00
            18         35 12121122212112222 122112121 0.000000e+00 0.000000e+00 8.843271e+00
            19         27 12121112212222211 122112212 1.834862e-02 1.834862e-02 8.843271e+00
            20         46 22121122212122211 122212122 3.669725e-02 3.669725e-02 8.949109e+00
            21         45 22121122212111212 122212121 9.174312e-03 9.174312e-03 9.124088e+00
            22          2 11112121122112112 211112121 9.174312e-03 9.174312e-03 9.186182e+00
            23         14 11212121121112111 211121121 9.174312e-03 9.174312e-03 9.248275e+00
            24         11 11212121121111112 211121121 9.174312e-03 9.174312e-03 9.310368e+00
            24         15 11212121121112112 211121121 2.752294e-02 3.669725e-02 9.310368e+00
            24         18 11212221121112112 211121121 9.174312e-03 4.587156e-02 9.310368e+00
            25         17 11212121122112112 211122121 9.174312e-02 9.174312e-02 9.514325e+00
            26         38 21121121122112112 211212121 4.587156e-02 4.587156e-02 9.830495e+00
            27         42 21221121122112112 211222121 9.174312e-03 9.174312e-03 1.003445e+01
            28         29 12121121121112112 212111121 9.174312e-03 9.174312e-03 1.009655e+01
            29         25 12111121122112112 212112121 9.174312e-03 9.174312e-03 1.015864e+01
            29         31 12121121122112112 212112121 7.339450e-02 8.256881e-02 1.015864e+01
            29         32 12121121122122112 212112122 9.174312e-03 9.174312e-02 1.015864e+01
            30         48 22221211121112112 212221111 9.174312e-03 9.174312e-03 1.047481e+01
            31          1 11112121121111222 221111121 9.174312e-03 9.174312e-03 1.053690e+01
            32          8 11212111121111212 221121111 1.834862e-02 1.834862e-02 1.059900e+01
            32          5 11211111121111222 221121111 9.174312e-03 2.752294e-02 1.059900e+01
            33         12 11212121121111212 221121121 9.174312e-03 9.174312e-03 1.074165e+01
            33         13 11212121121111222 221121121 3.669725e-02 4.587156e-02 1.074165e+01
            34         16 11212121122111222 221122121 0.000000e+00 0.000000e+00 1.094561e+01
            35         24 11221221122212212 221122221 9.174312e-03 9.174312e-03 1.094561e+01
            36         43 21221121122112212 221222121 1.834862e-02 1.834862e-02 1.100770e+01
            37         28 12121121121111212 222111121 9.174312e-03 9.174312e-03 1.111354e+01
            37         30 12121121121112222 222111121 2.752294e-02 3.669725e-02 1.111354e+01
            38         47 22221211121111212 222221111 4.587156e-02 4.587156e-02 1.149248e+01
          **************************
          total_haps = 48 len = 17 done = 1 used = 10 n_haps = 40 entropy 4.808085e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 
                          ***---*--****-**-
             1          3 11121122212112112 1111121211 9.174312e-03 9.174312e-03 6.913974e+00
             2          6 11211112212112112 1111221111 9.174312e-03 9.174312e-03 6.976067e+00
             3         10 11212121112112112 1111221211 1.834862e-02 1.834862e-02 7.038161e+00
             4         39 21121122212112112 1112121211 3.669725e-02 3.669725e-02 7.143999e+00
             5         41 21121122212212112 1112122211 1.834862e-02 1.834862e-02 7.318978e+00
             6         33 12121122212112112 1121121211 0.000000e+00 0.000000e+00 7.424816e+00
             6         34 12121122212112112 1121121211 1.009174e-01 1.009174e-01 7.424816e+00
             7         36 12121122212212112 1121122211 4.587156e-02 4.587156e-02 7.758727e+00
             8         37 12221122212212112 1121222211 9.174312e-03 9.174312e-03 7.962684e+00
             9          4 11121122212122211 1211121221 3.669725e-02 3.669725e-02 8.024777e+00
            10          9 11212121111111222 1211211212 9.174312e-03 9.174312e-03 8.199757e+00
            11         19 11221112212112212 1211221111 9.174312e-03 9.174312e-03 8.261850e+00
            12         20 11221112212122211 1211221121 9.174312e-03 9.174312e-03 8.323944e+00
            13         21 11221122212112212 1211221211 2.752294e-02 2.752294e-02 8.386037e+00
            14          7 11211122212122211 1211221221 9.174312e-03 9.174312e-03 8.528694e+00
            14         22 11221122212122211 1211221221 9.174312e-03 1.834862e-02 8.528694e+00
            15         23 11221122212122212 1211221221 1.834862e-02 1.834862e-02 8.634533e+00
            16         40 21121122212122211 1212121221 1.834862e-02 1.834862e-02 8.740371e+00
            17         44 21221122212122211 1212221221 1.834862e-02 1.834862e-02 8.846209e+00
            18         35 12121122212112222 1221121212 0.000000e+00 0.000000e+00 8.952047e+00
            19         26 12111122212122211 1221121221 9.174312e-03 9.174312e-03 8.952047e+00
            20         27 12121112212222211 1221122121 1.834862e-02 1.834862e-02 9.014141e+00
            21         45 22121122212111212 1222121211 9.174312e-03 9.174312e-03 9.119979e+00
            22         46 22121122212122211 1222121221 3.669725e-02 3.669725e-02 9.182073e+00
            23          2 11112121122112112 2111121211 9.174312e-03 9.174312e-03 9.357052e+00
            24         14 11212121121112111 2111211211 9.174312e-03 9.174312e-03 9.419145e+00
            25         11 11212121121111112 2111211211 9.174312e-03 9.174312e-03 9.481239e+00
            25         15 11212121121112112 2111211211 2.752294e-02 3.669725e-02 9.481239e+00
            25         18 11212221121112112 2111211211 9.174312e-03 4.587156e-02 9.481239e+00
            26         17 11212121122112112 2111221211 9.174312e-02 9.174312e-02 9.685195e+00
            27         38 21121121122112112 2112121211 4.587156e-02 4.587156e-02 1.000137e+01
            28         42 21221121122112112 2112221211 9.174312e-03 9.174312e-03 1.020532e+01
            29         29 12121121121112112 2121111211 9.174312e-03 9.174312e-03 1.026742e+01
            30         25 12111121122112112 2121121211 9.174312e-03 9.174312e-03 1.032951e+01
            30         31 12121121122112112 2121121211 7.339450e-02 8.256881e-02 1.032951e+01
            31         32 12121121122122112 2121121221 9.174312e-03 9.174312e-03 1.062661e+01
            32         48 22221211121112112 2122211111 9.174312e-03 9.174312e-03 1.068871e+01
            33          1 11112121121111222 2211111212 9.174312e-03 9.174312e-03 1.075080e+01
            34          8 11212111121111212 2211211111 1.834862e-02 1.834862e-02 1.081289e+01
            34          5 11211111121111222 2211211112 9.174312e-03 2.752294e-02 1.081289e+01
            35         12 11212121121111212 2211211211 9.174312e-03 9.174312e-03 1.095555e+01
            35         13 11212121121111222 2211211212 3.669725e-02 4.587156e-02 1.095555e+01
            36         16 11212121122111222 2211221212 0.000000e+00 0.000000e+00 1.115951e+01
            37         24 11221221122212212 2211222211 9.174312e-03 9.174312e-03 1.115951e+01
            38         43 21221121122112212 2212221211 1.834862e-02 1.834862e-02 1.122160e+01
            39         28 12121121121111212 2221111211 9.174312e-03 9.174312e-03 1.132744e+01
            39         30 12121121121112222 2221111212 2.752294e-02 3.669725e-02 1.132744e+01
            40         47 22221211121111212 2222211111 4.587156e-02 4.587156e-02 1.170638e+01
          **************************
          total_haps = 48 len = 17 done = 1 used = 11 n_haps = 42 entropy 4.880563e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 
                          ***---*--****-***
             1          3 11121122212112112 11111212112 9.174312e-03 9.174312e-03 7.041164e+00
             2          6 11211112212112112 11112211112 9.174312e-03 9.174312e-03 7.103258e+00
             3         10 11212121112112112 11112212112 1.834862e-02 1.834862e-02 7.165351e+00
             4         39 21121122212112112 11121212112 3.669725e-02 3.669725e-02 7.271190e+00
             5         41 21121122212212112 11121222112 1.834862e-02 1.834862e-02 7.446169e+00
             6         33 12121122212112112 11211212112 0.000000e+00 0.000000e+00 7.552007e+00
             6         34 12121122212112112 11211212112 1.009174e-01 1.009174e-01 7.552007e+00
             7         36 12121122212212112 11211222112 4.587156e-02 4.587156e-02 7.885918e+00
             8         37 12221122212212112 11212222112 9.174312e-03 9.174312e-03 8.089875e+00
             9          4 11121122212122211 12111212211 3.669725e-02 3.669725e-02 8.151968e+00
            10          9 11212121111111222 12112112122 9.174312e-03 9.174312e-03 8.326947e+00
            11         19 11221112212112212 12112211112 9.174312e-03 9.174312e-03 8.389041e+00
            12         20 11221112212122211 12112211211 9.174312e-03 9.174312e-03 8.451134e+00
            13         21 11221122212112212 12112212112 2.752294e-02 2.752294e-02 8.513228e+00
            14          7 11211122212122211 12112212211 9.174312e-03 9.174312e-03 8.655885e+00
            14         22 11221122212122211 12112212211 9.174312e-03 1.834862e-02 8.655885e+00
            15         23 11221122212122212 12112212212 1.834862e-02 1.834862e-02 8.761723e+00
            16         40 21121122212122211 12121212211 1.834862e-02 1.834862e-02 8.867562e+00
            17         44 21221122212122211 12122212211 1.834862e-02 1.834862e-02 8.973400e+00
            18         35 12121122212112222 12211212122 0.000000e+00 0.000000e+00 9.079238e+00
            19         26 12111122212122211 12211212211 9.174312e-03 9.174312e-03 9.079238e+00
            20         27 12121112212222211 12211221211 1.834862e-02 1.834862e-02 9.141331e+00
            21         45 22121122212111212 12221212112 9.174312e-03 9.174312e-03 9.247170e+00
            22         46 22121122212122211 12221212211 3.669725e-02 3.669725e-02 9.309263e+00
            23          2 11112121122112112 21111212112 9.174312e-03 9.174312e-03 9.484242e+00
            24         14 11212121121112111 21112112111 9.174312e-03 9.174312e-03 9.546336e+00
            25         11 11212121121111112 21112112112 9.174312e-03 9.174312e-03 9.608429e+00
            25         15 11212121121112112 21112112112 2.752294e-02 3.669725e-02 9.608429e+00
            25         18 11212221121112112 21112112112 9.174312e-03 4.587156e-02 9.608429e+00
            26         17 11212121122112112 21112212112 9.174312e-02 9.174312e-02 9.812386e+00
            27         38 21121121122112112 21121212112 4.587156e-02 4.587156e-02 1.012856e+01
            28         42 21221121122112112 21122212112 9.174312e-03 9.174312e-03 1.033251e+01
            29         29 12121121121112112 21211112112 9.174312e-03 9.174312e-03 1.039461e+01
            30         25 12111121122112112 21211212112 9.174312e-03 9.174312e-03 1.045670e+01
            30         31 12121121122112112 21211212112 7.339450e-02 8.256881e-02 1.045670e+01
            31         32 12121121122122112 21211212212 9.174312e-03 9.174312e-03 1.075380e+01
            32         48 22221211121112112 21222111112 9.174312e-03 9.174312e-03 1.081590e+01
            33          1 11112121121111222 22111112122 9.174312e-03 9.174312e-03 1.087799e+01
            34          8 11212111121111212 22112111112 1.834862e-02 1.834862e-02 1.094008e+01
            35          5 11211111121111222 22112111122 9.174312e-03 9.174312e-03 1.104592e+01
            36         12 11212121121111212 22112112112 9.174312e-03 9.174312e-03 1.110802e+01
            37         13 11212121121111222 22112112122 3.669725e-02 3.669725e-02 1.117011e+01
            38         16 11212121122111222 22112212122 0.000000e+00 0.000000e+00 1.134509e+01
            39         24 11221221122212212 22112222112 9.174312e-03 9.174312e-03 1.134509e+01
            40         43 21221121122112212 22122212112 1.834862e-02 1.834862e-02 1.140718e+01
            41         28 12121121121111212 22211112112 9.174312e-03 9.174312e-03 1.151302e+01
            42         30 12121121121112222 22211112122 2.752294e-02 2.752294e-02 1.157511e+01
            43         47 22221211121111212 22222111112 4.587156e-02 4.587156e-02 1.192173e+01
          **************************
          total_haps = 48 len = 17 done = 1 used = 12 n_haps = 44 entropy 4.940465e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 4 
                          ****--*--****-***
             1          3 11121122212112112 111112121122 9.174312e-03 9.174312e-03 0.000000e+00
             2          6 11211112212112112 111122111121 9.174312e-03 9.174312e-03 6.209343e-02
             3         10 11212121112112112 111122121121 1.834862e-02 1.834862e-02 1.241869e-01
             4         39 21121122212112112 111212121122 3.669725e-02 3.669725e-02 2.300251e-01
             5         41 21121122212212112 111212221122 1.834862e-02 1.834862e-02 4.050044e-01
             6         33 12121122212112112 112112121122 0.000000e+00 0.000000e+00 5.108426e-01
             6         34 12121122212112112 112112121122 1.009174e-01 1.009174e-01 5.108426e-01
             7         36 12121122212212112 112112221122 4.587156e-02 4.587156e-02 8.447534e-01
             8         37 12221122212212112 112122221122 9.174312e-03 9.174312e-03 1.048710e+00
             9          4 11121122212122211 121112122112 3.669725e-02 3.669725e-02 1.110804e+00
            10          9 11212121111111222 121121121221 9.174312e-03 9.174312e-03 1.285783e+00
            11         19 11221112212112212 121122111122 9.174312e-03 9.174312e-03 1.347876e+00
            12         20 11221112212122211 121122112112 9.174312e-03 9.174312e-03 1.409970e+00
            13         21 11221122212112212 121122121122 2.752294e-02 2.752294e-02 1.472063e+00
            14          7 11211122212122211 121122122111 9.174312e-03 9.174312e-03 1.614721e+00
            14         22 11221122212122211 121122122112 9.174312e-03 1.834862e-02 1.614721e+00
            15         23 11221122212122212 121122122122 1.834862e-02 1.834862e-02 1.720559e+00
            16         40 21121122212122211 121212122112 1.834862e-02 1.834862e-02 1.826397e+00
            17         44 21221122212122211 121222122112 1.834862e-02 1.834862e-02 1.932235e+00
            18         35 12121122212112222 122112121222 0.000000e+00 0.000000e+00 2.038074e+00
            19         26 12111122212122211 122112122111 9.174312e-03 9.174312e-03 2.038074e+00
            20         27 12121112212222211 122112212112 1.834862e-02 1.834862e-02 2.100167e+00
            21         45 22121122212111212 122212121122 9.174312e-03 9.174312e-03 2.206005e+00
            22         46 22121122212122211 122212122112 3.669725e-02 3.669725e-02 2.268099e+00
            23          2 11112121122112112 211112121121 9.174312e-03 9.174312e-03 2.443078e+00
            24         14 11212121121112111 211121121111 9.174312e-03 9.174312e-03 2.505171e+00
            25         11 11212121121111112 211121121121 9.174312e-03 9.174312e-03 2.567265e+00
            26         15 11212121121112112 211121121121 2.752294e-02 2.752294e-02 2.629358e+00
            26         18 11212221121112112 211121121121 9.174312e-03 3.669725e-02 2.629358e+00
            27         17 11212121122112112 211122121121 9.174312e-02 9.174312e-02 2.804337e+00
            28         38 21121121122112112 211212121122 4.587156e-02 4.587156e-02 3.120508e+00
            29         42 21221121122112112 211222121122 9.174312e-03 9.174312e-03 3.324464e+00
            30         29 12121121121112112 212111121122 9.174312e-03 9.174312e-03 3.386558e+00
            31         25 12111121122112112 212112121121 9.174312e-03 9.174312e-03 3.448651e+00
            31         31 12121121122112112 212112121122 7.339450e-02 8.256881e-02 3.448651e+00
            32         32 12121121122122112 212112122122 9.174312e-03 9.174312e-03 3.745755e+00
            33         48 22221211121112112 212221111122 9.174312e-03 9.174312e-03 3.807849e+00
            34          1 11112121121111222 221111121221 9.174312e-03 9.174312e-03 3.869942e+00
            35          8 11212111121111212 221121111121 1.834862e-02 1.834862e-02 3.932036e+00
            36          5 11211111121111222 221121111221 9.174312e-03 9.174312e-03 4.037874e+00
            37         12 11212121121111212 221121121121 9.174312e-03 9.174312e-03 4.099967e+00
            38         13 11212121121111222 221121121221 3.669725e-02 3.669725e-02 4.162061e+00
            39         16 11212121122111222 221122121221 0.000000e+00 0.000000e+00 4.337040e+00
            40         24 11221221122212212 221122221122 9.174312e-03 9.174312e-03 4.337040e+00
            41         43 21221121122112212 221222121122 1.834862e-02 1.834862e-02 4.399133e+00
            42         28 12121121121111212 222111121122 9.174312e-03 9.174312e-03 4.504972e+00
            43         30 12121121121112222 222111121222 2.752294e-02 2.752294e-02 4.567065e+00
            44         47 22221211121111212 222221111122 4.587156e-02 4.587156e-02 4.913679e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 13 n_haps = 45 entropy 4.973581e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 4 6 
                          ****-**--****-***
             1          3 11121122212112112 1111121211221 9.174312e-03 9.174312e-03 0.000000e+00
             2          6 11211112212112112 1111221111211 9.174312e-03 9.174312e-03 6.209343e-02
             3         10 11212121112112112 1111221211211 1.834862e-02 1.834862e-02 1.241869e-01
             4         39 21121122212112112 1112121211221 3.669725e-02 3.669725e-02 2.300251e-01
             5         41 21121122212212112 1112122211221 1.834862e-02 1.834862e-02 4.050044e-01
             6         33 12121122212112112 1121121211221 0.000000e+00 0.000000e+00 5.108426e-01
             6         34 12121122212112112 1121121211221 1.009174e-01 1.009174e-01 5.108426e-01
             7         36 12121122212212112 1121122211221 4.587156e-02 4.587156e-02 8.447534e-01
             8         37 12221122212212112 1121222211221 9.174312e-03 9.174312e-03 1.048710e+00
             9          4 11121122212122211 1211121221121 3.669725e-02 3.669725e-02 1.110804e+00
            10          9 11212121111111222 1211211212211 9.174312e-03 9.174312e-03 1.285783e+00
            11         19 11221112212112212 1211221111221 9.174312e-03 9.174312e-03 1.347876e+00
            12         20 11221112212122211 1211221121121 9.174312e-03 9.174312e-03 1.409970e+00
            13         21 11221122212112212 1211221211221 2.752294e-02 2.752294e-02 1.472063e+00
            14          7 11211122212122211 1211221221111 9.174312e-03 9.174312e-03 1.614721e+00
            15         22 11221122212122211 1211221221121 9.174312e-03 9.174312e-03 1.676814e+00
            16         23 11221122212122212 1211221221221 1.834862e-02 1.834862e-02 1.738907e+00
            17         40 21121122212122211 1212121221121 1.834862e-02 1.834862e-02 1.844746e+00
            18         44 21221122212122211 1212221221121 1.834862e-02 1.834862e-02 1.950584e+00
            19         35 12121122212112222 1221121212221 0.000000e+00 0.000000e+00 2.056422e+00
            20         26 12111122212122211 1221121221111 9.174312e-03 9.174312e-03 2.056422e+00
            21         27 12121112212222211 1221122121121 1.834862e-02 1.834862e-02 2.118516e+00
            22         45 22121122212111212 1222121211221 9.174312e-03 9.174312e-03 2.224354e+00
            23         46 22121122212122211 1222121221121 3.669725e-02 3.669725e-02 2.286447e+00
            24          2 11112121122112112 2111121211211 9.174312e-03 9.174312e-03 2.461427e+00
            25         14 11212121121112111 2111211211111 9.174312e-03 9.174312e-03 2.523520e+00
            26         11 11212121121111112 2111211211211 9.174312e-03 9.174312e-03 2.585613e+00
            27         15 11212121121112112 2111211211211 2.752294e-02 2.752294e-02 2.647707e+00
            27         18 11212221121112112 2111211211212 9.174312e-03 3.669725e-02 2.647707e+00
            28         17 11212121122112112 2111221211211 9.174312e-02 9.174312e-02 2.822686e+00
            29         38 21121121122112112 2112121211221 4.587156e-02 4.587156e-02 3.138856e+00
            30         42 21221121122112112 2112221211221 9.174312e-03 9.174312e-03 3.342813e+00
            31         29 12121121121112112 2121111211221 9.174312e-03 9.174312e-03 3.404907e+00
            32         25 12111121122112112 2121121211211 9.174312e-03 9.174312e-03 3.467000e+00
            33         31 12121121122112112 2121121211221 7.339450e-02 7.339450e-02 3.529093e+00
            34         32 12121121122122112 2121121221221 9.174312e-03 9.174312e-03 3.805657e+00
            35         48 22221211121112112 2122211111222 9.174312e-03 9.174312e-03 3.867751e+00
            36          1 11112121121111222 2211111212211 9.174312e-03 9.174312e-03 3.929844e+00
            37          8 11212111121111212 2211211111211 1.834862e-02 1.834862e-02 3.991938e+00
            38          5 11211111121111222 2211211112211 9.174312e-03 9.174312e-03 4.097776e+00
            39         12 11212121121111212 2211211211211 9.174312e-03 9.174312e-03 4.159869e+00
            40         13 11212121121111222 2211211212211 3.669725e-02 3.669725e-02 4.221963e+00
            41         16 11212121122111222 2211221212211 0.000000e+00 0.000000e+00 4.396942e+00
            42         24 11221221122212212 2211222211222 9.174312e-03 9.174312e-03 4.396942e+00
            43         43 21221121122112212 2212221211221 1.834862e-02 1.834862e-02 4.459035e+00
            44         28 12121121121111212 2221111211221 9.174312e-03 9.174312e-03 4.564874e+00
            45         30 12121121121112222 2221111212221 2.752294e-02 2.752294e-02 4.626967e+00
            46         47 22221211121111212 2222211111222 4.587156e-02 4.587156e-02 4.973581e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 14 n_haps = 46 entropy 5.003353e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 4 6 14 
                          ****-**--********
             1          3 11121122212112112 11111212112212 9.174312e-03 9.174312e-03 0.000000e+00
             2          6 11211112212112112 11112211112112 9.174312e-03 9.174312e-03 6.209343e-02
             3         10 11212121112112112 11112212112112 1.834862e-02 1.834862e-02 1.241869e-01
             4         39 21121122212112112 11121212112212 3.669725e-02 3.669725e-02 2.300251e-01
             5         41 21121122212212112 11121222112212 1.834862e-02 1.834862e-02 4.050044e-01
             6         33 12121122212112112 11211212112212 0.000000e+00 0.000000e+00 5.108426e-01
             6         34 12121122212112112 11211212112212 1.009174e-01 1.009174e-01 5.108426e-01
             7         36 12121122212212112 11211222112212 4.587156e-02 4.587156e-02 8.447534e-01
             8         37 12221122212212112 11212222112212 9.174312e-03 9.174312e-03 1.048710e+00
             9          4 11121122212122211 12111212211212 3.669725e-02 3.669725e-02 1.110804e+00
            10          9 11212121111111222 12112112122111 9.174312e-03 9.174312e-03 1.285783e+00
            11         19 11221112212112212 12112211112212 9.174312e-03 9.174312e-03 1.347876e+00
            12         20 11221112212122211 12112211211212 9.174312e-03 9.174312e-03 1.409970e+00
            13         21 11221122212112212 12112212112212 2.752294e-02 2.752294e-02 1.472063e+00
            14          7 11211122212122211 12112212211112 9.174312e-03 9.174312e-03 1.614721e+00
            15         22 11221122212122211 12112212211212 9.174312e-03 9.174312e-03 1.676814e+00
            16         23 11221122212122212 12112212212212 1.834862e-02 1.834862e-02 1.738907e+00
            17         40 21121122212122211 12121212211212 1.834862e-02 1.834862e-02 1.844746e+00
            18         44 21221122212122211 12122212211212 1.834862e-02 1.834862e-02 1.950584e+00
            19         35 12121122212112222 12211212122212 0.000000e+00 0.000000e+00 2.056422e+00
            20         26 12111122212122211 12211212211112 9.174312e-03 9.174312e-03 2.056422e+00
            21         27 12121112212222211 12211221211212 1.834862e-02 1.834862e-02 2.118516e+00
            22         45 22121122212111212 12221212112211 9.174312e-03 9.174312e-03 2.224354e+00
            23         46 22121122212122211 12221212211212 3.669725e-02 3.669725e-02 2.286447e+00
            24          2 11112121122112112 21111212112112 9.174312e-03 9.174312e-03 2.461427e+00
            25         14 11212121121112111 21112112111112 9.174312e-03 9.174312e-03 2.523520e+00
            26         11 11212121121111112 21112112112111 9.174312e-03 9.174312e-03 2.585613e+00
            27         15 11212121121112112 21112112112112 2.752294e-02 2.752294e-02 2.647707e+00
            28         18 11212221121112112 21112112112122 9.174312e-03 9.174312e-03 2.790364e+00
            29         17 11212121122112112 21112212112112 9.174312e-02 9.174312e-02 2.852458e+00
            30         38 21121121122112112 21121212112212 4.587156e-02 4.587156e-02 3.168628e+00
            31         42 21221121122112112 21122212112212 9.174312e-03 9.174312e-03 3.372585e+00
            32         29 12121121121112112 21211112112212 9.174312e-03 9.174312e-03 3.434678e+00
            33         25 12111121122112112 21211212112112 9.174312e-03 9.174312e-03 3.496772e+00
            34         31 12121121122112112 21211212112212 7.339450e-02 7.339450e-02 3.558865e+00
            35         32 12121121122122112 21211212212212 9.174312e-03 9.174312e-03 3.835429e+00
            36         48 22221211121112112 21222111112222 9.174312e-03 9.174312e-03 3.897522e+00
            37          1 11112121121111222 22111112122111 9.174312e-03 9.174312e-03 3.959616e+00
            38          8 11212111121111212 22112111112111 1.834862e-02 1.834862e-02 4.021709e+00
            39          5 11211111121111222 22112111122111 9.174312e-03 9.174312e-03 4.127548e+00
            40         12 11212121121111212 22112112112111 9.174312e-03 9.174312e-03 4.189641e+00
            41         13 11212121121111222 22112112122111 3.669725e-02 3.669725e-02 4.251734e+00
            42         16 11212121122111222 22112212122111 0.000000e+00 0.000000e+00 4.426714e+00
            43         24 11221221122212212 22112222112222 9.174312e-03 9.174312e-03 4.426714e+00
            44         43 21221121122112212 22122212112212 1.834862e-02 1.834862e-02 4.488807e+00
            45         28 12121121121111212 22211112112211 9.174312e-03 9.174312e-03 4.594645e+00
            46         30 12121121121112222 22211112122212 2.752294e-02 2.752294e-02 4.656739e+00
            47         47 22221211121111212 22222111112221 4.587156e-02 4.587156e-02 5.003353e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 15 n_haps = 46 entropy 5.003353e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 4 6 14 5 
                          *******--********
             1          3 11121122212112112 111112121122121 9.174312e-03 9.174312e-03 0.000000e+00
             2          6 11211112212112112 111122111121121 9.174312e-03 9.174312e-03 6.209343e-02
             3         10 11212121112112112 111122121121122 1.834862e-02 1.834862e-02 1.241869e-01
             4         39 21121122212112112 111212121122121 3.669725e-02 3.669725e-02 2.300251e-01
             5         41 21121122212212112 111212221122121 1.834862e-02 1.834862e-02 4.050044e-01
             6         33 12121122212112112 112112121122121 0.000000e+00 0.000000e+00 5.108426e-01
             6         34 12121122212112112 112112121122121 1.009174e-01 1.009174e-01 5.108426e-01
             7         36 12121122212212112 112112221122121 4.587156e-02 4.587156e-02 8.447534e-01
             8         37 12221122212212112 112122221122121 9.174312e-03 9.174312e-03 1.048710e+00
             9          4 11121122212122211 121112122112121 3.669725e-02 3.669725e-02 1.110804e+00
            10          9 11212121111111222 121121121221112 9.174312e-03 9.174312e-03 1.285783e+00
            11         19 11221112212112212 121122111122121 9.174312e-03 9.174312e-03 1.347876e+00
            12         20 11221112212122211 121122112112121 9.174312e-03 9.174312e-03 1.409970e+00
            13         21 11221122212112212 121122121122121 2.752294e-02 2.752294e-02 1.472063e+00
            14          7 11211122212122211 121122122111121 9.174312e-03 9.174312e-03 1.614721e+00
            15         22 11221122212122211 121122122112121 9.174312e-03 9.174312e-03 1.676814e+00
            16         23 11221122212122212 121122122122121 1.834862e-02 1.834862e-02 1.738907e+00
            17         40 21121122212122211 121212122112121 1.834862e-02 1.834862e-02 1.844746e+00
            18         44 21221122212122211 121222122112121 1.834862e-02 1.834862e-02 1.950584e+00
            19         35 12121122212112222 122112121222121 0.000000e+00 0.000000e+00 2.056422e+00
            20         26 12111122212122211 122112122111121 9.174312e-03 9.174312e-03 2.056422e+00
            21         27 12121112212222211 122112212112121 1.834862e-02 1.834862e-02 2.118516e+00
            22         45 22121122212111212 122212121122111 9.174312e-03 9.174312e-03 2.224354e+00
            23         46 22121122212122211 122212122112121 3.669725e-02 3.669725e-02 2.286447e+00
            24          2 11112121122112112 211112121121122 9.174312e-03 9.174312e-03 2.461427e+00
            25         14 11212121121112111 211121121111122 9.174312e-03 9.174312e-03 2.523520e+00
            26         11 11212121121111112 211121121121112 9.174312e-03 9.174312e-03 2.585613e+00
            27         15 11212121121112112 211121121121122 2.752294e-02 2.752294e-02 2.647707e+00
            28         18 11212221121112112 211121121121222 9.174312e-03 9.174312e-03 2.790364e+00
            29         17 11212121122112112 211122121121122 9.174312e-02 9.174312e-02 2.852458e+00
            30         38 21121121122112112 211212121122121 4.587156e-02 4.587156e-02 3.168628e+00
            31         42 21221121122112112 211222121122121 9.174312e-03 9.174312e-03 3.372585e+00
            32         29 12121121121112112 212111121122121 9.174312e-03 9.174312e-03 3.434678e+00
            33         25 12111121122112112 212112121121121 9.174312e-03 9.174312e-03 3.496772e+00
            34         31 12121121122112112 212112121122121 7.339450e-02 7.339450e-02 3.558865e+00
            35         32 12121121122122112 212112122122121 9.174312e-03 9.174312e-03 3.835429e+00
            36         48 22221211121112112 212221111122221 9.174312e-03 9.174312e-03 3.897522e+00
            37          1 11112121121111222 221111121221112 9.174312e-03 9.174312e-03 3.959616e+00
            38          8 11212111121111212 221121111121112 1.834862e-02 1.834862e-02 4.021709e+00
            39          5 11211111121111222 221121111221111 9.174312e-03 9.174312e-03 4.127548e+00
            40         12 11212121121111212 221121121121112 9.174312e-03 9.174312e-03 4.189641e+00
            41         13 11212121121111222 221121121221112 3.669725e-02 3.669725e-02 4.251734e+00
            42         16 11212121122111222 221122121221112 0.000000e+00 0.000000e+00 4.426714e+00
            43         24 11221221122212212 221122221122221 9.174312e-03 9.174312e-03 4.426714e+00
            44         43 21221121122112212 221222121122121 1.834862e-02 1.834862e-02 4.488807e+00
            45         28 12121121121111212 222111121122111 9.174312e-03 9.174312e-03 4.594645e+00
            46         30 12121121121112222 222111121222121 2.752294e-02 2.752294e-02 4.656739e+00
            47         47 22221211121111212 222221111122211 4.587156e-02 4.587156e-02 5.003353e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 16 n_haps = 46 entropy 5.003353e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 4 6 14 5 8 
                          ********-********
             1          3 11121122212112112 1111121211221212 9.174312e-03 9.174312e-03 0.000000e+00
             2          6 11211112212112112 1111221111211212 9.174312e-03 9.174312e-03 6.209343e-02
             3         10 11212121112112112 1111221211211221 1.834862e-02 1.834862e-02 1.241869e-01
             4         39 21121122212112112 1112121211221212 3.669725e-02 3.669725e-02 2.300251e-01
             5         41 21121122212212112 1112122211221212 1.834862e-02 1.834862e-02 4.050044e-01
             6         33 12121122212112112 1121121211221212 0.000000e+00 0.000000e+00 5.108426e-01
             6         34 12121122212112112 1121121211221212 1.009174e-01 1.009174e-01 5.108426e-01
             7         36 12121122212212112 1121122211221212 4.587156e-02 4.587156e-02 8.447534e-01
             8         37 12221122212212112 1121222211221212 9.174312e-03 9.174312e-03 1.048710e+00
             9          4 11121122212122211 1211121221121212 3.669725e-02 3.669725e-02 1.110804e+00
            10          9 11212121111111222 1211211212211121 9.174312e-03 9.174312e-03 1.285783e+00
            11         19 11221112212112212 1211221111221212 9.174312e-03 9.174312e-03 1.347876e+00
            12         20 11221112212122211 1211221121121212 9.174312e-03 9.174312e-03 1.409970e+00
            13         21 11221122212112212 1211221211221212 2.752294e-02 2.752294e-02 1.472063e+00
            14          7 11211122212122211 1211221221111212 9.174312e-03 9.174312e-03 1.614721e+00
            15         22 11221122212122211 1211221221121212 9.174312e-03 9.174312e-03 1.676814e+00
            16         23 11221122212122212 1211221221221212 1.834862e-02 1.834862e-02 1.738907e+00
            17         40 21121122212122211 1212121221121212 1.834862e-02 1.834862e-02 1.844746e+00
            18         44 21221122212122211 1212221221121212 1.834862e-02 1.834862e-02 1.950584e+00
            19         35 12121122212112222 1221121212221212 0.000000e+00 0.000000e+00 2.056422e+00
            20         26 12111122212122211 1221121221111212 9.174312e-03 9.174312e-03 2.056422e+00
            21         27 12121112212222211 1221122121121212 1.834862e-02 1.834862e-02 2.118516e+00
            22         45 22121122212111212 1222121211221112 9.174312e-03 9.174312e-03 2.224354e+00
            23         46 22121122212122211 1222121221121212 3.669725e-02 3.669725e-02 2.286447e+00
            24          2 11112121122112112 2111121211211221 9.174312e-03 9.174312e-03 2.461427e+00
            25         14 11212121121112111 2111211211111221 9.174312e-03 9.174312e-03 2.523520e+00
            26         11 11212121121111112 2111211211211121 9.174312e-03 9.174312e-03 2.585613e+00
            27         15 11212121121112112 2111211211211221 2.752294e-02 2.752294e-02 2.647707e+00
            28         18 11212221121112112 2111211211212221 9.174312e-03 9.174312e-03 2.790364e+00
            29         17 11212121122112112 2111221211211221 9.174312e-02 9.174312e-02 2.852458e+00
            30         38 21121121122112112 2112121211221211 4.587156e-02 4.587156e-02 3.168628e+00
            31         42 21221121122112112 2112221211221211 9.174312e-03 9.174312e-03 3.372585e+00
            32         29 12121121121112112 2121111211221211 9.174312e-03 9.174312e-03 3.434678e+00
            33         25 12111121122112112 2121121211211211 9.174312e-03 9.174312e-03 3.496772e+00
            34         31 12121121122112112 2121121211221211 7.339450e-02 7.339450e-02 3.558865e+00
            35         32 12121121122122112 2121121221221211 9.174312e-03 9.174312e-03 3.835429e+00
            36         48 22221211121112112 2122211111222211 9.174312e-03 9.174312e-03 3.897522e+00
            37          1 11112121121111222 2211111212211121 9.174312e-03 9.174312e-03 3.959616e+00
            38          8 11212111121111212 2211211111211121 1.834862e-02 1.834862e-02 4.021709e+00
            39          5 11211111121111222 2211211112211111 9.174312e-03 9.174312e-03 4.127548e+00
            40         12 11212121121111212 2211211211211121 9.174312e-03 9.174312e-03 4.189641e+00
            41         13 11212121121111222 2211211212211121 3.669725e-02 3.669725e-02 4.251734e+00
            42         16 11212121122111222 2211221212211121 0.000000e+00 0.000000e+00 4.426714e+00
            43         24 11221221122212212 2211222211222211 9.174312e-03 9.174312e-03 4.426714e+00
            44         43 21221121122112212 2212221211221211 1.834862e-02 1.834862e-02 4.488807e+00
            45         28 12121121121111212 2221111211221111 9.174312e-03 9.174312e-03 4.594645e+00
            46         30 12121121121112222 2221111212221211 2.752294e-02 2.752294e-02 4.656739e+00
            47         47 22221211121111212 2222211111222111 4.587156e-02 4.587156e-02 5.003353e+00
          **************************
          total_haps = 48 len = 17 done = 1 used = 17 n_haps = 46 entropy 5.003353e+00 5.003353e+00 SNPS used: 10 15 2 1 3 11 12 7 13 16 17 4 6 14 5 8 9 
                          *****************
             1          3 11121122212112112 11111212112212122 9.174312e-03 9.174312e-03 0.000000e+00
             2          6 11211112212112112 11112211112112122 9.174312e-03 9.174312e-03 6.209343e-02
             3         10 11212121112112112 11112212112112211 1.834862e-02 1.834862e-02 1.241869e-01
             4         39 21121122212112112 11121212112212122 3.669725e-02 3.669725e-02 2.300251e-01
             5         41 21121122212212112 11121222112212122 1.834862e-02 1.834862e-02 4.050044e-01
             6         33 12121122212112112 11211212112212122 0.000000e+00 0.000000e+00 5.108426e-01
             6         34 12121122212112112 11211212112212122 1.009174e-01 1.009174e-01 5.108426e-01
             7         36 12121122212212112 11211222112212122 4.587156e-02 4.587156e-02 8.447534e-01
             8         37 12221122212212112 11212222112212122 9.174312e-03 9.174312e-03 1.048710e+00
             9          4 11121122212122211 12111212211212122 3.669725e-02 3.669725e-02 1.110804e+00
            10          9 11212121111111222 12112112122111211 9.174312e-03 9.174312e-03 1.285783e+00
            11         19 11221112212112212 12112211112212122 9.174312e-03 9.174312e-03 1.347876e+00
            12         20 11221112212122211 12112211211212122 9.174312e-03 9.174312e-03 1.409970e+00
            13         21 11221122212112212 12112212112212122 2.752294e-02 2.752294e-02 1.472063e+00
            14          7 11211122212122211 12112212211112122 9.174312e-03 9.174312e-03 1.614721e+00
            15         22 11221122212122211 12112212211212122 9.174312e-03 9.174312e-03 1.676814e+00
            16         23 11221122212122212 12112212212212122 1.834862e-02 1.834862e-02 1.738907e+00
            17         40 21121122212122211 12121212211212122 1.834862e-02 1.834862e-02 1.844746e+00
            18         44 21221122212122211 12122212211212122 1.834862e-02 1.834862e-02 1.950584e+00
            19         35 12121122212112222 12211212122212122 0.000000e+00 0.000000e+00 2.056422e+00
            20         26 12111122212122211 12211212211112122 9.174312e-03 9.174312e-03 2.056422e+00
            21         27 12121112212222211 12211221211212122 1.834862e-02 1.834862e-02 2.118516e+00
            22         45 22121122212111212 12221212112211122 9.174312e-03 9.174312e-03 2.224354e+00
            23         46 22121122212122211 12221212211212122 3.669725e-02 3.669725e-02 2.286447e+00
            24          2 11112121122112112 21111212112112211 9.174312e-03 9.174312e-03 2.461427e+00
            25         14 11212121121112111 21112112111112211 9.174312e-03 9.174312e-03 2.523520e+00
            26         11 11212121121111112 21112112112111211 9.174312e-03 9.174312e-03 2.585613e+00
            27         15 11212121121112112 21112112112112211 2.752294e-02 2.752294e-02 2.647707e+00
            28         18 11212221121112112 21112112112122211 9.174312e-03 9.174312e-03 2.790364e+00
            29         17 11212121122112112 21112212112112211 9.174312e-02 9.174312e-02 2.852458e+00
            30         38 21121121122112112 21121212112212111 4.587156e-02 4.587156e-02 3.168628e+00
            31         42 21221121122112112 21122212112212111 9.174312e-03 9.174312e-03 3.372585e+00
            32         29 12121121121112112 21211112112212111 9.174312e-03 9.174312e-03 3.434678e+00
            33         25 12111121122112112 21211212112112111 9.174312e-03 9.174312e-03 3.496772e+00
            34         31 12121121122112112 21211212112212111 7.339450e-02 7.339450e-02 3.558865e+00
            35         32 12121121122122112 21211212212212111 9.174312e-03 9.174312e-03 3.835429e+00
            36         48 22221211121112112 21222111112222111 9.174312e-03 9.174312e-03 3.897522e+00
            37          1 11112121121111222 22111112122111211 9.174312e-03 9.174312e-03 3.959616e+00
            38          8 11212111121111212 22112111112111211 1.834862e-02 1.834862e-02 4.021709e+00
            39          5 11211111121111222 22112111122111111 9.174312e-03 9.174312e-03 4.127548e+00
            40         12 11212121121111212 22112112112111211 9.174312e-03 9.174312e-03 4.189641e+00
            41         13 11212121121111222 22112112122111211 3.669725e-02 3.669725e-02 4.251734e+00
            42         16 11212121122111222 22112212122111211 0.000000e+00 0.000000e+00 4.426714e+00
            43         24 11221221122212212 22112222112222111 9.174312e-03 9.174312e-03 4.426714e+00
            44         43 21221121122112212 22122212112212111 1.834862e-02 1.834862e-02 4.488807e+00
            45         28 12121121121111212 22211112112211111 9.174312e-03 9.174312e-03 4.594645e+00
            46         30 12121121121112222 22211112122212111 2.752294e-02 2.752294e-02 4.656739e+00
            47         47 22221211121111212 22222111112221111 4.587156e-02 4.587156e-02 5.003353e+00
                  
          

big_haplotype

produces a text output file containing three main sections:

       # Reading haplotypes from test.haps.txt
       # min acceptable haplotype freq 0.000000e+00
       # 80 haplotypes, 899 markers read from test.haps.txt
       
       # SNP search in each window of length 20
  1. The first (very long) section describes the SNPs selected in each window. For example, with a window-size of 20, the first two windows (with SNP-ID 0 to 19 and 1 to 20) might produce output like this:
           window 0 to 19 entropy 5.478379e+00
           index   SNP-ID      delta-E        cum-E   fraction-E       
               0        3 9.995491e-01 9.995491e-01 1.824534e-01       
               1       10 9.977877e-01 1.997337e+00 3.645854e-01       
               2        8 8.968164e-01 2.894153e+00 5.282864e-01       
               3        0 7.233388e-01 3.617492e+00 6.603216e-01       
               4        4 6.635405e-01 4.281032e+00 7.814415e-01       
               5       19 4.579576e-01 4.738990e+00 8.650351e-01       
               6       15 2.256496e-01 4.964640e+00 9.062242e-01       
               7       16 1.987517e-01 5.163391e+00 9.425035e-01       
               8        7 1.545566e-01 5.317948e+00 9.707156e-01       
               9       17 6.043084e-02 5.378379e+00 9.817464e-01       
              10        9 5.000000e-02 5.428379e+00 9.908732e-01       
              11       13 5.000000e-02 5.478379e+00 1.000000e-00       
              12        1 0.000000e+00 5.478379e+00 1.000000e-00       
              13        2 0.000000e+00 5.478379e+00 1.000000e-00       
              14        5 0.000000e+00 5.478379e+00 1.000000e-00       
              15        6 0.000000e+00 5.478379e+00 1.000000e-00       
              16       11 0.000000e+00 5.478379e+00 1.000000e-00       
              17       12 0.000000e+00 5.478379e+00 1.000000e-00       
              18       14 0.000000e+00 5.478379e+00 1.000000e-00       
              19       18 0.000000e+00 5.478379e+00 1.000000e-00       
    
    
    • The total entropy for all 20 markers in this window is 5.4739.
    • The markers are added one at a time using a greedy algorithm that miximises delta-E, the increase in entropy. the ids of the first few markers added are 3, 10, 8, 0,...
    • cum-E is the cumulative total entropy.
    • fraction-E is the ratio cum-E/entropy
    • We see that only 10 markers are required to explain the haplotypic diversity in this window
    
           window 1 to 20 entropy 5.359184e+00
           index   SNP-ID      delta-E        cum-E   fraction-E
               0        3 9.995491e-01 9.995491e-01 1.865114e-01              
               1       10 9.977877e-01 1.997337e+00 3.726942e-01
               2       20 9.183086e-01 2.915645e+00 5.440466e-01       
               3        8 8.468226e-01 3.762468e+00 7.020599e-01       
               4       19 5.206436e-01 4.283112e+00 7.992097e-01       
               5        4 4.114076e-01 4.694519e+00 8.759765e-01       
               6       16 2.606844e-01 4.955204e+00 9.246191e-01       
               7       15 1.784914e-01 5.133695e+00 9.579248e-01       
               8        7 1.438722e-01 5.277567e+00 9.847707e-01       
               9       17 5.661656e-02 5.334184e+00 9.953351e-01       
              10        9 2.500000e-02 5.359184e+00 1.000000e-00       
              11        1 0.000000e+00 5.359184e+00 1.000000e-00       
              12        2 0.000000e+00 5.359184e+00 1.000000e-00       
              13        5 0.000000e+00 5.359184e+00 1.000000e-00       
              14        6 0.000000e+00 5.359184e+00 1.000000e-00       
              15       11 0.000000e+00 5.359184e+00 1.000000e-00       
              16       12 0.000000e+00 5.359184e+00 1.000000e-00       
              17       13 0.000000e+00 5.359184e+00 1.000000e-00       
              18       14 0.000000e+00 5.359184e+00 1.000000e-00       
              19       18 0.000000e+00 5.359184e+00 1.000000e-00       
           
    
    • For the second window (markers 1..20) the markers are added in the order 3,10,20,8,...
    • Note that marker 0 is outside of the window and so is no longer available
  2. The second section lists all the selected markers with their scores S:
    
           #score threshold for identifying top 449 SNPs = 1.255623e-01
           # Selected SNPs:
           
             SNP-ID        score
                  0 7.233388e-01
                  3 9.995491e-01
                  4 4.390538e-01
                  8 8.798721e-01
                 10 9.973898e-01
                 16 1.858399e-01
                 18 1.255623e-01
                 19 6.669792e-01
                 20 6.232454e-01
                 21 2.354860e-01
                 .....
                 .....
    
    
  3. The final section lists summary statistics for each window (plotted in this figure)
           #Evalution of SNP Selection
             window    #SNPs      Total-E    Partial-E Fractional-E
                  0        8 5.076319e+00 5.478379e+00 9.266097e-01
                  1        8 5.075071e+00 5.359184e+00 9.469857e-01
                  2        9 5.268943e+00 5.493620e+00 9.591022e-01
                  3       10 5.678056e+00 5.842371e+00 9.718752e-01
                  4        9 5.440312e+00 5.731687e+00 9.491641e-01
                  5        8 5.043420e+00 5.643943e+00 8.935987e-01
                  6        9 5.496439e+00 5.848499e+00 9.398034e-01
    
    • window is the id of the first (left-most) SNP in the window
    • #SNPs is the number of SNPs in the window which are selected
    • Total-E is the total entropy for the window
    • Partial-E is the entropy explained by just the selected SNPs in the window
    • Fractional-E is the fraction of the total entropy explained in the window by the selected SNPs

Drop me an email for more details.

 
 
spacer