# HG changeset patch # User mheinzl # Date 1682535352 0 # Node ID 24f166c1dba77152dadffa7890c5fef85f5e9c60 # Parent 1a5974404d4fb3b95ba1f1c7b6f0e22c60b06756 planemo upload for repository https://github.com/Single-Molecule-Genetics/VariantAnalyzerGalaxy/tree/master/tools/variant_analyzer commit ee4a8e6cf290e6c8a4d55f9cd2839d60ab3b11c8-dirty diff -r 1a5974404d4f -r 24f166c1dba7 read2mut.py --- 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())[0] # get index of max HD diff -r 1a5974404d4f -r 24f166c1dba7 read2mut.xml --- a/read2mut.xml Tue Apr 25 17:06:38 2023 +0000 +++ b/read2mut.xml Wed Apr 26 18:55:52 2023 +0000 @@ -1,5 +1,5 @@ - + Looks for reads with a mutation at known positions and calculates frequencies and stats. va_macros.xml