File: loc_pic.py

package info (click to toggle)
spades 3.13.1+dfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, sid
  • size: 22,172 kB
  • sloc: cpp: 136,213; ansic: 48,218; python: 16,809; perl: 4,252; sh: 2,115; java: 890; makefile: 507; pascal: 348; xml: 303
file content (51 lines) | stat: -rw-r--r-- 1,466 bytes parent folder | download | duplicates (2)
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
#!/usr/bin/python3

############################################################################
# Copyright (c) 2015 Saint Petersburg State University
# Copyright (c) 2011-2014 Saint Petersburg Academic University
# All Rights Reserved
# See file LICENSE for details.
############################################################################


import sys
import os
import string
import re
import subprocess
import datetime

init_dir = "before_simplification/pos_loc"
dirs = ["tip_clipping_", "err_con_removal_", "bulge_removal_"]
add_dirs = ["final_tip_clipping", "final_bulge_removal", "final_err_con_removal", "final_simplified"]
try:
	g_pos = subprocess.check_output('ls -1 ' + init_dir, shell=True)
except:
	print "Not found:", directory
	exit(0)
g_pos = g_pos.strip().split('\n')

output_root = "interest_pos/"
os.mkdir(output_root)

for pos in g_pos:
	os.mkdir(output_root + pos)

cnt = 0;
for i in range(0, 10):
	for dir_pref in dirs:
		if (os.path.exists(dir_pref + str(i))):
			print "found dir " + dir_pref + str(i) 
			for pos in g_pos:
				file_from = dir_pref + str(i) + "/pos_loc/" + pos + "/kmer1_.dot"
				print "looking for file " + file_from
				if (os.path.exists(file_from)):
					print "found file " + file_from
					file_to = output_root + pos + "/" + str(cnt) + "_" + dir_pref + str(i) + ".dot"
					print "copying to " + file_to
					os.system('cp ' + file_from + " " + file_to) 
			cnt = cnt + 1	
		else:
			print "Finished"
			exit(0)