Find the best model of evolution for Trpc2 using modelomatic:
#For this analysis, all stop codons need to be removed.
#For troubleshooting modelOmatic, I wrote a more detailed post here.
#Navigate to the correct directory and run this command:

modelomatic_OSX Trpc2_alignment_noStops.phy bionj Trpc2_modelomatic_normal.txt 0 normal

Interestingly REV+4dG is the best model, despite it being a protein coding gene. This is essentially "GTR+G". The 4dG is showing the four rate categories (this is an arbitrary number) for the gamma parameter.

Find the "best" (ML) tree in GALI using search replicates:
#run this command
garli Trpc2_garli.conf
#Note how to implement GTR+G in GARLI. See here for more examples of other models.
#Here is my configuration file (Trpc2_garli.conf):
[general]
datafname = Trpc2_garli.nex
constraintfile = none
streefname = stepwise
attachmentspertaxon = 50
ofprefix = Trpc2.GTR4DG
randseed = -1
availablememory = 512
logevery = 10
saveevery = 100
refinestart = 1
outputeachbettertopology = 0
outputcurrentbesttopology = 0
enforcetermconditions = 1
genthreshfortopoterm = 20000
scorethreshforterm = 0.05
significanttopochange = 0.01
outputphyliptree = 0
outputmostlyuselessfiles = 0
writecheckpoints = 0
restart = 0
outgroup = 1
resampleproportion = 1.0
inferinternalstateprobs = 0
outputsitelikelihoods = 0
optimizeinputonly = 0
collapsebranches = 1

searchreps = 8
bootstrapreps = 0

[model1]
datatype = nucleotide
ratematrix = 6rate
statefrequencies = estimate
ratehetmodel = gamma
numratecats = 4
invariantsites = none

[master]
nindivs = 4
holdover = 1
selectionintensity = 0.5
holdoverpenalty = 0
stopgen = 5000000
stoptime = 5000000

startoptprec = 0.5
minoptprec = 0.01
numberofprecreductions = 10
treerejectionthreshold = 50.0
topoweight = 1.0
modweight = 0.05
brlenweight = 0.2
randnniweight = 0.1
randsprweight = 0.3
limsprweight =  0.6
intervallength = 100
intervalstostore = 5
limsprrange = 6
meanbrlenmuts = 5
gammashapebrlen = 1000
gammashapemodel = 1000
uniqueswapbias = 0.1
distanceswapbias = 1.0

Putting bootstrap values on GARLI:
You can bootstrap in GARLI, but to "put" the bootstraps on the nodes, you have to summarize the bootstrap results using an outside method. SumTrees.py in the Dendropy python package works well. I basically followed the exact instructions somebody wrote here.

Running bootstrap iterations in GARLI:
#I ran four independent runs simultaneously of 250 bootstrap iterations that will combine to = 1000.
#Here is my GARLI configuration file:
[general]
datafname = Trpc2_garli.nex
constraintfile = none
streefname = stepwise
attachmentspertaxon = 50
ofprefix = Trpc2.GTR4DG.boot
randseed = -1
availablememory = 512
logevery = 10
saveevery = 100
refinestart = 1
outputeachbettertopology = 0
outputcurrentbesttopology = 0
enforcetermconditions = 1
genthreshfortopoterm = 20000
scorethreshforterm = 0.05
significanttopochange = 0.01
outputphyliptree = 0
outputmostlyuselessfiles = 0
writecheckpoints = 0
restart = 0
outgroup = 1
resampleproportion = 1.0
inferinternalstateprobs = 0
outputsitelikelihoods = 0
optimizeinputonly = 0
collapsebranches = 1

searchreps = 1
bootstrapreps = 250

[model1]
datatype = nucleotide
ratematrix = 6rate
statefrequencies = estimate
ratehetmodel = gamma
numratecats = 4
invariantsites = none

