# HG changeset patch # User greg # Date 1675793794 0 # Node ID fcf01dad2ceb73f0202c970c97fd0b66b68ec5b7 # Parent 7d47800ee5ac4d28c6774864b5dc270cdf48eb9d Uploaded diff -r 7d47800ee5ac -r fcf01dad2ceb draw_features.xml --- 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'