Run_stats.tpl 13.9 KB
Newer Older
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
{*
Copyright (C) 2009 INRA
 
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program.  If not, see <http://www.gnu.org/licenses/>.
*}

{extends file='AnalysisTemplate.tpl'}
{block name=params}
	{assign var="params" value=" "|explode:$analyse.params}
	{foreach from=$analyse_results key=sample item=sample_results}
	{assign var="sample" value=$sample}
	{/foreach}
	<ul>
		{if in_array("minSubReadLength", $params)}
		{assign var="minSubReadLength" value=$params|@array_keys:"minSubReadLength"}
		<li class="parameter">Subreads shorter than {$params[$minSubReadLength[0]+1]} (in base pairs) are filtered out and excluded from analysis.</li>
		{/if}
		{if in_array("readScore", $params)}
		{assign var="polymerase_read_qual" value=$params|@array_keys:"readScore"}
		<li class="parameter">Polymerase reads with lower quality than {$params[$polymerase_read_qual[0]+1]} are filtered out and excluded from analysis. </li>
		{/if}
		{if in_array("minLength", $params)}
		{assign var="polymerase_read_length" value=$params|@array_keys:"minLength"}
		<li class="parameter">Polymerase reads shorter than {$params[$polymerase_read_length[0]+1]} (in base pairs) are filtered out and excluded from analysis.</li>
		{/if}
		{if in_array("barcode_file", $params)}
		{assign var="barcode_file" value=$params|@array_keys:"barcode_file"}
		<li class="parameter">Input barcode file : {$params[$barcode_file[0]+1]}.</li>
		{/if}
		{if in_array("barcode_score", $params)}
		{assign var="barcode_score" value=$params|@array_keys:"barcode_score"}
		<li class="parameter">Min identical base for barcode : {$params[$barcode_score[0]+1]}.</li>
		{/if}
	</ul>
{/block}


{block name=results_title} Reports {/block}
{block name=results}
51
	{assign var="analyse_results_sorted" value=$analyse_results["ONT_sample"]|@ksort}
52
53
54
55
	{assign var="metrics" value=$analyse_results["metrics"]}
	{assign var="basic" value=$analyse_results["ONT_sample"].basic}
	{assign var="quality" value=$analyse_results["ONT_sample"].quality}
	{assign var="plots" value=$analyse_results["ONT_sample"].plots}
56
	{assign var="plots_barcode" value=$analyse_results["ONT_sample"].plots_barcode}
57
58
59
60
61
62
63
64
	
	{assign var='basic_headers' value=','|explode:$metrics['basic'].headers}
	{assign var='basic_count' value=$basic_headers|@count}
    
	{assign var='quality_headers' value=','|explode:$metrics['quality'].headers}
	{assign var='quality_count' value=$quality_headers|@count}
	
	{assign var='plots_headers' value=','|explode:$metrics['plots'].headers}
