Mercurial > repos > davidvanzessen > argalaxy_tools
diff igblast/igblast.sh @ 53:0ef7f80ea061 draft
Uploaded
author | davidvanzessen |
---|---|
date | Thu, 25 Feb 2021 13:39:50 +0000 |
parents | 124b7fd92a3e |
children |
line wrap: on
line diff
--- a/igblast/igblast.sh Thu Feb 25 13:36:15 2021 +0000 +++ b/igblast/igblast.sh Thu Feb 25 13:39:50 2021 +0000 @@ -1,28 +1,28 @@ -set -e - -dir="$(cd "$(dirname "$0")" && pwd)" - -input=$1 -species=$2 -locus=$3 -output=$4 - -declare -A speciesdict - -speciesdict=(["Rattus norvegicus functional"]="rat" ["Rattus norvegicus non-functional"]="rat" ["Oryctolagus cuniculus functional"]="rabbit" ["Oryctolagus cuniculus non-functional"]="rabbit" ["Mus musculus functional"]="mouse" ["Mus musculus non-functional"]="mouse" ["Homo sapiens functional"]="human" ["Homo sapiens non-functional"]="human" ["Macaca mulatta non-functional"]="rhesus_monkey" ["Macaca mulatta functional"]="rhesus_monkey") - -echo "Species: $species ${speciesdict[$species]}" - -species="${speciesdict[$species]}" - -if [ "$species" == "" ] -then - >&2 echo "Species not possible with igBLASTn, use IMGT" - exit 1 -fi - -echo "$input $species $locus $output" - -java -Xmx16G -jar $IGBLASTWRP/igblastwrp.jar -p 4 -S $species -R $locus ${input} $PWD/blasted_output 2>&1 - -Rscript --verbose $dir/igblast.r "$PWD/blasted_output.L2.txt" "$output" 2>&1 +set -e + +dir="$(cd "$(dirname "$0")" && pwd)" + +input=$1 +species=$2 +locus=$3 +output=$4 + +declare -A speciesdict + +speciesdict=(["Rattus norvegicus functional"]="rat" ["Rattus norvegicus non-functional"]="rat" ["Oryctolagus cuniculus functional"]="rabbit" ["Oryctolagus cuniculus non-functional"]="rabbit" ["Mus musculus functional"]="mouse" ["Mus musculus non-functional"]="mouse" ["Homo sapiens functional"]="human" ["Homo sapiens non-functional"]="human" ["Macaca mulatta non-functional"]="rhesus_monkey" ["Macaca mulatta functional"]="rhesus_monkey") + +echo "Species: $species ${speciesdict[$species]}" + +species="${speciesdict[$species]}" + +if [ "$species" == "" ] +then + >&2 echo "Species not possible with igBLASTn, use IMGT" + exit 1 +fi + +echo "$input $species $locus $output" + +java -Xmx16G -jar $IGBLASTWRP/igblastwrp.jar -p 4 -S $species -R $locus ${input} $PWD/blasted_output 2>&1 + +Rscript --verbose $dir/igblast.r "$PWD/blasted_output.L2.txt" "$output" 2>&1