
seed used = 140026221

K12                                  ATG AAT AAA ATC CTG TTA GTT GAT GAT GAC CGA GAG CTG ACT TCC CTA TTA AAG GAG CTG CTC GAG ATG GAA GGC TTC AAC GTG ATT GTT GCC CAC GAT GGG GAA CAG GCG CTT GAT CTT CTG GAC GAC AGC ATT GAT TTA CTT TTG CTT GAC GTA ATG ATG CCG AAG AAA AAT GGT ATC GAC ACA TTA AAA GCA CTT CGC CAG ACA CAC CAG ACG CCT GTC ATT ATG TTG ACG GCG CGC GGC AGT GAA CTT GAT CGC GTT CTC GGC CTT GAG CTG GGC GCA GAT GAC TAT CTC CCG AAA CCG TTT AAT GAT CGT GAG CTG GTG GCA CGT ATT CGC GCG ATC CTG CGC CGT TCG CAC TGG AGC GAG CAA CAG CAA AAC AAC GAC AAC GGT TCA CCG ACA CTG GAA GTT GAT GCC TTA GTG CTG AAT CCA GGC CGT CAG GAA GCC AGC TTC GAC GGG CAA ACG CTG GAG TTA ACC GGT ACT GAG TTT ACC CTG CTC TAT TTG CTG GCA CAG CAT CTG GGT CAG GTG GTT TCC CGT GAA CAT TTA AGC CAG GAA GTG TTG GGC AAA CGC CTG ACG CCT TTC GAC CGC GCT ATT GAT ATG CAC ATT TCC AAC CTG CGT CGT AAA CTG CCG GAT CGT AAA GAT GGT CAC CCG TGG TTT AAA ACC TTG CGT GGT CGC GGC TAT CTG ATG GTT TCT GCT TCA 
S2457T                               ATG AAT AAA ATC CTG TTA GTT GAT GAT GAC CGA GAG CTG ACT TCC CTA TTA AAG GAG CTG CTC GAG ATG GAA GGC TTC AAC GTG ATT GTT GCC CAC GAT GGG GAA CAG GCG CTT GAT CTT CTG GAC GAC AGC ATT GAT TTA CTT TTG CTT GAC GTA ATG ATG CCG AAG AAA AAT GGT ATC GAC ACA TTA AAA GCA CTT CGC CAG ACA CAC CAG ACG CCT GTC ATT ATG TTG ACG GCG CGC GGC AGT GAA CTT GAT CGC GTT CTC GGC CTT GAG CTG GGC GCA GAT GAC TAT CTC CCG AAA CCG TTT AAT GAT CGT GAG CTG GTG GCA CGT ATT CGC GCA ATC CTG CGC CGT TCG CAC TGG AGC GAG CAA CAG CAA AAC AAC GAC AAC GGT TCA CCG ACA CTG GAA GTT GAT GCC TTA GTG CTG AAT CCA GGC CGT CAG GAA GCC AGC TTC GAC GGG CAA ACG CTG GAG TTA ACC GGT ACT GAG TTT ACC CTG CTC TAT TTG CTG GCA CAG CAT CTG GGT CAG GTG GTT TCC CGT GAA CAT TTA AGC CAG GAA GTG TTG GGC AAA CGC CTG ACG CCT TTC GAC CGC GCT ATT GAT ATG CAC ATT TCC AAC CTG CGT CGT AAA CTG CCG GAT CGT AAA GAT GGT CAC CCG TGG TTT AAA ACC TTG CGT GGT CGC GGC TAT CTG ATG GTT TCT GCT TCA 
Shig4337                             ATG AAT AAA ATC CTG TTA GTT GAT GAT GAC CGA GAG CTG ACT TCC CTA TTA AAG GAG CTG CTC GAG ATG GAA GGC TTC AAC GTG ATT GTT GCC CAC GAT GGG GAA CAG GCG CTT GAT CTT CTG GAC GAC AGC ATT GAT TTA CTT TTG CTT GAC GTA ATG ATG CCG AAG AAA AAT GGT ATC GAC ACA TTA AAA GCA CTT CGC CAG ACA CAC CAG ACG CCT GTC ATT ATG TTG ACG GCG CGC GGC AGT GAA CTT GAT CGC GTT CTC GGC CTT GAG CTG GGC GCA GAT GAC TAT CTC CCG AAA CCG TTT AAT GAT CGT GAG CTG GTG GCA CGT ATT CGC GCA ATC CTG CGC CGT TCG CAC TGG AGC GAG CAA CAG CAA AAC AAC GAC AAC GGT TCA CCG ACA CTG GAA GTT GAT GCC TTA GTG CTG AAT CCA GGC CGT CAG GAA GCC AGC TTC GAC GGG CAA ACG CTG GAG TTA ACC GGT ACT GAG TTT ACC CTG CTC TAT TTG CTG GCA CAG CAT CTG GGT CAG GTG GTT TCC CGT GAA CAT TTA AGC CAG GAA GTG TTG GGC AAA CGC CTG ACG CCT TTC GAC CGC GCT ATT GAT ATG CAC ATT TCC AAC CTG CGT CGT AAA CTG CCG GAT CGT AAA GAT GGT CAC CCG TGG TTT AAA ACC TTG CGT GGT CGC GGC TAT CTG ATG GTT TCT GCT TCA 
CFT073                               ATG AAT AAA ATC CTG TTA GTT GAT GAT GAC CGA GAG CTG ACT TCC CTA TTA AAG GAG CTG CTC GAG ATG GAA GGC TTC AAC GTG ATT GTT GCC CAC GAT GGG GAA CAG GCG CTT GAT CTT CTG GAC GAC AGC ATT GAT TTA CTT TTG CTT GAT GTA ATG ATG CCG AAG AAA AAT GGT ATC GAC ACA TTA AAA GCA CTT CGC CAG ACA CAC CAG ACG CCT GTC ATT ATG TTG ACG GCG CGC GGC AGT GAA CTT GAT CGC GTT CTC GGG CTT GAG CTG GGC GCA GAT GAC TAT CTC CCG AAA CCG TTT AAT GAT CGT GAG CTG GTG GCA CGT ATT CGC GCG ATC CTG CGC CGT TCG CAC TGG AGC GAA CAA CAG CAA AAC AAC GAC AAC GGT TCA CCG ACA CTG GAA GTT GAT GCC TTA AGT CTC AAT CCT GGT CGC CAG GAA GCC AGC TTC GAC GGG CAA ACG CTG GAG TTA ACC GGT ACT GAG TTT ACC CTG CTC TAT TTG CTG GCA CAG CAT CTG GGT CAG GTG GTT TCC CGT GAG CAT TTA AGC CAG GAA GTG CTG GGC AAG CGC CTG ACG CCT TTC GAC CGC GCT ATC GAT ATG CAC ATT TCC AAC CTG CGT CGT AAA CTG CCG GAT CGT AAA GAT GGT CAC CCG TGG TTT AAA ACC CTG CGT GGT CGC GGC TAT CTG ATG GTT TCT GCT TCA 
O157                                 ATG AAT AAA ATC CTG TTA GTT GAT GAT GAC CGA GAG CTG ACT TCC CTA TTA AAG GAG CTG CTC GAG ATG GAA GGC TTC AAC GTG ATT GTT GCC CAC GAT GGG GAA CAG GCG CTT GAT CTT CTG GAC GAC AGC ATT GAT TTA CTT TTG CTT GAC GTA ATG ATG CCG AAG AAA AAT GGT ATC GAC ACA TTA AAA GCA CTT CGC CAG ACA CAC CAG ACG CCT GTC ATT ATG TTG ACG GCG CGC GGC AGT GAA CTT GAT CGC GTT CTC GGC CTT GAG CTG GGC GCA GAT GAC TAT CTC CCG AAA CCG TTT AAT GAT CGT GAG CTG GTG GCA CGT ATT CGC GCG ATC CTG CGC CGT TCG CAC TGG AGC GAG CAA CAG CAA AAC AAC GAC AAC GGT TCA CCG ACA CTG GAA GTT GAT GCC TTA GTA CTG AAT CCA GGC CGT CAG GAA GCC AGC TTC GAC GGG CAA ACG CTG GAG TTA ACC GGT ACT GAG TTT ACC CTG CTC TAT TTG CTG GCA CAG CAT CTG GGT CAG GTG GTT TCC CGT GAA CAT TTA AGC CAG GAA GTG CTG GGC AAA CGC CTG ACG CCT TTC GAC CGC GCT ATC GAT ATG CAC ATT TCC AAC CTG CGT CGT AAA CTG CCG GAT CGT AAA GAT GGT CAC CCG TGG TTT AAA ACC TTG CGT GGT CGC GGC TAC CTG ATG GTT TCT GCT TCA 
EDL933                               ATG AAT AAA ATC CTG TTA GTT GAT GAT GAC CGA GAG CTG ACT TCC CTA TTA AAG GAG CTG CTC GAG ATG GAA GGC TTC AAC GTG ATT GTT GCC CAC GAT GGG GAA CAG GCG CTT GAT CTT CTG GAC GAC AGC ATT GAT TTA CTT TTG CTT GAC GTA ATG ATG CCG AAG AAA AAT GGT ATC GAC ACA TTA AAA GCA CTT CGC CAG ACA CAC CAG ACG CCT GTC ATT ATG TTG ACG GCG CGC GGC AGT GAA CTT GAT CGC GTT CTC GGC CTT GAG CTG GGC GCA GAT GAC TAT CTC CCG AAA CCG TTT AAT GAT CGT GAG CTG GTG GCA CGT ATT CGC GCG ATC CTG CGC CGT TCG CAC TGG AGC GAG CAA CAG CAA AAC AAC GAC AAC GGT TCA CCG ACA CTG GAA GTT GAT GCC TTA GTA CTG AAT CCA GGC CGT CAG GAA GCC AGC TTC GAC GGG CAA ACG CTG GAG TTA ACC GGT ACT GAG TTT ACC CTG CTC TAT TTG CTG GCA CAG CAT CTG GGT CAG GTG GTT TCC CGT GAA CAT TTA AGC CAG GAA GTG CTG GGC AAA CGC CTG ACG CCT TTC GAC CGC GCT ATC GAT ATG CAC ATT TCC AAC CTG CGT CGT AAA CTG CCG GAT CGT AAA GAT GGT CAC CCG TGG TTT AAA ACC TTG CGT GGT CGC GGC TAC CTG ATG GTT TCT GCT TCA 


