comparison agat.xml @ 2:e009d8260be2 draft

planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/agat commit 8eb20f601bc1d2a50c8877b7d0ade057e8f86eae
author bgruening
date Thu, 07 Sep 2023 05:29:24 +0000
parents f7c0a0030254
children a6318f87f2cd
comparison
equal deleted inserted replaced
1:6d8444408ff1 2:e009d8260be2
7 <expand macro="requirements" /> 7 <expand macro="requirements" />
8 <version_command>agat_sq_stat_basic.pl --version</version_command> 8 <version_command>agat_sq_stat_basic.pl --version</version_command>
9 <command detect_errors="exit_code"><![CDATA[ 9 <command detect_errors="exit_code"><![CDATA[
10 #if $tool.selector == 'fix' 10 #if $tool.selector == 'fix'
11 @input_annotation_single@ 11 @input_annotation_single@
12 agat_convert_sp_gxf2gxf.pl -gff $input_annotation --output 'output.gff' && 12 agat_convert_sp_gxf2gxf.pl
13 cat 'output.gff' > '${annotation_gff}' 13 --gxf $input_annotation
14 --config $agat_configfile
15 --output 'output' &&
16 cat 'output' > '${annotation}'
14 #else if $tool.selector == 'convert_GFF2GTF' 17 #else if $tool.selector == 'convert_GFF2GTF'
15 @input_annotation_single@ 18 @input_annotation_single@
16 agat_convert_sp_gff2gtf.pl --gff $input_annotation --gtf_version $tool.gtf_version --output 'output.gtf' && 19 agat_convert_sp_gff2gtf.pl
20 --gff $input_annotation
21 --gtf_version $tool.gtf_version
22 --output 'output.gtf' &&
17 cat 'output.gtf' > '${annotation_gtf}' 23 cat 'output.gtf' > '${annotation_gtf}'
18 #else if $tool.selector == 'convert_GTF2GFF' 24 #else if $tool.selector == 'convert_GTF2GFF'
19 @input_annotation_single@ 25 @input_annotation_single@
20 agat_convert_sp_gxf2gxf.pl --gff $input_annotation --output 'output.gff' && 26 agat_convert_sp_gxf2gxf.pl
27 --gff $input_annotation
28 --output 'output.gff' &&
21 cat 'output.gff' > '${annotation_gff}' 29 cat 'output.gff' > '${annotation_gff}'
22 #else if $tool.selector == 'compare' 30 #else if $tool.selector == 'compare'
23 @input_annotation_double@ 31 @input_annotation_double@
24 agat_sp_compare_two_annotations.pl --gff1 $input1 --gff2 $input2 --output 'temp_output' && 32 agat_sp_compare_two_annotations.pl
33 --gff1 $input1
34 --gff2 $input2
35 --output 'temp_output' &&
25 cat 'temp_output' > '${stats_output}' 36 cat 'temp_output' > '${stats_output}'
26 #else if $tool.selector == 'extract' 37 #else if $tool.selector == 'extract'
27 @input_annotation_single@ 38 @input_annotation_single@
28 @input_reference@ 39 @input_reference@
29 agat_sp_extract_sequences.pl 40 agat_sp_extract_sequences.pl
54 --output '${sequence_output}' 65 --output '${sequence_output}'
55 #else if $tool.selector == 'functional_analysis' 66 #else if $tool.selector == 'functional_analysis'
56 @input_annotation_single@ 67 @input_annotation_single@
57 @input_reference@ 68 @input_reference@
58 mkdir -p './statistics' && 69 mkdir -p './statistics' &&
59 agat_sp_statistics.pl 70 agat_sp_functional_statistics.pl
60 --gff $input_annotation 71 --gff $input_annotation
61 --gs $ref_genome 72 --gs $ref_genome
62 --output 'temp_output' && 73 --output 'temp_output' &&
63 cat 'temp_output' > '$stats_output' 74 cat 'temp_output/gene@transcript/table_per_feature_type.txt' > '$stats_output'
75
64 #else if $tool.selector == 'merge_annotations' 76 #else if $tool.selector == 'merge_annotations'
65 @input_annotation_double@ 77 @input_annotation_double@
66 agat_sp_merge_annotations.pl -gff $input1 --gff $input2 --output 'temp_output' && 78 agat_sp_merge_annotations.pl
67 cat 'temp_output' > '${annotation_gff}' 79 --gff $input1
80 --gff $input2
81 --config $agat_configfile
82 --output 'output' &&
83 cat 'output' > '${annotation}'
68 #else if $tool.selector == 'annotation_statistics' 84 #else if $tool.selector == 'annotation_statistics'
69 @input_annotation_single@ 85 @input_annotation_single@
70 @input_reference@ 86 @input_reference@
71 agat_sp_statistics.pl --gff $input_annotation --gs $ref_genome -d --output 'temp_output' && 87 agat_sp_statistics.pl
88 --gff $input_annotation
89 --gs $ref_genome
90 -d
91 --output 'temp_output' &&
72 cat 'temp_output' > '$stats_output' 92 cat 'temp_output' > '$stats_output'
73 #else if $tool.selector == 'filter_feature_fasta' 93 #else if $tool.selector == 'filter_feature_fasta'
74 @input_annotation_single@ 94 @input_annotation_single@
75 @input_reference@ 95 @input_reference@
76 agat_sq_filter_feature_from_fasta.pl --gff $input_annotation --fasta $ref_genome --output 'temp_output' && 96 agat_sq_filter_feature_from_fasta.pl
77 cat 'temp_output' > '${features_filtered}' 97 --gff $input_annotation
98 --fasta $ref_genome
99 --config $agat_configfile
100 --output 'output' &&
101 cat 'output' > '${annotation}'
78 #else if $tool.selector == 'complement' 102 #else if $tool.selector == 'complement'
79 @input_annotation_double@ 103 @input_annotation_double@
80 agat_sp_complement_annotations.pl --ref $input1 --add $input2 --size_min $tool.size_min --output 'temp_output' && 104 agat_sp_complement_annotations.pl
81 cat 'temp_output' > '${annotation_gff}' 105 --ref $input1
106 --add $input2
107 --size_min $tool.size_min
108 --config $agat_configfile
109 --output 'temp_output' &&
110 cat 'temp_output' > '${annotation}'
111 #else if $tool.selector == 'splice_sites'
112 @input_annotation_single@
113 agat_sp_add_splice_sites.pl
114 --gff $input_annotation
115 --config $agat_configfile
116 --output 'output' &&
117 cat 'output' > '${annotation}'
82 #end if 118 #end if
83 ]]> 119 ]]>
84 </command> 120 </command>
121 <configfiles>
122 <configfile name="agat_configfile"><![CDATA[
123 #if $tool.selector in ['fix','merge_annotations','complement','splice_sites','filter_feature_fasta']
124 ---
125 output_format: $tool.output_format.selector
126 #if $tool.output_format.selector == "GFF"
127 gff_output_version: $tool.output_format.version
128 gtf_output_version: relax
129 #else
130 gff_output_version: 3
131 gtf_output_version: $tool.output_format.version
132 #end if
133 verbose: 1
134 progress_bar: true
135 log: true
136 debug: false
137 tabix: false
138 merge_loci: $tool.merge_loci
139 throw_fasta: false
140 force_gff_input_version: 0
141 create_l3_for_l2_orphan: $tool.create_exon
142 locus_tag:
143 - locus_tag
144 - gene_id
145 prefix_new_id: nbis
146 check_sequential: true
147 check_l2_linked_to_l3: true
148 check_l1_linked_to_l2: true
149 remove_orphan_l1: true
150 check_all_level3_locations: true
151 check_cds: true
152 check_exons: true
153 check_utrs: true
154 check_all_level2_locations: true
155 check_all_level1_locations: true
156 check_identical_isoforms: true
157 #end if
158 ]]></configfile>
159 </configfiles>
85 <inputs> 160 <inputs>
86 <conditional name="tool"> 161 <conditional name="tool">
87 <param name="selector" type="select" label="AGAT tool selector" help="As AGAT is a toolkit, it contains a lot of tools. If any of them is missing, please contact the server admin."> 162 <param name="selector" type="select" label="AGAT tool selector" help="As AGAT is a toolkit, it contains a lot of tools. If any of them is missing, please contact the server admin.">
163 <option value="splice_sites">Add splice sites</option>
88 <option value="annotation_statistics">Annotation statistics (agat_sp_statistics.pl)</option> 164 <option value="annotation_statistics">Annotation statistics (agat_sp_statistics.pl)</option>
89 <option value="compare">Compare annotation files (agat_sp_compare_two_annotations.pl)</option> 165 <option value="compare">Compare annotation files (agat_sp_compare_two_annotations.pl)</option>
90 <option value="complement">Complement annotation file (agat_sp_complement_annotations.pl)</option> 166 <option value="complement">Complement annotation file (agat_sp_complement_annotations.pl)</option>
91 <option value="extract">Extract sequences (agat_sp_extract_sequences.pl)</option> 167 <option value="extract">Extract sequences (agat_sp_extract_sequences.pl)</option>
92 <option value="convert_GFF2GTF">GFF to GTF format conversion (agat_convert_sp_gff2gtf.pl)</option> 168 <option value="convert_GFF2GTF">GFF to GTF format conversion (agat_convert_sp_gff2gtf.pl)</option>
111 <option value="gene">Gene</option> 187 <option value="gene">Gene</option>
112 <option value="transcript">Transcript</option> 188 <option value="transcript">Transcript</option>
113 <option value="exon">Exon</option> 189 <option value="exon">Exon</option>
114 <option value="cds">CDS</option> 190 <option value="cds">CDS</option>
115 <option value="trna">tRNA</option> 191 <option value="trna">tRNA</option>
116 <option value="three_prime_utr">3' UTR</option> 192 <option value="three_prime_utr">3 UTR</option>
117 <option value="five_prime_utr">5' UTR</option> 193 <option value="five_prime_utr">5 UTR</option>
118 </param> 194 </param>
119 <param argument="--mrna" type="boolean" truevalue="--mrna" falsevalue="" checked="false" label="Extract mRNA sequences" help=" This extract the mrna 195 <param argument="--mrna" type="boolean" truevalue="--mrna" falsevalue="" checked="false" label="Extract mRNA sequences" help=" This extract the mrna
120 sequence (i.e transcribed sequence (devoid of introns, but containing untranslated exons))." /> 196 sequence (i.e transcribed sequence (devoid of introns, but containing untranslated exons))." />
121 <param argument="--cdna" type="boolean" truevalue="--cdna" falsevalue="" checked="false" label="Extract the cDNA sequence" 197 <param argument="--cdna" type="boolean" truevalue="--cdna" falsevalue="" checked="false" label="Extract the cDNA sequence"
122 help=" This extract the cdna sequence (i.e reverse complement of the mRNA: transcribed sequence (devoid of introns, but 198 help=" This extract the cdna sequence (i.e reverse complement of the mRNA: transcribed sequence (devoid of introns, but
125 help=" This option allows removing the translation of the final stop codons that is represented by the '*' character. This character can be 201 help=" This option allows removing the translation of the final stop codons that is represented by the '*' character. This character can be
126 disturbing for many programs (e.g interproscan)" /> 202 disturbing for many programs (e.g interproscan)" />
127 <param argument="--clean_internal_stop" type="boolean" truevalue="--clean_internal_stop" falsevalue="" checked="false" label="Clean internal 203 <param argument="--clean_internal_stop" type="boolean" truevalue="--clean_internal_stop" falsevalue="" checked="false" label="Clean internal
128 stop codons" help="The Clean Internal Stop option allows replacing the translation of the stop codons present among the sequence that is 204 stop codons" help="The Clean Internal Stop option allows replacing the translation of the stop codons present among the sequence that is
129 represented by the '*' character by . This character can be disturbing for many programs (e.g interproscan)" /> 205 represented by the '*' character by . This character can be disturbing for many programs (e.g interproscan)" />
130 <param argument="--upstream" type="integer" min="0" value="" optional="true" label="Upstream nucleotides" help="It will take that number of nucleotide in more at the 5' extremity." /> 206 <param argument="--upstream" type="integer" min="0" value="" optional="true" label="Upstream nucleotides" help="It will take that number of nucleotide in more at the 5 extremity." />
131 <param argument="--downstream" type="integer" min="0" value="" optional="true" label="Downstream nucleotides" help="It will take that number of downstream nucleotides." /> 207 <param argument="--downstream" type="integer" min="0" value="" optional="true" label="Downstream nucleotides" help="It will take that number of downstream nucleotides." />
132 <param argument="--full" type="boolean" truevalue="--full" falsevalue="" checked="false" label="Full" help="This option allows dealing 208 <param argument="--full" type="boolean" truevalue="--full" falsevalue="" checked="false" label="Full" help="This option allows dealing
133 with feature that may span over several locations like CDS or exon, in order to extract the full sequence from the start extremity 209 with feature that may span over several locations like CDS or exon, in order to extract the full sequence from the start extremity
134 of the first chunck to the end extremity of the last chunk. The use of that option with '--type exon' will extract the pre-mRNA 210 of the first chunck to the end extremity of the last chunk. The use of that option with '--type exon' will extract the pre-mRNA
135 sequence (i.e with introns). Use of that option on CDS will give the pre-mRNA without the untraslated regions (UTRs). " /> 211 sequence (i.e with introns). Use of that option on CDS will give the pre-mRNA without the untraslated regions (UTRs). " />
169 <expand macro="ANNOTATION_INPUT" format="gtf"/> 245 <expand macro="ANNOTATION_INPUT" format="gtf"/>
170 </when> 246 </when>
171 <when value="filter_feature_fasta"> 247 <when value="filter_feature_fasta">
172 <expand macro="ANNOTATION_INPUT" /> 248 <expand macro="ANNOTATION_INPUT" />
173 <expand macro="REFERENCE_FASTA"/> 249 <expand macro="REFERENCE_FASTA"/>
250 <expand macro="AGAT_CONFIG"/>
174 </when> 251 </when>
175 <when value="fix"> 252 <when value="fix">
176 <expand macro="ANNOTATION_INPUT" format="gff,gff3,gff3.gz"/> 253 <expand macro="ANNOTATION_INPUT" format="gff,gff3,gff3.gz"/>
254 <expand macro="AGAT_CONFIG"/>
177 </when> 255 </when>
178 <when value="functional_analysis"> 256 <when value="functional_analysis">
179 <expand macro="ANNOTATION_INPUT" format="gff,gtf,gff3,gff3.gz"/> 257 <expand macro="ANNOTATION_INPUT" format="gff,gtf,gff3,gff3.gz"/>
180 <expand macro="REFERENCE_FASTA"/> 258 <expand macro="REFERENCE_FASTA"/>
181 </when> 259 </when>
182 <when value="merge_annotations"> 260 <when value="merge_annotations">
183 <param argument="--gff1" name="input_annotation1" type="data" format="gff,gtf,gff3,gff3.gz" label="Annotation file 1" help="Input GTF/GFF file" /> 261 <param argument="--gff1" name="input_annotation1" type="data" format="gff,gtf,gff3,gff3.gz" label="Annotation file 1" help="Input GTF/GFF file" />
184 <param argument="--gff2" name="input_annotation2" type="data" format="gff,gtf,gff3,gff3.gz" label="Annotation file 2" help="Input GTF/GFF file" /> 262 <param argument="--gff2" name="input_annotation2" type="data" format="gff,gtf,gff3,gff3.gz" label="Annotation file 2" help="Input GTF/GFF file" />
263 <expand macro="AGAT_CONFIG"/>
185 </when> 264 </when>
186 <when value="complement"> 265 <when value="complement">
187 <param argument="--ref" name="input_annotation1" type="data" format="gff,gtf,gff3,gff3.gz" label="Reference annotaiton" help="Reference GTF/GFF file" /> 266 <param argument="--ref" name="input_annotation1" type="data" format="gff,gtf,gff3,gff3.gz" label="Reference annotaiton" help="Reference GTF/GFF file" />
188 <param argument="--add" name="input_annotation2" type="data" format="gff,gtf,gff3,gff3.gz" label="Annotation to complement" help="Annotation file you would like to use to complement the reference annotation." /> 267 <param argument="--add" name="input_annotation2" type="data" format="gff,gtf,gff3,gff3.gz" label="Annotation to complement" help="Annotation file you would like to use to complement the reference annotation." />
189 <param argument="--size_min" type="integer" min="0" value="0" label="Minimun CDS size" help="Option to keep the non-overlping gene only if the CDS size (in nucleotide) is over the minimum 268 <param argument="--size_min" type="integer" min="0" value="0" label="Minimun CDS size" help="Option to keep the non-overlping gene only if the CDS size (in nucleotide) is over the minimum
190 size defined. Default = 0 that means all of them are kept." /> 269 size defined. Default = 0 that means all of them are kept." />
270 <expand macro="AGAT_CONFIG"/>
271 </when>
272 <when value="splice_sites">
273 <expand macro="ANNOTATION_INPUT" format="gff,gff3,gff3.gz"/>
274 <expand macro="AGAT_CONFIG"/>
191 </when> 275 </when>
192 </conditional> 276 </conditional>
193 </inputs> 277 </inputs>
194 <outputs> 278 <outputs>
195 <data name="annotation_gff" format="gff" label="${tool.name} on ${on_string}: annotation file (GFF)"> 279 <data name="annotation_gff" format="gff" label="${tool.name} on ${on_string}: annotation file (GFF)">
196 <filter>tool['selector'] not in ['annotation_statistics','extract','functional_analysis','compare','convert_GFF2GTF','filter_feature_fasta']</filter> 280 <filter>tool['selector'] == 'convert_GTF2GFF'</filter>
197 </data> 281 </data>
198 <data name="annotation_gtf" format="gtf" label="${tool.name} on ${on_string}: annotation file (GTF)"> 282 <data name="annotation_gtf" format="gtf" label="${tool.name} on ${on_string}: annotation file (GTF)">
199 <filter>tool['selector'] == 'convert_GFF2GTF'</filter> 283 <filter>tool['selector'] == 'convert_GFF2GTF'</filter>
200 </data> 284 </data>
201 <data name="features_filtered" format="tabular" label="${tool.name} on ${on_string}: filtered results"> 285 <data name="annotation" format="gff" label="${tool.name} on ${on_string}: annotation file">
202 <filter>tool['selector'] == 'filter_feature_fasta'</filter> 286 <filter>tool['selector'] in ['fix','merge_annotations','complement','filter_feature_fasta','splice_sites','bam2gff']</filter>
287 <change_format>
288 <when input="output_format.selector" value="GTF" format="gtf" />
289 </change_format>
203 </data> 290 </data>
204 <data name="sequence_output" format="fasta" label="${tool.name} on ${on_string}: FASTA file"> 291 <data name="sequence_output" format="fasta" label="${tool.name} on ${on_string}: FASTA file">
205 <filter>tool['selector'] =='extract'</filter> 292 <filter>tool['selector'] =='extract'</filter>
206 </data> 293 </data>
207 <data name="stats_output" format="txt" label="${tool.name} on ${on_string}: stats file"> 294 <data name="stats_output" format="txt" label="${tool.name} on ${on_string}: stats file">
226 <param name="source" value="history"/> 313 <param name="source" value="history"/>
227 <param name="history_item" value="genome.fasta.gz"/> 314 <param name="history_item" value="genome.fasta.gz"/>
228 </conditional> 315 </conditional>
229 </conditional> 316 </conditional>
230 <output name="stats_output" file="test01_stats.txt" ftype="txt"/> 317 <output name="stats_output" file="test01_stats.txt" ftype="txt"/>
231 <output_collection name="distribution_plots_woiso" type="list" count="4">
232 <element name="transcriptClass_cds" file="test01_plot2.pdf" ftype="pdf" compare="sim_size" delta="100"/>
233 </output_collection>
234 <output_collection name="distribution_plots_wiso" type="list" count="4"> 318 <output_collection name="distribution_plots_wiso" type="list" count="4">
235 <element name="transcriptClass_cds" file="test01_plot1.pdf" ftype="pdf" compare="sim_size" delta="100"/> 319 <element name="transcriptClass_cds" file="test01_plot1.pdf" ftype="pdf" compare="sim_size" delta="100"/>
236 </output_collection> 320 </output_collection>
237 </test> 321 </test>
238 <!-- Test 02: extract sequences --> 322 <!-- Test 02: extract sequences -->
257 <param name="input_annotation1" value="annotation.gtf"/> 341 <param name="input_annotation1" value="annotation.gtf"/>
258 <param name="input_annotation2" value="annotation_small.gtf"/> 342 <param name="input_annotation2" value="annotation_small.gtf"/>
259 </conditional> 343 </conditional>
260 <output name="stats_output" file="test03.txt" ftype="txt" lines_diff="2"/> 344 <output name="stats_output" file="test03.txt" ftype="txt" lines_diff="2"/>
261 </test> 345 </test>
262 <!-- Test 04: comlement annotation --> 346 <!-- Test 04: complement annotation -->
263 <test expect_num_outputs="1"> 347 <test expect_num_outputs="1">
264 <conditional name="tool"> 348 <conditional name="tool">
265 <param name="selector" value="complement"/> 349 <param name="selector" value="complement"/>
266 <param name="input_annotation1" value="annotation_small.gtf" ftype="gtf"/> 350 <param name="input_annotation1" value="annotation_small.gtf" ftype="gtf"/>
267 <param name="input_annotation2" value="annotation_unique.gtf" ftype="gtf"/> 351 <param name="input_annotation2" value="annotation_unique.gtf" ftype="gtf"/>
268 <param name="size_min" value="10"/> 352 <param name="size_min" value="10"/>
353 <conditional name="output_format">
354 <param name="selector" value="gff"/>
355 <param name="version" value="3"/>
356 </conditional>
269 </conditional> 357 </conditional>
270 <output name="annotation_gff" file="test04.gff" ftype="gff"/> 358 <output name="annotation_gff" file="test04.gff" ftype="gff"/>
271 </test> 359 </test>
272 <!-- Test 05: Convert GFF2GTF --> 360 <!-- Test 05: Convert GFF2GTF -->
273 <test expect_num_outputs="1"> 361 <test expect_num_outputs="1">
294 <conditional name="reference_genome"> 382 <conditional name="reference_genome">
295 <param name="source" value="history"/> 383 <param name="source" value="history"/>
296 <param name="history_item" value="genome.fasta.gz"/> 384 <param name="history_item" value="genome.fasta.gz"/>
297 </conditional> 385 </conditional>
298 </conditional> 386 </conditional>
299 <output name="features_filtered" file="test07.tabular" ftype="tabular"/> 387 <output name="annotation" file="test07.gff" ftype="gff"/>
300 </test> 388 </test>
301 <!-- Test 08: Fix annotation file --> 389 <!-- Test 08: Fix annotation file -->
302 <test expect_num_outputs="1"> 390 <test expect_num_outputs="1">
303 <conditional name="tool"> 391 <conditional name="tool">
304 <param name="selector" value="fix"/> 392 <param name="selector" value="fix"/>
326 <conditional name="tool"> 414 <conditional name="tool">
327 <param name="selector" value="merge_annotations"/> 415 <param name="selector" value="merge_annotations"/>
328 <param name="input_annotation1" value="annotation_small.gtf"/> 416 <param name="input_annotation1" value="annotation_small.gtf"/>
329 <param name="input_annotation2" value="annotation_unique.gtf"/> 417 <param name="input_annotation2" value="annotation_unique.gtf"/>
330 </conditional> 418 </conditional>
419 <conditional name="output_format">
420 <param name="selector" value="gff"/>
421 <param name="version" value="3"/>
422 </conditional>
331 <output name="annotation_gff" file="test10.gff" ftype="gff"/> 423 <output name="annotation_gff" file="test10.gff" ftype="gff"/>
332 </test> 424 </test>
333 <!-- Test 11: Test compressed files --> 425 <!-- Test 11: Test compressed files -->
334 <test expect_num_outputs="1"> 426 <test expect_num_outputs="1">
335 <conditional name="tool"> 427 <conditional name="tool">
354 </conditional> 446 </conditional>
355 <assert_stdout> 447 <assert_stdout>
356 <has_text text="Job done" /> 448 <has_text text="Job done" />
357 </assert_stdout> 449 </assert_stdout>
358 </test> 450 </test>
451 <!-- Test 13: Add splicing sites -->
452 <test expect_num_outputs="1">
453 <conditional name="tool">
454 <param name="selector" value="splice_sites"/>
455 <param name="gff" value="test04.gff" ftype="gff"/>
456 </conditional>
457 <output name="annotation" file="test13.gff" ftype="gff"/>
458 </test>
359 </tests> 459 </tests>
360 <help><![CDATA[ 460 <help><![CDATA[
361 461
362 .. class:: infomark 462 .. class:: infomark
363 463