Description
Counts the number of extra lineages contributed by a species tree and a set of gene trees. The trees must be specified in the Rich Newick Format. Specifying more than one species tree will run the command multiple times with one execution per species tree.
Usage
| Code Block |
|---|
|
DeepCoalCount_tree {species_tree_ident1 [, species_tree_ident2...]} {gene_tree_ident1 [, gene_tree_ident2...]} [-u] [-b threshold] [-a taxa map] [result output file]
|
...
By default, it is assumed that only one individual is sampled per species in gene trees. However, the option [-a taxa map] allows multiple alleles to be sampled.
Examples
| Code Block |
|---|
|
#NEXUS
BEGIN TREES;
Tree speceiesTree = ((((a,b),c),d),e);
Tree geneTree1 = ((((a,b),c),d),e);
Tree geneTree2 = ((a,b),((c,e),d));
Tree geneTree3 = ((a,c),((b,e),d));
END;
BEGIN PHYLONET;
DeepCoalCount {speceiesTree} {geneTree1, geneTree2, geneTree3};
END;
|
| Code Block |
|---|
|
#NEXUS
BEGIN TREES;
Tree speceiesTree = ((((x,y),w),m),n);
Tree geneTree1 = ((((a::.5,b::.5)::.5,c::.5)::.5,d::.5)::.5,e::.5)::.5;
Tree geneTree2 = ((a::.5,b::.5)::.5,((c::.5,e::.5)::.5,d::.5)::.5)::.5;
Tree geneTree3 = ((a::.5,c::.5)::.5,((b::.5,e::.5)::.5,d::.5)::.5)::.5;
END;
BEGIN PHYLONET;
DeepCoalCount {speciesTree} {geneTree1, geneTree2, geneTree3} -u -b .5 -a <x:a; y:b; w:c; m:d; n:e>;
END;
|
See Also