CODONML (in paml 3.14, January 2004)    H:\users\mt269\runpaml\NEWdna80\b3912.phy   Model: One dN/dS ratio 
Codon frequencies: F3x4
Site-class models:  PositiveSelection

ns = 6  	ls = 232
# site patterns = 66
    7    4    7    3   16    7    6   12    8    1    7    2    3    1    2
    4    6    5    3    5    4    5    3    5    2    8    2    7    4    3
    1    1    6    6    3    4    8    4    2    1    1    1    2    3    8
    1    1    2    1    3    2    1    1    1    1    1    3    2    1    1
    1    2    1    1    1    1

1       
K12                   ATG AAT AAA ATC CTG TTA GTT GAT GAC CGA GAG ACT TCC CTA AAG CTC GAA GGC TTC AAC GTG ATT GCC CAC GGG CAG GCG CTT AGC TTG GAC GTA CCG GGT ACA GCA CGC ACG CCT GTC AGT GGC TAT TTT CGT GCG TCG TGG GAG CAA TCA GTG CTG CCA GGC CGT ACC CAT GAA TTG AAA GCT ATT TTG TAT TCT 
S2457T                ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..A ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 
Shig4337              ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..A ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 
CFT073                ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..T ... ... ... ... ... ... ... ... ... ... ..G ... ... ... ... ... ... ..A ... ... AGT ..C ..T ..T ..C ... ... ..G C.. ..G ... ..C C.. ... ... 
O157                  ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..A ... ... ... ... ... ... ... C.. ... ... ..C ... ..C ... 
EDL933                ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..A ... ... ... ... ... ... ... C.. ... ... ..C ... ..C ... 

