planemo upload for repository https://github.com/Single-Molecule-Genetics/VariantAnalyzerGalaxy/tree/master/tools/variant_analyzer commit ee4a8e6cf290e6c8a4d55f9cd2839d60ab3b11c8-dirty
author mheinzl Wed, 26 Apr 2023 18:55:52 +0000 1a5974404d4f
line wrap: on
line diff
```--- a/read2mut.py	Tue Apr 25 17:06:38 2023 +0000
+++ b/read2mut.py	Wed Apr 26 18:55:52 2023 +0000
@@ -24,7 +24,6 @@

import argparse
import csv
-import itertools
import json
import operator
import os
@@ -1295,14 +1294,14 @@
half1_mate2 = array2_half2
half2_mate2 = array2_half
# calculate HD of "a" in the tag to all "a's" or "b" in the tag to all "b's"
-                                    dist = np.array([sum(itertools.map(operator.ne, half1_mate1, c)) for c in half1_mate2])
+                                    dist = np.array([sum(map(operator.ne, half1_mate1, c)) for c in half1_mate2])
min_index = np.where(dist == dist.min())  # get index of min HD
# get all "b's" of the tag or all "a's" of the tag with minimum HD
min_tag_half2 = half2_mate2[min_index]
min_tag_array2 = array2[min_index]  # get whole tag with min HD
min_value = dist.min()
# calculate HD of "b" to all "b's" or "a" to all "a's"
-                                    dist_second_half = np.array([sum(itertools.map(operator.ne, half2_mate1, e))
+                                    dist_second_half = np.array([sum(map(operator.ne, half2_mate1, e))
for e in min_tag_half2])
dist2 = dist_second_half.max()
max_index = np.where(dist_second_half == dist_second_half.max())  # get index of max HD```