execute ctree_with_refseq.nex;

begin mrbayes;
    set nowarn=yes seed=42 swapseed=1337;
    lset nst=6 rates=invgamma ngammacat=4;
    prset statefreqpr=fixed(equal) topologypr=uniform brlenspr=unconstrained:exponential(29900);
    mcmcp nruns=1 nchains=4 ngen=50000000 samplefreq=1000 printfreq=1000 diagnfreq=1000 file=test_output append=no;
    mcmc;
    sump burninfrac=0.5;
    sumt burninfrac=0.5;
end;