Here is an example mrbayes block begin mrbayes; lset coding = all Nucmodel= 4by4 Nst= 1 Rates= equal Nbetacat= 5; prset ratepr=variable; mcmc ngen=1500 printfreq=100 samplefreq=1000 nruns=2 nchains=8 temp=0.02; sump burnin=0 nruns=2 printtofile=yes outputname=sumpoutput.out plot=Yes marglike=Yes table=Yes; sumt burnin=0 nruns=2 ntrees=1 displaygeq=0.05 contype=Halfcompat; end;