Codon usage in sequences
--------------------------------------------------------------------------------------------------------------
Phe TTT  3  3  3  3  3  3 | Ser TCT  1  1  1  1  1  1 | Tyr TAT  3  3  3  3  2  2 | Cys TGT  0  0  0  0  0  0
    TTC  3  3  3  3  3  3 |     TCC  3  3  3  3  3  3 |     TAC  0  0  0  0  1  1 |     TGC  0  0  0  0  0  0
Leu TTA  7  7  7  7  7  7 |     TCA  2  2  2  2  2  2 | *** TAA  0  0  0  0  0  0 | *** TGA  0  0  0  0  0  0
    TTG  5  5  5  3  4  4 |     TCG  1  1  1  1  1  1 |     TAG  0  0  0  0  0  0 | Trp TGG  2  2  2  2  2  2
--------------------------------------------------------------------------------------------------------------
Leu CTT  7  7  7  7  7  7 | Pro CCT  2  2  2  3  2  2 | His CAT  2  2  2  2  2  2 | Arg CGT  9  9  9  8  9  9
    CTC  4  4  4  5  4  4 |     CCC  0  0  0  0  0  0 |     CAC  5  5  5  5  5  5 |     CGC  8  8  8  9  8  8
    CTA  1  1  1  1  1  1 |     CCA  1  1  1  0  1  1 | Gln CAA  3  3  3  3  3  3 |     CGA  1  1  1  1  1  1
    CTG 17 17 17 18 18 18 |     CCG  6  6  6  6  6  6 |     CAG  8  8  8  8  8  8 |     CGG  0  0  0  0  0  0
