annotate mega_galaxy_input.py @ 0:76d433f51d11 draft default tip

Uploaded
author jaredgk
date Wed, 03 Jul 2019 15:12:04 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
1 import sys
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
2 import shutil
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
3 import os
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
4
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
5 def check_data_format(fn):
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
6 f = open(fn,'r')
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
7 l = f.readline()
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
8 if l[0:5] == '#mega':
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
9 return 'meg'
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
10 os.environ['MDTYPE'] = 1
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
11 return 'fasta'
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
12
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
13 data_name = str(sys.argv[1])
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
14 data_mod_name = 'mega_data.'+check_data_format(data_name)
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
15 os.symlink(data_name,data_mod_name)
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
16 if len(sys.argv) == 3:
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
17 tree_name = str(sys.argv[2])
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
18 tree_mod_name = 'mega_tree.nwk'
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
19 os.symlink(tree_name,tree_mod_name)
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
20 lf = open('mega_data.txt','w')
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
21 lf.write(data_mod_name+'\n')
76d433f51d11 Uploaded
jaredgk
parents:
diff changeset
22