
seed used = 139991425

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


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

ns = 4  	ls = 344
# site patterns = 122
    4    6    7    6    7    2    2    1    1    4   14    5    1   13    2
   12   12    8    1    7    1   17    3    3    3    5    2    6   10    1
    5    1    5   10    3    7    4   18    5    1    3    5    1    1    4
    2    6    1    5    1    1    1    1    1    3    1    3    2    3    1
    1    2    5    1    1    1    1    2    1    1    1    2    1    1    1
    1    1    1    1    1    1    1    1    1    1    2    1    2    1    1
    1    1    1    1    1    1    2    1    1    1    1    2    1    1    1
    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
    1    1

K12                                  ATG TCT GCT CAG GCA GCT GAA ATT TAT AAT AAA GAC AGT AAC AAG CTG GAT TAC GGG GTT CAC TTC TCC AAT GGT GAT ACT ACT TAT GCC CGT CTT GGC GAA ACC CAA ATC GGT TGG CGC TCC ACC CTT GCC TAC GTA GAC GCG ACT GTC CCA ACC GTG CAA CGT GCA TTT CAA GGC AAC AGC GAT --- TTC AAC TAC GGT TCT GCT GGA ATC GCG GAC ACT AAT GCA GGG AAA GTT CTT CCT GAA GTA TTT AAA GCA GTT GCC GCT ATT GCC GTA GCT GAT CAC TTC GTT GCT GCT CCG GGT AAG CTT GGA GTA TTA GTC AAA TTC CTT AAA AAT ACT GAA GGT --- --- AAC AAG TCC CTG ATA 
S2457T                               ... ... ... ... ... ..A ..G ... ..C ... ... ... ... ... ... ... ... ... ..A ... ... ... ..T ..C ..C ..C ..A ... ... .TT ... ..G ... ... ..T ... ... ... ... ..T ... ... ... ... ..T ... ..T ..A ..C ..T ... ..A ... GGT ..C A.C ... ..G ..T ..T .CT --- --- --- G.. GTA .C. ..C A.. ... G.A ..T A.T GA. .TC TAC ..T ..C AAC AGC --- --- C.G AA. C.. ..T ..A ..A ..G ..C ..T ..T .G. A.. A.. CAT A.. ..C .TA ... ..A ..A T.G ..T .CG C.G A.T G.. ..T A.. ..G .GC ..G .C. TC. GGC GTT .C. G.T GAT A.C G.T 
Shig4337                             ... ... ... ... ... ..A ..G ... ..C ... ... ... ... ... ... ... ... ... ..A ... ... ... ..T ..C ..C ..C ..A ... ... .TT ... ..G ... ... ..T ... ... ... ... ..T ... ... ... ... ..T ... ..T ..A ..C ..T ... ..A ... GGT ..C A.C ... ..G ..T ..T .CT --- --- --- G.. GTA .C. ..C A.. ... G.A ..T A.T GA. .TC TAC ..T ..C AAC AGC --- --- C.G AA. C.. ..T ..A ..A ..G ..C ..T ..T .G. A.. A.. CAT A.. ..C .TA ... ..A ..A T.G ..T .CG C.G A.T G.. ..T A.. ..G .GC ..G .C. TC. GGC GTT .C. G.T GAT A.C G.T 
CFT073                               ... ... ... ... ... ..A ..G ... ..C ... ... ... ... ... ... ... ... ... ..A ... ... ... ..T ..C ..C ..C ..A ... ... .TT ... ..G ... ... ..T ... ... ... ... ..T ... ... ... ... ..T ... ..T ..A ..C ..T ... ..A ... GGT ..C A.C ... ..G ..T ..T .CT --- --- --- G.. GTA .C. ..C A.. ... G.A ..T A.T GA. .TC TAC ..T ..C AAC AGC --- --- C.G AA. C.. ..T ..A ..A ..G ..C ..T ..T .G. A.. A.. CAT A.. ..C .TA ... ..A ..A T.G ..T .CG C.G ..T G.. ..T A.. ..G .GC ..G .C. TC. GGC GTT .C. G.T GAT A.C G.T 



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

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

#1: K12            
position  1:    T:0.20475    C:0.14243    A:0.25519    G:0.39763
position  2:    T:0.21958    C:0.22255    A:0.40356    G:0.15430
position  3:    T:0.39466    C:0.29970    A:0.18694    G:0.11869

#2: S2457T         
position  1:    T:0.19760    C:0.13473    A:0.27545    G:0.39222
position  2:    T:0.21856    C:0.22156    A:0.39820    G:0.16168
position  3:    T:0.39521    C:0.29341    A:0.18263    G:0.12874

#3: Shig4337       
position  1:    T:0.19760    C:0.13473    A:0.27545    G:0.39222
position  2:    T:0.21856    C:0.22156    A:0.39820    G:0.16168
position  3:    T:0.39521    C:0.29341    A:0.18263    G:0.12874