[master]
nindivs = 4
holdover = 1
selectionintensity = 0.5
holdoverpenalty = 0
stopgen = 5000000
stoptime = 5000000

startoptprec = 0.5
minoptprec = 0.01
numberofprecreductions = 10
treerejectionthreshold = 20.0
topoweight = 1.0
modweight = 0.05
brlenweight = 0.2
randnniweight = 0.1
randsprweight = 0.3
limsprweight =  0.6
intervallength = 100
intervalstostore = 5
limsprrange = 6
meanbrlenmuts = 5
gammashapebrlen = 1000
gammashapemodel = 1000
uniqueswapbias = 0.1
distanceswapbias = 1.0

Installing dendropy
#We have "pip" set up on our server, making the python package easy to install:

sudo pip install -U dendropy

Running SumTrees
#place all of your bootstrap outputs and your best tree from GARLI in its own folder

cd /Volumes/Yango/Trpc2/garli/

mkdir sumtrees

cp Trpc2.GTR4DG.best.tre sumtrees/

cp Trpc2.GTR4DG.boot.boot.tre sumtrees/

cp Trpc2.GTR4DG.boot2.boot.tre sumtrees/

cp Trpc2.GTR4DG.boot3.boot.tre sumtrees/

cp Trpc2.GTR4DG.boot4.boot.tre sumtrees/

#from within the folder, run sumtrees.py

#your target tree is the "best" tree from your GARLI output

sumtrees.py --decimals=0 --percentages --no-analysis-metainformation --no-taxa-block Trpc2.GTR4DG.boot.boot.tre Trpc2.GTR4DG.boot2.boot.tre Trpc2.GTR4DG.boot3.boot.tre Trpc2.GTR4DG.boot4.boot.tre --target=Trpc2.GTR4DG.best.tre --output=supportOnBest.phy

Viewing tree with bootstraps

Now, open supportOnBest.phy in FigTree. Click on "node labels" and select "support" as the node label. These are your bootstrap values.

I am going to make a gene tree of all the Trpc2 sequences to demonstrate how different the pseudogenes are from the conserved gene--the hypothesis being that the pseudogenes will have much longer branch lengths.

Running modelo

modelomatic_OSX Trpc2_modelomatic.phy bionj Trpc2_modelomatic_normal.txt 0 normal

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: <Trpc2_modelomatic.phy>: 118 sequences of length 489 (DataMatrix: 118 x 380)

Checking for sparse (>85% gaps) sequences

Starting with 118 ... after removal there are 118 sequences

Creating start tree ... 

Assertion failed: (0), function Error, file tools.cxx, line 235.

Unrecognised codon: TAG for data codon position 153 in sequence 52...Abort trap: 6

105-238:modelomatic loloyohe$ modelomatic_OSX Trpc2_modelomatic.phy bionj Trpc2_modelomatic_normal.txt 0 normal

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: <Trpc2_modelomatic.phy>: 118 sequences of length 489 (DataMatrix: 118 x 380)

Checking for sparse (>85% gaps) sequences

Starting with 118 ... after removal there are 118 sequences

Creating start tree ...  estimated using bionj (0.813652s)

Optimisation settings: normal 

Scanning for modelomatic.ini file ...

Amino acid models skipping rtREV, HIVb, HIVw

Codon models skipping F0, F1X4, F3X4, F64 done

Working with genetic code: Universal

>>> Doing model analysis <<< 

RY Done  (2.84812s)

NT Done  (65.4292s)

AA Done  (267.533s)...

Codon Done  (0.255969s)

Outputting results to <Trpc2_modelomatic_normal.txt>

Successful exit

For once, something worked. The model with the best fit is the K2P+4dG. This is the Kimura two-parameter + gamma with four site rates. A good summary of the nucleotide models can be found here. K2P means that the transitions (alpha) and transversions (beta) will have different rates, but that all nucleotides occur at the same frequency. This model can easily be implemented in GARLI. The K2P model is just a simpler HKY model--HKY allows for different transition and transversion rates, but also allows 4 parameters for nucleotide substitutions. I am using the GARLI on CIPRES. 

