Hi! How could I define the following evolution model (ProtTest output) to analyze 158 aminoacid sequences (lenght=219aa) using MrBayes, please? What would be the command lines?
Model................................ : MtMam+G
Number of parameters............... : 314 (1 + 313 branch length estimates)
gamma shape (4 rate categories).. = 0.123
-lnL................................ = 987.60
(seconds))
I'm trying:
> lset rates=gamma ngammacat=4
> prset aamodelpr=fixed(MtMam) shapepr=fixed(0.123)
> mcmc ngen=2500000 nchain=4
But the estimated analysis time is over 100 hours!!! Is it correct?
Thanks for some help!