|News and Announcements »|
Many downstream analyses require that the phylogenetic tree relating the OTUs in a study be present. The script make_phylogeny.py produces this tree from a multiple sequence alignment. Trees are constructed with a set of sequences representative of the OTUs, by default using FastTree (Price, Dehal, & Arkin, 2009).
Usage: make_phylogeny.py [options]
The result of make_phylogeny.py consists of a newick formatted tree file (.tre) and optionally a log file. The tree file is formatted using the Newick format and this file can be viewed using most tree visualization tools, such as TopiaryTool, FigTree, etc.
The tips of the tree are the first word from the input sequences from the fasta file, e.g.: ‘>101 PC.481_71 RC:1..220’ is represented in the tree as ‘101’.
A simple example of make_phylogeny.py is shown by the following command, where we use the default tree building method (fasttree) and write the file to the current working directory without a log file:
make_phylogeny.py -i $PWD/aligned.fasta -o $PWD/rep_phylo.tre
Alternatively, if the user would prefer using another tree building method (i.e. clearcut (Sheneman, Evans, & Foster, 2006)), then they could use the following command:
make_phylogeny.py -i $PWD/aligned.fasta -t clearcut