#everything is default except
Maximum hours to run=48
-searchreps=8
-bootstrapreps=0
-datatype=nucleotide
#to make it K2P make sure to change base statefrequencies
-ratematrix=HKY(2rate)
-statefrequencies=fixed
-invariant sites=none
#this is where the 4dG comes in from modelomatic
-ratehetmodel=gamma
-numratecats=4

ERROR: state frequencies specified as fixed, but no

        Garli block found in Trpc2_garli.phy!!

Oops, it choked. -statefrequencies should be set to equal:

#everything is default except
Maximum hours to run=48
-searchreps=8
-bootstrapreps=0
-datatype=nucleotide
#to make it K2P make sure to change base statefrequencies
-ratematrix=HKY(2rate)
-statefrequencies=equal
-invariant sites=none
#this is where the 4dG comes in from modelomatic
-ratehetmodel=gamma

-numratecats=4
From the best tree from GARLI:

Now onto testing if some branch lengths are significantly longer...

Running GARLI from our Model-O-Matic code.

If we remember from our model-o-matic output and sort by the lowest AIC, we can see the best model for OR6 is CodonM0.0.EQU+4dG. It is interesting how all codon models are preferred to the nucleotide models! This turned out to be the case in every set of alignments.

Things are looking good. I'll keep you posted on the results.

Here is the following best-fit models for each of the OR gene families (see summary files with AIC and likelihoods in /Volumes/Yango/Hayden_batORs/model_o_matic/excel_summaries):

OR137      CodonM0.0.EQU+4dG

OR213      CodonM0.0F64+4dG

OR4          CodonM0.0.F64+4dG

OR589      CodonM0.0.F64+4dG

OR6          CodonM0.0.EQU + 4dG

OR10        CodonM0.0.F64+4dG

OR11        CodonM0.0F64+4dG

OR51        CodonM0.0F64+4dG

OR52        CodonM0.0F64+4dG

Those highlighted in pink were run in model-O-matic as the -100normal parameter, as the program would crash if it was run as -normal. My thoughts are there are too many taxa, but other than than that I am at a loss.

Now we can apply this model to GARLI using our tree output. There are two options for running GARLI (at least that I use). One is on the CIPRES server, the other can be from command line. I am currently having CPU usage problems on CIPRES so I am going to run it from command.

Installing GARLI is easy if you have homebrew.

105-238:~ loloyohe$ brew install garli

Should be ready to go.

Let's take a second to interpret the output from

105-238:OR6 loloyohe$ garli garli.conf

Running GARLI Version 2.01.1067 (18 May 2012)

->Single processor version for 64-bit OS<-

##############################################################

 This is GARLI 2.0, the first "official" release including 

          partitioned models.  It is a merging of

   official release 1.0 and beta version GARLI-PART 0.97

  Briefly, it includes models for nucleotides, amino acids,

 codons, and morphology-like characters, any of which can be 

  mixed together and applied to different subsets of data.

    General program usage is extensively documented here:

            http://www.nescent.org/wg/garli/

      see this page for details on partitioned usage:

  http://www.nescent.org/wg/garli/Partition_testing_version

   and this page for details on Mkv mophology model usage:

    http://www.nescent.org/wg/garli/Mkv_morphology_model

         PLEASE LET ME KNOW OF ANY PROBLEMS AT:

                garli.support@gmail.com

##############################################################

Compiled Apr 27 2015 16:06:23 using GNU gcc compiler version 4.9.2

Using NCL version 2.1.17

#######################################################

Reading config file garli.conf

###################################################

READING OF DATA

Attempting to read data file in Nexus format (using NCL):

Phyllos_OR6_Funct_tAlign_noStops.nex ...

Reading DATA block...storing implied block: TAXA