--------------------------------------------------------------------------------------------------------------
Ile ATT  6  6  6  5  5  5 | Thr ACT  2  2  2  2  2  2 | Asn AAT  4  4  4  4  4  4 | Ser AGT  1  1  1  2  1  1
    ATC  3  3  3  4  4  4 |     ACC  3  3  3  3  3  3 |     AAC  5  5  5  5  5  5 |     AGC  4  4  4  4  4  4
    ATA  0  0  0  0  0  0 |     ACA  3  3  3  3  3  3 | Lys AAA  8  8  8  7  8  8 | Arg AGA  0  0  0  0  0  0
Met ATG  7  7  7  7  7  7 |     ACG  4  4  4  4  4  4 |     AAG  2  2  2  3  2  2 |     AGG  0  0  0  0  0  0
--------------------------------------------------------------------------------------------------------------
Val GTT  6  6  6  6  6  6 | Ala GCT  2  2  2  2  2  2 | Asp GAT 12 12 12 13 12 12 | Gly GGT  6  6  6  7  6  6
    GTC  1  1  1  1  1  1 |     GCC  3  3  3  3  3  3 |     GAC  9  9  9  8  9  9 |     GGC  7  7  7  5  7  7
    GTA  1  1  1  1  2  2 |     GCA  4  5  5  4  4  4 | Glu GAA  7  7  7  7  7  7 |     GGA  0  0  0  0  0  0
    GTG  5  5  5  4  4  4 |     GCG  3  2  2  3  3  3 |     GAG  8  8  8  8  8  8 |     GGG  2  2  2  3  2  2
--------------------------------------------------------------------------------------------------------------

Codon position x base (3x4) table for each sequence.

#1: K12            
position  1:    T:0.12931    C:0.31897    A:0.22414    G:0.32759
position  2:    T:0.32759    C:0.17241    A:0.32759    G:0.17241
position  3:    T:0.28448    C:0.25000    A:0.16379    G:0.30172

#2: S2457T         
position  1:    T:0.12931    C:0.31897    A:0.22414    G:0.32759
position  2:    T:0.32759    C:0.17241    A:0.32759    G:0.17241
position  3:    T:0.28448    C:0.25000    A:0.16810    G:0.29741

