Skip to content

A simple Python script that annotates an input rooted tree with counts of potentially decisive and congruent genes in a large set of gene trees (some of which are missing taxa).

License

Notifications You must be signed in to change notification settings

claumer/parse_gene_trees.py

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

4 Commits
 
 
 
 
 
 

Repository files navigation

Parse_gene_trees.py v1.0

Requires the Environment for Tree Exploration v2 (ETE2), argparse 1.1+.

The purpose of this script is to parse a large set of gene trees with missing taxa and, for each rooted node in a single user-defined input tree, annotate the proportion of genes in the large set that are a.) potentially decisive for that node, that is, containing at least one member of each descendant branch plus two distinct outgroups, and b.) the number within the above set that are actually congruent with that node.

Usage is as follows (by default prints to stdout):

python parse_gene_trees.py -i [newick-formatted input gene tree containing all taxa, whose nodes are to be annotated]
-t [large set of gene trees with missing taxa, newick-formatted, one per line] \ -o [String representing the name of some outgroup with which to root the input gene tree]

About

A simple Python script that annotates an input rooted tree with counts of potentially decisive and congruent genes in a large set of gene trees (some of which are missing taxa).

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published