storing read block: DATA

 successful

###################################################

PARTITIONING OF DATA AND MODELS

GARLI data subset 1

CHARACTERS block #1 ("Untitled DATA Block 1")

Data read as Nucleotide data,

modeled as Codon data

Gaps or ambiguity codes found within codon for taxon Lsi_29_OR6.

Codons coded as missing for that taxon: 6 120 235 236 

Gaps or ambiguity codes found within codon for taxon Phe_43_OR6.

Codons coded as missing for that taxon: 6 

Gaps or ambiguity codes found within codon for taxon Sti_17_OR6.

Codons coded as missing for that taxon: 185 

Summary of data:

  35 sequences.

  6 constant characters.

  232 parsimony-informative characters.

  5 uninformative variable characters.

  243 total characters.

  243 unique patterns in compressed data matrix.

Pattern processing required < 1 second

###################################################

NOTE: Unlike many programs, the amount of system memory that Garli will

use can be controlled by the user.

(This comes from the availablememory setting in the configuration file.

Availablememory should NOT be set to more than the actual amount of 

physical memory that your computer has installed)

For this dataset:

 Mem level availablememory setting

  great     >= 150 MB

  good approx 149 MB to 99 MB

  low approx 98 MB to 42 MB

  very low approx 42 MB to 31 MB

the minimum required availablememory is 31 MB

You specified that Garli should use at most 512.0 MB of memory.

Garli will actually use approx. 224.4 MB of memory

