annotate falco-report.sh @ 34:b51c36920d7b draft

Uploaded
author stef
date Fri, 10 Oct 2014 10:13:55 -0400
parents 2210b89487c6
children f1b32c09deaf
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
25
2210b89487c6 Uploaded
stef
parents:
diff changeset
1 #!/bin/bash
2210b89487c6 Uploaded
stef
parents:
diff changeset
2 TOOLDIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"
2210b89487c6 Uploaded
stef
parents:
diff changeset
3
2210b89487c6 Uploaded
stef
parents:
diff changeset
4 echo "Start falco sh"
2210b89487c6 Uploaded
stef
parents:
diff changeset
5
2210b89487c6 Uploaded
stef
parents:
diff changeset
6 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
7 ## Variables setup
2210b89487c6 Uploaded
stef
parents:
diff changeset
8 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
9 #bam_file=$1
2210b89487c6 Uploaded
stef
parents:
diff changeset
10 #bam_name=$2
2210b89487c6 Uploaded
stef
parents:
diff changeset
11 #html_out=$3
2210b89487c6 Uploaded
stef
parents:
diff changeset
12 #vcf_out=$4
2210b89487c6 Uploaded
stef
parents:
diff changeset
13 #out_path=$5
2210b89487c6 Uploaded
stef
parents:
diff changeset
14 source $1
2210b89487c6 Uploaded
stef
parents:
diff changeset
15 JAR_PATH=$2
2210b89487c6 Uploaded
stef
parents:
diff changeset
16 #JAR_PATH=$3
2210b89487c6 Uploaded
stef
parents:
diff changeset
17
2210b89487c6 Uploaded
stef
parents:
diff changeset
18 bam_base=`echo $bam_name | sed 's#.bam$##' - `
2210b89487c6 Uploaded
stef
parents:
diff changeset
19
2210b89487c6 Uploaded
stef
parents:
diff changeset
20 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
21 ## Create extra init script
2210b89487c6 Uploaded
stef
parents:
diff changeset
22 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
23 #if [ $design_file != 'None' && $design_file != '' ]; then BED_PATH=$design_file; fi
2210b89487c6 Uploaded
stef
parents:
diff changeset
24 #if [ $manifest_file != 'None' && $manifest_file != '' ]; then MAN_PATH=$manifest_file; fi
2210b89487c6 Uploaded
stef
parents:
diff changeset
25 if [ $design_file != 'None' ] # Galaxy default is "None"
2210b89487c6 Uploaded
stef
parents:
diff changeset
26 then
2210b89487c6 Uploaded
stef
parents:
diff changeset
27 BED_PATH=$design_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
28 else
2210b89487c6 Uploaded
stef
parents:
diff changeset
29 BED_PATH="$TOOLDIR/tool-data/GALAXY_TSACP.bed"
2210b89487c6 Uploaded
stef
parents:
diff changeset
30 fi
2210b89487c6 Uploaded
stef
parents:
diff changeset
31
2210b89487c6 Uploaded
stef
parents:
diff changeset
32 if [ $manifest_file != 'None' ] # Galaxy default is "None"
2210b89487c6 Uploaded
stef
parents:
diff changeset
33 then
2210b89487c6 Uploaded
stef
parents:
diff changeset
34 MAN_PATH=$manifest_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
35 else
2210b89487c6 Uploaded
stef
parents:
diff changeset
36 MAN_PATH="$TOOLDIR/tool-data/GALAXY_TruSeq_Amplicon_Cancer_Panel_Manifest_AFP1_PN15032433.txt"
2210b89487c6 Uploaded
stef
parents:
diff changeset
37 fi
2210b89487c6 Uploaded
stef
parents:
diff changeset
38
2210b89487c6 Uploaded
stef
parents:
diff changeset
39 echo "BED:"$design_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
40 echo "MAN:"$manifest_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
41 echo "BED chosen:"$BED_PATH
2210b89487c6 Uploaded
stef
parents:
diff changeset
42 echo "MAN chosen:"$MAN_PATH
2210b89487c6 Uploaded
stef
parents:
diff changeset
43 #echo "VARSCAN_PATH1:"$SNPEFF_JAR_PATH
2210b89487c6 Uploaded
stef
parents:
diff changeset
44 #echo "VARSCAN_PATH2:"\$SNPEFF_JAR_PATH
2210b89487c6 Uploaded
stef
parents:
diff changeset
45 echo "JAR_PATH1:"$JAR_PATH
2210b89487c6 Uploaded
stef
parents:
diff changeset
46 #echo "JAR_PATH2:"$jar_path2
2210b89487c6 Uploaded
stef
parents:
diff changeset
47
2210b89487c6 Uploaded
stef
parents:
diff changeset
48 init_file="galaxy_init.sh"
2210b89487c6 Uploaded
stef
parents:
diff changeset
49 echo '#!/bin/bash' > $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
50
2210b89487c6 Uploaded
stef
parents:
diff changeset
51 ## executables
2210b89487c6 Uploaded
stef
parents:
diff changeset
52 #echo "snpEff=/ccagc/lib/snpEff/snpEff_v3_4" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
53 echo "snpEff=\$SNPEFF_JAR_PATH/" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
54 #echo "snpSift=/ccagc/lib/SnpSift" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
55 echo "snpSift=\$SNPEFF_JAR_PATH/" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
56 #echo "varscan=/ccagc/lib/VarScan/VarScan.v2.2.11.jar" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
57 echo "varscan=\$JAVA_JAR_PATH/VarScan.v2.3.5.jar" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
58
2210b89487c6 Uploaded
stef
parents:
diff changeset
59 ## files
2210b89487c6 Uploaded
stef
parents:
diff changeset
60 echo "targetBed=$BED_PATH" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
61 echo "manifest=$MAN_PATH" >> $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
62
2210b89487c6 Uploaded
stef
parents:
diff changeset
63 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
64 ## Status / debug
2210b89487c6 Uploaded
stef
parents:
diff changeset
65 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
66 echo "BAM: "$bam_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
67 echo "OUT: "$bam_name
2210b89487c6 Uploaded
stef
parents:
diff changeset
68 echo "BAMBASE: "$bam_base
2210b89487c6 Uploaded
stef
parents:
diff changeset
69 echo "OUTPATH: "$out_path
2210b89487c6 Uploaded
stef
parents:
diff changeset
70 echo "HTMLOUT: "$html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
71 echo "VCFOUT : "$vcf_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
72 echo "INIT : "$init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
73
2210b89487c6 Uploaded
stef
parents:
diff changeset
74 source $init_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
75
2210b89487c6 Uploaded
stef
parents:
diff changeset
76 echo "--- sam ---"
2210b89487c6 Uploaded
stef
parents:
diff changeset
77 samtools
2210b89487c6 Uploaded
stef
parents:
diff changeset
78 echo "--- bcf ---"
2210b89487c6 Uploaded
stef
parents:
diff changeset
79 bcftools
2210b89487c6 Uploaded
stef
parents:
diff changeset
80 echo "--- snpsift ---"
2210b89487c6 Uploaded
stef
parents:
diff changeset
81 java -jar $snpSift/SnpSift.jar
2210b89487c6 Uploaded
stef
parents:
diff changeset
82 echo "--- snpeff ---"
2210b89487c6 Uploaded
stef
parents:
diff changeset
83 java -jar $snpEff/snpeff.jar
2210b89487c6 Uploaded
stef
parents:
diff changeset
84 echo "--- varscan ---"
2210b89487c6 Uploaded
stef
parents:
diff changeset
85 java -jar $varScan
2210b89487c6 Uploaded
stef
parents:
diff changeset
86
2210b89487c6 Uploaded
stef
parents:
diff changeset
87 exit 0
2210b89487c6 Uploaded
stef
parents:
diff changeset
88
2210b89487c6 Uploaded
stef
parents:
diff changeset
89 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
90 ## create output files dir
2210b89487c6 Uploaded
stef
parents:
diff changeset
91 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
92 mkdir $out_path
2210b89487c6 Uploaded
stef
parents:
diff changeset
93 mkdir $out_path/$bam_base
2210b89487c6 Uploaded
stef
parents:
diff changeset
94
2210b89487c6 Uploaded
stef
parents:
diff changeset
95 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
96 ## running analysis
2210b89487c6 Uploaded
stef
parents:
diff changeset
97 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
98 echo "Starting VcAn"
2210b89487c6 Uploaded
stef
parents:
diff changeset
99 $TOOLDIR/falco/amplicon-VcAn-TSACP.sh $bam_file $bam_base
2210b89487c6 Uploaded
stef
parents:
diff changeset
100 #$TOOLDIR/falco/amplicon-AlVcAn-TSACP.sh ## this contains more
2210b89487c6 Uploaded
stef
parents:
diff changeset
101
2210b89487c6 Uploaded
stef
parents:
diff changeset
102 echo "Starting amplicon-Report"
2210b89487c6 Uploaded
stef
parents:
diff changeset
103 $TOOLDIR/falco/amplicon-Report.sh $bam_base
2210b89487c6 Uploaded
stef
parents:
diff changeset
104
2210b89487c6 Uploaded
stef
parents:
diff changeset
105 echo "Starting QC"
2210b89487c6 Uploaded
stef
parents:
diff changeset
106 $TOOLDIR/falco/amplicon-QC.sh $bam_file
2210b89487c6 Uploaded
stef
parents:
diff changeset
107
2210b89487c6 Uploaded
stef
parents:
diff changeset
108 #echo "Starting plots.R"
2210b89487c6 Uploaded
stef
parents:
diff changeset
109 #$TOOLDIR/falco/misc/plots.R $bam_base.qc.ann.txt $bam_base.qc.targets.txt $bam_base.pdf
2210b89487c6 Uploaded
stef
parents:
diff changeset
110
2210b89487c6 Uploaded
stef
parents:
diff changeset
111 echo "Starting mkHtmlReport"
2210b89487c6 Uploaded
stef
parents:
diff changeset
112 $TOOLDIR/falco/misc/mkHtmlReport.pl ./ ./
2210b89487c6 Uploaded
stef
parents:
diff changeset
113
2210b89487c6 Uploaded
stef
parents:
diff changeset
114
2210b89487c6 Uploaded
stef
parents:
diff changeset
115
2210b89487c6 Uploaded
stef
parents:
diff changeset
116 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
117 ## create index html for main galaxy output
2210b89487c6 Uploaded
stef
parents:
diff changeset
118 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
119 echo "<!DOCTYPE html>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
120 echo "<html>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
121 echo "<head>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
122 echo "<style>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
123 echo " body{ padding: 0px 20px; }" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
124 echo " h1{ color: red; }" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
125 echo " table{ border: 1px solid black; padding: 5px }" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
126 echo "</style>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
127 echo "</head>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
128 echo "<body>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
129 echo " <h1>FALCO</h1>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
130 echo " <a href=\"index.html\">HTML</a>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
131 echo " <table><tbody>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
132 #for file in *.vcf *.txt *.tsv
2210b89487c6 Uploaded
stef
parents:
diff changeset
133 for file in *
2210b89487c6 Uploaded
stef
parents:
diff changeset
134 do
2210b89487c6 Uploaded
stef
parents:
diff changeset
135 lineCount=`wc -l $file | cut -f 1 -d " "`
2210b89487c6 Uploaded
stef
parents:
diff changeset
136 echo " <tr><td><a href=\"$file\">$file</a> has $lineCount lines</td></tr>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
137 #echo " <tr><td> HEAD: " `head $file` "</td></tr>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
138 done
2210b89487c6 Uploaded
stef
parents:
diff changeset
139 echo " </tbody></table>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
140 echo "</body>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
141 echo "</html>" >> $html_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
142
2210b89487c6 Uploaded
stef
parents:
diff changeset
143 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
144 ## creating galaxy history outputs
2210b89487c6 Uploaded
stef
parents:
diff changeset
145 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
146 #cp 'index.html' $html_out # this is the overview of samples html
2210b89487c6 Uploaded
stef
parents:
diff changeset
147 cp $bam_base'.html' $out_path/'out.html' # this is the sample html
2210b89487c6 Uploaded
stef
parents:
diff changeset
148 cp $bam_base'.vcf' $vcf_out
2210b89487c6 Uploaded
stef
parents:
diff changeset
149
2210b89487c6 Uploaded
stef
parents:
diff changeset
150 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
151 ## copy files to keep to output path
2210b89487c6 Uploaded
stef
parents:
diff changeset
152 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
153 cp -r ./$bam_base/*png $out_path/$bam_base/
2210b89487c6 Uploaded
stef
parents:
diff changeset
154 cp -r ./* $out_path
2210b89487c6 Uploaded
stef
parents:
diff changeset
155 #cp *.vcf $out_path; cp *.tsv $out_path; cp *.txt $out_path; cp *.html $out_path
2210b89487c6 Uploaded
stef
parents:
diff changeset
156
2210b89487c6 Uploaded
stef
parents:
diff changeset
157 ## ----------
2210b89487c6 Uploaded
stef
parents:
diff changeset
158 echo "END falco sh"
2210b89487c6 Uploaded
stef
parents:
diff changeset
159 exit 0