annotate structurefold/reactivity_cal/parse_dis_react.py @ 113:aedb21527abd draft

Uploaded
author tyty
date Tue, 14 Apr 2015 14:09:42 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
113
aedb21527abd Uploaded
tyty
parents:
diff changeset
1 #!/usr/bin/env python
aedb21527abd Uploaded
tyty
parents:
diff changeset
2 # -*- coding: utf-8 -*-
aedb21527abd Uploaded
tyty
parents:
diff changeset
3 import sys
aedb21527abd Uploaded
tyty
parents:
diff changeset
4
aedb21527abd Uploaded
tyty
parents:
diff changeset
5 def parse_dist(in_file):
aedb21527abd Uploaded
tyty
parents:
diff changeset
6 result = []
aedb21527abd Uploaded
tyty
parents:
diff changeset
7 distribution = {}
aedb21527abd Uploaded
tyty
parents:
diff changeset
8 name = []
aedb21527abd Uploaded
tyty
parents:
diff changeset
9 f = open(in_file)
aedb21527abd Uploaded
tyty
parents:
diff changeset
10 flag = 0
aedb21527abd Uploaded
tyty
parents:
diff changeset
11 for aline in f.readlines():
aedb21527abd Uploaded
tyty
parents:
diff changeset
12 line = aline.strip()
aedb21527abd Uploaded
tyty
parents:
diff changeset
13 dis = line.strip()
aedb21527abd Uploaded
tyty
parents:
diff changeset
14 dist = dis.split('\t')
aedb21527abd Uploaded
tyty
parents:
diff changeset
15 if len(dist) > 0:
aedb21527abd Uploaded
tyty
parents:
diff changeset
16 if len(dist) == 1:
aedb21527abd Uploaded
tyty
parents:
diff changeset
17 if dist[0].strip().find('coverage')==-1:
aedb21527abd Uploaded
tyty
parents:
diff changeset
18 if flag == 0:
aedb21527abd Uploaded
tyty
parents:
diff changeset
19 name.append(line)
aedb21527abd Uploaded
tyty
parents:
diff changeset
20 flag = 1
aedb21527abd Uploaded
tyty
parents:
diff changeset
21 t_name = line
aedb21527abd Uploaded
tyty
parents:
diff changeset
22 else:
aedb21527abd Uploaded
tyty
parents:
diff changeset
23 distribution[t_name] = 'null'
aedb21527abd Uploaded
tyty
parents:
diff changeset
24 name.append(line)
aedb21527abd Uploaded
tyty
parents:
diff changeset
25 flag = 1
aedb21527abd Uploaded
tyty
parents:
diff changeset
26 t_name = line
aedb21527abd Uploaded
tyty
parents:
diff changeset
27 else:
aedb21527abd Uploaded
tyty
parents:
diff changeset
28 distri = []
aedb21527abd Uploaded
tyty
parents:
diff changeset
29 for i in range(0, len(dist)):
aedb21527abd Uploaded
tyty
parents:
diff changeset
30 distri.append(dist[i].strip())
aedb21527abd Uploaded
tyty
parents:
diff changeset
31 distribution[t_name] = distri
aedb21527abd Uploaded
tyty
parents:
diff changeset
32 flag = 0
aedb21527abd Uploaded
tyty
parents:
diff changeset
33 result.append(name)
aedb21527abd Uploaded
tyty
parents:
diff changeset
34 result.append(distribution)
aedb21527abd Uploaded
tyty
parents:
diff changeset
35 f.close()
aedb21527abd Uploaded
tyty
parents:
diff changeset
36 return result
aedb21527abd Uploaded
tyty
parents:
diff changeset
37
aedb21527abd Uploaded
tyty
parents:
diff changeset
38
aedb21527abd Uploaded
tyty
parents:
diff changeset
39
aedb21527abd Uploaded
tyty
parents:
diff changeset
40
aedb21527abd Uploaded
tyty
parents:
diff changeset
41
aedb21527abd Uploaded
tyty
parents:
diff changeset
42
aedb21527abd Uploaded
tyty
parents:
diff changeset
43
aedb21527abd Uploaded
tyty
parents:
diff changeset
44
aedb21527abd Uploaded
tyty
parents:
diff changeset
45
aedb21527abd Uploaded
tyty
parents:
diff changeset
46
aedb21527abd Uploaded
tyty
parents:
diff changeset
47
aedb21527abd Uploaded
tyty
parents:
diff changeset
48
aedb21527abd Uploaded
tyty
parents:
diff changeset
49
aedb21527abd Uploaded
tyty
parents:
diff changeset
50
aedb21527abd Uploaded
tyty
parents:
diff changeset
51