#4: CFT073         
position  1:    T:0.19760    C:0.13473    A:0.27246    G:0.39521
position  2:    T:0.21856    C:0.22156    A:0.39820    G:0.16168
position  3:    T:0.39521    C:0.29341    A:0.18263    G:0.12874

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

(Ambiguity data are not used in the counts.)


Codon position x base (3x4) table, overall

position  1:    T:0.19940    C:0.13667    A:0.26960    G:0.39432
position  2:    T:0.21882    C:0.22181    A:0.39955    G:0.15982
position  3:    T:0.39507    C:0.29500    A:0.18372    G:0.12621

Codon frequencies under model, for use in evolver:
  0.01778142  0.01327724  0.00826886  0.00568064
  0.01802417  0.01345850  0.00838175  0.00575819
  0.03246778  0.02424343  0.00000000  0.00000000
  0.01298711  0.00969737  0.00000000  0.00414900
  0.01218726  0.00910013  0.00566742  0.00389347
  0.01235364  0.00922436  0.00574479  0.00394663
  0.02225320  0.01661628  0.01034837  0.00710924
  0.00890128  0.00664651  0.00413935  0.00284370
  0.02404154  0.01795162  0.01118000  0.00768057
  0.02436976  0.01819670  0.01133263  0.00778542
  0.04389838  0.03277856  0.02041399  0.01402425
  0.01755935  0.01311143  0.00816560  0.00560970
  0.03516325  0.02625612  0.01635191  0.01123363
  0.03564330  0.02661456  0.01657514  0.01138699
  0.06420594  0.04794205  0.02985758  0.02051192
  0.02568238  0.01917682  0.01194303  0.00820477



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

K12                 
S2457T               0.1115 (0.0706 0.6328)
Shig4337             0.1115 (0.0706 0.6328)-1.0000 (0.0000 0.0000)
CFT073               0.1094 (0.0691 0.6320)-1.0000 (0.0013 0.0000)-1.0000 (0.0013 0.0000)


TREE #  1:  ((1, (2, 3)), 4);   MP score: -1
This is a rooted tree.  Please check!
lnL(ntime:  6  np: 11):  -1787.068032     +0.000000
   5..6     6..1     6..7     7..2     7..3     5..4  
  0.00000  9.32771  0.00277  0.00000  0.00000  0.00000  2.30933  0.93462  0.00000  0.04745 105.79200
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   9.33049

((1: 9.327713, (2: 0.000004, 3: 0.000004): 0.002769): 0.000000, 4: 0.000004);

((K12: 9.327713, (S2457T: 0.000004, Shig4337: 0.000004): 0.002769): 0.000000, CFT073: 0.000004);

Detailed output identifying parameters

kappa (ts/tv) =  2.30933


dN/dS for site classes (K=3)

p:   0.93462  0.00000  0.06538
w:   0.04745  1.00000105.79200

dN & dS for each branch

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

   5..6       0.000    262.5    769.5   6.9610   0.0000   0.0000    0.0    0.0
   6..1       9.328    262.5    769.5   6.9610   3.9751   0.5710  149.9 3058.8
   6..7       0.003    262.5    769.5   6.9610   0.0012   0.0002    0.0    0.9
   7..2       0.000    262.5    769.5   6.9610   0.0000   0.0000    0.0    0.0
   7..3       0.000    262.5    769.5   6.9610   0.0000   0.0000    0.0    0.0
   5..4       0.000    262.5    769.5   6.9610   0.0000   0.0000    0.0    0.0


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

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

   131 Q      0.998**       105.602
   133 A      0.506         53.510
   164 S      0.539         57.031
   174 Y      0.998**       105.548
   195 I      0.510         53.940
   207 T      0.926         97.930
   211 A      0.998**       105.529
   214 V      0.963*        101.896
   215 L      0.997**       105.434
   218 V      0.718         75.996
   219 F      0.972*        102.805
   258 F      0.987*        104.445
   301 V      0.949         100.373
   336 G      0.977*        103.392
   341 K      0.912         96.504
   342 S      0.989*        104.652
   343 L      0.886         93.777


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

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

   131 Q      0.866         7.290 +- 3.118
   174 Y      0.906         7.599 +- 2.872
   207 T      0.585         5.160 +- 3.964
   211 A      0.903         7.579 +- 2.892
   214 V      0.645         5.640 +- 3.912
   215 L      0.883         7.431 +- 3.025
   219 F      0.722         6.205 +- 3.722
   258 F      0.768         6.581 +- 3.592
   336 G      0.745         6.376 +- 3.647
   341 K      0.519         4.666 +- 3.991
   342 S      0.798         6.800 +- 3.472
   343 L      0.589         5.114 +- 4.158



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:   1.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
w2:   0.002  0.010  0.026  0.047  0.072  0.100  0.131  0.166  0.203  0.243

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.016
 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.187 0.796

sum of density on p0-p1 =   1.000000
