# HG changeset patch
# User ziru-zhou
# Date 1355840412 18000
# Node ID e623e9c7bb23e67462fde6a5852d514b1e6c5af7
# Parent  5f4cb6ee3427e298db2beaa865ba49a0e31473a2
Uploaded

diff -r 5f4cb6ee3427 -r e623e9c7bb23 split.sh
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/split.sh	Tue Dec 18 09:20:12 2012 -0500
@@ -0,0 +1,25 @@
+#!bin/bash
+
+#purpose: shell script which performs the split option
+#author: Ziru Zhou
+#date: October, 2012
+
+main()
+{
+	output1="${1}"
+	output2="${2}"
+	input="${3}"
+
+	$(samtools view -H ${input} > tmp)
+	linecount=$(samtools view ${input} | wc -l)
+	half=$(((${linecount} + 1) / 2))
+	samtools view ${input} | shuf | split -d -l ${half}
+
+	$(cat tmp x00 > x00h)
+	$(cat tmp x01 > x01h)
+        samtools view -bSo ${output1} x00h
+        samtools view -bSo ${output2} x01h
+
+	rm tmp x00h x01h x00 x01
+}
+main "${@}"