Mercurial > repos > greg > draw_features
comparison draw_features.xml @ 6:d5a5d145cfce draft
Uploaded
author | greg |
---|---|
date | Fri, 10 Mar 2023 18:02:00 +0000 |
parents | c318e520dca1 |
children |
comparison
equal
deleted
inserted
replaced
5:c318e520dca1 | 6:d5a5d145cfce |
---|---|
15 | 15 |
16 mkdir feature_hits_dir && | 16 mkdir feature_hits_dir && |
17 mkdir output_dir && | 17 mkdir output_dir && |
18 | 18 |
19 #for $i in $feature_hits: | 19 #for $i in $feature_hits: |
20 #set file_name = $i.file_name | |
21 #set identifier = re.sub('[^\s\w\-\\.]', '_', str($i.element_identifier)) | 20 #set identifier = re.sub('[^\s\w\-\\.]', '_', str($i.element_identifier)) |
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' && | 21 ln -s $i 'feature_hits_dir/$identifier' && |
23 #end for | 22 #end for |
24 | 23 |
25 python '$__tool_directory__/draw_features.py' | 24 python '$__tool_directory__/draw_features.py' |
26 --contigs 'contigs.fasta' | 25 --contigs 'contigs.fasta' |
27 --feature_hits_dir 'feature_hits_dir' | 26 --feature_hits_dir 'feature_hits_dir' |
58 </collection> | 57 </collection> |
59 </param> | 58 </param> |
60 <output_collection name="features_bed" type="list" count="1"> | 59 <output_collection name="features_bed" type="list" count="1"> |
61 <element name="amr_cds.bed" ftype="bed"> | 60 <element name="amr_cds.bed" ftype="bed"> |
62 <assert_contents> | 61 <assert_contents> |
63 <has_size value="129" delta="10"/> | 62 <has_size value="610" delta="10"/> |
64 </assert_contents> | 63 </assert_contents> |
65 </element> | 64 </element> |
66 </output_collection> | 65 </output_collection> |
67 <output_collection name="features_png" type="list" count="1"> | 66 <output_collection name="features_png" type="list" count="1"> |
68 <element name="contig_1" ftype="png"> | 67 <element name="contig_1" ftype="png"> |