Mercurial > repos > mheinzl > variant_analyzer2

comparison read2mut.py @ 53:27b00e38e13d 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:07:30 +0000
parents 49ecab32c83f
children 95c27bcb1b7a
comparison
equal deleted inserted replaced
52:49ecab32c83f 53:27b00e38e13d
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(sorted(tier_dict[key1].keys()))) 1074 print(list(sorted(tier_dict[key1].keys())))
1075 print(list(sorted(tier_dict[key1].keys()))[:6]) 1075 print(list(sorted(tier_dict[key1].keys()))[:6])
1076 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(sorted(tier_dict[key1].keys()))[:6]])
1077 print(sum_highTiers) 1077 print(sum_highTiers)
1078 if tier_dict[key1]["tier 4"] > 0 and sum_highTiers > 0: 1078 if tier_dict[key1]["tier 4"] > 0 and sum_highTiers > 0:
1079 tier_dict[key1]["tier 2.5"] = tier_dict[key1]["tier 4"] 1079 tier_dict[key1]["tier 2.5"] = tier_dict[key1]["tier 4"]
1080 tier_dict[key1]["tier 4"] = 0 1080 tier_dict[key1]["tier 4"] = 0
1081 lines = change_tier_after_print[key1] 1081 lines = change_tier_after_print[key1]