diff nextalign.xml @ 14:4d93e708218c draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nextclade commit 2c63aeec74c260ccd9b2cdbbe1869bef2f3d5cda
author iuc
date Thu, 04 Aug 2022 06:59:23 +0000
parents 1c785a6b4d90
children aadc5d8099e2
line wrap: on
line diff
--- a/nextalign.xml	Mon Mar 21 22:09:49 2022 +0000
+++ b/nextalign.xml	Thu Aug 04 06:59:23 2022 +0000
@@ -10,10 +10,9 @@
     <version_command>nextalign --version-detailed</version_command>
     <command detect_errors="exit_code"><![CDATA[
 @REF_FASTA@
-ln -s '$sequences' sequences.fasta &&
-nextalign
-    --sequences sequences.fasta
-    --reference reference.fa
+@QUERY_FASTA@
+nextalign run
+    --input-ref reference.fa
     --output-fasta '$output_fasta'
 #if $output_insertions:
     --output-insertions '$output_csv'
@@ -21,23 +20,25 @@
 #if $translation.translation_select == "yes":
     --genes '${translation.genes}'
     --genemap '${translation.genemap}'
+    --output-translations 'sequences.gene.{gene}.fasta'
 #end if
-    --min-length '${min_length}'
-    --penalty-gap-extend '${penalty_gap_extend}'
-    --penalty-gap-open '${penalty_gap_open}'
-    --penalty-gap-open-in-frame '${penalty_gap_open_in_frame}'
-    --penalty-gap-open-out-of-frame '${penalty_gap_open_out_of_frame}'
-    --penalty-mismatch '${penalty_mismatch}'
-    --score-match '${score_match}'
-    --max-indel '${max_indel}'
-    --nuc-seed-length '${nuc_seed_length}'
-    --nuc-min-seeds '${nuc_min_seeds}'
-    --nuc-seed-spacing '${nuc_seed_spacing}'
-    --nuc-mismatches-allowed '${nuc_mismatches_allowed}'
+    --min-length ${min_length}
+    --penalty-gap-extend ${penalty_gap_extend}
+    --penalty-gap-open ${penalty_gap_open}
+    --penalty-gap-open-in-frame ${penalty_gap_open_in_frame}
+    --penalty-gap-open-out-of-frame ${penalty_gap_open_out_of_frame}
+    --penalty-mismatch ${penalty_mismatch}
+    --score-match ${score_match}
+    --max-indel ${max_indel}
+    --seed-length ${seed_length}
+    --min-seeds ${min_seeds}
+    --seed-spacing ${seed_spacing}
+    --mismatches-allowed ${mismatches_allowed}
+    $query
     ]]></command>
     <inputs>
         <expand macro="reference"/>
-        <param argument="--sequences" type="data" format="fasta" label="FASTA file with input sequences"/>
+        <param name="input_fasta" type="data" format="fasta,fasta.gz" label="FASTA file with input sequences"/>
         <param argument="--output-insertions" type="boolean" checked="true" label="Output insertion sequences?" help="Outputs stripped insertions relative to reference as CSV"/>
         <conditional name="translation">
             <param name="translation_select" type="select" label="Translate annotated genes based on GFF and gene list?">
@@ -54,7 +55,7 @@
                         </valid>
                     </sanitizer>
                 </param>
-                <param argument="--genemap" type="data" format="gtf" label="GTF file containing custom gene map"/> <!-- TODO - make sure they need GFF and not GTF or GFF3 -->
+                <param argument="--genemap" type="data" format="gff3" label="GFF3 file containing custom gene map"/>
             </when>
             <when value="no"/>
         </conditional>
@@ -82,15 +83,15 @@
         <param argument="--max-indel" type="integer" value="400" min="0"
             label="Maximum length of insertions or deletions allowed to proceed with alignment."
             help="Alignments with long indels are slow to compute and require substantial memory in the current implementation. Alignment of sequences with indels that are longer than this value will not be attempted and a warning will be emitted." />
