Merge pull request #33 from ammaraziz/master

Support for amino acids, correct ambiguous handling
aineniamh · Apr 12, 2024 · 863e518 · 863e518
2 parents 55aa31b + 258303a
commit 863e518
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/snipit/scripts/snp_functions.py b/snipit/scripts/snp_functions.py
@@ -261,7 +261,14 @@ def find_snps(reference_seq,input_seqs,show_indels):
 
     return snp_dict,record_snps,len(var_counter)
 
-def find_ambiguities(alignment, snp_dict):
+def find_ambiguities(alignment,snp_dict,sequence_type):
+
+
+    if sequence_type == "nt":
+        amb = NT_AMBIG
+    if sequence_type == "aa":
+        amb = AA_AMBIG
+
 
     snp_sites = collections.defaultdict(list)
     for seq in snp_dict: