Mercurial > repos > mheinzl > variant_analyzer2
diff read2mut.py @ 51:26e53b5b8bcb draft
planemo upload for repository https://github.com/Single-Molecule-Genetics/VariantAnalyzerGalaxy/tree/master/tools/variant_analyzer commit ee4a8e6cf290e6c8a4d55f9cd2839d60ab3b11c8
| author | mheinzl |
|---|---|
| date | Thu, 11 Mar 2021 14:58:29 +0000 |
| parents | b32c973fe8ab |
| children | 49ecab32c83f |
line wrap: on
line diff
--- 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