# HG changeset patch
# User wolma
# Date 1418509216 18000
# Node ID b2d889cd87f013ead81d7737215cd00cb8548790
Imported from capsule None
diff -r 000000000000 -r b2d889cd87f0 snap_caller.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/snap_caller.xml Sat Dec 13 17:20:16 2014 -0500
@@ -0,0 +1,235 @@
+
+ Map sequence reads to a reference genome using SNAP
+
+ mimodd
+
+ mimodd version -q
+
+ mimodd snap-batch -s
+ ## SNAP calls (considering different cases)
+
+ #for $i in $datasets
+ "snap ${i.mode_choose.mode} $ref_genome
+ #if $str($i.mode_choose.mode) == "paired" and $str($i.mode_choose.input.iformat) in ("fastq", "gz"):
+${i.mode_choose.input.ifile1} ${i.mode_choose.input.ifile2}
+ #else:
+${i.mode_choose.input.ifile}
+ #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
+ #if $i.mode_choose.input.header:
+--header ${i.mode_choose.input.header}
+ #end if
+ #if $str($i.mode_choose.mode) == "paired":
+--spacing $set.sp_min $set.sp_max
+ #end if
+ #if $str($set.selectivity) != "off":
+--selectivity $set.selectivity
+ #end if
+ #if $str($set.filter_output) != "off":
+--filter-output $set.filter_output
+ #end if
+ #if $str($set.sort) != "off":
+--sort $set.sort
+ #end if
+ #if $str($set.mmatch_notation) == "general":
+-M
+ #end if
+--max-mate-overlap $set.max_mate_overlap
+--verbose
+"
+ #end for
+
+
+
+ ## mandatory arguments (and mode-conditionals)
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ ## optional arguments
+
+
+
+
+
+
+
+ ## default settings
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ ## change settings
+
+
+
+
+
+ ## paired-end specific options
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+.. class:: infomark
+
+ **What it does**
+
+The tool aligns the sequenced reads in an arbitrary number of input datasets against a common reference genome and stores the results in a single, possibly multi-sample output file. It supports a variety of different sequenced reads input formats, i.e., SAM, BAM, fastq and gzipped fastq, and both single-end and paired-end data.
+
+Internally, the tool uses the ultrafast, hashtable-based aligner SNAP (http://snap.cs.berkeley.edu), hence its name.
+
+**Notes:**
+
+1) In its standard configuration Galaxy will decompress any .gz files during their upload, so the option to align gzipped fastq input is useful only with customized Galaxy instances or by using linked files as explained in our `recipe for using gzipped fastq files in Galaxy`_ from the `MiModD user guide`_.
+
+2) To use paired-end fastq data with the tool the read mate information needs to be split over two fastq files in corresponding order.
+
+ **TIP:** If your paired-end data is arranged differently, you may look into the *fastq splitter* and *fastq de-interlacer* tools for Galaxy from the `Fastq Manipulation category`_ of the Galaxy Tool Shed to see if they can convert your files to the expected format.
+
+3) The tool supports the alignment of reads from the same sequencing run, but distributed across several input files.
+
+ Generally, it expects the reads from each input dataset to belong to one read-group and will abort with an error message if any input dataset declares more than one read group or sample names in its header. Different datasets, however, are allowed to contain reads from the same read-group (as indicated by matching read-group IDs and sample names in their headers), in which case the reads will be combined into one group in the output.
+
+4) Read-group information is required for every input dataset!
+
+ We generally recommend to store NGS datasets in SAM/BAM format with run metadata stored in the file header. You can use the *NGS Run Annotation* and *Convert* tools to convert data in fastq format to SAM/BAM with added run information.
+
+ While it is not our recommended approach, you can, if you prefer it, align reads from fastq files or SAM/BAM files without header read-group information. To do so, you **must** specify a SAM file that provides the missing information in its header along with the input dataset. You can generate a SAM header file with the *NGS Run Annotation* tool.
+
+ Optionally, a SAM header file can also be used to replace existing read-group information in a headered SAM/BAM input file. This can be used to resolve read-group ID conflicts between multiple input files at tool runtime.
+
+4) Currently, you cannot configure aligner-specific options separately for specific input files from within this Galaxy tool. If you need this advanced level of control, you should use the command line tool ``mimodd snap-batch``.
+
+.. _Fastq Manipulation category: https://toolshed.g2.bx.psu.edu/repository/browse_repositories_in_category?id=310ff67d4caf6531
+.. _recipe for using gzipped fastq files in Galaxy: http://mimodd.readthedocs.org/en/latest/recipes.html#use-gzipped-fastq-files-in-galaxy
+.. _MiModD user guide: http://mimodd.readthedocs.org/en/latest
+
+
+
+
diff -r 000000000000 -r b2d889cd87f0 tool_dependencies.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml Sat Dec 13 17:20:16 2014 -0500
@@ -0,0 +1,6 @@
+
+
+
+
+
+