diff wrapper.sh @ 22:d84c9791d8c4 draft

Uploaded
author davidvanzessen
date Tue, 07 Apr 2015 03:52:34 -0400
parents c9f9623f1f76
children 58a62d2c0377
line wrap: on
line diff
--- a/wrapper.sh	Thu Apr 02 03:31:23 2015 -0400
+++ b/wrapper.sh	Tue Apr 07 03:52:34 2015 -0400
@@ -6,6 +6,7 @@
 output=$3
 outdir=$4
 title=$5
+include_fr1=$6
 mkdir $outdir
 
 unzip $input -d $PWD/files/ > $PWD/unziplog.log
@@ -14,7 +15,7 @@
 cat $PWD/files/*/8_* > $PWD/mutationstats.txt
 cat $PWD/files/*/10_* > $PWD/hotspots.txt
 
-#BLASTN_DIR="/home/galaxy/tmp/blast/ncbi-blast-2.2.30+/bin"
+BLASTN_DIR="/home/galaxy/tmp/blast/ncbi-blast-2.2.30+/bin"
 
 echo "${BLASTN_DIR}"
 
@@ -46,7 +47,7 @@
 
 genes="ca,ca1,ca2,cg,cg1,cg2,cg3,cg4,cm"
 echo "R mutation analysis"
-Rscript $dir/mutation_analysis.r $outdir/merged.txt $genes $outdir 2>&1
+Rscript $dir/mutation_analysis.r $outdir/merged.txt $genes $outdir ${include_fr1} 2>&1
 echo "python mutation analysis"
 python $dir/mutation_analysis.py --input $outdir/merged.txt --genes $genes --output $outdir/hotspot_analysis.txt
 
@@ -82,6 +83,10 @@
 then
 	echo "<img src='cg.png'/><br />" >> $output
 fi
+if [ -a $outdir/scatter.png ]
+then
+	echo "<img src='scatter.png'/><br />" >> $output
+fi
 
 for gene in ${genes[@]}
 do