Mercurial > repos > wolma > mimodd

changeset 6:85214e4428fd

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
upgrade to mimodd version 0.1.7.0
author Wolfgang Maier
date Fri, 11 Sep 2015 16:25:49 +0200
parents bdd1995c9e66
children 3619e85a5477
files annotate_variants.xml bamsort.xml cloudmap.xml convert.xml covstats.xml deletion_predictor.xml fileinfo.xml reheader.xml sam_header.xml snap_caller.xml snp_caller_caller.xml snpeff_genomes.xml tool_dependencies.xml toolshed_macros.xml varextract.xml vcf_filter.xml
diffstat 16 files changed, 318 insertions(+), 61 deletions(-) [+]
line wrap: on
line diff
--- a/annotate_variants.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/annotate_variants.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="annotate_variants" name="Variant Annotation" version="0.1.6.1">
+<tool id="annotate_variants" name="Variant Annotation" version="0.1.7.0">
   <description>Predict the effects of SNPs and indels on known genes in the reference genome using SnpEff</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/bamsort.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/bamsort.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="bamsort" name="Sort BAM file" version="0.1.6.1">
+<tool id="bamsort" name="Sort BAM file" version="0.1.7.0">
   <description>Sort a BAM file by coordinates (or names) of the mapped reads</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/cloudmap.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/cloudmap.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,54 +1,247 @@
-<tool id="cloudmap_prepare" name="Prepare variant data for mapping" version="0.1.6.1">
-  <description>with the CloudMap series of tools.</description>
+<tool id="nacreousmap" name="NacreousMap" version="0.1.7.0">
+  <description>Map causative mutations by multi-variant linkage analysis.</description>
   <macros>
     <import>toolshed_macros.xml</import>
   </macros>
   <expand macro="requirements"/>
   <version_command>mimodd version -q</version_command>
   <command>
-    mimodd cloudmap "$ifile" ${run.mode} "$sample"
-    
-    #if $str($run.related_parent_sample):
-        -r "${run.related_parent_sample}"
+    mimodd map "${opt.source.ifile}" ${opt.mode}
+    #if $str($opt.source.sample):
+      -m "${opt.source.sample}"
+    #end if
+    #if $str($opt.source.related_parent_sample):
+        -r "${opt.source.related_parent_sample}"
+    #end if
+    #if $str($opt.source.unrelated_parent_sample):
+        -u "${opt.source.unrelated_parent_sample}"
+    #end if
+    $opt.source.infer_missing
+    -o "$ofile"
+    #if $str($opt.source.seqdict_required.required) == "yes":
+        -s "${opt.source.seqdict_required.seqdict}"
+    #end if
+    $opt.source.norm
+    #if $len($opt.source.bin_sizes):
+      --bin-sizes
+      #for $size in $opt.source.bin_sizes:
+        "${size.bin_size}"
+      #end for
+    #end if
+    #if $str($opt.source.tabfile):
+      $str($opt.source.tabfile) $tfile
     #end if
-    #if $str($run.unrelated_parent_sample):
-        -u "${run.unrelated_parent_sample}"
+    #if $str($opt.source.plotopts.plots):
+      $str($opt.source.plotopts.plots) "$pfile"
+      $str($opt.source.plotopts.xlim)
+      #if $str($opt.source.plotopts.hylim):
+        --ylim-hist $str($opt.source.plotopts.hylim)
+      #end if
+      #if $str($opt.source.plotopts.hcols) and $len($opt.source.plotopts.hcols):
+        --hist-colors
+        #for $color in $opt.source.plotopts.hcols:
+          "${color.hcolor}"
+        #end for
+      #end if
+      #if $str($opt.source.plotopts.sylim):
+        --ylim-scatter $str($opt.source.plotopts.sylim)
+      #end if
+      #if $str($opt.source.plotopts.pcol):
+        --points-color "$str($opt.source.plotopts.pcol)"
+      #end if
+      #if $str($opt.source.plotopts.lcol):
+        --loess-color "$str($opt.source.plotopts.lcol)"
+      #end if
+      #if $str($opt.source.plotopts.span):
+        --loess-span "$str($opt.source.plotopts.span)"
+      #end if      
     #end if
     
