annotate reactivity_cal/map_ex.py @ 12:a47e0bfd8de8 draft

Uploaded
author tyty
date Mon, 20 Oct 2014 01:49:28 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
12
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
1 #!/usr/bin/env python
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
2 # -*- coding: utf-8 -*-
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
3
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
4 import sys
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
5 from read_file import *
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
6 from Bio import SeqIO
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
7
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
8 map_file = sys.argv[1]
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
9 result_file = sys.argv[2]
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
10
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
11
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
12 #reads = read_t_file(read_file);
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
13
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
14 f = open(map_file);
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
15 h = file(result_file, 'w')
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
16
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
17 for aline in f.readlines():
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
18 tline = aline.strip();
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
19 tl = tline.split('\t');
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
20 if len(tl)>4:
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
21 if int(tl[1].strip())== 0:
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
22 h.write(tline)
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
23 h.write('\n')
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
24
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
25
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
26 f.close();
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
27 h.close()
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
28
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
29
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
30
a47e0bfd8de8 Uploaded
tyty
parents:
diff changeset
31