#3: Shig4337       
position  1:    T:0.12931    C:0.31897    A:0.22414    G:0.32759
position  2:    T:0.32759    C:0.17241    A:0.32759    G:0.17241
position  3:    T:0.28448    C:0.25000    A:0.16810    G:0.29741

#4: CFT073         
position  1:    T:0.12069    C:0.32759    A:0.22845    G:0.32328
position  2:    T:0.32328    C:0.17241    A:0.32759    G:0.17672
position  3:    T:0.29310    C:0.25000    A:0.15517    G:0.30172

#5: O157           
position  1:    T:0.12500    C:0.32328    A:0.22414    G:0.32759
position  2:    T:0.32759    C:0.17241    A:0.32759    G:0.17241
position  3:    T:0.27586    C:0.25862    A:0.16810    G:0.29741

#6: EDL933         
position  1:    T:0.12500    C:0.32328    A:0.22414    G:0.32759
position  2:    T:0.32759    C:0.17241    A:0.32759    G:0.17241
position  3:    T:0.27586    C:0.25862    A:0.16810    G:0.29741

Sums of codon usage counts
------------------------------------------------------------------------------
Phe F TTT      18 | Ser S TCT       6 | Tyr Y TAT      16 | Cys C TGT       0
      TTC      18 |       TCC      18 |       TAC       2 |       TGC       0
Leu L TTA      42 |       TCA      12 | *** * TAA       0 | *** * TGA       0
      TTG      26 |       TCG       6 |       TAG       0 | Trp W TGG      12
------------------------------------------------------------------------------
Leu L CTT      42 | Pro P CCT      13 | His H CAT      12 | Arg R CGT      53
      CTC      25 |       CCC       0 |       CAC      30 |       CGC      49
      CTA       6 |       CCA       5 | Gln Q CAA      18 |       CGA       6
      CTG     105 |       CCG      36 |       CAG      48 |       CGG       0
------------------------------------------------------------------------------
Ile I ATT      33 | Thr T ACT      12 | Asn N AAT      24 | Ser S AGT       7
      ATC      21 |       ACC      18 |       AAC      30 |       AGC      24
      ATA       0 |       ACA      18 | Lys K AAA      47 | Arg R AGA       0
Met M ATG      42 |       ACG      24 |       AAG      13 |       AGG       0
------------------------------------------------------------------------------
Val V GTT      36 | Ala A GCT      12 | Asp D GAT      73 | Gly G GGT      37
      GTC       6 |       GCC      18 |       GAC      53 |       GGC      40
      GTA       8 |       GCA      26 | Glu E GAA      42 |       GGA       0
      GTG      27 |       GCG      16 |       GAG      48 |       GGG      13
------------------------------------------------------------------------------


Codon position x base (3x4) table, overall

position  1:    T:0.12644    C:0.32184    A:0.22486    G:0.32687
position  2:    T:0.32687    C:0.17241    A:0.32759    G:0.17313
position  3:    T:0.28305    C:0.25287    A:0.16523    G:0.29885

Codon frequencies under model, for use in evolver:
  0.01197116  0.01069505  0.00698824  0.01263960
  0.00631446  0.00564134  0.00368611  0.00666704
  0.01199747  0.01071855  0.00000000  0.00000000
  0.00634077  0.00566485  0.00000000  0.00669482
  0.03047205  0.02722376  0.01778825  0.03217353
  0.01607317  0.01435978  0.00938281  0.01697065
  0.03053902  0.02728359  0.01782735  0.03224424
  0.01614014  0.01441962  0.00942191  0.01704137
  0.02128962  0.01902017  0.01242795  0.02247838
  0.01122969  0.01003262  0.00655540  0.01185673
  0.02133641  0.01906197  0.01245527  0.02252779
  0.01127648  0.01007442  0.00658272  0.01190613
  0.03094817  0.02764913  0.01806619  0.03267624
  0.01632431  0.01458416  0.00952942  0.01723582
  0.03101619  0.02770990  0.01810590  0.03274806
  0.01639233  0.01464492  0.00956913  0.01730764



Nei & Gojobori 1986. dN/dS (dN, dS)
(Note: This matrix is not used in later m.l. analysis.
Use runmode = -2 for ML pairwise comparison.)