65
66
67
68
69
70
71
72
73
	{assign var='plots_count' value=$plots_headers|@count}
	
	{assign var='plots_barcode_headers' value=','|explode:$metrics['plots_barcode'].headers}
	{assign var='plots_barcode_count' value=$plots_barcode_headers|@count}
    
	{assign var='barcode_headers' value=','|explode:$metrics['barcode'].headers|@ksort}
	{assign var='barcode_headers_count' value=$barcode_headers|@count}
    {assign var='barcode_name_sample' value=','|explode:$metrics['barcode'].names|@ksort}
    {assign var='barcode_name_count' value=$barcode_name_sample|@count}
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
    
    {*debug*}
	<div>
		<legend>Analyse results report - Sample name : {$descriptions.sample_1}</legend>
		
		{assign var='data_col' value=2}
		<div class="row">
			{if $metrics['basic']}
				<div class="col-md-4 col-lg-4">
					<table id="ont_table_basic" class="table table-striped table-bordered dataTable analysis-result-table">
					<thead><tr><th colspan="2"><h4>Basic matrics report</h4></thead></tr></th>
						<tbody>
						{foreach from=$basic_headers key=k item=head}
							<tr>
								{if ($head == 'median_yield_per_sec')  or ($head == 'mean_yield_per_sec')}
									<td><label>{$head}</label></td>
									<td>{($basic.$head)|number_format:2:',':' '}</td>
								{else}
									<td><label>{$head}</label></td>
									<td>{($basic.$head)|number_format:0:',':' '}</td>
								{/if}
							</tr>
							{$data_col = $data_col + 1}
						{/foreach}
						</body>
					</table>
				</div>
			{/if}
			{if $metrics['quality']}
				<div class="col-md-4 col-lg-4">
					<table id="ont_table_quality" class="table table-striped table-bordered dataTable analysis-result-table">
						<thead><tr><th colspan="2"><h4>Quality report</h4></thead></tr></th>
					<tbody>
						{foreach from=$quality_headers key=k item=head}
							<tr>
								{if ($head == 'median_read_quality')  or ($head == 'mean_read_quality')}
									<td><label>{$head}</label></td>
									<td>{($quality.$head)|number_format:2:',':' '}</td>
								{elseif $head|strstr:"nb_read_Q"}
									<td><label>{$head}</label></td>
									<td>{($quality.$head)|number_format:0:',':' '} ({($quality.$head / $basic.nb_reads*100)|number_format:2:',':' '}%)</td>
								{else}
									<td><label>{$head}</label></td>
									<td>{($quality.$head)|number_format:0:',':' '} ({($quality.$head / $basic.total_bases*100)|number_format:2:',':' '}%)</td>
								{/if}
							</tr>
							{$data_col = $data_col + 1}
						{/foreach}
					</tbody>
					</table>
				</div>
			{/if}
			{if $metrics['plots']}
				<div class="col-md-3 col-lg-3">
					<table id="ont_table_plot" class="table table-striped table-bordered dataTable analysis-result-table">
						<thead><tr><th colspan="2"><h4>Plots report</h4></h4></thead></tr></th>
					<tbody>
						{foreach from=$plots_headers key=k item=head}
							<tr>
									<td><label>{$head}</label></td>
									<td><a class="imglink" href="{$plots.$head}" ><i class="glyphicon glyphicon-picture"></i></a></td>
							</tr>
							{$data_col = $data_col + 1}
						{/foreach}
138
139
140
141
142
143
144
145
146
						{if $barcode_name_count > 1 }
						{foreach from=$plots_barcode_headers key=k item=head}
							<tr>
									<td><label>{$head}</label></td>
									<td><a class="imglink" href="{$plots_barcode.$head}" ><i class="glyphicon glyphicon-picture"></i></a></td>
							</tr>
							{$data_col = $data_col + 1}
						{/foreach}
						{/if}
147
148
149
150
151
152
					</tbody>
					</table>
				</div>
			{/if}
		</div>
	</div>
153
154
155
156
	
	
	
	
Maxime Manno's avatar
Maxime Manno committed
157
	<br><br>
158
	{*If there are barcodes file in the run*}
159
	{if $barcode_name_count > 1 }	
160
	
161
	<table id="ont_stats_datatable" class="table table-striped table-bordered dataTable analysis-result-table">
162
163
		<thead>
			<tr>
164
				{assign var="nb_samples" value=$barcode_name_count}
165
				<th class="string-sort" rowspan="2" id="th_id_1"><center>Sample {if $nb_samples > 1 }({$nb_samples}){/if}</center></th>
166
        		<th colspan="{$barcode_headers_count}"><center>Barcode ONT results</center></th>
