--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/get.seqs.xml	Fri May 19 05:44:07 2017 -0400
@@ -0,0 +1,154 @@
+<tool profile="16.07" id="mothur_get_seqs" name="Get.seqs" version="@WRAPPER_VERSION@.0">
+    <description>Picks sequences by name</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements"/>
+    <expand macro="stdio"/>
+    <expand macro="version_command"/>
+    <command><![CDATA[
+        @SHELL_OPTIONS@
+
+        ## create symlinks to input datasets
+        ln -s "$accnos" accnos.dat &&
+        ln -s "$fasta_in" fasta_in.dat &&
+        ln -s "$fastq_in" fastq_in.dat &&
+        ln -s "$count_in" count_in.dat &&
+        ln -s "$qfile_in" qfile_in.dat &&
+        ln -s "$name_in" name_in.dat &&
+        ln -s "$group_in" group_in.dat &&
+        ln -s "$alignreport_in" alignreport_in.dat &&
+        ln -s "$list_in" list_in.dat &&
+        ln -s "$taxonomy_in" taxonomy_in.dat &&
+
+        echo 'get.seqs(
+            accnos=accnos.dat,
+            #if $fasta_in:
+                fasta=fasta_in.dat,
+            #end if
+            #if $fastq_in:
+                fastq=fastq_in.dat,
+            #end if
+            #if $count_in:
+                count=count_in.dat,
+            #end if
+            #if $qfile_in:
+                qfile=qfile_in.dat,
+            #end if
+            #if $name_in:
+                name=name_in.dat,
+            #end if
+            #if $group_in:
+                group=group_in.dat,
+            #end if
+            #if $alignreport_in:
+                alignreport=alignreport_in.dat,
+            #end if
+            #if $list_in:
+                list=list_in.dat,
+            #end if
+            #if $taxonomy_in:
+                taxonomy=taxonomy_in.dat,
+            #end if
+            dups=$dups
+        )'
+        | sed 's/ //g'  ## mothur trips over whitespace
+        | mothur
+        | tee mothur.out.log
+    ]]></command>
+    <inputs>
+        <param name="accnos" type="data" format="mothur.accnos" label="accnos - Accession Names"/>
+        <param name="fasta_in" type="data" format="fasta" optional="true" label="fasta - Fasta Sequences"/>
+        <param name="qfile_in" type="data" format="qual" optional="true" label="qfile - Fasta Quality"/>
+        <param name="fastq_in" type="data" format="fastq" optional="true" label="count - a count_table" help="fastq - allows you to select sequences from your fastq file"/>
+        <param name="count_in" type="data" format="mothur.count_table" optional="true" label="count - a count_table" help="generated by count.seqs"/>
+        <param name="name_in" type="data" format="mothur.names" optional="true" label="name - Sequences Name reference"/>
+        <param name="group_in" type="data" format="mothur.groups" optional="true" label="group - Sequences Groups"/>
+        <param name="alignreport_in" type="data" format="mothur.align.report" optional="true" label="alignreport - Align Report"/>
+        <param name="list_in" type="data" format="mothur.list" optional="true" label="list - OTU List"/>
+        <param name="taxonomy_in" type="data" format="mothur.seq.taxonomy" optional="true" label="taxonomy - Taxonomy"/>
+        <param name="dups" type="boolean" truevalue="dups" falsevalue="false" checked="true" label="dups - Apply to duplicates"/>
+    </inputs>
+    <outputs>
+        <expand macro="logfile-output"/>
+        <data name="fasta_out" format_source="fasta_in" from_work_dir="fasta_in*.pick.*" label="${tool.name} on ${on_string}: pick.fasta">
+            <filter>fasta_in</filter>
+        </data>
+        <data name="fastq_out" format_source="fastq_in" from_work_dir="fastq_in*.pick.*" label="${tool.name} on ${on_string}: pick.fastq">
+            <filter>fastq_in</filter>
+        </data>
+        <data name="count_out" format_source="count_in" from_work_dir="count_in*.pick.*" label="${tool.name} on ${on_string}: pick.count">
+            <filter>count_in</filter>
+        </data>
+        <data name="qfile_out" format_source="qfile_in" from_work_dir="qfile_in*.pick.*" label="${tool.name} on ${on_string}: pick.qfile">
+            <filter>qfile_in</filter>
+        </data>
+        <data name="name_out" format="mothur.names" from_work_dir="name_in*.pick.*" label="${tool.name} on ${on_string}: pick.names">
+            <filter>name_in</filter>
+        </data>
+        <data name="group_out" format="mothur.groups" from_work_dir="group_in*.pick.*" label="${tool.name} on ${on_string}: pick.groups">
+            <filter>group_in</filter>
+        </data>
+        <data name="alignreport_out" format="mothur.align.report" from_work_dir="alignreport_in*.pick.*" label="${tool.name} on ${on_string}: pick.align.report">
+            <filter>alignreport_in</filter>
+        </data>
+        <data name="list_out" format="mothur.list" from_work_dir="list_in*.pick.*" label="${tool.name} on ${on_string}: pick.list">
+            <filter>list_in</filter>
+        </data>
+        <data name="taxonomy_out" format="mothur.seq.taxonomy" from_work_dir="taxonomy_in*.pick.*" label="${tool.name} on ${on_string}: pick.taxonomy">
+            <filter>taxonomy_in</filter>
+        </data>
+    </outputs>
+    <tests>
+        <test>
+            <param name="accnos" value="Mock_S280_L001_R1_001_small.trim.contigs.bad.accnos"/>
+            <param name="fasta_in" value="Mock_S280_L001_R1_001_small.trim.contigs.fasta"/>
+            <param name="dups" value=""/>
+            <output name="fasta_out" md5="f9335b6b70dc639e420d33809e51431d"/>
+            <expand macro="logfile-test"/>
+        </test>
+        <test>
+            <param name="accnos" value="Mock_S280_L001_R1_001_small.trim.contigs.bad.accnos"/>
+            <param name="fastq_in" value="Mock_S280_L001_R1_001_small.fastq"/>
+            <output name="fastq_out" md5="e9aae546bf9c00cb7e884dff267d4ba1"/>
+            <expand macro="logfile-test"/>
+        </test>
+        <test>
+            <!-- test two input files -->
+            <param name="accnos" value="Mock_S280_L001_R1_001_small.trim.contigs.bad.accnos"/>
+            <param name="fasta_in" value="Mock_S280_L001_R1_001_small.trim.contigs.fasta"/>
+            <param name="fastq_in" value="Mock_S280_L001_R1_001_small.fastq"/>
+            <param name="dups" value="false"/>
+            <output name="fasta_out" md5="f9335b6b70dc639e420d33809e51431d"/>
+            <output name="fastq_out" md5="e9aae546bf9c00cb7e884dff267d4ba1"/>
+            <expand macro="logfile-test"/>
+        </test>
+        <test>
+            <param name="accnos" value="amazon.bad.accnos"/>
+            <param name="count_in" value="amazon.count_table"/>
+            <output name="count_out" md5="6892dd99850ce9e9f8f15e77b28f57e2"/>
+            <expand macro="logfile-test"/>
+        </test>
+    </tests>
+    <help>
+<![CDATA[
+
+@MOTHUR_OVERVIEW@
+
+**Command Documentation**
+
+The get.seqs_ command takes a list of sequence names and either a fasta, name_, group_, list_, align.report_ or taxonomy_ file to generate a new file that contains only the sequences in the list. This command may be used in conjunction with the list.seqs_ command to help screen a sequence collection.
+
+.. _name: https://www.mothur.org/wiki/Name_file
+.. _group: https://www.mothur.org/wiki/Group_file
+.. _list: https://www.mothur.org/wiki/List_file
+.. _align.report: https://www.mothur.org/wiki/Align.seqs
+.. _taxonomy: https://www.mothur.org/wiki/Taxonomy_outline
+.. _list.seqs: https://www.mothur.org/wiki/list.seqs
+.. _get.seqs: https://www.mothur.org/wiki/Get.seqs
+
+v.1.27.0 : Updated to Mothur 1.33, added count and fastq params
+]]>
+    </help>
+    <expand macro="citations"/>
+</tool>