Mercurial > repos > davidvanzessen > argalaxy_tools

view igblast/igblast.sh @ 6:e61d06912e86 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Thu, 30 Jul 2015 09:48:10 -0400
parents c7c7000de220
children c5cffbd8fe9d
line wrap: on
line source

set -e

dir="$(cd "$(dirname "$0")" && pwd)"

input=$1
species=$2
locus=$3
output=$4


echo "$input $species $locus $output"

java -jar igblastwrp.jar -cf -S $species -R $locus ${input} $PWD/blasted_output

Rscript --verbose $dir/igblast.r "$PWD/blasted_output.L2.txt" "$output" 2>&1