annotate vidjil.sh @ 6:7346cd210a43 draft

Uploaded
author davidvanzessen
date Thu, 12 Feb 2015 09:54:07 -0500
parents fe4001bdf8d4
children e55dfd9a8c00
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
1 #!/bin/bash
6
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
2 set -e
0
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
3
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
4 inputFile=$1
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
5 outputFile=$2
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
6 outputDir=$3
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
7 database=$4
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
8 limit_reads=$5
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
9 limit_ratio=$6
6
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
10 selection=$7
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
11 germline=$8
0
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
12
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
13 dir="$(cd "$(dirname "$0")" && pwd)"
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
14 mkdir $outputDir
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
15
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
16 unzip $dir/browser.zip -d $outputDir/ > $outputDir/unziplog.log
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
17
6
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
18 mkdir $outputDir/data
0
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
19
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
20 if [[ "$database" == "download" ]] ; then
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
21 echo "Downloading up to date database from IMGT..."
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
22 mkdir $outputDir/tmp
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
23 mkdir $outputDir/tmp/germline
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
24 cp $dir/germline/get-germline $outputDir/tmp/germline
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
25 cp $dir/germline/split-from-imgt.py $outputDir/tmp/germline
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
26 cp $dir/vidjil $outputDir/tmp
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
27 cd $outputDir/tmp/germline
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
28 ./get-germline 2>&1
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
29 cd $outputDir/tmp
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
30 else
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
31 cd $dir
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
32 fi
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
33
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
34
6
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
35 inputFiles=$(echo "$inputFile" | tr "| " "\n")
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
36 files_for_config=""
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
37 last_file=""
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
38 for id_file in $inputFiles
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
39 do
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
40 tmp=($(echo "${id_file}" | tr ";" "\n"))
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
41 id=${tmp[0]}
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
42 file=${tmp[1]}
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
43
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
44 ./vidjil -G germline/$germline -y all -r ${limit_reads} -% ${limit_ratio} -b $id -o $outputDir/data $file
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
45
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
46 if [[ "$files_for_config" != "" ]] ; then
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
47 files_for_config="${files_for_config},"
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
48 fi
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
49 files_for_config="${files_for_config}'${id}.vidjil'"
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
50 last_file="${id}.vidjil"
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
51 done
0
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
52
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
53 echo "<html><center><h1><a href='index.html'>Click here for the results</a></h1>Tip: Open it in a new tab (middle mouse button or right mouse button -> 'open in new tab' on the link above)</center></html>" > $outputFile
fe4001bdf8d4 Uploaded
davidvanzessen
parents:
diff changeset
54
6
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
55 echo -e \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
56 "var config = {\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
57 " 'file_menu' : {\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
58 " 'path' : 'data/',\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
59 " 'file' : [\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
60 " ${files_for_config}\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
61 " ]\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
62 " },\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
63 " 'autoload' : 'data/${last_file}'\n" \
7346cd210a43 Uploaded
davidvanzessen
parents: 0
diff changeset
64 "}" > $outputDir/js/conf.js;