--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/seqsero2.xml	Tue Nov 06 09:45:57 2018 -0500
@@ -0,0 +1,111 @@
+<tool id="seqsero_v2" name="SeqSero 2" version="2.0">
+    <description>Salmonella serotype prediction</description>
+    <requirements>
+      <requirement type="package" version="2.7.10">python</requirement>
+      <requirement type="package" version="1.60">biopython</requirement>
+      <requirement type="package" version="2.2.31">blast</requirement>
+      <requirement type="package" version="1.9">samtools</requirement>
+      <requirement type="package" version="2.6.2">sra-tools</requirement>
+      <requirement type="package" version="0.7.17">bwa</requirement>
+      <requirement type="package" version="3.9.0">spades</requirement>
+      <requirement type="package" version="2.27.1">bedtools</requirement>
+    </requirements>
+    <command detect_errors="exit_code"><![CDATA[
+      #if $reads.reads_select == 'history'
+      #set name = $reads.forward.name.split('_')[0]
+      #set forward = $reads.forward
+      #set reverse = $reads.reverse
+      #else
+      #set name = $reads.coll.name
+      #set forward = $reads.coll.forward
+      #set reverse = $reads.coll.reverse
+      #end if
+      #if forward.is_of_type('fastq.gz', 'fastqsanger.gz')
+      gunzip -c $forward > forward.fastq;
+      #set forward = './forward.fastq'
+      #end if
+      #if reverse.is_of_type('fastq.gz', 'fastqsanger.gz')
+      gunzip -c $reverse > reverse.fastq;
+      #set reverse = './reverse.fastq'
+      #end if
+      ln -s forward.fastq ${name}_1.fastq;
+      ln -s reverse.fastq ${name}_2.fastq;
+      $__tool_directory__/SeqSero2.py
+      -p \${GALAXY_SLOTS:-4}
+      -i ${name}_1.fastq ${name}_2.fastq
+      -b $maptype && 
+      cat  SeqSero_result*/Seqsero_result.tsv > results.tsv 
+
+    ]]></command>
+    <inputs>
+      
+        
+        <conditional name="reads">
+            <param name="reads_select" type="select" label="Paired-end collection, or two datasets from your history">
+                <option value="collection">Paired collection from your history</option>
+                <option value="history">Two FASTQ datasets from your history</option>
+            </param>
+            <when value="collection">
+                <param label="Paired reads" name="coll" type="data_collection" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz" collection_type="paired" />
+            </when>
+            <when value="history">
+                <param label="Forward reads" type="data" name="forward" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz" />
+                <param label="Reverse reads" type="data" name="reverse" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz" />
+            </when>
+        </conditional>
+     
+          <!-- <param name="fastq1" type="data" format="fastq" label="FASTQ paired end read 1" />
+          <param name="fastq2" type="data" format="fastq" label="FASTQ paired end read 2" /> -->
+          <!-- <param name="numofthr" type="select" label="Number of threads">
+          <option value="1">1</option>
+          <option value="2">2</option>
+          <option value="3">3</option>
+          <option value="4">4</option> -->
+        <!-- </param> -->
+        <param name="maptype" type="select" label="Algorithms for BWA mapping?">
+          <option value="mem">mem</option>
+          <option value="sam">sam</option>
+        </param>
+    
+
+    </inputs>
+    <outputs>
+      <data format="tsv" label="SeqSero Results" name="results" from_work_dir="results.tsv"/>
+    </outputs>
+    <tests>
+       <test>
+         <param name="reads_select" value="history" />
+         <param name="forward" value="forward_25k.fastq.gz" ftype="fastqsanger.gz" />
+         <param name="reverse" value="reverse_25k.fastq.gz" ftype="fastqsanger.gz" />
+       </test>
+       <test>
+        <param name="reads_select" value="collection" />
+        <param name="coll">
+            <collection type="paired">
+                <element name="forward" value="forward_25k.fastq.gz" ftype="fastqsanger.gz" />
+                <element name="reverse" value="reverse_25k.fastq.gz" ftype="fastqsanger.gz" />
+            </collection>
+        </param>
+       </test>
+    </tests>
+    <help><![CDATA[
+    
+**Usage: SeqSero2.py** 
+
+**Algorithms for BWA mapping**
+
+'mem' for mem, 'sam' for samse/sampe; default=mem; optional; for now SeqSero2 is only optimized for "mem" mode
+   
+    ]]></help>
+    <citations>
+       <citation type="bibtex">
+        @misc{zhang_yin_jones_zhang_deathrage_dinsmore_fitzgeral_fields_deng_2015,
+        title={Salmonella serotype determination utilizing high-throughput genome sequencing data.},
+        journal={J Clin Microbiol}, publisher={ASM},
+        author={Zhang S, Yin Y, Jones MB, Zhang Z, Deatherage Kaiser BL, Dinsmore BA, Fitzgerald C, Fields PI, Deng X.},
+        year={2015}, month={Max},
+        url={http://http://jcm.asm.org/content/early/2015/03/05/JCM.00323-15}},
+       }</citation>
+    </citations>
+
+</tool>