K12                 
S2457T               0.0000 (0.0000 0.0059)
Shig4337             0.0000 (0.0000 0.0059)-1.0000 (0.0000 0.0000)
CFT073               0.0617 (0.0048 0.0773) 0.0569 (0.0048 0.0838) 0.0569 (0.0048 0.0838)
O157                 0.0000 (0.0000 0.0238) 0.0000 (0.0000 0.0299) 0.0000 (0.0000 0.0299) 0.0621 (0.0045 0.0717)
EDL933               0.0000 (0.0000 0.0238) 0.0000 (0.0000 0.0299) 0.0000 (0.0000 0.0299) 0.0621 (0.0045 0.0717)-1.0000 (0.0000 0.0000)


TREE #  1:  ((1, (2, 3)), (4, (5, 6)));   MP score: 16
This is a rooted tree.  Please check!
lnL(ntime: 10  np: 15):  -1018.367740     +0.000000
   7..8     8..1     8..9     9..2     9..3     7..10   10..4    10..11   11..5    11..6  
  0.01542  0.00000  0.00482  0.00000  0.00000  0.00833  0.12025  0.01260  0.00000  0.00000  4.78439  0.99565  0.00000  0.00000 166.54204
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.16144

((1: 0.000004, (2: 0.000004, 3: 0.000004): 0.004822): 0.015420, (4: 0.120246, (5: 0.000004, 6: 0.000004): 0.012600): 0.008332);

((K12: 0.000004, (S2457T: 0.000004, Shig4337: 0.000004): 0.004822): 0.015420, (CFT073: 0.120246, (O157: 0.000004, EDL933: 0.000004): 0.012600): 0.008332);

Detailed output identifying parameters

kappa (ts/tv) =  4.78439


dN/dS for site classes (K=3)

p:   0.99565  0.00000  0.00435
w:   0.00000  1.00000166.54204

dN & dS for each branch

 branch           t        S        N    dN/dS       dN       dS   S*dS   N*dN

   7..8       0.015    221.5    474.5   0.7238   0.0046   0.0063    1.4    2.2
   8..1       0.000    221.5    474.5   0.7238   0.0000   0.0000    0.0    0.0
   8..9       0.005    221.5    474.5   0.7238   0.0014   0.0020    0.4    0.7
   9..2       0.000    221.5    474.5   0.7238   0.0000   0.0000    0.0    0.0
   9..3       0.000    221.5    474.5   0.7238   0.0000   0.0000    0.0    0.0
   7..10      0.008    221.5    474.5   0.7238   0.0025   0.0034    0.8    1.2
  10..4       0.120    221.5    474.5   0.7238   0.0357   0.0494   10.9   17.0
  10..11      0.013    221.5    474.5   0.7238   0.0037   0.0052    1.1    1.8
  11..5       0.000    221.5    474.5   0.7238   0.0000   0.0000    0.0    0.0
  11..6       0.000    221.5    474.5   0.7238   0.0000   0.0000    0.0    0.0


Positively selected sites (*: P>95%; **: P>99%)

            Pr(w>1)     post mean +- SE for w

   140 V      1.000**       166.542


Bayes Empirical Bayes (BEB) analysis
Positively selected sites (*: P>95%; **: P>99%)

            Pr(w>1)     post mean +- SE for w

   140 V      0.917         5.405 +- 2.985



The grid (see ternary graph for p0-p1)

w0:   0.050  0.150  0.250  0.350  0.450  0.550  0.650  0.750  0.850  0.950
w2:   1.500  2.500  3.500  4.500  5.500  6.500  7.500  8.500  9.500 10.500


Posterior on the grid

w0:   0.943  0.054  0.003  0.000  0.000  0.000  0.000  0.000  0.000  0.000
w2:   0.125  0.119  0.115  0.109  0.102  0.096  0.090  0.085  0.081  0.078

Posterior for p0-p1 (see the ternary graph)

 0.000
 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.005
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.003 0.066
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.005 0.040 0.880

sum of density on p0-p1 =   1.000000