-    $run.infer_missing
-    -o "$ofile"
-    
-    #if $seqdict:
-        -s "$dictfile"
-    #end if
   </command>
 
+  <macros>
+    <macro name="svd_unconditional">
+      <expand macro="hidden_vaf_algo_params" />
+      <expand macro="seqdict_param" />
+      <expand macro="bins" />
+      <param name="norm" type="boolean" label="normalize variant counts to bin-width" truevalue="" falsevalue="--no-normalize" checked="true" help="without normalization the tool will just report the number of nucleotides per bin; with normalization the results for different bin-widths will be comparable." />
+      <conditional name="plotopts">
+        <param name="plots" type="select" label="graphical output settings">
+          <option value="">Do not generate graphs.</option>
+          <option value="-p">Give me graphics.</option>
+        </param>
+        <when value="-p">
+          <expand macro="scatter_default" />
+          <param name="hylim" type="text" label="upper limit for the histogram y-axis (leave blank for automatic scaling)" />
+          <param name="xlim" type="select" label="x-axis scaling">
+            <option value="">preserve relative contig sizes</option>
+            <option value="--fit-width">scale each contig to fit the plot width</option>
+          </param>
+          <expand macro="hist_colors" />
+        </when>
+      </conditional>
+    </macro>
+    <macro name="vaf_unconditional">
+      <expand macro="bins" />
+      <param name="norm" type="boolean" label="normalize variant counts to bin-width" truevalue="" falsevalue="--no-normalize" checked="true" />
+      <conditional name="plotopts">
+        <param name="plots" type="select" label="graphical output settings">
+          <option value="">Do not generate graphs.</option>
+          <option value="--no-scatter -p">Generate only histograms</option>
+          <option value="--no-hist -p">Generate only scatter plots</option>
+          <option value="-p">Give me everything (scatter plots and histograms)</option>
+        </param>
+        <when value="--no-scatter -p">
+          <expand macro="scatter_default" />
+          <param name="hylim" type="text" label="upper limit for the histogram y-axis (leave blank for automatic scaling)" />
+          <param name="xlim" type="select" label="x-axis scaling">
+            <option value="">preserve relative contig sizes</option>
+            <option value="--fit-width">scale each contig to fit the plot width</option>
+          </param>
+          <expand macro="hist_colors" />
+        </when>
+        <when value="--no-hist -p">
+          <expand macro="hist_default" />
+          <param name="sylim" type="text" label="upper limit for the scatter plot y-axis (default: 1)" />
+          <param name="xlim" type="select" label="x-axis scaling">
+            <option value="">preserve relative contig sizes</option>
+            <option value="--fit-width">scale each contig to fit the plot width</option>
+          </param>
+          <param name="span" type="text" label="span value to be used in calculating the Loess regression line through the scatter data (default=0.1)" help="smaller values give a more responsive curve that often picks up local evidence for tight linkage better, but too small values lead to plotting failures (in that case just rerun the tool with a larger value)." />
+          <expand macro="scatter_colors" />
+        </when>
+        <when value="-p">
+          <expand macro="plot_all" />
+        </when>
+      </conditional>
+    </macro>
+    <macro name="hidden_vaf_algo_params">
+      <param name="sample" type="hidden" value="" />
+      <param name="related_parent_sample" type="hidden" value="" />
+      <param name="unrelated_parent_sample" type="hidden" value="" />
+      <param name="infer_missing" type="hidden" value="" />
+    </macro>
+    <macro name="bins">
+      <repeat name="bin_sizes" title="bin sizes to analyze variants in (defaults to: 1Mb and 500Kb" default="0" min="0" help="Values can be entered in bases (e.g., 1000000), kilobases (e.g., 500Kb) or megabases (e.g., 1Mb), but must be integral, i.e. no decimal numbers are allowed.">
+        <param name="bin_size" type="text" />
+      </repeat>
+    </macro>
+    <macro name="scatter_default">
+      <param name="sylim" type="hidden" value="" />
+      <param name="span" type="hidden" value="" />
+      <param name="pcol" type="hidden" value="" />
+      <param name="lcol" type="hidden" value="" />
+    </macro>
+    <macro name="hist_default">
+      <param name="hylim" type="hidden" value="" />
+      <param name="hcols" type="hidden" value="" />
+    </macro>
+    <macro name="hist_colors">
+      <repeat name="hcols" title="histogram colors" default="0" min="0" help="For each bin size chosen above a histogram will be generated with its color selected from the list provided here (defaults to alternating darkgrey, red).">
+        <param name="hcolor" type="color" value="darkgrey">
+          <sanitizer><valid><add value="#" /></valid></sanitizer>
+        </param>
+      </repeat>
+    </macro>
+    <macro name="scatter_colors">
+      <param name="pcol" type="color" value="#454545" label="color to be used for the scatter plot data points (default: gray27)">
+        <sanitizer><valid><add value="#" /></valid></sanitizer>
+      </param>
+      <param name="lcol" type="color" value="red" label="color to be used for the regression line (default: red)">
+        <sanitizer><valid><add value="#" /></valid></sanitizer>
+      </param>
+    </macro>
+    <macro name="plot_all">
+      <param name="hylim" type="text" label="upper limit for the histogram y-axis (leave blank for automatic scaling)" />
+      <param name="sylim" type="text" label="upper limit for the scatter plot y-axis (default: 1)" />
+      <param name="xlim" type="select" label="x-axis scaling">
+        <option value="">preserve relative contig sizes</option>
+        <option value="--fit-width">scale each contig to fit the plot width</option>
+      </param>
+      <param name="span" type="text" label="span value to be used in calculating the Loess regression line through the scatter data (default=0.1)" help="smaller values give a more responsive curve that often picks up local evidence for tight linkage better, but too small values lead to plotting failures (in that case just rerun the tool with a larger value)." />
+      <expand macro="hist_colors" />
+      <expand macro="scatter_colors" />
+    </macro>
+    <macro name="seqdict_param">
+      <conditional name="seqdict_required">
+        <param name="required" type="select" label="does this input file require a CloudMap-style sequence dictionary?" help="A sequence dictionary file is required ONLY if the input file does not provide information about the sizes of the chromosomes defined in it. It is NEVER needed for MiModD-generated input files.">
+          <option value="no">No</option>
+          <option value="yes">Yes</option>
+        </param>
+        <when value="yes">
+          <param name="seqdict" type="data" format="tabular" label="CloudMap-style sequence dictionary file" />
+        </when>
+      </conditional>
+    </macro>
+  </macros>
+  
   <inputs>
