
seed used = 140023641

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


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

ns = 4  	ls = 389
# site patterns = 102
   11    8    3   15   12    8    1    2   16   11    6    6    2    1    6
    9    3   22    3   12    2    6    2    9   12    9    9    1    9    8
    9    6    5    5    6    2    5    6    1    3    1    1   10    1    4
    1    2    1    7    4    4    2    2    1    1    3    3    1    1    1
    8    1    3    3    6    4    1    1    1    1    1    3    1    1    1
    1    1    1    2    2    1    1    1    1    1    1    1    2    1    1
    1    1    1    1    1    1    1    1    1    1    1    1

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

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

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

#1: K12            
position  1:    T:0.13111    C:0.22365    A:0.28278    G:0.36247
position  2:    T:0.29049    C:0.24936    A:0.30077    G:0.15938
position  3:    T:0.31105    C:0.24679    A:0.17995    G:0.26221

#2: S2457T         
position  1:    T:0.12853    C:0.22365    A:0.28535    G:0.36247
position  2:    T:0.29049    C:0.24936    A:0.29820    G:0.16195
position  3:    T:0.31362    C:0.24936    A:0.16967    G:0.26735

#3: Shig4337       
position  1:    T:0.12853    C:0.22365    A:0.28535    G:0.36247
position  2:    T:0.29049    C:0.24936    A:0.29820    G:0.16195
position  3:    T:0.31362    C:0.24936    A:0.16967    G:0.26735

#4: CFT073         
position  1:    T:0.13368    C:0.22108    A:0.28278    G:0.36247
position  2:    T:0.29306    C:0.24679    A:0.30077    G:0.15938
position  3:    T:0.32134    C:0.22879    A:0.18252    G:0.26735

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


Codon position x base (3x4) table, overall

position  1:    T:0.13046    C:0.22301    A:0.28406    G:0.36247
position  2:    T:0.29113    C:0.24871    A:0.29949    G:0.16067
position  3:    T:0.31491    C:0.24357    A:0.17545    G:0.26607

Codon frequencies under model, for use in evolver:
  0.01221651  0.00944910  0.00680634  0.01032170
  0.01043662  0.00807241  0.00581469  0.00881788
  0.01256709  0.00972026  0.00000000  0.00000000
  0.00674200  0.00521473  0.00000000  0.00569630
  0.02088241  0.01615190  0.01163448  0.01764350
  0.01783994  0.01379865  0.00993939  0.01507293
  0.02148168  0.01661542  0.01196837  0.01814983
  0.01152451  0.00891385  0.00642080  0.00973703
  0.02659949  0.02057389  0.01481972  0.02247386
  0.02272407  0.01757637  0.01266055  0.01919952
  0.02736283  0.02116431  0.01524501  0.02311880
  0.01467963  0.01135425  0.00817865  0.01240279
  0.03394143  0.02625266  0.01891023  0.02867705
  0.02899632  0.02242777  0.01615510  0.02449894
  0.03491547  0.02700605  0.01945291  0.02950001
  0.01873148  0.01448822  0.01043611  0.01582619



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.1686 (0.0088 0.0523)
Shig4337             0.1686 (0.0088 0.0523)-1.0000 (0.0000 0.0000)
CFT073               0.0522 (0.0080 0.1524) 0.0928 (0.0123 0.1321) 0.0928 (0.0123 0.1321)


TREE #  1:  ((1, (2, 3)), 4);   MP score: 53
This is a rooted tree.  Please check!
lnL(ntime:  6  np: 11):  -1863.541041     +0.000000
   5..6     6..1     6..7     7..2     7..3     5..4  
  0.05396  0.05533  0.03424  0.00000  0.00000  0.11350  2.77517  0.99736  0.00000  0.08263 128.49764
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.25704

((1: 0.055335, (2: 0.000004, 3: 0.000004): 0.034239): 0.053961, 4: 0.113502);

((K12: 0.055335, (S2457T: 0.000004, Shig4337: 0.000004): 0.034239): 0.053961, CFT073: 0.113502);

Detailed output identifying parameters

kappa (ts/tv) =  2.77517


dN/dS for site classes (K=3)

p:   0.99736  0.00000  0.00264
w:   0.08263  1.00000128.49764

dN & dS for each branch

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

   5..6       0.054    328.3    838.7   0.4214   0.0130   0.0308   10.1   10.9
   6..1       0.055    328.3    838.7   0.4214   0.0133   0.0316   10.4   11.2
   6..7       0.034    328.3    838.7   0.4214   0.0082   0.0195    6.4    6.9
   7..2       0.000    328.3    838.7   0.4214   0.0000   0.0000    0.0    0.0
   7..3       0.000    328.3    838.7   0.4214   0.0000   0.0000    0.0    0.0
   5..4       0.114    328.3    838.7   0.4214   0.0273   0.0648   21.3   22.9


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

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

   388 F      1.000**       128.473


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

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

   388 F      0.865         2.943 +- 1.970



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.994  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
w2:   0.363  0.262  0.154  0.086  0.050  0.031  0.020  0.014  0.011  0.009

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.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.010
 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.001 0.019 0.970

sum of density on p0-p1 =   1.000000
