annotate upload/Iterative_mapping/map_ex.py @ 34:d74ed492efdd draft

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