-    <param name="ifile" type="data" format="vcf" label="vcf input file" />
-    <conditional name="run">
-      <param name="mode" type="select" label="Type of CloudMap analysis to prepare data for" help="select Simple Variant Density (SVD) Mapping to generate output for the CloudMap EMS Variant Density Mapping tool or Variant Allele Frequency (VAF) Mapping for output for the CloudMap Variant Discovery or Hawaiian Variant Mapping tools.">
+    <conditional name="opt">
+      <param name="mode" type="select" label="type of mapping analysis to perform" help="Select Simple Variant Density (SVD) Mapping to map mutations based on linked inheritance in near isogenic populations, Variant Allele Frequency (VAF) Mapping for bulk segregant analysis. Select Reprocess for rapidly replotting the result of a previous VAF analysis.">
         <option value="SVD">Simple Variant Density Mapping</option>
         <option value="VAF">Variant Allele Frequency Mapping</option>
       </param>
       <when value="SVD">
-        <param name="related_parent_sample" type="hidden" value="" />
-        <param name="unrelated_parent_sample" type="hidden" value="" />
-        <param name="infer_missing" type="hidden" value="" />
+        <conditional name="source">
+          <param name="inputtype" type="select" label="data source to use">
+            <option value="vcf">VCF file of variants (for de-novo mapping)</option>
+            <option value="rep">per-variant report file (for remapping a previous analysis)</option>
+          </param>
+          <when value="vcf">
+            <param name="ifile" type="data" format="vcf" label="input file with variants to analyze" />
+            <expand macro="svd_unconditional" />
+            <param name="tabfile" type="select" label="additional per-variant output file" help="You can either choose to produce a tabular per-variant report, which is useful for fast replotting with different plot settings or a vcf-like CloudMap-compatibility file that can be used as input for the CloudMap EMS Variant Density Mapping tool as an alternative plotting tool.">
+              <option value="">Do not generate per-variant output</option>
+              <option value="-t">Tabular per-variant report</option>
+              <option value="--cloudmap -t">CloudMap compatibility file</option>
+            </param>
+          </when>
+          <when value="rep">
+            <param name="ifile" type="data" format="tabular" label="input file with variants to analyze" />
+            <param name="tabfile" type="hidden" value="" />
+            <expand macro="svd_unconditional" />
+          </when>
+        </conditional>
       </when>
       <when value="VAF">
