#NEXUS

begin trees;
translate
 1 Ch_maritimum_canescens_16021.assembled,
 2 Ch_maritimum_canescens_16412a.assembled,
 3 Ch_maritimum_canescens_16414a.assembled,
 4 Ch_maritimum_canescens_16414b.assembled;
tree bestREP2 = [&U][!GarliScore -555.847095][!GarliModel  M1 r  0.00125 0.00125 0.00125 0.02355 999.89936 1.00000 e 0.20247 0.21728 0.27654 0.30370 a 0.05025 p 1.00000 ](3:0.00000001,1:0.00000001,2:0.00000001,4:0.00000001);
end;
begin paup;
clear;
gett file=locus000ML.best.tre storebr;
lset userbr nst=6 rmat=(0.001250 0.001250 0.001250 0.023549 999.899356) base=( 0.202469 0.217284 0.276544) rates=gamma shape=0.050250 ncat=4 pinv=1.000000;
end;