Created

Embed URL

HTTPS clone URL

SSH clone URL

You can clone with HTTPS or SSH.

Download Gist

Substitutions in a phylogenetic tree file (http://silveiraneto.net/?p=4073)

View tree_substitution.py
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
#!/usr/bin/env python
 
import csv
import sys
 
# check parameters
if len(sys.argv) < 2:
sys.exit('Usage: %s TREE_FILE ID_FILE' % sys.argv[0])
 
# initiates the tree
tree = ''
f = open(sys.argv[1], 'rw')
for line in f:
tree = line
 
# open the ID keys and replace them in the tree
markers_file = open(sys.argv[2], 'rt')
reader = csv.DictReader(markers_file, delimiter=";")
for row in reader:
if row['FICD'] != '':
before = row['FICD']
after = row['Taxon Order']+'_'+row['Family']+'_'+row['Genus']+'_'+row['ID']
tree = tree.replace(before, after)
 
print tree
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Something went wrong with that request. Please try again.