Commit 38246e33 authored by LANORE Vincent's avatar LANORE Vincent
Browse files

Added ratio output to inferred tree script

parent 0e9c7c3b
......@@ -98,7 +98,7 @@ else:
MESSAGE("Done.")
#===================================================================================================
STEP("Computingt tree lengths")
STEP("Computing tree lengths")
from ete3 import Tree
def length(tf):
......@@ -111,4 +111,5 @@ def length(tf):
trees = [["tree", tree_filename], ["inferred_tree_aa", inferred_tree_aa_filename], ["inferred_tree_nt", inferred_tree_nt_filename]]
list(map(lambda x: x.append(length(x[1])), trees))
for t in trees:
MESSAGE("Tree " + data(t[0]) + " has length " + data(t[2]))
\ No newline at end of file
MESSAGE("Tree " + data(t[0]) + " has length " + data(t[2]))
MESSAGE("Ratio aa/nt is " + data(trees[1][2] / trees[2][2]))
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment