Mercurial > repos > iuc > virulign

comparison virulign.xml @ 0:0620360e441e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/main/tools/virulign/ commit b18f2f0794274e906cf75c8ed07e1a7eb1c9b7f8
author iuc
date Thu, 14 Aug 2025 17:13:56 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:0620360e441e
1 <tool id="virulign" name="VIRULIGN" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="23.2">
2 <description>Codon-correct pairwise alignments.</description>
3 <macros>
4 <token name="@TOOL_VERSION@">1.1.1</token>
5 <token name="@VERSION_SUFFIX@">0</token>
6 </macros>
7 <creator>
8 <person givenName="Lieven" familyName="Sterck" url="https://github.com/lsterck" />
9 <organization name="Sciensano-BioIT" url="https://github.com/BioinformaticsPlatformWIV-ISP" />
10 </creator>
11 <requirements>
12 <requirement type="package" version="@TOOL_VERSION@">virulign</requirement>
13 </requirements>
14 <version_command>echo @TOOL_VERSION@</version_command>
15 <command detect_errors="exit_code"><![CDATA[
16 ##link named files
17 ln -sf '$reference' reference.$reference.ext &&
18 ln -sf '$sequences' sequences.$sequences.ext &&
19
20 #if $nt_debug == 'yes'
21 mkdir Failed/ &&
22 #end if
23
24 ## run virulign
25 virulign
26 reference.$reference.ext
27 sequences.$sequences.ext
28 --exportKind $exportKind
29 --exportAlphabet $exportAlphabet
30 --exportWithInsertions $exportWithInsertions
31 --exportReferenceSequence $exportReferenceSequence
32 --gapExtensionPenalty $gapExtensionPenalty
33 --gapOpenPenalty $gapOpenPenalty
34 --maxFrameShifts $maxFrameShifts
35
36 #if $nt_debug == 'yes'
37 --nt-debug Failed
38 #end if
39
40 > '$output'
41
42 ]]></command>
43 <inputs>
44 <param name="reference" type="data" format="fasta,xml" label="Reference (FASTA or ORF‑description XML)" />
45 <param name="sequences" type="data" format="fasta" label="FASTA sequences to align" />
46
47 <param argument="--exportKind" type="select" label="Which type of output to generate:" help="Defines the output type of the alignment, either a FASTA alignment file or a CSV mutation file">
48 <option value="GlobalAlignment">GlobalAlignment</option>
49 <option value="PairwiseAlignments">PairwiseAlignments</option>
50 <option value="PositionTable">PositionTable</option>
51 <option value="MutationTable">MutationTable</option>
52 <option value="Mutations">Mutations</option>
53 </param>
54
55 <param argument="--exportAlphabet" type="select" label="Sequence alphabet" help="Defines the alphabet in which the alignment is generated.">
56 <option value="AminoAcids">Amino acids</option>
57 <option value="Nucleotides">Nucleotides</option>
58 </param>
59
60 <param argument="--exportWithInsertions" type="boolean" truevalue="yes" falsevalue="no" label="Include insertions?" checked="true" help="Determines whether insertions can be added to the reference sequence. The option 'yes' will insert gaps into the reference sequence to accommodate the identification of codon insertions in the target sequences. Option 'No' removes codon insertions in the reference sequence that were generated during the alignment procedure."/>
61
62 <param argument="--exportReferenceSequence" type="boolean" truevalue="yes" falsevalue="no" label="Add the reference sequence to output?" checked="false" help="Controls whether the reference sequence is to be added to the alignment."/>
63
64 <param argument="--gapOpenPenalty" type="float" value="10.0" label="Gap open penalty" help="Defines the value of the penalty to start a new gap."/>
65 <param argument="--gapExtensionPenalty" type="float" value="3.3" min="0" label="Gap extension penalty" help="Defines the value of the penalty to extend an existing gap."/>
66 <param argument="--maxFrameShifts" type="integer" value="3" min="0" label="Maximum allowed frame shifts" help="Defines the maximum number of frame-shifts allowed."/>
67 <param argument="--nt-debug" type="boolean" truevalue="yes" falsevalue="no" label="Report failed sequences?" checked="false" help="This parameter allows to visualise sequences that could not be aligned by VIRULIGN. When this option is used, pairwise sequence alignments of failed target sequences and the reference sequence will be reported. This feature allows to inspect each failed target sequence individually to understand why the target sequence did not pass the quality control of the alignment."/>
68 </inputs>
69
70 <outputs>
71 <data name="output" format="txt" label="${tool.name} on ${sequences.name} ($on_string): Output">
72 <change_format>
73 <when input="exportKind" value="GlobalAlignment" format="fasta" />
74 <when input="exportKind" value="Mutations" format="tsv" />
75 <when input="exportKind" value="MutationTable" format="csv" />
76 <when input="exportKind" value="PositionTable" format="csv" />
77 <when input="exportKind" value="PairwiseAlignments" format="fasta" />
78 </change_format>
79 </data>
80 <collection name="failed_outputs" type="list" label="${tool.name} on ${sequences.name}: Failed alignments">
81 <filter>nt_debug is True</filter>
82 <discover_datasets pattern="(?P&lt;designation&gt;.+)\.fasta" directory="Failed" format="fasta" visible="false"/>
83 </collection>
84 </outputs>
85 <tests>
86 <test expect_num_outputs="1">
87 <param name="reference" value="denv-NC_001477.fasta" ftype="fasta"/>
88 <param name="sequences" value="denv-1_small.fasta" ftype="fasta"/>
89 <param name="exportKind" value="GlobalAlignment" />
90 <param name="exportAlphabet" value="Nucleotides" />
91 <output name="output" file="denv-output.fasta" ftype="fasta"/>
92 </test>
93 <test expect_num_outputs="1">
94 <param name="reference" value="ZIKV-rega.xml" ftype="xml"/>
95 <param name="sequences" value="zikv-small.fasta" ftype="fasta"/>
96 <param name="exportKind" value="PositionTable" />
97 <param name="exportReferenceSequence" value="yes" />
98 <param name="exportAlphabet" value="AminoAcids" />
99 <output name="output" file="zikv-output.csv" ftype="csv"/>
100 </test>
101 <test expect_num_outputs="1">
102 <param name="reference" value="HIV-HXB2-gag.xml" ftype="xml"/>
103 <param name="sequences" value="HIV_small.fasta" ftype="fasta"/>
104 <param name="exportKind" value="GlobalAlignment" />
105 <param name="exportAlphabet" value="Nucleotides" />
106 <param name="exportReferenceSequence" value="yes" />
107 <param name="exportWithInsertions" value="no" />
108 <output name="output" file="HIV-output.fasta" ftype="fasta"/>
109 </test>
110 </tests>
111 <help><![CDATA[
112 **VIRULIGN**
113
114 A tool for codon-correct pairwise alignments, with an augmented functionality to annotate the alignment according the positions of the proteins. We present a `tutorial`_ demonstrating VIRULIGN's functionality for three different pathogens (i.e., HIV-1, Dengue virus, Zika virus).
115
116 VIRULIGN minimally requires a FASTA file with target sequences and a reference sequence in order to generate a codon-correct alignment in a predefined output format (see below). The reference sequence can be either provided to VIRULIGN in FASTA format or embedded in an XML file. In this XML file, also an annotation of the different proteins, regions or other structures can be given by the positions relative to the reference genome.
117
118 OUTPUT:
119
120 - The option GlobalAlignment will generate a FASTA file of the target sequences aligned against a single reference sequence and formatted as a MSA.
121 - The option PairwiseAlignments will generate a FASTA file of the target sequences, with each sequence aligned separately against the reference sequence.
122 - The option PositionTable will create a comma-separated value (CSV) file where each position of the alignment is given as a separate column. The CSV file is annotated according to the numerical positions in the protein.
123 - The option MutationTable will create a CSV file where each mutation present at a specific position is given as a separate column in Boolean representation. The CSV file is annotated according to the numerical position in the protein
124 - the option Mutations will output, for each sequence, a list of amino acids changes compared to the reference sequence
125
126
127 **More Information**
128
129 - **Official Repository**: https://github.com/rega-cev/virulign
130 - **User Guide & Tutorial**: https://github.com/rega-cev/virulign-tutorial
131
132 .. _tutorial: https://github.com/rega-cev/virulign-tutorial
133
134 ]]></help>
135 <citations>
136 <citation type="doi">10.1093/bioinformatics/bty851</citation>
137 </citations>
138 </tool>