Mercurial > repos > mheinzl > variant_analyzer2

comparison read2mut.py @ 52:49ecab32c83f draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 15:03:42 +0000
parents 26e53b5b8bcb
children 27b00e38e13d
comparison
equal deleted inserted replaced
51:26e53b5b8bcb 52:49ecab32c83f
1069 else: 1069 else:
1070 chimeric_dcs_high_tiers += high_tiers 1070 chimeric_dcs_high_tiers += high_tiers
1071 chimera_dict[key1] = (chimeric_dcs, chimeric_dcs_high_tiers) 1071 chimera_dict[key1] = (chimeric_dcs, chimeric_dcs_high_tiers)
1072 1072
1073 # move tier 4 counts to tier 2.5 if there other mutations with tier <= 2.4 1073 # move tier 4 counts to tier 2.5 if there other mutations with tier <= 2.4
1074 print(list(tier_dict[key1].keys())[:6]) 1074 print(list(sorted(tier_dict[key1].keys())))
1075 print(list(sorted(tier_dict[key1].keys()))[:6])
1075 sum_highTiers = sum([tier_dict[key1][ij] for ij in list(tier_dict[key1].keys())[:6]]) 1076 sum_highTiers = sum([tier_dict[key1][ij] for ij in list(tier_dict[key1].keys())[:6]])
1076 print(sum_highTiers) 1077 print(sum_highTiers)
1077 if tier_dict[key1]["tier 4"] > 0 and sum_highTiers > 0: 1078 if tier_dict[key1]["tier 4"] > 0 and sum_highTiers > 0:
1078 tier_dict[key1]["tier 2.5"] = tier_dict[key1]["tier 4"] 1079 tier_dict[key1]["tier 2.5"] = tier_dict[key1]["tier 4"]
1079 tier_dict[key1]["tier 4"] = 0 1080 tier_dict[key1]["tier 4"] = 0