Mercurial > repos > bcclaywell > microbiome_pplacer_suite
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}