# HG changeset patch
# User eugen
# Date 1345035194 14400
# Node ID 120a27b53be7ed7c4c8e21f929ed20dd1c019e8e
# Parent  ee3bd6ad3958a0042cc583ee19ce0fc986a93683
Uploaded

diff -r ee3bd6ad3958 -r 120a27b53be7 bismark_meth_caller.sh
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/bismark_meth_caller.sh	Wed Aug 15 08:53:14 2012 -0400
@@ -0,0 +1,35 @@
+#!/bin/bash
+#
+# Galaxy wrapper for Bismark Methylation Caller
+#
+
+set -e
+
+#get parameters
+
+until [ $# -eq 0 ]
+do
+	case $1 in
+		input=*)
+			input=${1#input=}
+			;;
+		method=*)
+			method=${1#method=}
+			;;
+		output=*)
+			output=${1#output=}
+			;;
+		tempdir=*)
+			tempdir=${1#tempdir=}
+			;;
+	esac
+	shift
+done
+(grep ^@ $input; grep -v ^@ $input | sort -k3,3 -k1,1) > $tempdir"bismark_sam_sorted.sam"
+methylation_extractor -$method -o $output $tempdir"bismark_sam_sorted.sam" 2>&1 > /dev/null
+
+
+
+
+
+