--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/median_size.py	Fri Mar 10 13:22:58 2023 +0000
@@ -0,0 +1,26 @@
+#!/usr/bin/env
+
+## Errol Strain (estrain@gmail.com)
+## calculate median insert size from sam file
+
+import numpy as np
+
+def get_data(infile):
+    lengths = []
+    for line in infile:
+        if line.startswith('@'):
+            pass
+        else:
+            line = line.rsplit()
+            length = int(line[8])
+            if length > 0:
+                lengths.append(length)
+            else:
+                pass
+    return lengths
+
+if __name__ == "__main__":
+    import sys
+    lengths = get_data(sys.stdin)
+    md = int(np.median(lengths))
+print(md)