/tools/rgenetics/rgWebLogo3.xml

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