-
Notifications
You must be signed in to change notification settings - Fork 0
/
complete.sh
72 lines (62 loc) · 2.06 KB
/
complete.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
#!/bin/bash
set -e
inputFiles=($1) # Input files can be multiple files
outputDir=$3
outputFile=$3/index.html #$1
clonalType=$4
species=$5
locus=$6
filterproductive=$7
clonality_method=$8
html=$2
dir="$(cd "$(dirname "$0")" && pwd)"
array=("$@")
echo "<html><h3>Progress</h3><table><tr><td>info</td></tr>" > "$html"
echo "<tr><td>-----------------------------------</td></tr>" >> "$html"
#mkdir $PWD/igblastdatabase
#unzip $dir/database.zip -d $PWD/igblastdatabase/
#export IGDATA=$PWD/igblastdatabase/
echo "python: $(which python)"
echo "R: $(which R)"
echo "Rscript: $(which Rscript)"
id=""
forwardSlash="/"
mergerInput=()
echo "Before loop"
count=1
for current in "${inputFiles[@]}"
do
if [[ "$current" != *"$forwardSlash"* ]]; then
id="$current"
mergerInput+=("$id")
count=1
continue
fi
echo "working on $current"
fileName=$(basename "$current")
fileName="${fileName%.*}"
parsedFileName="$PWD/$fileName.parsed"
f=$(file "$current")
zipType="Zip archive"
zxType="XZ compressed data"
echo "filetype of ${id}: $f"
if [[ "$f" == *"$zipType"* ]] || [[ "$f" == *"$zxType"* ]]
then
echo "<tr><td>Sample $count of patient $id is an archive file, using IMGT Loader</td></tr>" >> $html
fileName=$(basename $current)
bash ${dir}/imgt_loader/imgt_loader.sh $current $parsedFileName "${fileName}"
else
echo "<tr><td>Sample $count of patient $id is not a zip file so assuming fasta/fastq, using igBLASTn</td></tr>" >> $html
bash ${dir}/igblast/igblast.sh $current "$species" $locus $parsedFileName
fi
mergerInput+=("$parsedFileName")
count=$((count+1))
done
echo "<tr><td>-----------------------------------</td></tr>" >> $html
echo "<tr><td>merging</td></tr>" >> $html
bash $dir/experimental_design/experimental_design.sh ${mergerInput[*]} $PWD/merged.txt
echo "<tr><td>done</td></tr>" >> $html
echo "<tr><td>-----------------------------------</td></tr>" >> $html
echo "<tr><td>plotting</td></tr>" >> $html
echo "after ED"
bash $dir/report_clonality/r_wrapper.sh $PWD/merged.txt $2 $outputDir $clonalType "$species" "$locus" $filterproductive $clonality_method