begin mrbayes;
   set autoclose=yes nowarn=yes;
   execute data50/time0.50_birth06.0_death_02.0_dup06.0_loss02.0_glen02.0_p0.005/run000/Greal_seq.nex;
   lset nst=1;
   prset statefreqpr=fixed(equal);
   showmodel;
   mcmc nruns=1 nchains=2 ngen=500000 samplefreq=500 diagnfreq=5000 file=data50/time0.50_birth06.0_death_02.0_dup06.0_loss02.0_glen02.0_p0.005/run000/G.mb;
   sumt conformat=simple relburnin=yes burninfrac=0.25;
end;
