begin mrbayes;
log start filename=Gene1.log replace;
mcmcp ngen=50000000 printfreq=10000 samplefreq=100
mcmcdiagn=yes diagnfreq=500000 stoprule=yes stopval=0.01
nchains=4 savebrlens=yes filename=Gene1;
mcmc;
log stop;
end;