polishing gene cov feature in analyze_assembly

[biopieces.git] / bp_bin / analyze_assembly

diff --git a/bp_bin/analyze_assembly b/bp_bin/analyze_assembly
index d6b65a6d7bd4c4fc6e514df2126234724efc1d34..c6423f8e93faecefb05f899289f8bbbdf6bab9c1 100755 (executable)
--- a/bp_bin/analyze_assembly
+++ b/bp_bin/analyze_assembly
@@ -93,7 +93,7 @@ new_record[:COUNT]      = lengths.size
 
 if options[:gene_cov]
   new_record[:GENE_COV]   = gene_cov
-  new_record[:GENE_RATIO] = (gene_cov.to_f / (2 * total.to_f)).round(3)
+  new_record[:GENE_RATIO] = (gene_cov.to_f / total.to_f).round(4)
 end
 
 bp.puts new_record