comparison splitSHAPE.py @ 9:0690d59881b9 draft

planemo upload for repository https://github.com/eteriSokhoyan/galaxytools/tree/branchForIterations/tools/GraphClust commit 65d322f9ab2f24d65b307f3553589149a1d678d5
author rnateam
date Thu, 01 Jun 2017 08:09:49 -0400
parents a04e93fdb40a
children c0c9d19bc7b2
comparison
equal deleted inserted replaced
8:a04e93fdb40a 9:0690d59881b9
34 with open(shape_file, 'r') as shape: 34 with open(shape_file, 'r') as shape:
35 content = shape.read() 35 content = shape.read()
36 lines = content.split('\n') 36 lines = content.split('\n')
37 for line in lines: 37 for line in lines:
38 if pattern.match(line): 38 if pattern.match(line):
39 line = line.replace('>','').strip() 39 line = line.replace('>','').split()[0]
40 react_arr=[] 40 react_arr=[]
41 react_dict[line] = react_arr 41 react_dict[line] = react_arr
42 continue 42 continue
43 else: 43 else:
44 react_arr.append(line) 44 react_arr.append(line)
45 45
46 toWrite = "" 46 toWrite = ""
47 chunks = [] 47 chunks = []
48 for i in range(len(orig_id)): 48 for i in range(len(orig_id)):
49 if not orig_id[i] in react_dict: 49 if not orig_id[i] in react_dict:
50 raise RuntimeError('Error key {} not found'.format(orig_id)) 50 raise RuntimeError('Error key {} {} not found'.format(i, orig_id[i]))
51 51
52 react_val = react_dict[orig_id[i]] 52 react_val = react_dict[orig_id[i]]
53 toWrite += '>' + str(seq_id[i]) + " " + seq_string[i] + "\n" 53 toWrite += '>' + str(seq_id[i]) + " " + seq_string[i] + "\n"
54 chunks = re.findall(r'\d+', seq_string[i]) 54 chunks = re.findall(r'\d+', seq_string[i])
55 55