diff snippy.xml @ 0:0801bffdfcc8 draft

Initial upload

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/snippy.xml	Sun Jun 05 21:12:39 2016 -0400
@@ -0,0 +1,159 @@
+<tool id="snippy" name="snippy" version="0.2.0">
+    <requirements>
+      <requirement type="package" version="3.0">snippy</requirement>
+      <requirement type="package" version="1.2">samtools</requirement>
+      <requirement type="package" version="0.9.20">freebayes</requirement>
+    </requirements>
+    <stdio>
+        <exit_code range="1:" />
+    </stdio>
+
+    <command><![CDATA[
+      cp $ref foo.fna &&
+      snippy
+      --outdir out
+      --cpus "\${GALAXY_SLOTS:-1}"
+      --ref foo.fna
+      $cleanup
+      #if str( $advanced.is_advanced ) == "advanced"
+        --mapqual $advanced.mapqual
+        --mincov $advanced.mincov
+        --minfrac $advanced.minfrac
+        #if $advanced.rgid
+          --rgid $advanced.rgid
+        #end if
+        #if $advanced.bwaopt
+          --bwaopt $advanced.bwaopt
+        #end if
+      #end if
+      #if str( $fastq_input.fastq_input_selector ) == "paired"
+        --pe1 $fastq_input.fastq_input1
+        --pe2 $fastq_input.fastq_input2
+      #end if
+      #if str( $fastq_input.fastq_input_selector ) == "paired_collection"
+        --pe1 $fastq_input.fastq_input1.forward
+        --pe2 $fastq_input.fastq_input1.reverse
+      #end if
+      #if str( $fastq_input.fastq_input_selector ) == "single"
+        --se $fastq_input.fastq_input1
+      #end if
+      #if str( $fastq_input.fastq_input_selector ) == "paired_iv"
+        --peil $fastq_input.fastq_input1
+      #end if
+
+      &&
+
+      gunzip out/snps.depth.gz
+
+
+    ]]></command>
+    <inputs>
+      <param name="ref" type="data" format="fasta" label="Reference Fasta" help="Fasta file to use as the reference" />
+      <conditional name="fastq_input">
+        <param name="fastq_input_selector" type="select" label="Single or Paired-end reads" help="Select between paired and single end data">
+          <option value="paired">Paired</option>
+          <option value="single">Single</option>
+          <option value="paired_collection">Paired Collection</option>
+          <option value="paired_iv">Paired Interleaved</option>
+        </param>
+        <when value="paired">
+          <param name="fastq_input1" type="data" format="fastqsanger,fasta" label="Select first set of reads" help="Specify dataset with forward reads"/>
+          <param name="fastq_input2" type="data" format="fastqsanger,fasta" label="Select second set of reads" help="Specify dataset with reverse reads"/>
+        </when>
+        <when value="single">
+          <param name="fastq_input1" type="data" format="fastqsanger,fasta" label="Select fastq dataset" help="Specify dataset with single reads"/>
+        </when>
+        <when value="paired_collection">
+          <param name="fastq_input1" format="fastqsanger,fasta" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
+        </when>
+        <when value="paired_iv">
+          <param name="fastq_input1" type="data" format="fastqsanger" label="Select fastq dataset" help="Specify dataset with interleaved reads"/>
+        </when>
+      </conditional>
+      <param name="cleanup" type="boolean" checked="true" truevalue="--cleanup" falsevalue="" label="Cleanup the non-snp output files" help="Remove all non-SNP files: BAMs, indices etc" />
+      <conditional name="advanced">
+        <param name="is_advanced" type="select" label="Advanced parameters" help="unhide advanced parameter settings">
+          <option value="advanced">Show advanced settings</option>
+          <option value="simple" selected="true">Hide advanced settings</option>
+        </param>
+        <when value="advanced">
+          <param name="mapqual" type="float" value="60" label="Minimum mapping quality" help="Minimum mapping quality to allow" />
+          <param name="mincov" type="float" value="10" label="Minimum coverage" help="Minimum coverage to call a snp" />
+          <param name="minfrac" type="float" value="0.9" label="Minumum proportion for variant evidence" help="Minumum proportion for variant evidence" />
+          <param name="rgid" type="text" value="" label="Bam header @RG ID" help="Use this @RG ID: in the BAM header" />
+          <param name="bwaopt" type="text" value="" label="Extra BWA MEM options" help="Extra BWA MEM options, eg. -x pacbio" />
+        </when>
+        <when value="simple">
+
+        </when>
+      </conditional>
+    </inputs>
+    <outputs>
+      <data format="vcf" name="snpvcf" label="${tool.name} on ${on_string} snps vcf file" from_work_dir="out/snps.vcf"/>
+      <data format="gff3" name="snpgff" label="${tool.name} on ${on_string} snps gff file" from_work_dir="out/snps.gff"/>
+      <data format="tabular" name="snptab" label="${tool.name} on ${on_string} snps table" from_work_dir="out/snps.tab"/>
+      <data format="tabular" name="snpsum" label="${tool.name} on ${on_string} snps summary" from_work_dir="out/snps.txt"/>
+      <data format="text" name="snplog" label="${tool.name} on ${on_string} log file" from_work_dir="out/snps.log"/>
+      <data format="fasta" name="snpalign" label="${tool.name} on ${on_string} aligned fasta" from_work_dir="out/snps.aligned.fa"/>
+      <data format="fasta" name="snpconsensus" label="${tool.name} on ${on_string} consensus fasta" from_work_dir="out/snps.consensus.fa"/>
+      <data format="tabular" name="snpsdepth" label="${tool.name} on ${on_string} mapping depth" from_work_dir="out/snps.depth"/>
+      <data format="bam" name="snpsbam" label="${tool.name} on ${on_string} mapped reads (bam)" from_work_dir="out/snps.bam">
+        <filter>cleanup is False</filter>
+      </data>
+    </outputs>
+
+    <tests>
+      <test>
+        <param name="ref" value="Ecoli.fna" ftype="fasta" />
+        <param name="fastq_input_selector" value="paired" />
+        <param name="fastq_input1" ftype="fastq" value="reads_1.fq" />
+        <param name="fastq_input2" ftype="fastq" value="reads_2.fq" />
+        <output name="snpsum" ftype="tabular" file="test/snps.txt" lines-diff="5" />
+      </test>
+    </tests>
+
+
+    <help><![CDATA[
+      This is a change to force a reinstall
+        Synopsis:
+  snippy 3.0 - fast bacterial variant calling from NGS reads
+Author:
+  Torsten Seemann <[e-mail]>
+Usage:
+  snippy [options] --outdir <dir> --ref <ref> --pe1 <R1.fq.gz> --pe2 <R2.fq.gz>
+  snippy [options] --outdir <dir> --ref <ref> --se <454.fastq>
+  snippy [options] --outdir <dir> --ref <ref> --peil <velvet.fa.gz>
+Options:
+  --help            This help
+  --version         Print version and exit
+  --citation        Print citation for referencing snippy
+  --quiet           No screen output (default OFF)
+  --cpus [N]        Maximum number of CPU cores to use (default '8')
+  --reference [X]   Reference genome. Supports FASTA, GenBank, EMBL (not GFF) (default '')
+  --outdir [X]      Output folder (default '')
+  --prefix [X]      Prefix for output files (default 'snps')
+  --force           Force overwrite of existing output folder (default OFF)
+  --pe1|R1|left [X] Reads, paired-end R1 (left) (default '')
+  --pe2|R2|right [X] Reads, paired-end R2 (right) (default '')
+  --se|single [X]   Single-end reads (default '')
+  --peil [X]        Reads, paired-end R1/R2 interleaved (default '')
+  --mapqual [n.n]   Minimum mapping quality to allow (default '60')
+  --mincov [N]      Minimum coverage of variant site (default '10')
+  --minfrac [n.n]   Minumum proportion for variant evidence (default '0.9')
+  --report          Produce long report with visual alignment (slow) (default OFF)
+  --cleanup         Remove all non-SNP files: BAMs, indices etc (default OFF)
+  --rgid [X]        Use this @RG ID: in the BAM header (default '')
+  --bwaopt [X]      Extra BWA MEM options, eg. -x pacbio (default '')
+
+    ]]></help>
+
+    <citations>
+      <citation type="bibtex">@UNPUBLISHED{Seemann2013,
+        author = "Seemann T",
+        title = "snippy: fast bacterial variant calling from NGS reads",
+        year = "2015",
+        note = "https://github.com/tseemann/snippy"}
+      </citation>
+    </citations>
+
+</tool>