Question: Prottest output to Raxml input
0
gravatar for yp19
17 days ago by
yp1920
yp1920 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 • 123 views
ADD COMMENTlink modified 17 days ago by Mensur Dlakic210 • written 17 days ago by yp1920
3
gravatar for flogin
17 days ago by
flogin150
FioCruz/Brazil
flogin150 wrote:

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

ADD COMMENTlink written 17 days ago by flogin150
2
gravatar for Mensur Dlakic
17 days ago by
Mensur Dlakic210
USA
Mensur Dlakic210 wrote:

Your command correctly interprets the prottest output.

ADD COMMENTlink written 17 days ago by Mensur Dlakic210
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: 1976 users visited in the last hour