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!

More Jandui Amorim's questions See All
Similar questions and discussions