-        <param name="related_parent_sample" type="text" label="name of the related parent sample" help="the sample that provides variants present in your original mutant strain or in an ancestor (like the pre-mutagenesis strain); leave blank if not available" />
-        <param name="unrelated_parent_sample" type="text" label="name of the unrelated parent sample" help="the sample that provides variants present in the unrelated mapping strain (or in an ancestor of it) used in the mapping cross; leave blank if not available" />
-        <param name="infer_missing" type="boolean" checked="false" truevalue="--infer-missing" falsevalue="" label="Infer alleles for missing parent" help="if variant data for either the related or the unrelated parent strain is not available, the tool can try to infer the alleles present in that parent from the allele spectrum found in the mapping sample. This is an experimental option that will give a benefit only in certain situations. Enable at your own risk." />
+        <conditional name="source">
+          <param name="inputtype" type="select" label="data source to use">
+            <option value="vcf">VCF file of variants (for de-novo mapping)</option>
+            <option value="rep">per-variant report file (for remapping a previous analysis)</option>
+          </param>
+          <when value="vcf">
+            <param name="ifile" type="data" format="vcf" label="input file with variants to analyze" />
+            <expand macro="seqdict_param" />
+            <param name="sample" type="text" label="mapping sample name" help="the sample to perform mutation mapping for" />
+            <param name="related_parent_sample" type="text" label="name of the related parent sample" help="the sample that provides variants present in your original mutant strain or in an ancestor (like the pre-mutagenesis strain); leave blank if not available" />
+            <param name="unrelated_parent_sample" type="text" label="name of the unrelated parent sample" help="the sample that provides variants present in the unrelated mapping strain (or in an ancestor of it) used in the mapping cross; leave blank if not available" />
+            <param name="infer_missing" type="boolean" checked="false" truevalue="--infer-missing" falsevalue="" label="Infer alleles for missing parent" help="if variant data for either the related or the unrelated parent strain is not available, the tool can try to infer the alleles present in that parent from the allele spectrum found in the mapping sample. This is an EXPERIMENTAL option that will give a benefit only in certain situations. Enable at your own risk." />
+            <expand macro="vaf_unconditional" />
+            <param name="tabfile" type="select" label="additional per-variant output file" help="You can either choose to produce a tabular per-variant report, which is useful for fast replotting with different plot settings or a vcf-like CloudMap-compatibility file that can be used as input for the CloudMap Hawaiian Variant Mapping tool as an alternative plotting tool.">
+              <option value="">Do not generate per-variant output</option>
+              <option value="-t">Tabular per-variant report</option>
+              <option value="--cloudmap -t">CloudMap compatibility file</option>
+            </param>
+          </when>
+          <when value="rep">
+            <param name="ifile" type="data" format="tabular" label="input file with variants to analyze" />
+            <expand macro="seqdict_param" />
+            <param name="tabfile" type="hidden" value="" />
+            <expand macro="hidden_vaf_algo_params" />
+            <expand macro="vaf_unconditional" />
+          </when>
+        </conditional>
       </when>
     </conditional>
