Repository 'bed_window'
hg clone https://toolshed.g2.bx.psu.edu/repos/xuebing/bed_window

Changeset 1:0f849231e4e0 (2012-03-31)
Previous changeset 0:21b4cb068fab (2012-03-31)
Commit message:
Uploaded
added:
bed_window.py
b
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]))