Skip to content

Instantly share code, notes, and snippets.

Created March 9, 2012 02:44
What would you like to do?
Substitutions in a phylogenetic tree file (
#!/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