Mercurial > repos > mheinzl > variant_analyzer2

changeset 85:d1cd4cd9f18d 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-dirty
author mheinzl
date Wed, 24 Aug 2022 09:47:08 +0000
parents e46d5e377760
children 97bd9c7a1b44
files read2mut.py
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/read2mut.py	Fri Aug 19 11:23:37 2022 +0000
+++ b/read2mut.py	Wed Aug 24 09:47:08 2022 +0000
@@ -477,9 +477,9 @@
                     continue
 
                 if (((key2[:-5] in tag_dict.keys()) and (key2[:-5] in pure_tags_dict_short[key1].keys()) and (key1 in tag_dict[key2[:-5]].keys()) and (key2[:-5] not in used_keys)) or ((key2[:-5] in tag_dict_ref.keys()) and (key2[:-5] in pure_tags_dict_ref[key1].keys()) and (key1 in tag_dict_ref[key2[:-5]].keys()) and (key2[:-5] not in used_keys))):
-                    if key2[:-5] in tag_dict.keys():
+                    if key2[:-5] in tag_dict.keys() and key1 in tag_dict[key2[:-5]].keys():
                         variant_type = "alt"
-                    elif key2[:-5] in tag_dict_ref.keys():
+                    elif key2[:-5] in tag_dict_ref.keys() and key1 in tag_dict_ref[key2[:-5]].keys():
                         variant_type = "ref"
 
                     if refalttiers is False and variant_type == "ref":  # if we only want information about alt tiers, skip all refs