--- a/ribowaltz.xml	Thu Sep 22 20:30:54 2022 +0000
+++ b/ribowaltz.xml	Mon Oct 23 12:10:31 2023 +0000
@@ -1,16 +1,16 @@
-<tool id="ribowaltz_process" name="riboWaltz" version="@VERSION@" profile="@PROFILE@">
+<tool id="ribowaltz_process" name="riboWaltz" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
     <description>calculation of optimal P-site offsets and diagnostic analysis</description>
     <macros>
         <import>macros.xml</import>
     </macros>
-    <expand macro='requirements'/>
     <expand macro='edam_ontology' />
     <expand macro='xrefs'/>
+    <expand macro='requirements'/>
     <command detect_errors="exit_code"><![CDATA[
         #for $i, $sample in enumerate($rep_samples):
             ln -s $sample.bam_file $sample.sample_name'.bam' &&
         #end for
-        Rscript '${__tool_directory__}/ribowaltz.R' -b . -g '$gtf' --refseq_sep '$refseq_sep'
+        Rscript '${__tool_directory__}/ribowaltz.R' -b . -g '$gtf' --refseq_sep '$refseq_sep' --indel_threshold $indel_threshold
         #import json
         #if $filtering.duplicates.filter == 'yes':
             #set params_duplicate_filterting = []
@@ -59,6 +59,8 @@
         <param name="gtf" type="data" format="gtf,gff" label="Annotation in GTF format"/>
         <param name="refseq_sep" type="text" optional="true" label="Separator between reference sequences' name and additional information to discard"
             help=" All characters before the first occurrence of the specified separator are kept"/>
+        <param name="indel_threshold" type="integer" value="5" min="0"
+                label="Positive integer value specifying the maximum number of indels allowed for each read."/>
         <param name="save_rdata" type="boolean" truevalue="1" falsevalue="0" checked="true"
             label="Save p-site info RDATA file?"
             help="Useful for advanced plotting using riboWaltz-plot tool"/>
@@ -198,6 +200,46 @@
                 </assert_contents>
             </output>
         </test>
+        <test expect_num_outputs="5">
+            <param name="gtf" value="rep1_annot.gtf.gz"/>
+            <param name="refseq_sep" value="."/>
+            <param name="indel_threshold" value="3"/>
+            <param name="save_rdata" value="true"/>
+            <repeat name="rep_samples">
+                <param name="sample_name" value="Replicate1"/>
+                <param name="bam_file" value="rep1.bam"/>
+            </repeat>
+            <section name="filtering"/>
+            <section name="psite_additional">
+                <param name="flanking" value="3"/>
+                <param name="use_start" value="true"/>
+                <param name="psite_extrimity" value="auto"/>
+                <param name="cl" value="99"/>
+            </section>
+            <section name="coverage_additional">
+                <param name="psites_per_region" value="true"/>
+                <param name="min_overlap" value="1"/>
+                <param name="start_nts" value="0"/>
+                <param name="stop_nts" value="0"/>
+            </section>
+            <output name="codon_coverage_out">
+                <assert_contents>
+                    <has_text_matching expression="transcript\tstart\tend\tfrom_cds_start\tfrom_cds_stop\tregion\tReplicate1"/>
+                    <has_text_matching expression="ENSMUST00000000137\t2656\t2659\t865\t471\t3utr\t0"/>
+                </assert_contents>
+            </output>
+            <output name="cds_coverage_out">
+                <assert_contents>
+                    <has_text_matching expression="transcript\tlength_cds\tReplicate1"/>
+                    <has_text_matching expression="ENSMUST00000019109\t741\t5"/>
+                </assert_contents>
+            </output>
+            <output name="psite_rdata_out">
+                <assert_contents>
+                    <has_size value="169141" delta="1000"/>
+                </assert_contents>
+            </output>
+        </test>        
     </tests>
     <help><![CDATA[
 riboWaltz is an R package for calculation of optimal P-site offsets, diagnostic analysis and visual inspection of ribosome profiling data. Taking advantage of a two-step algorithm where offset information is passed through populations of reads with different length in order to maximize offset coherence, riboWaltz computes with high precision the P-site offset. riboWaltz also provides a variety of graphical representations, laying the foundations for further accurate RiboSeq analyses and improved interpretation of positional information. More information can be found here: https://github.com/LabTranslationalArchitectomics/riboWaltz