# HG changeset patch
# User xuebing
# Date 1333217149 14400
# Node ID 0f849231e4e0dab57c8254162ce2c27d001ee6b8
# Parent  21b4cb068fabf7b649e20d6ed028d83c3f9bde6c
Uploaded

diff -r 21b4cb068fab -r 0f849231e4e0 bed_window.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/bed_window.py	Sat Mar 31 14:05:49 2012 -0400
@@ -0,0 +1,19 @@
+def makeWindow(filename,outfile,window):
+    window = window/2
+    f=open(filename)
+    out = open(outfile,'w')
+    for line in f:
+        flds = line.strip().split()
+        #new position
+        center = (int(flds[1]) + int(flds[2]))/2
+        start = center - window
+        end = center + window
+        if start >= 0:
+            flds[1] = str(start)
+            flds[2] = str(end)
+            out.write('\t'.join(flds)+'\n')
+    f.close()
+    out.close()
+
+import sys
+makeWindow(sys.argv[1],sys.argv[2],int(sys.argv[3]))