**Your memory level is: great (you don't need to change anything)**

#######################################################

Found outgroup specification:  1

#######################################################

STARTING RUN

>>>Search rep 1 (of 2)<<<

MODEL REPORT - Parameters are at their INITIAL values (not yet optimized)

Model 1

  Number of states = 61 (codon data, standard code)

  Nucleotide Relative Rate Matrix Assumed by Codon Model:     2 rates (transition and transversion) K param = 4.0000

  Equilibrium State Frequencies: equal (1/61 = 0.01639, fixed)

  Rate Heterogeneity Model:

    4 nonsynonymous rate categories, rate and proportion of each estimated

     (this is effectively the M3 model of PAML)

      dN/dS Proportion

      0.2500 0.2500

      0.5000 0.2500

      0.7500 0.2500

      1.0000 0.2500

Starting with seed=287687

creating likelihood stepwise addition starting tree...

number of taxa added:

 4  5  6  7  8  9  10  11  12  13  14  15  16  17

We can set the exact same thing up in CIPRES. 

I am using modelOmatic to determine which models best fit the sequences, and it include comparisons of codon models. Installing and running modelOmatic. For awhile I thought I was having Yosemite issues but because they are all binaries, everything seems to checking out okay. Lots of "duh" moments were had.

To install on Mac OSX Yosemite:
[1] Download the binaries and move to where you store your programs.

mv ~/Downloads/modelomatic_OSX /usr/local/bin/

[2] Set modelOmatic into your path.

vi .bash_profile  #in your home directory, this opens up your bash profile

#press "i" to edit and then paste this into your .bash_profile

export PATH=/usr/local/bin/modelomatic_OSX:$PATH 

#press ESC, then type "wq!"; restart the terminal

[3] ModelOmatic requires a tree file to run. It can make a NJ tree if you have bioNJ set up. Download the binaries for bioNJ here.

[4] Move the bioNJ into where you store your programs. Set your path. You have to rename your path in lower case so that modelOmatic can read it.

mv ~/Downloads/BIONJ/ /usr/local/bin/ 

mv /usr/local/bin/BIONJ/BIONJ /usr/local/bin/BIONJ/bionj

#add to your path variable
vi .bash_profile

export PATH=/usr/local/bin/BIONJ/BIONJ:$PATH

export PATH=/usr/local/bin/BIONJ:$PATH

The format for input 

./ModelOMatic <data> bionj <output> <genetic_code> fast

I am running it as modelomatic_OSX because it is in my path. 

Using the bionj option tell modelOmatic to make its own tree using bioNJ.

Specify in <output> slot where to put your output folder.

The <genetic_code> has been tricky. It can choke and tell you that you have stop codons if you don't use the correct number:

     0:  Universal code
     1:  Vertebrate mt
     2:  Yeast mt
     3:  Mould mt

     4:  Invertebrate mt
     5:  Ciliate nuclear
     6:  Echinoderm mt
     7:  Euplotid mt

     8:  Alternative yeast nuclear
     9:  Ascidian mt

     10: Blepharisma nuclear

     11: Everything codes (64 character state-space)

My adviser was having trouble but then realized she was using cytb so it actually should have been "1" instead of "0". 

ModelOmatic uses .phylip file. Allegedly can use other formats but I have not tried. However, it actually seems ridiculously flexible in reading a .phylip file.

OR_alignments loloyohe$ modelomatic_OSX Phyllos_OR6_Funct_tAlign.phylip bionj test.txt 1

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: <Phyllos_OR6_Funct_tAlign.phylip>: 35 sequences of length 729 (DataMatrix: 35 x 570)

Checking for sparse (>85% gaps) sequences

Starting with 35 ... after removal there are 35 sequences

Creating start tree ... 

Assertion failed: (0), function Error, file tools.cxx, line 235.

Unrecognised codon: AGA for data codon position 241 in sequence 1...Abort trap: 6

We had stop codons. Okay, so I translated the sequences, looked for where the stop codons were, and then I replaced them in the original .phylip file using "---".  Note in the future that all .phylip sequence files in /Volumes/Yango/Hayden_ORs/OR_alignments have the stop codon replaced, while the .nex files have the original sequence alignments. Lets try to get it running...

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR4_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/Phyllos_OR4_Funct_tAlign_mOm_results.txt 0 fast

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR4_Funct_tAlign.phylip>: 127 sequences of length 753 (DataMatrix: 127 x 666)

Checking for sparse (>85% gaps) sequences

Starting with 127 ... after removal there are 127 sequences

Creating start tree ...  estimated using bionj (1.15646s)

Optimisation settings: fast 

Scanning for modelomatic.ini file ... done

Working with genetic code: Universal

>>> Doing model analysis <<< 

RY Done  (3.56514s)

NT Done  (10.6517s)

AA Done  (71.2184s).........

Codon Done  (232.733s)

Outputting results to <Phyllos_OR4_Funct_tAlign_mOm_results.txt>

Successful exit

####now for the others

modelomatic_OSX /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR6_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/Phyllos_OR6_Funct_tAlign_mOm_results.txt 0 fast

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR10_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/Phyllos_OR10_Funct_tAlign_mOm_results.txt 0 fast

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR51_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/Phyllos_OR51_Funct_tAlign_mOm_results.txt 0 fast

Had a hiccup with this data:

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR2_13_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/model_o_matic/Phyllos_OR2_13_Funct_tAlign_mOm_results.txt 0 normal

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR2_13_Funct_tAlign.phylip>: 264 sequences of length 1530 (DataMatrix: 264 x 943)

Checking for sparse (>85% gaps) sequences

Starting with 264 ... after removal there are 264 sequences

Creating start tree ...  estimated using bionj (4.68751s)

Optimisation settings: normal 

Scanning for modelomatic.ini file ... done

Working with genetic code: Universal

>>> Doing model analysis <<< 

Broken Prob(): -0.000104812

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR137_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/Phyllos_OR137_Funct_tAlign_mOm_results.txt 0 fast

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR137_Funct_tAlign.phylip>: 452 sequences of length 804 (DataMatrix: 452 x 766)

Checking for sparse (>85% gaps) sequences

Starting with 452 ... after removal there are 452 sequences

Creating start tree ...  estimated using bionj (4.7425s)

Optimisation settings: fast 

Scanning for modelomatic.ini file ... done

Working with genetic code: Universal

>>> Doing model analysis <<< 

RY Done  (8.05337s)

NT Done  (44.7082s)

Broken Prob(): -40.8111

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR52_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/Phyllos_OR52_Funct_tAlign_mOm_results.txt 0 fast

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR51_Funct_tAlign.phylip>: 152 sequences of length 750 (DataMatrix: 152 x 712)

Checking for sparse (>85% gaps) sequences

Starting with 152 ... after removal there are 152 sequences

Creating start tree ...  estimated using bionj (1.32059s)

Optimisation settings: fast 

Scanning for modelomatic.ini file ... done

Working with genetic code: Universal

>>> Doing model analysis <<< 

RY Done  (4.69988s)

NT Done  (13.8241s)

AA Done  (85.6409s).........

...

Broken Prob(): -2.64535e+41

I think the problem is in running the "fast" option with large datasets. When I use "trimfast" and "normal", they both work. 

modelomatic_OSX  /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR51_Funct_tAlign.phylip bionj /Volumes/Yango/Hayden_batORs/model_o_matic/Phyllos_OR51_Funct_tAlign_mOm_results.txt 0 trimfast

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR51_Funct_tAlign.phylip>: 152 sequences of length 750 (DataMatrix: 152 x 712)

Checking for sparse (>85% gaps) sequences

Starting with 152 ... after removal there are 152 sequences

Creating start tree ...  estimated using bionj (1.58694s)

Optimisation settings: fast trim=10

TRIMMING: Tree contains 152 sequences (>TrimTree=10) ...

          Loose optimisation of start tree under JC ... done

          Obtaining greedy start tree with 10 sequences ... done

          Reinitialising objects for trimmed data ... done

          New files available: data = </Volumes/Yango/Hayden_batORs/model_o_matic/Phyllos_OR51_Funct_tAlign_mOm_results.txt.trim.data>; tree = </Volumes/Yango/Hayden_batORs/model_o_matic/Phyllos_OR51_Funct_tAlign_mOm_results.txt.trim.tree>

Scanning for modelomatic.ini file ... done

Working with genetic code: Universal

>>> Doing model analysis <<< 

RY Done  (1.08313s)

NT Done  (1.18703s)

AA Done  (6.64464s).........

Codon Done  (24.4936s)

Outputting results to </Volumes/Yango/Hayden_batORs/model_o_matic/Phyllos_OR51_Funct_tAlign_mOm_results.txt>

Successful exit

However, word of caution. The "trimfast" option gives a vastly different result:

Normal it is! With fewer parameters, its best just to wait it out.

Run modelomatic:

1) Download the binary from Google Code website.

2) Move to local bin

mv ~/Downloads/modelomatic_OSX /usr/local/bin

3) Give yourself permission to run it.

sudo chmod +x /usr/local/bin/modelomatic_OSX

4) To run the binary, use the command "xcrun"

xcrun modelomatic_OSX /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR6_Funct_tAlign.fasta

But I get this error:

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR6_Funct_tAlign.fasta>

Segmentation fault: 11

5) Okay I guess I need to make it executable (also make sure it is added to the path)

chmod a+x /usr/local/bin/modelomatic_OSX

6) Now run it

modelomatic_OSX /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR10_Funct_tAlign.fasta

Nope, still get same error.

modelomatic_OSX /Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR6_Funct_tAlign.fasta

But I get this error:

---------------------------------------------------------------------------------------------

  ModelOMatic (v1.01 (release)).

A program for choosing substitutions models for phylogenetic inference.

Written by Simon Whelan.

Contributions from James Allen, Ben Blackburne and David Talavera.

---------------------------------------------------------------------------------------------

Data: </Volumes/Yango/Hayden_batORs/OR_alignments/Phyllos_OR6_Funct_tAlign.fasta>

Segmentation fault: 11