# HG changeset patch # User tyty # Date 1413784168 14400 # Node ID a47e0bfd8de8e4acb120e64aba15fbddfa6e9de9 # Parent a71c7f99c92e0298b76f64b98c2962d65aabb19f Uploaded diff -r a71c7f99c92e -r a47e0bfd8de8 reactivity_cal/map_ex.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/reactivity_cal/map_ex.py Mon Oct 20 01:49:28 2014 -0400 @@ -0,0 +1,31 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- + +import sys +from read_file import * +from Bio import SeqIO + +map_file = sys.argv[1] +result_file = sys.argv[2] + + +#reads = read_t_file(read_file); + +f = open(map_file); +h = file(result_file, 'w') + +for aline in f.readlines(): + tline = aline.strip(); + tl = tline.split('\t'); + if len(tl)>4: + if int(tl[1].strip())== 0: + h.write(tline) + h.write('\n') + + +f.close(); +h.close() + + + +