diff diamond.xml @ 10:1e3323a44643 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/diamond commit b9f434bf71431a3b9548eb17fcc2639b3b6f8b01"

--- a/diamond.xml	Sat Nov 27 09:48:10 2021 +0000
+++ b/diamond.xml	Fri Apr 22 13:51:34 2022 +0000
@@ -1,4 +1,4 @@
-<tool id="bg_diamond" name="Diamond" version="@TOOL_VERSION@.1" profile="19.01">
+<tool id="bg_diamond" name="Diamond" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="19.01">
     <description>alignment tool for short sequences against a protein database</description>
     <macros>
         <import>macros.xml</import>
@@ -37,9 +37,15 @@
 
         @OUTPUT_ARGS@
 
-        --compress '0'
+        #if $output_section.output.outfmt != '100'
+            --compress '0'
+        #end if
         $sens_cond.sensitivity
-
+        $iterate
+        --algo $algo
+        #if $global_ranking
+            --global-ranking $global_ranking
+        #end if
         #if str($gapopen) != "":
           --gapopen '$gapopen'
         #end if
@@ -62,8 +68,8 @@
         --query-cover '$query_cover'
         --subject-cover '$subject_cover'
         --block-size '$sens_cond.block_size'
-        #if $output_unal
-            #if "--un" in $output_unal
+        #if $output_section.output_unal
+            #if "--un" in $output_section.output_unal
                 --un '$unalqueries'
                 #if $query.ext.startswith("fasta"):
                     --unfmt fasta
@@ -71,7 +77,7 @@
                     --unfmt fastq
                 #end if
             #end if
-            #if "--al" in $output_unal
+            #if "--al" in $output_section.output_unal
                 --al '$alqueries'
                 #if $query.ext.startswith("fasta"):
                     --alfmt fasta
@@ -80,33 +86,41 @@
                 #end if
             #end if
         #end if
+        #if $output_section.max_hsps
+            --max-hsps $output_section.max_hsps
+        #end if
         #if $tax_cond.tax_select == 'file':
             --taxonlist `cat '$tax_cond.taxonlistfile' | grep -v "^#" | grep -v "^$" | tr "\n" "," | sed 's/,$//'`
         #else if  $tax_cond.tax_select == 'list':
             --taxonlist '$tax_cond.taxonlist'
         #end if
+        #if $advanced_section.seed_cut
+            --seed-cut $advanced_section.seed_cut
+        #end if
+        $advanced_section.freq_masking
+        --motif-masking $advanced_section.motif_masking
 ]]>
     </command>
     <inputs>
         <conditional name="method_cond">
-            <param name="method_select" type="select" label="What do you want to align?" help="(blastp/blastx)">
-                <option value="blastp">Align amino acid query sequences (blastp)</option>
-                <option value="blastx">Align DNA query sequences (blastx)</option>
+            <param name="method_select" type="select" label="Alignment mode" help="(blastp/blastx)">
+                <option value="blastp">Amino acid query sequences (blastp)</option>
+                <option value="blastx">DNA query sequences (blastx)</option>
             </param>
             <when value="blastx">
                 <param argument="--query-gencode" type="select" label="Genetic code used for translation of query in BLASTX mode" help="">
