Question: Prottest output to Raxml input
0
gravatar for yp19
12 months ago by
yp1950
yp1950 wrote:

Hi all!

I ran Prottest to determine the best model for a given protein sequence alignment, so I can eventually use that model in RAxML. The output of Prottest tells me that the best model according to AIC is: WAG+I+G+F.

Now I have to set the model in RAxML command. Based on the RAxML manual I think "WAG" is the model, "I" is for invariant sites, "F" is for using empirical AA frequency, but i'm not sure what the "G" is for? This is what I have come up with, does anyone know if this is the right model based on the Prottest output?

raxmlHPC-PTHREADS-SSE3 -m PROTGAMMAIWAGF

Thanks!

raxml prottest • 502 views
ADD COMMENTlink modified 12 months ago by Mensur Dlakic5.8k • written 12 months ago by yp1950
3
gravatar for flogin
12 months ago by
flogin250
Brazil
flogin250 wrote:

G represents a Gamma Distribution (see http://treethinkers.org/tutorials/the-gamma-distribution/ for a well detailed explanation).

ADD COMMENTlink written 12 months ago by flogin250
2
gravatar for Mensur Dlakic
12 months ago by
Mensur Dlakic5.8k
USA
Mensur Dlakic5.8k wrote:

Your command correctly interprets the prottest output.

ADD COMMENTlink written 12 months ago by Mensur Dlakic5.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1563 users visited in the last hour