-    <param name="sample" type="text" label="mapping sample name" help="the sample to perform mutation mapping for" />
-    <param name="seqdict" type="boolean" checked="true" label="Generate species configuration file for CloudMap" />      
   </inputs>
 
   <outputs>
-    <data name="ofile" format="vcf" label="CloudMap-ready ${run.mode} File from ${on_string}" />
-    <data name="dictfile" format="tabular" label="Species Configuration File for CloudMap from ${on_string}">
-      <filter>seqdict</filter>
+    <data name="ofile" format="tabular" label="MiModD ${opt.mode} Mapping - binned variant counts for ${on_string}" />
+    <data name="tfile" format="tabular" label="MiModD ${opt.mode} Mapping - per-variant report for ${on_string}">
+      <filter>(opt['source']['tabfile'])</filter>
+    </data>
+    <data name="pfile" format="pdf" label="MiModD ${opt.mode} Mapping - linkage plots for ${on_string}">
+      <filter>(opt['source']['plotopts']['plots'])</filter>
     </data>
   </outputs>
 
@@ -57,29 +250,93 @@
 
    **What it does**
 
-The purpose of this tool is to provide compatibility of the MiModD analysis workflow with the external `CloudMap`_ *EMS Variant Density Mapping*, *Variant Discovery Mapping* and *Hawaiian Variant Mapping* tools. These tools complement MiModD by providing easily interpreted visualizations of mapping-by-sequencing analysis workflows.
+This tool is a complete rewrite of and improves the EMS Variant Density and Hawaiian Variant Mapping tools of `CloudMap`_. It is the most downstream tool in `mapping-by-sequencing analysis workflows in MiModD`_.
+
+It can be used to analyze and visualize the inheritance pattern of variants detected and selected by other MiModD tools or as an alternative (and more versatile) plotting engine for data generated with `CloudMap`_.
+
+-------------
+
+**Usage Modes:**
 
-The tool converts a VCF file as generated by the *Extract Variant Sites* or *VCF Filter* tools to the format expected by the *CloudMap* series of tools.
+This tool can be run in one of two different modes depending on the type of mapping analysis that should be performed:
+
+1) *Simple Variant Density (SVD) Mapping* mode analyzes the density of variants along the reference genome by dividing each chromosome into regions of user-defined size (bins) and counting the variants found in each bin. 
+    
+   All variants listed in the input file are analyzed in this mode, which means that as input you will typically want to use filtered lists of variants (as produced by the VCF Filter tool). 
+    
+   The aim of SVD analysis is to identify clusters of variants in an outcrossed strain carrying a selectable unknown mutation, which is interpreted as linkage between the corresponding genomic region and the unknown mutation.
+    
+   This mode corresponds roughly to EMS Variant Density Mapping in CloudMap.
 
-Optionally, it also extracts the chromosome names and sizes and reports them in the *CloudMap* *species configuration file* format.
-Such a file is required as input to the current versions of the *CloudMap* *Hawaiian* and *Variant Density* mapping tools, if you are working with a species other than the natively supported ones (i.e., other than *C. elegans*, *A. thaliana* or *Brachypodium distachyon*).
+2) *Variant Allele Frequency (VAF) Mapping** mode analyzes the inheritance pattern in cross-progeny at sites, at which the parents are homozygous for different alleles. 
+    
+   The aim of VAF analysis is to identify clusters of variants with (near) homozygous inheritance in a F2 (or later generation) population obtained from a cross between a strain carrying a selectable unknown mutation and an unrelated mapping strain. Such a cluster is interpreted as linkage between the corresponding genomic region and the unknown mutation selected for in the F2 generation. 
+    
+   This mode corresponds roughly to Hawaiian Variant Mapping in CloudMap, but can simultaneously take into account non-reference alleles found in either parent strain (CloudMap users may think of this as a combined Hawaiian Variant and Variant Discovery Mapping analysis).
+
+-------------
+
+**Input:**
 
