Mercurial > repos > mheinzl > hd

diff hd.py @ 18:a8581bf627fd draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/hd commit b09f6822cb5acf314c5034323501d6300fda25cc
author mheinzl
date Wed, 23 May 2018 14:47:43 -0400
parents 4a30328b1af9
children 2e9f7ea7ae93
line wrap: on
line diff
--- a/hd.py	Wed May 23 14:41:47 2018 -0400
+++ b/hd.py	Wed May 23 14:47:43 2018 -0400
@@ -977,7 +977,7 @@
             # chimeric reads
             if len(minHD_tags_zeros) != 0:
                 # absolute difference and tags where at least one half has HD=0
-                summary15, sumCol15 = createTableHD(listDifference1_zeros, "diff=")
+                summary15, sumCol15 = createTableHD(listDifference1_zeros, "HD=")
                 overallSum15 = sum(sumCol15)
                 # absolute difference and tags where at least one half has HD=0
            #     summary23, sumCol23 = createTableFSD2(familySizeList1_diff_zeros, diff=False)
@@ -996,7 +996,8 @@
                            diff=False)
 
             count = numpy.bincount(quant)
-            output_file.write("{}{}\n".format(sep, f))
+            #output_file.write("{}{}\n".format(sep, name_file))
+            output_file.write("\n")
             output_file.write("max. family size:{}{}\n".format(sep, max(quant)))
             output_file.write("absolute frequency:{}{}\n".format(sep, count[len(count) - 1]))
             output_file.write(