Mercurial > repos > xuebing > sharplab_seq_motif
annotate mytools/makewindow.py @ 0:39217fa39ff2
Uploaded
| author | xuebing |
|---|---|
| date | Tue, 13 Mar 2012 23:34:52 -0400 |
| parents | |
| children |
| rev | line source |
|---|---|
| 0 | 1 def makeWindow(filename,outfile,window): |
| 2 window = window/2 | |
| 3 f=open(filename) | |
| 4 out = open(outfile,'w') | |
| 5 for line in f: | |
| 6 flds = line.strip().split() | |
| 7 #new position | |
| 8 center = (int(flds[1]) + int(flds[2]))/2 | |
| 9 start = center - window | |
| 10 end = center + window | |
| 11 if start >= 0: | |
| 12 flds[1] = str(start) | |
| 13 flds[2] = str(end) | |
| 14 out.write('\t'.join(flds)+'\n') | |
| 15 f.close() | |
| 16 out.close() | |
| 17 | |
| 18 import sys | |
| 19 makeWindow(sys.argv[1],sys.argv[2],int(sys.argv[3])) |