-To use the output datasets of the tool with *CloudMap*, you only have to upload them to any public Galaxy server that hosts *CloudMap* like, e.g., the main Galaxy server at https://usegalaxy.org .
+Valid input for this tool are VCF files (any VCF file in SVD mode, a MiModD-generated multi-sample VCF file in VAF mode) or a CloudMap tabular report file as generated by the Hawaiian Variant Mapping tool. Alternatively, the tool can generate (in both modes) its own tabular report file, which can be used as input instead of the original VCF file when rerunning the tool with different plotting parameters to reduce analysis time. 
+
+.. class:: infomark
+
+   CloudMap-generated tabular input files require, as additional input, a CloudMap-style sequence dictionary (even if the original CloudMap analysis was possible without one) as described in the original CloudMap paper. This file has a simple two-column tab-delimited format, in which each line lists the chromosome name (as it appears in the input VCF file) and the up-rounded length of the chromosome in megabases.
+
+-------------
+
+**Output:**
+
+The tool produces up to three output files:
 
-**Notes:**
+1) a default tabular file of binned variant counts that can be used to plot the data with external software such as Excel,
+  
+  
+2) an optional pdf containing linkage plots, which should look just like the plots produced by CloudMap, but are optimized for file size and display speed and offer more user-configurable parameters and
+  
+  
+3) an optional tabular per-variant report file, which can be configured to be either a valid input file for the corresponding original CloudMap tool (for users who really, really want to continue using CloudMap for plotting) or to be reusable in fast reruns of the tool (which can be useful to experiment with different plotting parameters).
 
-1) Simple Variant Density (SVD) Mapping mode generates output for use with the CloudMap EMS Variant Density Mapping tool. The aim of SVD analysis is to identify clusters of variants that appear linked to a mutant phenotype selected for during several rounds of outcrossing or backcrossing to a non-mutagenized strain. The "mapping sample" is the out-/backcrossed strain and only its variants are taken into account for the analysis.
+-------------
 
-.. class:: warningmark
+**Settings:**
+
+1) Analysis settings
 
-   EMS Variant Density Mapping is currently limited to *C. elegans* and other species with six chromosomes on the *CloudMap* side.
+   *bin size to analyze variants in* - determines the width of the regions along each chromosome, in which variants are counted and analyzed together. 
+     
+   Several bin sizes can be specified and for each size you will get a corresponding report section in the binned variant counts file and a histogram plot in the linkage plots file.
+   
+   *normalize variant counts to bin-width* - if selected (as per default) the variant counts for different bin sizes are not absolute, but normalized to the bin width
+   
+   *sample names (in VAF mode only)* - to analyze inheritance patterns, VAF mode needs information about the relationship between the samples defined in the input VCF file:
+      
+   The *mapping sample name* should be set to the name of the sample for which the inheritance pattern is to be analyzed (the pooled progeny population). 
+     
+   The *name of the related sample* should be that of the parent sample that carried and brought in the unknown mutation to be mapped (or, alternatively, that of a closely related ancestor). 
+     
+   Finally, the *name of the unrelated sample* should be that of the other parent strain used in the cross. 
+     
+   At least one of the parent samples MUST be specified, but if the input file contains variant information for both parents, they can be analyzed together for higher mapping accuracy. If you are reanalyzing a tabular report file from a previous tool run or from CloudMap, the association between variants and samples is already incorporated into the input file and cannot be specified again.
 
-2) Variant Allele Frequency (VAF) Mapping mode generates output for use with the CloudMap Variant Discovery or Hawaiian Variant Mapping tools. The aim of VAF analysis is to identify clusters of variants with (near) homozygous inheritance in a F2 population obtained from a cross between a mutant strain of interest and an unrelated mapping strain. Here, the "mapping sample" is the pooled F2 population. To analyze inheritance patterns this mode **requires either** a list of variants that could have been inherited through the mapping strain, i.e. the "unrelated parent strain", or through the mutant parent, i.e. through the "related parent strain". If variants are available for both parents, they can be analyzed together for higher mapping accuracy.
+2) Graphical output settings
+
+   .. class:: warningmark
+  
+      To be able to generate plots the system running MiModD needs to have the statistical programming environment R and its Python interface rpy2 installed.
+
 
