# HG changeset patch # User mheinzl # Date 1615474709 0 # Node ID 26e53b5b8bcbbdcd03e507af72478deced42bf07 # Parent b32c973fe8abc4d48cc91867b9655ad23b727782 planemo upload for repository https://github.com/Single-Molecule-Genetics/VariantAnalyzerGalaxy/tree/master/tools/variant_analyzer commit ee4a8e6cf290e6c8a4d55f9cd2839d60ab3b11c8 diff -r b32c973fe8ab -r 26e53b5b8bcb read2mut.py --- a/read2mut.py Thu Mar 11 11:04:05 2021 +0000 +++ b/read2mut.py Thu Mar 11 14:58:29 2021 +0000 @@ -324,8 +324,8 @@ whole_array = pure_tags_dict_short[key1].keys() tier_dict[key1] = {} - values_tier_dict = [("tier 1.1", 0), ("tier 1.2", 0), ("tier 2.1", 0), ("tier 2.2", 0), ("tier 2.3", 0), ("tier 2.4", 0), ("tier 3.1", 0), - ("tier 3.2", 0), ("tier 2.5", 0), ("tier 4", 0), ("tier 5.1", 0), ("tier 5.2", 0), ("tier 5.3", 0), ("tier 5.4", 0), ("tier 5.5", 0), + values_tier_dict = [("tier 1.1", 0), ("tier 1.2", 0), ("tier 2.1", 0), ("tier 2.2", 0), ("tier 2.3", 0), ("tier 2.4", 0),("tier 2.5", 0), + ("tier 3.1", 0), ("tier 3.2", 0), ("tier 4", 0), ("tier 5.1", 0), ("tier 5.2", 0), ("tier 5.3", 0), ("tier 5.4", 0), ("tier 5.5", 0), ("tier 6", 0), ("tier 7", 0)] for k, v in values_tier_dict: tier_dict[key1][k] = v @@ -1071,7 +1071,9 @@ chimera_dict[key1] = (chimeric_dcs, chimeric_dcs_high_tiers) # move tier 4 counts to tier 2.5 if there other mutations with tier <= 2.4 - sum_highTiers = sum([tier_dict[key1][ij] for ij in tier_dict[key1].keys()[:6]]) + print(list(tier_dict[key1].keys())[:6]) + sum_highTiers = sum([tier_dict[key1][ij] for ij in list(tier_dict[key1].keys())[:6]]) + print(sum_highTiers) if tier_dict[key1]["tier 4"] > 0 and sum_highTiers > 0: tier_dict[key1]["tier 2.5"] = tier_dict[key1]["tier 4"] tier_dict[key1]["tier 4"] = 0