comparison predict/parse_dis_pac.py @ 27:564795252d1a draft

Uploaded
author tyty
date Mon, 20 Oct 2014 14:42:09 -0400
parents a292aaf51735
children
comparison
equal deleted inserted replaced
26:8bf6c855353b 27:564795252d1a
1 #parse reactivity file into a dictionary
2
3 import sys
4
5 def parse_dist(in_file):
6 result = []
7 distribution = {}
8 name = []
9 f = open(in_file)
10 for aline in f.readlines():
11 line = aline.strip()
12 dis = line.strip()
13 dist = dis.split('\t') #split the line and the reactivites or reads are in a list
14 if len(dist) > 0:
15 if len(dist) == 1:
16 if dist[0].strip().find('coverage')==-1:
17 name.append(line) #add the name in the name list
18 flag = 1
19 t_name = line
20 else:
21 distri = []
22 for i in range(0, len(dist)):
23 distri.append(dist[i].strip())
24 distribution[t_name] = distri #add the list of reactivities into a dictionary
25 result.append(name)
26 result.append(distribution) #Output the dictionary
27 f.close()
28 return result
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43