-3) More information on combining MiModD and CloudMap in mapping-by-sequencing analyses can be found in the `corresponding section of the MiModD User Guide`_.
+   *y-axes scaling* - if you want to override the defaults
+
+   *x-axis scaling* - choose *preserve relative contig sizes* if you want the largest chromosome to fit the page width and smaller chromosomes to appear according to their relative size or choose *scale each contig to fit the plot width* if all chromosomes should exploit the available space
+
+   *span value to be used in calculating the Loess regression line* - this value determines the degree of smoothing of the regression line through the scatterplot data. Information on loess regression and the loess span parameter can be found at http://en.wikipedia.org/wiki/Local_regression. The default is 0.1 as in CloudMap.
+
+   *colors used for plotting* - can be selected freely from the offered palette. For histogram colors, the list of selected colors will be used to provide the colors for the different histograms plotted. If less colors than histograms (determined by the number of bin sizes selected) are specified, colors from the list will be recycled.
+
 
 .. _CloudMap: https://usegalaxy.org/u/gm2123/p/cloudmap
-.. _corresponding section of the MiModD User Guide: http://mimodd.readthedocs.org/en/latest/cloudmap.html
-
+.. _mapping-by-sequencing analysis workflows in MiModD: http://mimodd.readthedocs.org/en/latest/cloudmap.html
   </help>
 </tool>
--- a/convert.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/convert.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="convert" name="Convert" version="0.1.6.1">
+<tool id="convert" name="Convert" version="0.1.7.0">
   <description>between different sequence data formats</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/covstats.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/covstats.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="coverage_stats" name="Coverage Statistics" version="0.1.6.1">
+<tool id="coverage_stats" name="Coverage Statistics" version="0.1.7.0">
   <description>Calculate coverage statistics for a BCF file as generated by the Variant Calling tool</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/deletion_predictor.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/deletion_predictor.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="deletion_prediction" name="Deletion Prediction for paired-end data" version="0.1.6.1">
+<tool id="deletion_prediction" name="Deletion Prediction for paired-end data" version="0.1.7.0">
   <description>Predicts deletions in one or more aligned read samples based on coverage of the reference genome and on insert sizes</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/fileinfo.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/fileinfo.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="fileinfo" name="Retrieve File Information" version="0.1.6.1">
+<tool id="fileinfo" name="Retrieve File Information" version="0.1.7.0">
   <description>for supported data formats.</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/reheader.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/reheader.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="reheader" name="Reheader BAM file" version="0.1.6.1">
+<tool id="reheader" name="Reheader BAM file" version="0.1.7.0">
   <description>From a BAM file generate a new file with the original header (if any) replaced or modified by that found in a second SAM file</description>
   <expand macro="requirements"/>
   <version_command>mimodd version -q</version_command>
--- a/sam_header.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/sam_header.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="ngs_run_annotation" name="NGS Run Annotation" version="0.1.6.1">
+<tool id="ngs_run_annotation" name="NGS Run Annotation" version="0.1.7.0">
   <description>Create a SAM format header from run metadata for sample annotation.</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/snap_caller.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/snap_caller.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="read_alignment" name="SNAP Read Alignment" version="0.1.6.1">
+<tool id="read_alignment" name="SNAP Read Alignment" version="0.1.7.0">
   <description>Map sequence reads to a reference genome using SNAP</description>
   <macros>
     <import>toolshed_macros.xml</import>
@@ -18,7 +18,7 @@
 		#end if
 --ofile '$outputfile' --iformat ${i.mode_choose.input.iformat} --oformat $oformat
 --idx-seedsize '$set.seedsize'
---idx-slack '$set.slack' --maxseeds '$set.maxseeds' --maxhits '$set.maxhits' --clipping=$set.clipping --maxdist '$set.maxdist' --confdiff '$set.confdiff' --confadapt '$set.confadpt'
+--idx-slack '$set.slack' --maxseeds '$set.maxseeds' --maxhits '$set.maxhits' --clipping $set.clipping --maxdist '$set.maxdist' --confdiff '$set.confdiff' --confadapt '$set.confadpt'
 		#if $i.mode_choose.input.header:
 --header '${i.mode_choose.input.header}'
 		#end if
