comparison wrapper.sh @ 1:75853bceec00 draft

Uploaded
author davidvanzessen
date Tue, 17 Jan 2017 07:24:44 -0500
parents ed6885c85660
children
comparison
equal deleted inserted replaced
0:ed6885c85660 1:75853bceec00
44 singles+=( "$patient" ) 44 singles+=( "$patient" )
45 fi 45 fi
46 oldLocus="" 46 oldLocus=""
47 sample1="$(echo ${sample1} | tr -d '\r' | tr -d '\n')" 47 sample1="$(echo ${sample1} | tr -d '\r' | tr -d '\n')"
48 sample2="$(echo ${sample2} | tr -d '\r' | tr -d '\n')" 48 sample2="$(echo ${sample2} | tr -d '\r' | tr -d '\n')"
49 tail -n+2 ${patient}_freq.txt | sed "s/>//" > tmp.txt 49 tail -n+2 "${patient}_freq.txt" | sed "s/>//" > tmp.txt
50 echo "<div class='tabber'>" >> "$html" 50 echo "<div class='tabber'>" >> "$html"
51 echo "<div class='tabbertab' title='Data frequency'>" >> "$html" 51 echo "<div class='tabbertab' title='Data frequency'>" >> "$html"
52 echo "<table><tr><td style='vertical-align:top;'>" >> "$html" 52 echo "<table><tr><td style='vertical-align:top;'>" >> "$html"
53 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_freq'>" >> "$html" 53 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_freq'>" >> "$html"
54 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_Both</th><th>Number of sequences_$sample1</th><th>Read Count $sample1</th><th>Number of sequences_$sample2</th><th>Read Count $sample2</th><th>Sum number of sequences $patient</th><th>Percentage of sequences ${patient}_both</th></thead>" >> "$html" 54 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_Both</th><th>Number of sequences_$sample1</th><th>Read Count $sample1</th><th>Number of sequences_$sample2</th><th>Read Count $sample2</th><th>Sum number of sequences $patient</th><th>Percentage of sequences ${patient}_both</th></thead>" >> "$html"
97 echo "<a href='${patient}_freq.png'><img src='${patient}_freq.png' width='1280' height='720' /></a><br />" >> "$html" 97 echo "<a href='${patient}_freq.png'><img src='${patient}_freq.png' width='1280' height='720' /></a><br />" >> "$html"
98 echo "<a href='${patient}_freq_both.png'><img src='${patient}_freq_both.png' width='1280' height='720' /></a><br />" >> "$html" 98 echo "<a href='${patient}_freq_both.png'><img src='${patient}_freq_both.png' width='1280' height='720' /></a><br />" >> "$html"
99 echo "<a href='${patient}_percent_freq.png'><img src='${patient}_percent_freq.png' width='1280' height='720' /></a></div>" >> "$html" 99 echo "<a href='${patient}_percent_freq.png'><img src='${patient}_percent_freq.png' width='1280' height='720' /></a></div>" >> "$html"
100 echo "${scatterplot_tab}</tr></table></div>" >> "$html" 100 echo "${scatterplot_tab}</tr></table></div>" >> "$html"
101 101
102 tail -n+2 ${patient}_reads.txt | sed "s/>//" > tmp.txt 102 tail -n+2 "${patient}_reads.txt" | sed "s/>//" > tmp.txt
103 echo "<div class='tabbertab' title='Data reads'>" >> "$html" 103 echo "<div class='tabbertab' title='Data reads'>" >> "$html"
104 echo "<table><tr><td style='vertical-align:top;'>" >> "$html" 104 echo "<table><tr><td style='vertical-align:top;'>" >> "$html"
105 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_reads'>" >> "$html" 105 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_reads'>" >> "$html"
106 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_Both</th><th>Number of sequences_$sample1</th><th>Read Count $sample1</th><th>Number of sequences_$sample2</th><th>Read Count $sample2</th><th>Sum number of sequences $patient</th><th>Percentage of sequences ${patient}_both</th></thead>" >> "$html" 106 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_Both</th><th>Number of sequences_$sample1</th><th>Read Count $sample1</th><th>Number of sequences_$sample2</th><th>Read Count $sample2</th><th>Sum number of sequences $patient</th><th>Percentage of sequences ${patient}_both</th></thead>" >> "$html"
107 echo "<tbody>" >> "$html" 107 echo "<tbody>" >> "$html"
187 sample3="$(echo ${sample3} | tr -d '\r' | tr -d '\n')" 187 sample3="$(echo ${sample3} | tr -d '\r' | tr -d '\n')"
188 patient="${sample1}_${sample2}_${sample3}" 188 patient="${sample1}_${sample2}_${sample3}"
189 echo "$patient" 189 echo "$patient"
190 html="${patient}.html" 190 html="${patient}.html"
191 echo "<tr><td><a href='${patient}.html'>$patient</a></td></tr>" >> "index.html" 191 echo "<tr><td><a href='${patient}.html'>$patient</a></td></tr>" >> "index.html"
192 echo "$header" > $html 192 echo "$header" > "$html"
193 oldLocus="" 193 oldLocus=""
194 tail -n+2 ${patient}_freq.txt | sed "s/>//" > tmp.txt 194 tail -n+2 "${patient}_freq.txt" | sed "s/>//" > tmp.txt
195 echo "<div class='tabber'>" >> "$html" 195 echo "<div class='tabber'>" >> "$html"
196 echo "<div class='tabbertab' title='Data frequency'>" >> "$html" 196 echo "<div class='tabbertab' title='Data frequency'>" >> "$html"
197 echo "<table><tr><td style='vertical-align:top;'>" >> "$html" 197 echo "<table><tr><td style='vertical-align:top;'>" >> "$html"
198 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_freq'>" >> "$html" 198 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_freq'>" >> "$html"
199 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_All</th><th>Number of sequences_$sample1</th><th>Number of sequences_$sample2</th><th>Number of sequences_$sample3</th><th>Number of sequences_${sample1}_${sample2}</th><th>Number of sequences_${sample1}_${sample3}</th><th>Number of sequences_${sample2}_${sample3}</th></thead>" >> "$html" 199 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_All</th><th>Number of sequences_$sample1</th><th>Number of sequences_$sample2</th><th>Number of sequences_$sample3</th><th>Number of sequences_${sample1}_${sample2}</th><th>Number of sequences_${sample1}_${sample3}</th><th>Number of sequences_${sample2}_${sample3}</th></thead>" >> "$html"
259 echo "<div class='tabbertab' title='Graphs frequency'>" >> "$html" 259 echo "<div class='tabbertab' title='Graphs frequency'>" >> "$html"
260 echo "<a href='${patient}_freq_total_all.png'><img src='${patient}_freq_total_all.png' width='1280' height='720' /></a><br />" >> "$html" 260 echo "<a href='${patient}_freq_total_all.png'><img src='${patient}_freq_total_all.png' width='1280' height='720' /></a><br />" >> "$html"
261 echo "<a href='${patient}_freq_indiv_all.png'><img src='${patient}_freq_indiv_all.png' width='1280' height='720' /></a><br /></div>" >> "$html" 261 echo "<a href='${patient}_freq_indiv_all.png'><img src='${patient}_freq_indiv_all.png' width='1280' height='720' /></a><br /></div>" >> "$html"
262 echo "${scatterplot_tab}</tr></table></div>" >> "$html" 262 echo "${scatterplot_tab}</tr></table></div>" >> "$html"
263 263
264 tail -n+2 ${patient}_reads.txt | sed "s/>//" > tmp.txt 264 tail -n+2 "${patient}_reads.txt" | sed "s/>//" > tmp.txt
265 echo "<div class='tabbertab' title='Data reads'>" >> "$html" 265 echo "<div class='tabbertab' title='Data reads'>" >> "$html"
266 echo "<table><tr><td style='vertical-align:top;'>" >> "$html" 266 echo "<table><tr><td style='vertical-align:top;'>" >> "$html"
267 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_reads'>" >> "$html" 267 echo "<table border = 1 class='result_table summary_table' id='summary_table_${patient}_reads'>" >> "$html"
268 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_All</th><th>Number of sequences_$sample1</th><th>Number of sequences_$sample2</th><th>Number of sequences_$sample3</th><th>Number of sequences_${sample1}_${sample2}</th><th>Number of sequences_${sample1}_${sample3}</th><th>Number of sequences_${sample2}_${sample3}</th></thead>" >> "$html" 268 echo "<thead><th>Ig/TCR gene rearrangement type</th><th>Proximal gene segment</th><th>Distal gene segment</th><th>Cut off value</th><th>Number of sequences ${patient}_All</th><th>Number of sequences_$sample1</th><th>Number of sequences_$sample2</th><th>Number of sequences_$sample3</th><th>Number of sequences_${sample1}_${sample2}</th><th>Number of sequences_${sample1}_${sample3}</th><th>Number of sequences_${sample2}_${sample3}</th></thead>" >> "$html"
269 echo "<tbody>" >> "$html" 269 echo "<tbody>" >> "$html"