--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/MUMmer/mummer_tool.sh	Tue Jun 07 17:22:27 2011 -0400
@@ -0,0 +1,107 @@
+#!/bin/bash
+## use #!/bin/bash -x for debugging
+
+## Galaxy wrapper for MUMmer (nucmer/promer)
+## Alex Bossers, CVI of Wageningen UR, NL
+## alex_dot_bossers_at_wur_dot_nl
+##
+## Sep 2010
+##
+## Wrapper runs MUMmer nucmer/promer and additional args
+## Calculates the comparison scores (delta and optional coords file)
+## Generates the optional STATIC comparison mummerplot to png (from delta file)
+##
+## finally the script renames (optional) output files to outfiles expected by Galaxy
+##
+##
+## INPUT args:
+## nucmer_tool.sh $input_ref $input_query $out_delta $out_coords $out_png $logfile 
+##                    @0          @1          @2          @3        @4       @5   
+##                $algorithm $keep_delta $make_coords $keep_log $make_image $cmd_extra 
+##                     @6        @7           @8          @9        @10         @11
+##
+
+# path to where mummer suite is installed
+# adjust this for your machine
+# this is the only hard coded path in the scripts
+mum_path="/opt/MUMmer/MUMmer"
+
+# since we have more than 9 arguments we need to shift the sections or use own array
+args=("$@")
+# to keep things readible assign vars
+input_ref="${args[0]}"
+input_query="${args[1]}"
+out_delta="${args[2]}"
+out_coords="${args[3]}"
+out_png="${args[4]}"
+logfile="${args[5]}"
+algorithm="${args[6]}"
+keep_delta="${args[7]}"
+make_coords="${args[8]}"
+keep_log="${args[9]}"
+make_image="${args[10]}"
+cmd_extra="${args[11]}"
+
+# enable/disable the STDOUT log file
+if [ "$keep_log" == "yes" ]; then
+	logfile_c="2>$logfile"
+	logfile_a="2>>$logfile"
+else
+	#dump to dev/null
+	logfile_c="2>&-"
+	logfile_a="2>&-"
+fi
+
+# extra mummer cmd line options
+
+## generate coords file on the fly?
+if [ "$make_coords" == "yes" ]; then
+	options=" --coords"
+fi
+## extra cmd line args to be concatenated in options? We need to prevent extra spaces!
+if [ "$cmd_extra" != "" ]; then
+	if [ "$options" == "" ]; then
+		options=" $cmd_extra"
+	else
+		options="$options $cmd_extra"
+	fi
+fi
+
+# run nucmer/promer
+eval "$mum_path/$algorithm$options $input_ref $input_query $logfile_c"
+
+## generate large png if option make_image = yes
+## suppress error from mummerplot since some is deprecated but not a real error
+## error can be easily avoided by modifying the source of mummerplot... just in case
+## however we need to check if a valid png was generated. This is not the case is alignment is none
+## 1 is stderr and 2 stdout. redirect to dev/null
+
+if [ "${make_image}" == "yes" ]; then
+	eval "$mum_path/mummerplot --large --png out.delta 1>&- $logfile_a"
+	if [ -f "out.png" ]; then
+		mv out.png $out_png
+		#cleanup temp gnuplot file
+		rm out.gp
+	else
+		echo "not exist the req png file!"
+	fi
+	
+	## clean up remaining files
+	rm out.fplot
+	rm out.rplot
+
+fi
+
+# keep/rename or delete delta file
+if [ "$keep_delta" == "yes" ]; then
+	mv out.delta "$out_delta" 
+else
+	rm out.delta
+fi
+
+# keep/rename coords file if it was created
+if [ "$make_coords" == "yes" ]; then
+	mv out.coords "$out_coords"
+fi
+
+# end script