--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/kma_map.xml	Sat Feb 15 15:32:58 2020 -0500
@@ -0,0 +1,179 @@
+<tool id="kma_map" name="Map with KMA" version="@TOOL_VERSION@+galaxy0">
+    <description></description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <requirements>
+        <requirement type="package" version="@TOOL_VERSION@">kma</requirement>
+    </requirements>
+    <version_command>kma -v</version_command>
+    <command detect_errors="exit_code">
+        <![CDATA[
+        kma
+            -t \${GALAXY_SLOTS:-1}
+	    -t_db '${kma_index.fields.path}'
+            #if $single_paired.single_paired_selector == 'paired'
+                -ipe '${single_paired.forward_input}' '${single_paired.reverse_input}'
+            #elif $single_paired.single_paired_selector == "paired_collection":
+                -ipe '${single_paired.input_pair.forward}' '${single_paired.input_pair.reverse}'
+            #else:
+                -i '${single_paired.input_sequences}'
+            #end if
+            #if str($settings.advanced) == "advanced"
+              #if str($settings.kmer_size)
+                -k '${settings.kmer_size}'
+              #end if
+              #if str($settings.p_value)
+                -p '${settings.p_value}'
+              #end if
+              ${settings.decontaminate}
+              ${settings.dense}
+              ${settings.ref_fsa}
+              ${settings.matrix}
+              ${settings.all_best_mappings}
+              #if str($settings.minimum_phred_score)
+                -mp '${settings.minimum_phred_score}'
+              #end if
+              #if str($settings.cut_5_prime)
+                -5p '${settings.cut_5_prime}'
+              #end if
+              ${settings.only_count_kmers}
+              #if str($settings.min_id)
+                -ID '${settings.min_id}'
+              #end if
+              #if str($settings.base_call_depth)
+                -bcd '${settings.base_call_depth}'
+              #end if
+              #if str($settings.minimum_mapping_quality)
+                -mq '${settings.minimum_mapping_quality}'
+              #end if
+              #if str($settings.reward)
+                -reward '${settings.reward}'
+              #end if
+              #if str($settings.penalty)
+                -penalty '${settings.penalty}'
+              #end if
+              #if str($settings.gapopen)
+                -gapopen '${settings.gapopen}'
+              #end if
+              #if str($settings.gapextend)
+                -gapextend '${settings.gapextend}'
+              #end if
+              ${settings.force_end_to_end}
+              ${settings.set_cge_penalties_and_rewards}
+            #end if
+            -o output
+	 
+        #if str($settings.advanced) == "advanced" and $settings.matrix
+            && gunzip output.mat.gz
+        #end if
+        && gunzip output.frag.gz
+	]]>
+    </command>
+    <inputs>
+        <conditional name="single_paired">
+            <param name="single_paired_selector" type="select" label="Single or paired reads" help="--paired">
+                <option value="paired_collection">Paired collection</option>
+                <option value="paired">Paired-end data</option>
+                <option selected="True" value="single">Single-end data</option>
+            </param>
+            <when value="paired_collection">
+                <param format="@INTYPES@" name="input_pair" type="data_collection" collection_type="paired" label="Collection of paired reads" help="FASTQ datasets" />
+            </when>
+            <when value="paired">
+                <param format="@INTYPES@" name="forward_input" type="data" label="Forward strand" help="FASTQ dataset"/>
+                <param format="@INTYPES@" name="reverse_input" type="data" label="Reverse strand" help="FASTQ dataset"/>
+            </when>
+            <when value="single">
+                <param format="@INTYPES@" label="Input sequences" name="input_sequences" type="data" help="FASTQ datasets"/>
+            </when>
+        </conditional>
+        <param name="kma_index" type="select">
+            <options from_data_table="kma_index">
+                <validator type="no_options" message="No KMA index available" />
+            </options>
+        </param>
+        <conditional name="settings">
+            <param name="advanced" type="select" label="Specify advanced parameters">
+                <option value="simple" selected="true">No, use program defaults</option>
+                <option value="advanced">Yes, see full parameter list.</option>
+            </param>
+            <when value="simple">
+            </when>
+            <when value="advanced">
+                <param name="kmer_size" type="integer" min="4" value="16" max="32" label="Kmer Size" />
+                <param name="p_value" type="float" min="0.0" value="0.05" max="1.0" label="p-value"/>
+                <param name="exhaustive_mode" type="boolean" truevalue="-ex_mode" falsevalue="" label="Exhaustive Mode" />
+                <param name="decontaminate" type="boolean" truevalue="-deCon" falsevalue="" label="Decontaminate" />
+                <param name="dense" type="boolean" truevalue="-dense" falsevalue="" label="Do not allow insertions in assembly" />
+                <param name="ref_fsa" type="boolean" truevalue="-ref_fsa" falsevalue="" label="Consensus sequence has 'n' instead of gaps" />
+                <param name="matrix" type="boolean" truevalue="-matrix" falsevalue="" label="Output assembly matrix" />
+                <param name="all_best_mappings" type="boolean" truevalue="-a" falsevalue="" label="Print all best mappings" />
+                <param name="minimum_phred_score" type="integer" min="0" value="20" max="60" label="Minimum phred score" />
+                <param name="cut_5_prime" type="integer" min="0" value="0" max="64" label="Cut a constant number of nucleotides from the 5 prime" />
+                <param name="only_count_kmers" type="boolean" truevalue="-Sparse" falsevalue="" label="Only count kmers" />
+                <param name="min_id" type="float" min="0.0" value="1.0" max="100.0" label="Minimum percent identity"/>
+                <param name="force_end_to_end" type="boolean" truevalue="-1t1" falsevalue="" label="Force end to end mapping" />
+                <param name="base_call_depth" type="integer" min="1" value="1" max="1000" label="Minimum depth at base" />
+                <param name="minimum_mapping_quality" type="integer" min="0" value="0" max="100" label="Minimum mapping quality" />
+                <param name="reward" type="integer" min="1" value="1" max="100" label="Score for match" />
+                <param name="penalty" type="integer" min="-100" value="-2" max="0" label="Penalty for mismatch" />
+                <param name="gapopen" type="integer" min="-100" value="-3" max="0" label="Penalty for gap opening" />
+                <param name="gapextend" type="integer" min="-100" value="-1" max="0" label="Penalty for gap extension" />
+                <param name="pairing_reward" type="integer" min="1" value="7" max="100" label="Reward for pairing reads" />
+                <param name="set_cge_penalties_and_rewards" type="boolean" truevalue="-cge" falsevalue="" label="Set CGE penalties and rewards" />
+            </when>
+        </conditional>
+    </inputs>
+    <outputs>
+        <data name="result_overview" label="Result overview" format="tabular" from_work_dir="output.res" />
+        <data name="consensus_alignment" label="Consensus alignment" format="fasta" from_work_dir="output.aln" />
+        <data name="consensus_sequences" label="Consensus sequences" format="fasta" from_work_dir="output.fsa" />
+        <data name="read_mapping" label="Read mapping info" format="tabular" from_work_dir="output.frag" />
+        <data name="assembly_matrix" label="Assembly matrix" format="txt" from_work_dir="output.mat">
+            <filter>settings['matrix']</filter>
+        </data>
+    </outputs>
+    <tests>
+        <test>
+            <param name="single_paired_selector" value="single"/>
+            <param name="input_sequences" value="ERR884056_ecoli_b0842.mapped_R1.fastq" ftype="fastq"/>
+            <param name="advanced" value="advanced"/>
+            <param name="kmer_size" value="8"/>
+            <param name="kma_index" value="test_index"/>
+            <output name="result_overview" file="ERR884056.res" ftype="tabular"/>
+            <output name="consensus_alignment" file="ERR884056.aln" ftype="fasta"/>
+            <output name="consensus_sequences" file="ERR884056.fsa" ftype="fasta"/>
+            <output name="read_mapping" file="ERR884056.frag" ftype="tabular"/>
+        </test>
+        <test>
+            <param name="single_paired_selector" value="single"/>
+            <param name="input_sequences" value="ERR884056_ecoli_b0842.mapped_R1.fastq" ftype="fastq"/>
+            <param name="advanced" value="advanced"/>
+            <param name="kmer_size" value="8"/>
+            <param name="matrix" value="true"/>
+            <param name="kma_index" value="test_index"/>
+            <output name="assembly_matrix" file="ERR884056.mat" ftype="txt"/>
+        </test>
+    </tests>
+    <help>
+      <![CDATA[
+
+When the mapping is done KMA will produce the following files:
+
+    *.res A result overview giving the most common statistics for each mapped template.
+    *.fsa The consensus sequences drawn from the alignments.
+    *.aln The consensus alignment of the reads against their template.
+    *.frag Mapping information on each mapped read, columns are: 
+        1. read
+        2. number of equally well mapping templates
+        3. mapping score
+        4. start position
+        5. end position (w.r.t. template)
+        6. the choosen template.
+    *.mat Base counts on each position in each template, (only if “-matrix” is enabled)
+
+      ]]>
+    </help>
+    <expand macro="citations" />
+</tool>