0
|
1 <tool id="rgweblogo3" name="Sequence Logo" version="0.4">
|
|
2 <description>generator for fasta (eg Clustal alignments)</description>
|
|
3 <requirements>
|
|
4 <requirement type="package" version="3.3">weblogo</requirement>
|
|
5 <requirement type="package" version="9.10">ghostscript</requirement>
|
|
6 </requirements>
|
|
7 <command interpreter="python">
|
|
8 rgWebLogo3.py -F $outformat -s $size -i $input -o $output -t "$logoname" -c "$colours" -U "$units"
|
|
9 #if $range.mode == 'part'
|
|
10 -l "$range.seqstart" -u "$range.seqend"
|
|
11 #end if
|
|
12 </command>
|
|
13 <inputs>
|
|
14 <page>
|
|
15 <param format="fasta" name="input" type="data" label="Fasta File" />
|
|
16 <param name="logoname" label="Title for output Sequence Logo" type="text" size="50" value="Galaxy-Rgenetics Sequence Logo" />
|
|
17 <param name="outformat" type="select" label="Output format for image (or text report)" >
|
|
18 <option value="png" selected="True">PNG screen quality</option>
|
|
19 <option value="png_print">High quality printable PNG</option>
|
|
20 <option value="pdf">PDF</option>
|
|
21 <option value="jpeg">JPG</option>
|
|
22 <option value="eps">EPS</option>
|
|
23 <option value="txt">Text (shows the detailed calculations for each position - no image)</option>
|
|
24 </param>
|
|
25 <param name="units" type="select" label="Display Units"
|
|
26 help="What the height of each logo element depicts - eg bits of entropy (default)">
|
|
27 <option value="bits" selected="True">Entropy (bits)</option>
|
|
28 <option value="probability">Probability</option>
|
|
29 <option value="nats">Nats</option>
|
|
30 <option value="kT">kT</option>
|
|
31 <option value="kJ/mol">kJ/mol</option>
|
|
32 <option value="kcal/mol">kcal/mol</option>
|
|
33 </param>
|
|
34 <param name="colours" type="select" label="Colour scheme for output Sequence Logo"
|
|
35 help="Note that some of these only make sense for protein sequences!">
|
|
36 <option value="auto" selected="True">Default automatic colour selection</option>
|
|
37 <option value="base pairing">Base pairing</option>
|
|
38 <option value="charge">Charge colours</option>
|
|
39 <option value="chemistry">Chemistry colours</option>
|
|
40 <option value="classic">Classical colours</option>
|
|
41 <option value="hydrophobicity">Hydrophobicity</option>
|
|
42 <option value="monochrome">monochrome</option>
|
|
43 </param>
|
|
44
|
|
45
|
|
46 <conditional name="range">
|
|
47 <param name="mode" type="select" label="Include entire sequence (default) or specify a subsequence range to use">
|
|
48 <option value="complete" selected="true">complete sequence</option>
|
|
49 <option value="part">Only use a part of the sequence</option>
|
|
50 </param>
|
|
51 <when value="complete">
|
|
52 </when>
|
|
53 <when value="part">
|
|
54 <param name="seqstart" size="5" type="integer" value="1" help="WARNING: Specifying indexes outside the sequence lengths will cause unpredictable but bad consequences!"
|
|
55 label="Index (eg 1=first letter) of the start of the sequence range to include in the logo">
|
|
56 </param>
|
|
57 <param name="seqend" size="5" type="integer" value="99999" label="Index (eg 75=75th letter) of the end of the sequence range to include in the logo" >
|
|
58 </param>
|
|
59 </when>
|
|
60 </conditional>
|
|
61 <param name="size" type="select" label="Output weblogo size" >
|
|
62 <option value="large" selected="True">Large</option>
|
|
63 <option value="medium">Medium</option>
|
|
64 <option value="small">Small</option>
|
|
65 </param>
|
|
66 </page>
|
|
67 </inputs>
|
|
68 <outputs>
|
|
69 <data format="pdf" name="output" label="${logoname}_output.${outformat}">
|
|
70 <change_format>
|
|
71 <when input="outformat" value="png_print" format="png" />
|
|
72 <when input="outformat" value="png" format="png" />
|
|
73 <when input="outformat" value="jpeg" format="jpg" />
|
|
74 <when input="outformat" value="eps" format="eps" />
|
|
75 <when input="outformat" value="txt" format="txt" />
|
|
76 </change_format>
|
|
77 </data>
|
|
78 </outputs>
|
|
79 <tests>
|
|
80 <test>
|
|
81 <param name="input" value="rgClustal_testout.fasta" />
|
|
82 <param name="logoname" value="Galaxy/Rgenetics weblogo" />
|
|
83 <param name="outformat" value="jpeg" />
|
|
84 <param name="mode" value="complete" />
|
|
85 <param name="size" value="medium" />
|
|
86 <param name="colours" value="auto" />
|
|
87 <param name="units" value="bits" />
|
|
88 <output name="output" file="rgWebLogo3_test.jpg" ftype="jpg" compare="sim_size" delta="10000" />
|
|
89 </test>
|
|
90 <test>
|
|
91 <param name="input" value="rgClustal_testout.fasta" />
|
|
92 <param name="logoname" value="Galaxy/Rgenetics weblogo" />
|
|
93 <param name="outformat" value="png" />
|
|
94 <param name="mode" value="complete" />
|
|
95 <param name="size" value="medium" />
|
|
96 <param name="colours" value="auto" />
|
|
97 <param name="units" value="probability" />
|
|
98 <output name="output" file="rgWebLogo3_test2.png" ftype="png" compare="sim_size" delta="10000" />
|
|
99 </test>
|
|
100 </tests>
|
|
101 <help>
|
|
102
|
|
103 **Note**
|
|
104
|
|
105 This tool uses Weblogo3_ in Galaxy to generate a sequence logo. The input file must be a fasta file in your current history.
|
|
106
|
|
107 It is recommended for (eg) viewing multiple sequence alignments output from the clustalw tool - set the output to fasta and feed
|
|
108 it in to this tool.
|
|
109
|
|
110 A typical output looks like this
|
|
111
|
|
112 .. image:: ${static_path}/images/rgWebLogo3_test.jpg
|
|
113
|
|
114 ----
|
|
115
|
|
116 **Warning about input Fasta format files**
|
|
117
|
|
118 The Weblogo3 program used by this tool will fail if your fasta sequences are not all EXACTLY the same length. The tool will provide a warning
|
|
119 and refuse to call the weblogo3 executable if irregular length sequences are detected.
|
|
120
|
|
121 Fasta alignments from the companion ClustalW Galaxy tool will work but many other fasta files may cause this tool to fail - please do not file
|
|
122 a Galaxy bug report - this is a feature of the tool and a problem with your source data - not a tool error - please make certain all your fasta
|
|
123 sequences are the same length!
|
|
124
|
|
125 ----
|
|
126
|
|
127 **Attribution**
|
|
128
|
|
129 Weblogo attribution and associated documentation are available at Weblogo3_
|
|
130
|
|
131 This Galaxy wrapper was written by Ross Lazarus for the rgenetics project and the source code is licensed under the LGPL_ like other rgenetics artefacts
|
|
132
|
|
133 .. _Weblogo3: http://weblogo.berkeley.edu/
|
|
134
|
|
135 .. _LGPL: http://www.gnu.org/copyleft/lesser.html
|
|
136
|
|
137 </help>
|
|
138
|
|
139 </tool>
|
|
140
|
|
141
|