annotate Iterative_mapping/unmap.py @ 80:5a13abe09699 draft

Uploaded

1 #!/usr/bin/env python
2 # -*- coding: utf-8 -*-
3 
4 import sys
5 from read_file import *
6 from Bio import SeqIO
7 
8 map_file = sys.argv[1]
9 result_file = sys.argv[2]
10 
11 
12 #reads = read_t_file(read_file);
13 
14 f = open(map_file);
15 h = file(result_file, 'w')
16 
17 for aline in f.readlines():
18     tline = aline.strip();
19     tl = tline.split('\t');
20     if len(tl)>4:
21         if int(tl[1].strip()) != 0:
22             h.write(tl[0].strip());
23             h.write('\n');
24 
25 
26 f.close();
27 h.close()
28 
29 
30 
31