Skip to content

Commit

Permalink
Manual
Browse files Browse the repository at this point in the history
  • Loading branch information
Cedric Chauve committed Apr 26, 2023
1 parent 19c1393 commit a042a62
Showing 1 changed file with 6 additions and 2 deletions.
8 changes: 6 additions & 2 deletions doc/manual.md
Original file line number Diff line number Diff line change
Expand Up @@ -100,16 +100,18 @@ This pipeline does not compute explicitly gene trees for the gene families, that

In the second pipeline (the `ALE` pipeline), the following steps are implemented:
- for each gene family, an MSA is computed using `MACSE`;
- for each gene family, a sample of gene trees is computed from its MSA using `IQ-TREE`;
- for each gene family, a sample of gene trees is computed from its MSA using `IQ-TREE` ultrafast bootstrap;
- for each gene family, a reconciled gene tree is computed from the sampled gene trees using `ALE`;
- for each ancestral species, a set of candidate ancestral gene adjacencies, that might not be compatible with a linear ancestral gene order, is computed using `DeCoSTAR`;
- for each ancestral species, a subset of the candidate ancestral adjacencies is computed, using `spp_dcj`, that is compatible with a linear gene order.

In the case of a gene family of size 3, `IQ-TREE` bootstrap method is not relevant, and a single gene tree is computed.

### Pipeline 3, based on ecceTERA

The third pipeline (the `ecceTERA` pipeline) takes advantage of the fact that the parsimony reconciliation algorithm ecceTERA is implemented within DeCoSTAR:
- for each gene family, an MSA is computed using MACSE;
- for each gene family, a sample of gene trees is computed from its MSA using IQ-TREE;
- for each gene family, a sample of gene trees is computed from its MSA using IQ-TREE ultrafast bootstrap;
- for each gene family, a reconciled gene tree is computed from the sampled gene trees using ecceTERA;
- for each ancestral species, a set of candidate ancestral gene adjacencies, that might not be compatible with a linear ancestral gene order, is computed using `DeCoSTAR`;
- for each ancestral species, a subset of the candidate ancestral adjacencies is computed, using `spp_dcj`, that is compatible with a linear gene order.
Expand Down Expand Up @@ -194,6 +196,8 @@ example.

**Any AGO pipeline requires a gene family input file.**

**AGO requires gene families having each at least three genes.**

### Gene orders

A *gene order* file is a tabulated file where each line contains the
Expand Down

0 comments on commit a042a62

Please sign in to comment.