-        <param argument="--nuc-seed-length" type="integer" value="21" min="1"
+        <param argument="--seed-length" type="integer" value="21" min="1"
             label="Seed length for nucleotide alignment."
             help="Seeds should be long enough to be unique, but short enough to match with high probability." />
-        <param argument="--nuc-min-seeds" type="integer" value="10" min="1"
+        <param argument="--min-seeds" type="integer" value="10" min="1"
             label="Minimum number of seeds to search for during nucleotide alignment."
             help="Relevant for short sequences. In long sequences, the number of seeds is determined by nucleotide seed spacing. This should be a positive integer." />
-        <param argument="--nuc-seed-spacing" type="integer" value="100" min="0"
+        <param argument="--seed-spacing" type="integer" value="100" min="0"
             label="Spacing between seeds during nucleotide alignment." />
-        <param argument="--nuc-mismatches-allowed" type="integer" value="3" min="0"
+        <param argument="--mismatches-allowed" type="integer" value="3" min="0"
             label="Maximum number of mismatching nucleotides."
             help="Maximum number of mismatching nucleotides allowed for a seed to be considered a match." />
     </inputs>
@@ -114,20 +115,18 @@
                 <param name="ref_file" value="reference.fasta"/>
             </conditional>
             <param name="output_insertions" value="true"/>
-            <param name="sequences" value="subsampled.fasta"/>
+            <param name="input_fasta" value="subsampled.fasta" ftype="fasta" />
             <conditional name="translation">
                 <param name="translation_select" value="yes"/>
-                <param name="genes" value="E,M,N,ORF10,ORF14,ORF1a,ORF1b,ORF3a,ORF6,ORF7a,ORF7b,ORF8,ORF9b,S"/>
-                <param name="genemap" value="genemap.gtf" ftype="gtf"/>
+                <param name="genes" value="E,M,N,ORF1a,ORF1b,ORF3a,ORF6,ORF7a,ORF7b,ORF8,ORF9b,ORF10,S"/>
+                <param name="genemap" value="genemap.gff3" ftype="gff3"/>
             </conditional>
             <output name="output_fasta" file="output.fasta" sort="true"/>
             <output name="output_csv" file="insertions.csv" sort="true"/>
-            <output_collection name="translations" type="list" count="14">
+            <output_collection name="translations" type="list" count="13">
                 <element name="E" file="subsampled.gene.E.fasta" ftype="fasta" sort="true"/>
                 <element name="M" file="subsampled.gene.M.fasta" ftype="fasta" sort="true"/>
                 <element name="N" file="subsampled.gene.N.fasta" ftype="fasta" sort="true"/>
-                <element name="ORF10" file="subsampled.gene.ORF10.fasta" ftype="fasta" sort="true"/>
-                <element name="ORF14" file="subsampled.gene.ORF14.fasta" ftype="fasta" sort="true"/>
                 <element name="ORF1a" file="subsampled.gene.ORF1a.fasta" ftype="fasta" sort="true"/>
                 <element name="ORF1b" file="subsampled.gene.ORF1b.fasta" ftype="fasta" sort="true"/>
                 <element name="ORF3a" file="subsampled.gene.ORF3a.fasta" ftype="fasta" sort="true"/>
@@ -148,7 +147,7 @@
                 <param name="ref_file" value="reference.fasta"/>
             </conditional>
             <param name="output_insertions" value="false"/>
-            <param name="sequences" value="subsampled.fasta"/>
+            <param name="input_fasta" value="subsampled.fasta" ftype="fasta" />
             <conditional name="translation">
                 <param name="translation_select" value="no"/>
             </conditional>
@@ -163,7 +162,7 @@
                 <param name="ref_file" value="reference_fasta"/>
             </conditional>
             <param name="output_insertions" value="false"/>
-            <param name="sequences" value="subsampled.fasta"/>
+            <param name="input_fasta" value="subsampled.fasta" ftype="fasta"/>
             <conditional name="translation">
                 <param name="translation_select" value="no"/>
             </conditional>