-                    <option value="1">The Standard Code</option>
-                    <option value="2">The Vertebrate Mitochondrial Code</option>
-                    <option value="3">The Yeast Mitochondrial Code</option>
-                    <option value="4">The Mold, Protozoan, and Coelenterate Mitochondrial Code and the Mycoplasma/Spiroplasma Code</option>
-                    <option value="5">The Invertebrate Mitochondrial Code</option>
-                    <option value="6">The Ciliate, Dasycladacean and Hexamita Nuclear Code</option>
-                    <option value="9">The Echinoderm and Flatworm Mitochondrial Code</option>
-                    <option value="10">The Euplotid Nuclear Code</option>
-                    <option value="11">The Bacterial, Archaeal and Plant Plastid Code</option>
-                    <option value="12">The Alternative Yeast Nuclear Code</option>
-                    <option value="13">The Ascidian Mitochondrial Code</option>
-                    <option value="14">The Alternative Flatworm Mitochondrial Code</option>
+                    <option value="1">Standard Code</option>
+                    <option value="2">Vertebrate Mitochondrial Code</option>
+                    <option value="3">Yeast Mitochondrial Code</option>
+                    <option value="4">Mold, Protozoan, and Coelenterate Mitochondrial Code and the Mycoplasma/Spiroplasma Code</option>
+                    <option value="5">Invertebrate Mitochondrial Code</option>
+                    <option value="6">Ciliate, Dasycladacean and Hexamita Nuclear Code</option>
+                    <option value="9">Echinoderm and Flatworm Mitochondrial Code</option>
+                    <option value="10">Euplotid Nuclear Code</option>
+                    <option value="11">Bacterial, Archaeal and Plant Plastid Code</option>
+                    <option value="12">Alternative Yeast Nuclear Code</option>
+                    <option value="13">Ascidian Mitochondrial Code</option>
+                    <option value="14">Alternative Flatworm Mitochondrial Code</option>
                     <option value="16">Chlorophycean Mitochondrial Code</option>
                     <option value="21">Trematode Mitochondrial Code</option>
                     <option value="22">Scenedesmus obliquus Mitochondrial Code</option>
@@ -116,7 +130,7 @@
                     <option value="26">Pachysolen tannophilus Nuclear Code</option>
                 </param>
                 <param argument="--min-orf" type="integer" value="1" label="ignore translated sequences without an open reading frame of at least this length" help="By default this feature is disabled for sequences of length below 30, set to 20 for sequences of length below 100, and set to 40 otherwise. Setting this option to 1 will disable this feature" />
-
+                
                 <param name="query_strand" argument="--strand" type="select" label="query strands to search" help="">
                     <option value="both" selected="True">Both</option>
                     <option value="plus">Plus</option>
@@ -140,7 +154,9 @@
                 </param>
             </when>
             <when value="blastp">
-                <param argument="--no-self-hits" type="boolean" truevalue="--no-self-hits" falsevalue="" checked="true" label="suppress reporting of identical self hits?" help=""/>
+                <param argument="--no-self-hits" type="boolean" truevalue="--no-self-hits" falsevalue="" checked="true" 
+                    label="Suppress reporting of identical self-hits between sequences" 
+                    help="The FASTA sequence identifiers as well as the sequences of query and target need to be identical for a hit to be deleted"/>
 
                 <param argument="--comp-based-stats" type="select" label="Composition based statistics" help="Compositionally biased sequences often cause false positive matches, which are effectively filtered by this algorithm in a way similar to the composition based statistics used by BLAST">
                     <option value="0">Disable</option>
@@ -169,16 +185,25 @@
                 <param name="reference_database" argument="--db" type="data" format="dmnd" label="Select the reference database" />
             </when>
         </conditional>
