annotate draw_features.xml @ 5:c318e520dca1 draft

Uploaded
author greg
date Thu, 02 Mar 2023 16:11:22 +0000
parents 3ba596014675
children d5a5d145cfce
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
1 <tool id="draw_features" name="PIMA: draw features" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
2 <description>of annotated assemblies</description>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
3 <macros>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
4 <import>macros.xml</import>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
5 </macros>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
6 <expand macro="requirements"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
7 <command detect_errors="exit_code"><![CDATA[
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
8 #import re
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
9
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
10 #if $contigs.is_of_type('fasta.gz'):
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
11 gunzip -c '$contigs' > 'contigs.fasta' &&
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
12 #else:
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
13 ln -s '$contigs' 'contigs.fasta' &&
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
14 #end if
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
15
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
16 mkdir feature_hits_dir &&
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
17 mkdir output_dir &&
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
18
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
19 #for $i in $feature_hits:
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
20 #set file_name = $i.file_name
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
21 #set identifier = re.sub('[^\s\w\-\\.]', '_', str($i.element_identifier))
2
fcf01dad2ceb Uploaded
greg
parents: 0
diff changeset
22 cat '$i' | awk '$7 != "."' | awk '{OFS="\t";locus=$7"\t"$8"\t"$9; if($5 > s[locus]){s[locus]=$5;id = sprintf("%.3f", $5); b[locus] = $1"\t"$2"\t"$3"\t"$4"\t"id"\t"$6}} END{for(i in b){print b[i]}}' | sort -k 1,1 -k2,2n > 'feature_hits_dir/$identifier' &&
0
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
23 #end for
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
24
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
25 python '$__tool_directory__/draw_features.py'
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
26 --contigs 'contigs.fasta'
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
27 --feature_hits_dir 'feature_hits_dir'
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
28 --output_dir 'output_dir'
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
29 #if str($output_process_log) == 'yes':
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
30 && mv 'process_log' '$process_log'
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
31 #end if
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
32 ]]></command>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
33 <inputs>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
34 <param name="contigs" type="data" format="fasta,fasta.gz" label="Fasta file of assembled contigs"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
35 <param name="feature_hits" format="bed" type="data_collection" collection_type="list" label="Collection of feature hits BED files"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
36 <param name="output_process_log" type="select" display="radio" label="Output process log file?">
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
37 <option value="no" selected="true">No</option>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
38 <option value="yes">Yes</option>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
39 </param>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
40 </inputs>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
41 <outputs>
4
3ba596014675 Uploaded
greg
parents: 2
diff changeset
42 <data name="process_log" format="txt" label="${tool.name} on ${on_string} (process log)">
0
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
43 <filter>output_process_log == 'yes'</filter>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
44 </data>
4
3ba596014675 Uploaded
greg
parents: 2
diff changeset
45 <collection name="features_bed" type="list" format="bed" label="${tool.name} on ${on_string} (best bed)">
5
c318e520dca1 Uploaded
greg
parents: 4
diff changeset
46 <discover_datasets pattern="__name__" directory="feature_hits_dir" format="bed"/>
4
3ba596014675 Uploaded
greg
parents: 2
diff changeset
47 </collection>
3ba596014675 Uploaded
greg
parents: 2
diff changeset
48 <collection name="features_png" type="list" format="png" label="${tool.name} on ${on_string} (png)">
0
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
49 <discover_datasets pattern="(?P&lt;designation&gt;.+)\.(?P&lt;ext&gt;png)" directory="output_dir"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
50 </collection>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
51 </outputs>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
52 <tests>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
53 <test>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
54 <param name="contigs" value="PS01519_contigs.fasta.gz" ftype="fasta.gz"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
55 <param name="feature_hits">
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
56 <collection type="list">
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
57 <element name="amr_cds.bed" value="amr_cds.bed"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
58 </collection>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
59 </param>
4
3ba596014675 Uploaded
greg
parents: 2
diff changeset
60 <output_collection name="features_bed" type="list" count="1">
5
c318e520dca1 Uploaded
greg
parents: 4
diff changeset
61 <element name="amr_cds.bed" ftype="bed">
4
3ba596014675 Uploaded
greg
parents: 2
diff changeset
62 <assert_contents>
3ba596014675 Uploaded
greg
parents: 2
diff changeset
63 <has_size value="129" delta="10"/>
3ba596014675 Uploaded
greg
parents: 2
diff changeset
64 </assert_contents>
3ba596014675 Uploaded
greg
parents: 2
diff changeset
65 </element>
3ba596014675 Uploaded
greg
parents: 2
diff changeset
66 </output_collection>
0
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
67 <output_collection name="features_png" type="list" count="1">
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
68 <element name="contig_1" ftype="png">
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
69 <assert_contents>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
70 <has_size value="25383" delta="100"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
71 </assert_contents>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
72 </element>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
73 </output_collection>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
74 </test>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
75 </tests>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
76 <help>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
77 **What it does**
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
78
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
79 Accepts a FASTA file of assembled contigs and a collection of 2 BED files containing the feature hits and plots
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
80 the features.
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
81 </help>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
82 <expand macro="citations"/>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
83 </tool>
636eeb4fb9ac Uploaded
greg
parents:
diff changeset
84