Mercurial > repos > davidvanzessen > prisca
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" |