#nexus

begin trees;
translate
  1  a,
  2  b,
  3  c,
  4  d,
  5  e,
  6  f,
  7  g,
  8  u,
  9  v,
  10  x,
  11  y,
  12  t;
tree rep1 = [&U][!GarliScore -2156.2393][!GarliModel  r 2.0670852 2.7917325 1.6309412 1.1469023 4.6921498 e 0.20110863 0.26326438 0.29601849 0.2396085  a 0.43811563 p 0.030165083 ](5:0.02228479,((((10:0.07237377,(((6:0.00000001,7:0.03261324):0.02558115,8:0.02901330):0.04798779,9:0.11530013):0.02312390):0.02637448,1:0.06274992):0.01305858,2:0.08196595):0.01324317,((12:0.06939769,11:0.04177658):0.14672066,3:0.04479000):0.00760926):0.05970612,4:0.00728599);
tree rep2BEST = [&U][!GarliScore -2156.22744757][!GarliModel  r 2.07168250 2.79864596 1.63652493 1.14801378 4.69942978 e 0.20085183 0.26350061 0.29600857 0.23963898  a 0.42708533 p 0.02131334 ]((7:0.03259502,6:0.00000001):0.02555910,((((2:0.08195637,((4:0.00728792,5:0.02226838):0.05972067,(3:0.04476716,(11:0.04178774,12:0.06937929):0.14672236):0.00761114):0.01322876):0.01301120,1:0.06271076):0.02639598,10:0.07236609):0.02314845,9:0.11522505):0.04800809,8:0.02900760);
end;
begin paup;
clear;
gett file=garli_run.best.all.tre storebr;
lset userbr nst=6 rmat=(2.07168250 2.79864596 1.63652493 1.14801378 4.69942978) base=(0.20085183 0.26350061 0.29600857) rates=gamma shape= 0.42708533 ncat=4 pinv= 0.02131334;
end;
[!****NOTE: The model parameters loaded are the final model estimates****
****from GARLI for the best scoring search replicate (#2).****
****The best model parameters for other trees may vary.****]