-        <expand macro="output_type_macro">
-            <!-- Taxonomy features are not supported for the DAA format (i.e.
-                 can't be used in diamond view) -->
-            <option value="staxids">unique Subject Taxonomy ID(s), separated by a ';' (in numerical order)</option>
-            <option value="sskingdoms">Subject super kingdoms</option>
-            <option value="skingdoms">Subject kingdoms</option>
-            <option value="sphylums">Subject phylums</option>
-        </expand>
+        <conditional name="tax_cond">
+            <param name="tax_select" type="select" label="Restrict search taxonomically?" help="Any taxonomic rank can be used, and only reference sequences matching one of the specified taxon ids will be searched against.">
+                <option value="no" selected="True">No</option>
+                <option value="list">List of taxids entered manually</option>
+                <option value="file">List of taxids from single column tabular file</option>
+            </param>
+            <when value="no"/>
+            <when value="list">
+                <param name="taxonlist" argument="--taxonlist" type="text" value="" label="Comma separated list of taxon ids" help="">
+                    <validator type="regex" message="Taxonlist needs to be a comma separated list of integers">[0-9,]*</validator>
+                </param>
+            </when>
+            <when value="file">
+                <param name="taxonlistfile" argument="--taxonlist" type="data" format="tabular" label="Keep alignments within the given percentage range of the top alignment score for a quer" help="" />
+            </when>
+        </conditional>
         <conditional name="sens_cond">
             <param name='sensitivity' type="select" label="Sensitivity Mode" help="Choose one of the sensitivity modes. The default mode is mainly designed for short read alignment, i.e. finding significant matches of >50 bits on 30-40aa fragments. The sensitive mode is a lot more sensitive than the default and generally recommended for aligning longer sequences. The more sensitive mode provides even more sensitivity. More sensitivity may increase computation time.">
+                <option value="--fast">Fast (--fast)</option>
                 <option value="" selected="True">Default</option>
                 <option value="--mid-sensitive">Mid Sensitive (--mid-sensitive)</option>
                 <option value="--sensitive">Sensitive (--sensitive)</option>
@@ -186,6 +211,9 @@
                 <option value="--very-sensitive">Very Sensitive (--very-sensitive)</option>
                 <option value="--ultra-sensitive">Ultra Sensitive (--ultra-sensitive)</option>
             </param>
+            <when value="--fast">
+                <expand macro="block_size_low_sens"/>
+            </when>
             <when value="">
                 <expand macro="block_size_low_sens"/>
             </when>
@@ -215,26 +243,14 @@
             <option value="PAM70">PAM70 ((6-8)/2; (9-11)/1) [10/1]</option>
             <option value="PAM30">PAM30 ((5-7)/2; (8-10)/1) [9/1]</option>
         </param>
-        <param argument="--gapopen" type="integer" optional="True" value="" label="Gap open penalty" help="leave empty for default (see scoring matrix)" />
-        <param argument="--gapextend" type="integer" optional="True" value="" label="Gap extension penalty" help="leave empty for default (see scoring matrix)" />
-        <param argument="--masking" type="boolean" truevalue="1" falsevalue="0" checked="true" label="enable masking of low complexity regions?" help="Masked residues appear in the output as X"/>
-        <conditional name="tax_cond">
-            <param name="tax_select" type="select" label="Restrict search taxonomically?" help="Any taxonomic rank can be used, and only reference sequences matching one of the specified taxon ids will be searched against">
-                <option value="no" selected="True">No</option>
-                <option value="list">list of taxids entered manually</option>
-                <option value="file">list of taxids from single column tabular file</option>
-            </param>
-            <when value="no"/>
-            <when value="list">
-                <param name="taxonlist" argument="--taxonlist" type="text" value="" label="comma separated list of taxon ids" help="">
-                    <validator type="regex" message="taxonlist needs to be a comma separated list of integers">[0-9,]*</validator>
-                    <validator type="expression" message="taxon ids 0 and 1 are not allowed">not ("0" in value.split(",") or "1" in value.split(","))</validator>
-                </param>
-            </when>
-            <when value="file">
-                <param name="taxonlistfile" argument="--taxonlist" type="data" format="tabular" label="Keep alignments within the given percentage range of the top alignment score for a quer" help="" />
-            </when>
-        </conditional>
+        <param argument="--gapopen" type="integer" optional="True" value="" label="Gap open penalty" help="Leave empty for default (see scoring matrix)" />
+        <param argument="--gapextend" type="integer" optional="True" value="" label="Gap extension penalty" help="Leave empty for default (see scoring matrix)" />
+        <param argument="--masking" type="select" label="Masking algorithm" help="DIAMOND by default applies the tantan repeat masking algorithm to the query and target sequences as described in (Frith, 2011). 
+            This masking procedure increases the specificity of alignments and serves to filter out spurious hits. Note that when using --comp-based-stats (2,3,4), tantan masking is disabled by default.">
+            <option value="0">Disabled</option>
+            <option value="1" selected="true">Tantan</option>
+            <option value="seg">SEG</option>
+        </param>
         <conditional name="filter_score">
             <param name="filter_score_select" type="select" label="Method to filter?" help="(--evalue/--min-score)">
                 <option value="evalue" selected="True">Maximum e-value to report alignments</option>
@@ -247,25 +263,60 @@
                 <param name="min_score" argument="--min-score" type="integer" value="0" label="Minimum bit score to keep an alignment" help="(--min-score)" />
             </when>
         </conditional>
+        <param argument="--iterate" type="boolean" truevalue="--iterate" falsevalue="" checked="false" 
+            label="Run multiple rounds of searches with increasing sensitivity" help="he query dataset will first be searched at a lower sensitivity setting, only searching those query sequences at 
+                the target sensitivity that fail to produce a significant alignment at a lower sensitivity." />
+        <param argument="--algo" type="select" label="Algorithm for seed search" help="Double-indexed is the main algorithm of the program, designed for large input files but less efficient for small 
+            query files. Query-indexed and improves performance for small query files. This mode will be automatically triggered based on the input. Contiguous-seed mode and further improves performance
+            for small query files. The modes differ slightly in their sensitivity, so results are not guaranteed to be 100% identical for different settings of this option.">
+            <option value="0">Doble-indexed (0)</option>
+            <option value="1">Query-indexed (1)</option>
+            <option value="ctg">Contiguous-seed mode (ctg)</option>
+        </param>
         <expand macro="hit_filter_macro" />
-        <param argument="--id" type="integer" value="0" label="Minimum identity percentage to report an alignment" help="" />
-        <param argument="--query-cover" type="integer" value="0" label="Minimum query cover percentage to report an alignment" help="" />
-        <param argument="--subject-cover" type="integer" value="0" label="Minimum subject cover percentage to report an alignment" help="" />
-        <param name="output_unal" type="select" optional="true" multiple="true" label="Output aligned/unaligned queries to separate file" help="">
-            <option value="--un">Output unaligned queries (--un)</option>
-            <option value="--al">Output alaligned queries (--al)</option>
-        </param>
+        <param argument="--global-ranking" type="integer" min="0" value="" optional="true" 
+            label="Limit on the number of Smith Waterman extensions" help="Target sequences will be ranked according to their ungapped extension scores at seed hits, and gapped extensions will only 
+                be computed for the best N targets for each query. Note that this option increases memory use." />
+        <param argument="--id" type="integer" value="0" label="Minimum identity percentage to report an alignment" help="Report only alignments above the given percentage of sequence identity" />
+        <param argument="--query-cover" type="integer" value="0" label="Minimum query cover percentage to report an alignment" help="Report only alignments above the given percentage of query cover" />
+        <param argument="--subject-cover" type="integer" value="0" label="Minimum subject cover percentage to report an alignment" help="Report only alignments above the given percentage of subject cover"/>
+        <section name="output_section" title="Output options">
+            <param argument="--max-hsps" type="integer" min="0" optional="true" label="Maximum number of HSPs" 
+                help="The maximum number of HSPs (High-Scoring Segment Pairs) per target sequence to report for each query. The default policy is to report only the highest-scoring 
+                    HSP for each target, while disregarding alternative, lower-scoring HSPs that are contained in the same target." />
+            <expand macro="output_type_macro">
+                <!-- Taxonomy features are not supported for the DAA format (i.e.
+                        can't be used in diamond view) -->
+                <option value="staxids">unique Subject Taxonomy ID(s), separated by a ';' (in numerical order)</option>
+                <option value="sskingdoms">Subject super kingdoms</option>
+                <option value="skingdoms">Subject kingdoms</option>
+                <option value="sphylums">Subject phylums</option>
+            </expand>
+            <param name="output_unal" type="select" optional="true" multiple="true" label="Output aligned/unaligned queries to separate file" help="">
+                <option value="--un">Output unaligned queries (--un)</option>
+                <option value="--al">Output alaligned queries (--al)</option>
+            </param>
+        </section>
+        <section name="advanced_section" title="Advanced options" expanded="false">
+            <param argument="--seed-cut" type="float" min="0" optional="true" label="Set a complexity cutoff for indexed seeds"/>
+            <param argument="--freq-masking" type="boolean" truevalue="--freq-masking" falsevalue="" checked="false" label="Enable masking seeds based on frequency" help="This option is incompatible with --sed-cut" />
+            <param argument="--motif-masking" type="select" label="Softmask abundant motifs" help="Enable or disable motif masking">
+                <option value="0">Disabled</option>
+                <option value="1">Enabled</option>
+            </param>
+        </section>
     </inputs>
     <outputs>
         <expand macro="output_macro" />
         <data format_source="query" name="unalqueries" label="${tool.name} on ${on_string}: unaligned queries">
-            <filter>output_unal and "--un" in output_unal</filter>
+            <filter>output_section['output_unal'] and "--un" in output_section['output_unal']</filter>
         </data>
         <data format_source="query" name="alqueries" label="${tool.name} on ${on_string}: aligned queries">
-            <filter>output_unal and "--un" in output_unal</filter>
+            <filter>output_section['output_unal'] and "--un" in output_section['output_unal']</filter>
         </data>
     </outputs>
     <tests>
+        <!--Test 01-->
         <test expect_num_outputs="3">
             <conditional name="method_cond">
                 <param name="method_select" value="blastp" />
@@ -275,12 +326,14 @@
                 <param name="db_source" value="history"/>
                 <param name="reference_database" value="db-wtax.dmnd"/>
             </conditional>
-            <conditional name="output">
-                <param name="outfmt" value="6"/>
-                <!-- removed ,cigar from test: https://github.com/bbuchfink/diamond/issues/532 -->
-                <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore,scovhsp,sskingdoms,skingdoms,sphylums"/>
-                <param name="unal" value="true"/>
-            </conditional>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <!-- removed ,cigar from test: https://github.com/bbuchfink/diamond/issues/532 -->
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore,scovhsp,sskingdoms,skingdoms,sphylums"/>
+                    <param name="unal" value="true"/>
+                </conditional>
+            </section>
             <conditional name="sens_cond">
                 <param name="sensitivity" value=""/>
             </conditional>
@@ -313,6 +366,7 @@
             </output>
             <output name="blast_tabular" file="diamond_results.tabular"/>
         </test>
+        <!--Test 02-->
         <test expect_num_outputs="1">
             <conditional name="method_cond">
                 <param name="method_select" value="blastp" />
@@ -326,10 +380,12 @@
                 <param name="tax_select" value="list"/>
                 <param name="taxonlist" value="2" />
             </conditional>
-            <conditional name="output">
-                <param name="outfmt" value="6"/>
-                <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
-            </conditional>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
             <conditional name="sens_cond">
                 <param name="sensitivity" value=""/>
             </conditional>
@@ -351,6 +407,7 @@
             </conditional>
             <output name="blast_tabular" file="diamond_results.wtax.tabular"/>
         </test>
+        <!--Test 03-->
         <test expect_num_outputs="1">
             <conditional name="method_cond">
                 <param name="method_select" value="blastx" />
@@ -363,9 +420,69 @@
                 <param name="db_source" value="history"/>
                 <param name="reference_database" value="db.dmnd"/>
             </conditional>
-            <conditional name="output">
-                <param name="outfmt" value="0"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="0"/>
+                </conditional>
+            </section>
+            <conditional name="sens_cond">
+                <param name="sensitivity" value=""/>
+            </conditional>
+            <param name="matrix" value="BLOSUM62"/>
+            <param name="comp_based_stats" value="1"/>
+            <param name="masking" value="1"/>
+            <conditional name="hit_filter">
+                <param name="hit_filter_select" value="top"/>
+                <param name="top" value="10" />
+            </conditional>
+            <conditional name="filter_score">
+                <param name="filter_score_select" value="min-score"/>
+                <param name="min_score" value="1" />
+            </conditional>
+            <param name="id" value="0"/>
+            <param name="query_cover" value="0"/>
+            <conditional name="sens_cond">
+                <param name="block_size" value="2"/>
+            </conditional>
+            <output name="blast_tabular" file="diamond_results.pairwise"/>
+        </test>
+        <!--Test 04-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastp" />
             </conditional>
+            <param name="query" value="protein.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="history"/>
+                <param name="reference_database" value="db-wtax.dmnd"/>
+            </conditional>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="100"/>
+                    <param name="salltitles" value="false"/>
+                    <param name="sallseqid" value="false"/>
+                </conditional>
+            </section>
+            <output name="daa_output" file="diamond_results.daa" compare="sim_size" delta="10"/>
+        </test>
+        <!--Test 05-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastx" />
+                <conditional name="frameshift_cond">
+                    <param name="frameshift_select" value="yes"/>
+                </conditional>
+            </conditional>
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
+            </conditional>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="0"/>
+                </conditional>
+            </section>
             <conditional name="sens_cond">
                 <param name="sensitivity" value=""/>
             </conditional>
@@ -387,57 +504,140 @@
             </conditional>
             <output name="blast_tabular" file="diamond_results.pairwise"/>
         </test>
+        <!-- Test 06 iterate option-->
         <test expect_num_outputs="1">
             <conditional name="method_cond">
-                <param name="method_select" value="blastp" />
+                <param name="method_select" value="blastx" />
             </conditional>
-            <param name="query" value="protein.fasta" ftype="fasta"/>
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
             <conditional name="ref_db_source">
-                <param name="db_source" value="history"/>
-                <param name="reference_database" value="db-wtax.dmnd"/>
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
             </conditional>
-            <conditional name="output">
-                <param name="outfmt" value="100"/>
-                <param name="salltitles" value="false"/>
-                <param name="sallseqid" value="false"/>
-            </conditional>
-            <output name="daa_output" file="diamond_results.daa" compare="sim_size" delta="10"/>
+            <param name="iterate" value="true"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_iterate.tabular"/>
         </test>
+        <!--Test 07 algo option-->
         <test expect_num_outputs="1">
             <conditional name="method_cond">
                 <param name="method_select" value="blastx" />
-                <conditional name="frameshift_cond">
-                    <param name="frameshift_select" value="yes"/>
+            </conditional>
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
+            </conditional>
+            <param name="algo" value="1"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
                 </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_algorithm.tabular"/>
+        </test>
+        <!--Test 08 global-ranking option-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastx" />
+            </conditional>
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
+            </conditional>
+            <param name="global_ranking" value="10"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_global_ranking.tabular"/>
+        </test>
+        <!--Test 09 max-hsps option-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastx" />
             </conditional>
             <param name="query" value="nucleotide.fasta" ftype="fasta"/>
             <conditional name="ref_db_source">
                 <param name="db_source" value="indexed"/>
                 <param name="index" value="testDb"/>
             </conditional>
-            <conditional name="output">
-                <param name="outfmt" value="0"/>
+            <param name="max_hsps" value="10"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_max_hsps.tabular"/>
+        </test>
+        <!--Test 10 seed-cut option-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastx" />
             </conditional>
-            <conditional name="sens_cond">
-                <param name="sensitivity" value=""/>
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
+            </conditional>
+            <param name="seed_cut" value="100"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_seed_cut.tabular"/>
+        </test>
+        <!--Test 11 freq-masking option-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastx" />
             </conditional>
-            <param name="matrix" value="BLOSUM62"/>
-            <param name="comp_based_stats" value="1"/>
-            <param name="masking" value="1"/>
-            <conditional name="hit_filter">
-                <param name="hit_filter_select" value="top"/>
-                <param name="top" value="10" />
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
             </conditional>
-            <conditional name="filter_score">
-                <param name="filter_score_select" value="min-score"/>
-                <param name="min_score" value="1" />
+            <param name="freq_masking" value="true"/>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_freq_masking.tabular"/>
+        </test>
+        <!--Test 12 motif-masking option-->
+        <test expect_num_outputs="1">
+            <conditional name="method_cond">
+                <param name="method_select" value="blastx" />
             </conditional>
-            <param name="id" value="0"/>
-            <param name="query_cover" value="0"/>
-            <conditional name="sens_cond">
-                <param name="block_size" value="2"/>
+            <param name="query" value="nucleotide.fasta" ftype="fasta"/>
+            <conditional name="ref_db_source">
+                <param name="db_source" value="indexed"/>
+                <param name="index" value="testDb"/>
             </conditional>
-            <output name="blast_tabular" file="diamond_results.pairwise"/>
+            <section name="advanced_section">
+                <param name="motif_masking" value="1"/>
+            </section>
+            <section name="output_section">
+                <conditional name="output">
+                    <param name="outfmt" value="6"/>
+                    <param name="fields" value="qseqid,sseqid,pident,length,mismatch,gapopen,qstart,qend,sstart,send,evalue,bitscore"/>
+                </conditional>
+            </section>
+            <output name="blast_tabular" file="diamond_results_motif_masking.tabular"/>
         </test>
     </tests>
     <help>