annotate igblast/igblast.sh @ 2:4f8a1d121384 draft

Uploaded
author davidvanzessen
date Mon, 05 Sep 2016 11:15:22 -0400
parents afe85eb6572e
children d001d0c05dbe
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
1 set -e
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
2
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
3 dir="$(cd "$(dirname "$0")" && pwd)"
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
4
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
5 input=$1
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
6 species=$2
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
7 locus=$3
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
8 output=$4
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
9
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
10
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
11 echo "$input $species $locus $output"
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
12
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
13 java -Xmx64G -jar $IGBLASTWRP/igblastwrp.jar -p 4 -S $species -R $locus ${input} $PWD/blasted_output 2>&1
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
14
afe85eb6572e Uploaded
davidvanzessen
parents:
diff changeset
15 Rscript --verbose $dir/igblast.r "$PWD/blasted_output.L2.txt" "$output" 2>&1