|
0
|
1 <tool id="sam_header" name="NGS Run Annotation">
|
|
|
2 <description>Create a SAM format header from run metadata for sample annotation.</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="0.1.3">mimodd</requirement>
|
|
|
5 </requirements>
|
|
|
6 <command>
|
|
|
7 mimodd header
|
|
|
8
|
|
|
9 --rg_id "$rg_id"
|
|
|
10 --rg_sm "$rg_sm"
|
|
|
11
|
|
|
12 #if $str($rg_cn):
|
|
|
13 --rg_cn "$rg_cn"
|
|
|
14 #end if
|
|
|
15 #if $str($rg_ds):
|
|
|
16 --rg_ds "$rg_ds"
|
|
|
17 #end if
|
|
|
18 #if $str($anno) and $str($month) and $str($day):
|
|
|
19 --rg_dt "$anno-$month-$day"
|
|
|
20 #end if
|
|
|
21 #if $str($rg_lb):
|
|
|
22 --rg_lb "$rg_lb"
|
|
|
23 #end if
|
|
|
24 #if $str($rg_pl):
|
|
|
25 --rg_pl "$rg_pl"
|
|
|
26 #end if
|
|
|
27 #if $str($rg_ds):
|
|
|
28 --rg_pi "$rg_pi"
|
|
|
29 #end if
|
|
|
30 #if $str($rg_pu):
|
|
|
31 --rg_pu "$rg_pu"
|
|
|
32 #end if
|
|
|
33
|
|
|
34 --outputfile $outputfile
|
|
|
35
|
|
|
36 </command>
|
|
|
37
|
|
|
38 <inputs>
|
|
|
39 <param name="rg_id" type="text" size="80" label="read-group ID (mandatory)">
|
|
|
40 <sanitizer invalid_char="">
|
|
|
41 <valid initial="string.printable">
|
|
|
42 <remove value=""" />
|
|
|
43 </valid>
|
|
|
44 <mapping initial="none">
|
|
|
45 <add source=""" target="\""/>
|
|
|
46 </mapping>
|
|
|
47 </sanitizer>
|
|
|
48 </param>
|
|
|
49 <param name="rg_sm" type="text" size="80" label="sample name (mandatory)">
|
|
|
50 <sanitizer invalid_char="">
|
|
|
51 <valid initial="string.printable">
|
|
|
52 <remove value=""" />
|
|
|
53 </valid>
|
|
|
54 <mapping initial="none">
|
|
|
55 <add source=""" target="\""/>
|
|
|
56 </mapping>
|
|
|
57 </sanitizer>
|
|
|
58 </param>
|
|
|
59 <param name="rg_cn" type="text" size="80" label="name of sequencing center">
|
|
|
60 <sanitizer invalid_char="">
|
|
|
61 <valid initial="string.printable">
|
|
|
62 <remove value=""" />
|
|
|
63 </valid>
|
|
|
64 <mapping initial="none">
|
|
|
65 <add source=""" target="\""/>
|
|
|
66 </mapping>
|
|
|
67 </sanitizer>
|
|
|
68 </param>
|
|
|
69 <param name="rg_ds" type="text" size="80" label="description">
|
|
|
70 <sanitizer invalid_char="">
|
|
|
71 <valid initial="string.printable">
|
|
|
72 <remove value=""" />
|
|
|
73 </valid>
|
|
|
74 <mapping initial="none">
|
|
|
75 <add source=""" target="\""/>
|
|
|
76 </mapping>
|
|
|
77 </sanitizer>
|
|
|
78 </param>
|
|
|
79 <param name="anno" type="text" label="year (YYYY) the run was produced" />
|
|
|
80 <param name="month" type="text" label="month (MM) the run was produced" />
|
|
|
81 <param name="day" type="text" label="day (DD) the run was produced" />
|
|
|
82 <param name="rg_lb" type="text" size="80" label="read-group library">
|
|
|
83 <sanitizer invalid_char="">
|
|
|
84 <valid initial="string.printable">
|
|
|
85 <remove value=""" />
|
|
|
86 </valid>
|
|
|
87 <mapping initial="none">
|
|
|
88 <add source=""" target="\""/>
|
|
|
89 </mapping>
|
|
|
90 </sanitizer>
|
|
|
91 </param>
|
|
|
92 <param name="rg_pl" type="text" label="platform/technology used to produce the reads" />
|
|
|
93 <param name="rg_pi" type="text" label="predicted median insert size" />
|
|
|
94 <param name="rg_pu" type="text" size="80" label="platform unit; unique identifier">
|
|
|
95 <sanitizer invalid_char="">
|
|
|
96 <valid initial="string.printable">
|
|
|
97 <remove value=""" />
|
|
|
98 </valid>
|
|
|
99 <mapping initial="none">
|
|
|
100 <add source=""" target="\""/>
|
|
|
101 </mapping>
|
|
|
102 </sanitizer>
|
|
|
103 </param>
|
|
|
104 </inputs>
|
|
|
105
|
|
|
106 <outputs>
|
|
|
107 <data name="outputfile" format="sam" label="${rg_sm} (${rg_id}) header information from MiModd ${tool.name} on ${on_string}"/>
|
|
|
108 </outputs>
|
|
|
109
|
|
|
110 <help>
|
|
|
111 .. class:: infomark
|
|
|
112
|
|
|
113 **What it does**
|
|
|
114
|
|
|
115 This tool takes the user-provided information about a next-generation sequencing run and constructs a valid header in the SAM file format from it.
|
|
|
116
|
|
|
117 The result file can be used by the tools *Convert* and *Reheader* or in the *SNAP Read Alignment* step to add run metadata to sequenced reads files (or to overwrite pre-existing information).
|
|
|
118
|
|
|
119 **Note:**
|
|
|
120
|
|
|
121 **MiModD requires run metadata for every input file at the Alignment step !**
|
|
|
122
|
|
|
123 **Tip:**
|
|
|
124
|
|
|
125 While you can do Alignments from fastq file format by providing a custom header file directly to the *SNAP Read Alignment* tool, the **recommended approach** is to first convert all input files to and archive all datasets in SAM/BAM format with appropriate header information prior to any downstream analysis. Although a bit more time-consuming this practice protects against information loss and ensures that the input datasets will remain useful for others in the future.
|
|
|
126
|
|
|
127 </help>
|
|
|
128 </tool>
|
|
|
129
|