annotate draw_amr_matrix.xml @ 10:03240ffe969a draft

Uploaded
author greg
date Tue, 21 Mar 2023 18:46:42 +0000
parents 70073df30a06
children da1c9c1be421
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
1 <tool id="draw_amr_matrix" name="PIMA: draw AMR matrix" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
2 <description>of genes and drugs</description>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
3 <macros>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
4 <import>macros.xml</import>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
5 </macros>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
6 <expand macro="requirements"/>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
7 <command detect_errors="exit_code"><![CDATA[
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
8 #import re
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
9
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
10 mkdir amr_feature_hits_dir &&
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
11 mkdir mutation_regions_dir &&
389c98d344ce Uploaded
greg
parents: 4
diff changeset
12 mkdir amr_matrix_png_dir &&
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
13
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
14 #if $reference_source.reference_source_selector == 'history':
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
15 ln -f -s '$reference_source.ref_file' reference.fa &&
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
16 #else:
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
17 ln -f -s '$reference_source.ref_file.fields.path' reference.fa &&
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
18 #end if
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
19
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
20 #for $i in $amr_feature_hits:
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
21 #set file_name = $i.file_name
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
22 #set identifier = re.sub('[^\s\w\-\\.]', '_', str($i.element_identifier))
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
23 ln -s '$file_name' 'amr_feature_hits_dir/$identifier' &&
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
24 #end for
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
25
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
26 python '$__tool_directory__/draw_amr_matrix.py'
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
27 --amr_feature_hits_dir 'amr_feature_hits_dir'
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
28 #if str($amr_deletions_file) != 'None':
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
29 --amr_deletions_file '$amr_deletions_file'
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
30 #end if
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
31 #if str($varscan_vcf_file) != 'None':
389c98d344ce Uploaded
greg
parents: 4
diff changeset
32 --varscan_vcf_file '$varscan_vcf_file'
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
33 #end if
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
34 #if str($amr_mutation_regions_bed_file) != 'None':
389c98d344ce Uploaded
greg
parents: 4
diff changeset
35 --amr_mutation_regions_bed_file '$amr_mutation_regions_bed_file'
389c98d344ce Uploaded
greg
parents: 4
diff changeset
36 --mutation_regions_dir 'mutation_regions_dir'
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
37 #end if
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
38 --amr_gene_drug_file '$amr_gene_drug_file'
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
39 --reference_genome reference.fa
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
40 --amr_matrix_png_dir 'amr_matrix_png_dir'
9
70073df30a06 Uploaded
greg
parents: 7
diff changeset
41 --errors '$errors'
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
42 #if str($output_process_log) == 'yes':
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
43 && mv 'process_log' '$process_log'
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
44 #end if
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
45 ]]></command>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
46 <inputs>
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
47 <conditional name="reference_source">
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
48 <param name="reference_source_selector" type="select" label="Select a reference genome from your history or use a cached genome index?">
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
49 <option value="cached">Use a cached genome index</option>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
50 <option value="history">Select a genome from the history and build the index</option>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
51 </param>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
52 <when value="cached">
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
53 <param name="ref_file" type="select" label="Using reference genome" help="Select reference genome">
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
54 <options from_data_table="all_fasta">
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
55 <filter type="sort_by" column="2"/>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
56 <validator type="no_options" message="No reference genomes are available"/>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
57 </options>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
58 </param>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
59 </when>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
60 <when value="history">
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
61 <param name="ref_file" type="data" format="fasta,fastq" label="Select the reference sequence" help="You can upload a FASTA file and use it as reference"/>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
62 </when>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
63 </conditional>
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
64 <param argument="--amr_feature_hits" format="bed" type="data_collection" collection_type="list" label="Collection of feature hits BED files"/>
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
65 <param argument="--amr_deletions_file" type="data" format="bed" optional="true" label="AMR deletions file" help="Optional, leave blank to ignore"/>
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
66 <param argument="--varscan_vcf_file" type="data" format="vcf" optional="true" label="Varscan VCF file" help="Optional, leave blank to ignore"/>
389c98d344ce Uploaded
greg
parents: 4
diff changeset
67 <param argument="--amr_mutation_regions_bed_file" type="data" format="bed" optional="true" label="AMR mutation regions BED file" help="Optional, leave blank to ignore"/>
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
68 <param argument="--amr_gene_drug_file" type="data" format="tabular,tsv" label="AMR gene drugs file"/>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
69 <param name="output_process_log" type="select" display="radio" label="Output process log file?">
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
70 <option value="no" selected="true">No</option>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
71 <option value="yes">Yes</option>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
72 </param>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
73 </inputs>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
74 <outputs>
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
75 <data name="process_log" format="txt" label="${tool.name} on ${on_string} (process log)">
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
76 <filter>output_process_log == 'yes'</filter>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
77 </data>
9
70073df30a06 Uploaded
greg
parents: 7
diff changeset
78 <data name="errors" format="txt" label="${tool.name} on ${on_string} (errors)"/>
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
79 <collection name="mutation_regions_tsv" type="list" format="tsv" label="${tool.name} on ${on_string} (mutation regions)">
389c98d344ce Uploaded
greg
parents: 4
diff changeset
80 <filter>amr_mutation_regions_bed_file not in [None, 'None']</filter>
389c98d344ce Uploaded
greg
parents: 4
diff changeset
81 <discover_datasets pattern="(?P&lt;designation&gt;.+)\.(?P&lt;ext&gt;tsv)" directory="mutation_regions_dir"/>
389c98d344ce Uploaded
greg
parents: 4
diff changeset
82 </collection>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
83 <collection name="amr_matrix_png" type="list" format="png">
7
389c98d344ce Uploaded
greg
parents: 4
diff changeset
84 <discover_datasets pattern="(?P&lt;designation&gt;.+)\.(?P&lt;ext&gt;png)" directory="amr_matrix_png_dir"/>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
85 </collection>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
86 </outputs>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
87 <tests>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
88 <test>
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
89 <param name="reference_source_selector" value="history"/>
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
90 <param name="ref_file" ftype="fasta" value="ref_genome.fasta"/>
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
91 <param name="amr_feature_hits">
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
92 <collection type="list">
4
33a0ea992043 Uploaded
greg
parents: 2
diff changeset
93 <element name="amr_pima_md.bed" value="amr_pima_md.bed"/>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
94 </collection>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
95 </param>
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
96 <param name="amr_gene_drug_file" value="amr_gene_drug.tsv" ftype="tsv"/>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
97 <output_collection name="amr_matrix_png" type="list" count="1">
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
98 <element name="contig_1" ftype="png">
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
99 <assert_contents>
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
100 <has_size value="25383" delta="100"/>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
101 </assert_contents>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
102 </element>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
103 </output_collection>
9
70073df30a06 Uploaded
greg
parents: 7
diff changeset
104 <output name="errors" value="errors.txt" ftype="txt"/>
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
105 </test>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
106 </tests>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
107 <help>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
108 **What it does**
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
109
2
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
110 Accepts a collection of best AMR feature hits, an optional AMR deletions BED file, an optional AMR mutations TSV file and a AMR gene
9fcc1ffd7526 Uploaded
greg
parents: 0
diff changeset
111 drug mappings file and draws an AMR matrix.
0
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
112 </help>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
113 <expand macro="citations"/>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
114 </tool>
6044ca44e9f6 Uploaded
greg
parents:
diff changeset
115