diff filter-wrapper.sh @ 2:ce6db18f5fd3 draft

Uploaded
author bcclaywell
date Thu, 26 Feb 2015 19:31:20 -0500
parents
children 2d023c621bd0
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/filter-wrapper.sh	Thu Feb 26 19:31:20 2015 -0500
@@ -0,0 +1,32 @@
+#!/bin/bash
+
+source $(dirname $0)/util.sh
+source $1
+
+INPUT_QUAL=$(extify qual ${INPUT_QUAL})
+BARCODES=$(extify csv ${BARCODES})
+RAW_SEQS=$(extify fasta ${RAW_SEQS})
+
+seqmagick quality-filter \
+    --input-qual ${INPUT_QUAL} \
+    --barcode-file ${BARCODES} \
+    --primer "${PRIMER}" \
+    --report-out ${FILTER_REPORT} \
+    --details-out ${FILTER_DETAILS} \
+    --map-out ${SPLIT_MAP} \
+    --barcode-header \
+    --min-length ${MIN_LENGTH} \
+    --min-mean-quality ${MIN_QUALITY} \
+    --quality-window 30 \
+    --quality-window-prop 0.9 \
+    --quality-window-mean-qual 15 \
+    ${RAW_SEQS} \
+    filtered.fasta
+
+if [[ ${REVERSE_COMPLEMENT} == "TRUE" ]]; then
+    seqmagick mogrify \
+        --reverse-complement \
+        filtered.fasta
+fi
+
+mv filtered.fasta ${FILTERED_SEQS}