diff draw_features.xml @ 2:fcf01dad2ceb draft

Uploaded
author greg
date Tue, 07 Feb 2023 18:16:34 +0000
parents 636eeb4fb9ac
children 3ba596014675
line wrap: on
line diff
--- a/draw_features.xml	Thu Jan 26 16:00:37 2023 +0000
+++ b/draw_features.xml	Tue Feb 07 18:16:34 2023 +0000
@@ -19,7 +19,7 @@
 #for $i in $feature_hits:
     #set file_name = $i.file_name
     #set identifier = re.sub('[^\s\w\-\\.]', '_', str($i.element_identifier))
-    ln -s '$file_name' 'feature_hits_dir/$identifier' &&
+    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' &&
 #end for
 
 python '$__tool_directory__/draw_features.py'