@@ -168,9 +168,9 @@
     	<param name="maxseeds" type="integer" value="25" label="maximum seeds per read (default: 25)" help="Number of seeds to use per read (SNAP option -n) when trying to match it to the reference genome; higher numbers will increase the rate of aligned reads and reduce the rate of misalignments, but will reduce performance."/>
 	    <param name="clipping" type="select" label="read clipping (default: from back and front)" help="Specifies from which end of a read low-quality bases should be clipped (SNAP option -Cxx)">
 	      <option value="++">from back and front</option>
-	      <option value="-+">from back only</option>
-	      <option value="+-">from front only</option>
-	      <option value="--">no clipping</option>
+	      <option value="x+">from back only</option>
+	      <option value="+x">from front only</option>
+	      <option value="xx">no clipping</option>
 	    </param>
 	    <param name="selectivity" type="integer" value="1" label="selectivity (default: 1)" help="randomly choose 1/selectivity of the reads to score (SNAP option -S). The tool uses the default of 1 (or a 0 setting) to indicate that all reads should be worked with." />
 	    <param name="filter_output" type="select" label="filter output (default: no filtering)" help="filter output (SNAP option -F for certain classes of reads.">
--- a/snp_caller_caller.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/snp_caller_caller.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="variant_calling" name="Variant Calling" version="0.1.6.1">
+<tool id="variant_calling" name="Variant Calling" version="0.1.7.0">
   <description>From a reference and aligned reads generate a BCF file with position-specific variant likelihoods and coverage information</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/snpeff_genomes.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/snpeff_genomes.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="snpeff_genomes" name="List Installed SnpEff Genomes" version="0.1.6.1">
+<tool id="snpeff_genomes" name="List Installed SnpEff Genomes" version="0.1.7.0">
   <description>Checks the local SnpEff installation to compile a list of currently installed genomes</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/tool_dependencies.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/tool_dependencies.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -7,10 +7,10 @@
     <repository changeset_revision="1c337560fa56" name="package_python3_zlib_dependent_1_0" owner="wolma" prior_installation_required="True" toolshed="https://toolshed.g2.bx.psu.edu" />
   </package>
     
-  <package name="mimodd" version="0.1.6.1">
+  <package name="mimodd" version="0.1.7.0">
     <install version="1.0">
       <actions>
-        <action type="download_by_url">http://sourceforge.net/projects/mimodd/files/MiModD-0.1.6.1.tar.gz</action>
+        <action type="download_by_url">http://sourceforge.net/projects/mimodd/files/MiModD-0.1.7.0.tar.gz</action>
         <action type="set_environment_for_install">
           <repository changeset_revision="1c337560fa56" name="package_python3_zlib_dependent_1_0" owner="wolma" toolshed="https://toolshed.g2.bx.psu.edu">
             <package name="python3" version="3.4.1" />
--- a/toolshed_macros.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/toolshed_macros.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,7 +1,7 @@
 <macros>
   <xml name="requirements">
     <requirements>
-      <requirement type="package" version="0.1.6.1">mimodd</requirement>
+      <requirement type="package" version="0.1.7.0">mimodd</requirement>
     </requirements>
   </xml>
 </macros>
--- a/varextract.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/varextract.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="extract_variants" name="Extract Variant Sites" version="0.1.6.1">
+<tool id="extract_variants" name="Extract Variant Sites" version="0.1.7.0">
   <description>from a BCF file</description>
   <macros>
     <import>toolshed_macros.xml</import>
--- a/vcf_filter.xml	Tue Jul 28 23:21:11 2015 +0200
+++ b/vcf_filter.xml	Fri Sep 11 16:25:49 2015 +0200
@@ -1,4 +1,4 @@
-<tool id="vcf_filter" name="VCF Filter" version="0.1.6.1">
+<tool id="vcf_filter" name="VCF Filter" version="0.1.7.0">
   <description>Extracts lines from a vcf variant file based on field-specific filters</description>
   <macros>
     <import>toolshed_macros.xml</import>