diff single_fastq_filtering_wrapper.sh @ 0:a4cd8608ef6b draft

Uploaded

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/single_fastq_filtering_wrapper.sh	Mon Apr 01 07:56:36 2019 -0400
@@ -0,0 +1,59 @@
+#!/bin/bash
+
+#set -euo pipefail # this is something like string in perl
+IFS=$'\n\t'
+
+# run filtering
+WD="`dirname $0`"
+ORIDIR=$PWD
+cd $WD
+WD=$PWD  # absolute path to this script
+cd $ORIDIR
+SAMPLING=""
+TRIM_END=""
+TRIM_START=""
+PERCENT_ABOVE="95"
+CUTADAPT=""
+FILTER_SEQ=""
+while getopts "a:o:n:c:p:e:s:N:C:G:F:" OPTION
+do
+    case $OPTION in
+        a)
+            FASTAA=$OPTARG;;
+        o)
+            OUTPUT=$OPTARG;;
+        n)
+            SAMPLING=( -n ${OPTARG} );;
+        c)
+            CUT_OFF=$OPTARG;;
+        p)
+            PERCENT_ABOVE=$OPTARG;;
+        e)
+            TRIM_END=( -e ${OPTARG} );;
+        s)
+            TRIM_START=( -s ${OPTARG} );;
+        N)
+            MAX_N=${OPTARG};;
+        C)
+            CUTADAPT=(-C " "${OPTARG}" " );;
+        G)
+            PNG_OUTPUT=${OPTARG};;
+        F)
+            FILTER_SEQ=( -F ${OPTARG} );;
+    esac
+done
+
+
+
+if [ -z "$CUTADAPT" ] # test if $CUTADAPT is empty
+then
+    ${WD}/single_fastq_filtering.R -a $FASTAA -x $OUTPUT  ${SAMPLING[@]} -c $CUT_OFF\
+         -p $PERCENT_ABOVE  ${TRIM_START[@]}  ${TRIM_END[@]} -N $MAX_N -G $PNG_OUTPUT ${FILTER_SEQ[@]}
+else
+    ${WD}/single_fastq_filtering.R -a $FASTAA -x $OUTPUT  ${SAMPLING[@]} -c $CUT_OFF -G $PNG_OUTPUT\
+         -p $PERCENT_ABOVE  ${TRIM_START[@]}  ${TRIM_END[@]} -N $MAX_N "${CUTADAPT[@]}" ${FILTER_SEQ[@]}
+fi
+
+
+
+