167
168
169
170
			</tr>
			<tr>
				{assign var="th_id" value=2}
				{foreach from=$barcode_headers key=k item=head}
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
					{if $head == 'total_bases'}
						<th class = "numeric-sort" id="th_id_{$th_id}">Total bases</th>
					{elseif $head == 'nb_reads'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Nb reads</th>
					{elseif $head == 'median_read_length'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Median read length</th>
					{elseif $head == 'median_read_quality'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Median read quality</th>
					{elseif $head == 'N50_read_length'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">N50 read length</th>
					{elseif $head == 'L50_read_length'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">L50 read length</th>
					{elseif $head == 'barcode_score'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Barcode score</th>
					{elseif $head == 'nb_read_Q>10'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Nb reads Q>10</th>
					{elseif $head == 'total_bases_Q>10'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Total bases Q>10</th>
					{elseif $head == 'median_yield_per_sec'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Median yield per sec</th>
					{elseif $head == 'nb_actif_channel'}
						<th class = "numeric-sort"  id="th_id_{$th_id}">Nb actif channel</th>
193
194
195
196
197
198
199
200
201
					{/if}
					{$th_id = $th_id +1}
				{/foreach}
			</tr>
		</thead>
		<tbody>
			{$i = 1}
			{assign var="totalReads" value=0}
			{assign var="totalBases" value=0}
202
203
204
205
206
207
			
            {foreach from=$analyse_results_sorted key=sample item=sample_results}
            	{if $sample|in_array:$barcode_name_sample}
            		<tr>
						<td id='sample_{$i}_col_1' class="sample_name">{$sample}</td>
						{$col_id = 2}
208
						{foreach from=$barcode_headers key=k item=head}
209
210
						
							{if $head == 'total_bases'}
Maxime Manno's avatar
Maxime Manno committed
211
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
212
							{elseif $head == 'nb_reads'}
Maxime Manno's avatar
Maxime Manno committed
213
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
214
							{elseif $head == 'median_read_length'}
Maxime Manno's avatar
Maxime Manno committed
215
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
216
							{elseif $head == 'median_read_quality'}
Maxime Manno's avatar
Maxime Manno committed
217
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:1:',':' '}</th>
218
							{elseif $head == 'N50_read_length'}
Maxime Manno's avatar
Maxime Manno committed
219
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
220
							{elseif $head == 'L50_read_length'}
Maxime Manno's avatar
Maxime Manno committed
221
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
222
							{elseif $head == 'barcode_score'}
Maxime Manno's avatar
Maxime Manno committed
223
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:1:',':' '}</th>
224
							{elseif $head == 'nb_read_Q>10'}
Maxime Manno's avatar
Maxime Manno committed
225
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
226
							{elseif $head == 'total_bases_Q>10'}
Maxime Manno's avatar
Maxime Manno committed
227
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
228
							{elseif $head == 'median_yield_per_sec'}
Maxime Manno's avatar
Maxime Manno committed
229
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
230
							{elseif $head == 'nb_actif_channel'}
Maxime Manno's avatar
Maxime Manno committed
231
								<th id="sample_{$i}_col_{$col_id}">{($sample_results.$head)|number_format:0:',':' '}</th>
232
233
234
							{/if}
							
							{$col_id = $col_id + 1}
235
						{/foreach}
236
237
238
239
240
					</tr>
					{$i = $i + 1}
          		{/if}
            {/foreach}

241
		</tbody>
Maxime Manno's avatar
Maxime Manno committed
242
	   {*	
243
	   	<tfoot>
244
			<tr>	
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
				<th>Total : </th>
				<th>&nbsp;</th>
				{foreach from=$barcode_headers key=k item=head}
					{if $head == 'Reads'}
						<th>{$totalReads|number_format:0:',':' '}</th>
					{elseif $head == 'Bases'}
						<th>{$totalBases|number_format:0:',':' '}</th>
					{/if}
				{/foreach}
			</tr>
			<tr>
				<th>Mean : </th>
				<th>&nbsp;</th>
				{foreach from=$barcode_headers key=k item=head}
					{if $head == 'Reads'}
						<th>{($totalReads/$nb_samples)|number_format:0:' ':' '}</th>
					{elseif $head == 'Bases'}
						<th>{($totalBases/$nb_samples)|number_format:0:' ':' '}</th>
					{/if}
				{/foreach}
			</tr>
			<tr>
				<th>All metrics :  <input type="checkbox" id="check_all_metrics"></th>
				<th>&nbsp;</th>
				{$th_id = 2}
				{foreach from=$barcode_headers key=k item=head}
					<th><center><input type="checkbox" id="chk_col_{$th_id}"></center></th>
					{$th_id = $th_id +1}
				{/foreach}
			</tr>
			<tr>
				<th align="left" colspan="{$th_id}">
					With selection :
					<button type="button" class="btn btn-default multiple-selection-btn" id="create_graph"><i class="glyphicon glyphicon-signal"></i> Compare</button>
				</th>
			</tr>
		</tfoot>
Maxime Manno's avatar
Maxime Manno committed
282
		*}
283
	</table>
Maxime Manno's avatar
Maxime Manno committed
284
	<br><br>
285
286
	{/if}

287
288
289
290




291
292
293
294
295
296
	<div class="tx-nG6-pi1-help">
		<img src="" alt="" class="img" />
		<p>Help for Basic metrics report :</p>
		<span class="meta">
			<ul>
				<li><strong>nb_reads</strong> : 
maxime mano's avatar
maxime mano committed
297
					The total number of reads for this Flow Cell (including the reads that are filtered during the basecalling).
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
				</li>
				<li><strong>total_bases</strong> : 
					The number of bases for this Flow Cell.
				</li>
				<li><strong>median_read_length</strong> : 
					The median length of the total reads (bases).
				</li>
				<li><strong>N50_read_length</strong> : 
					50% of all bases come from reads longer than this value.
				</li>
				<li><strong>median_yield_per_sec</strong> : 
					The median number of bases per second for the pores of this Flow Cell.
				</li>
				<li><strong>nb_actif_channel</strong> : 
					The number of activ channels during the sequencing.
				</li>
			</ul>
		</span>
	</div>
	<div class="tx-nG6-pi1-help">
		<img src="" alt="" class="img" />
		<p>Help for Quality report :</p>
		<span class="meta">
			<ul>
				<li><strong>median_read_quality</strong> : 
					The median quality (qscore ONT) for this Flow Cell.
				</li>
				<li><strong>nb_read_Q>N</strong> : 
					The total number of reads that have a quality > N.
				</li>
				<li><strong>total_bases_Q>N</strong> : 
					The total number of bases for the reads that have a quality > N.
				</li>
			</ul>
		</span>
	</div>
	
	<div class="tx-nG6-pi1-help">
		<img src="" alt="" class="img" />
		<p>Help for Plots report :</p>
		<span class="meta">
			<ul>
				<li><strong>cumulyieldperhour</strong> : 
					The cumulate number of bases per hour for this Flow Cell.
				</li>
				<li><strong>outrm_</strong> : 
					The outsiders have been removed.
				</li>
				<li><strong>outrm_distriblength</strong> : 
					The length distribution of reads for this Flow Cell.
				</li>
				<li><strong>outrm_distribqscore</strong> : 
					The quality (qscore ONT) distribution of reads for this Flow Cell.
				</li>
				<li><strong>outrm_lengthvsqscore_density</strong> : 
					The length and qscore distribution of reads with density for this Flow Cell.
				</li>
355
356
357
358
359
360
361
362
				{if $barcode_name_count > 1 }
				<li><strong>qscoreboxplot</strong> : 
					Boxplots for each barcode of the mean qscore of reads for this Flow Cell.
				</li>
				<li><strong>qscorepertimeintervalsboxplot</strong> : 
					Boxplots for each barcode of the mean qscore of reads per interval of time for this Flow Cell.
				</li>
				{/if}
363
364
365
366
367
			</ul>
		</span>
	</div>
{/block}