Question: Prottest output to Raxml input
0
gravatar for yp19
4 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 • 227 views
ADD COMMENTlink modified 4 months ago by Mensur Dlakic2.1k • written 4 months ago by yp1950
3
gravatar for flogin
4 months 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 4 months ago by flogin150
2
gravatar for Mensur Dlakic
4 months ago by
Mensur Dlakic2.1k
USA
Mensur Dlakic2.1k wrote:

Your command correctly interprets the prottest output.

ADD COMMENTlink written 4 months ago by Mensur Dlakic2.1k
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: 1771 users visited in the last hour