Add test for homopolymer indel

michael-ford · Mar 15, 2024 · ef3cfda · ef3cfda
1 parent d6421d7
commit ef3cfda
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/tests/test_alignment_variants/test_call_variants.py b/tests/test_alignment_variants/test_call_variants.py
@@ -18,6 +18,15 @@ def test_call_variants_insertion():
     variants = call_variants(assembly_seq, cigar_list, wildtype_seq)
     assert variants == {(3, "insG")}
 
+def test_call_variants_homo_insertion():
+    assembly_seq = "AGCGGGTA"
+    wildtype_seq = "AGCTA"  # Insertion at position 3
+    cigar_list = [(3, "="), (3, "D"), (2, "=")]
+
+    variants = call_variants(assembly_seq, cigar_list, wildtype_seq)
+    assert variants == {(3, "insGGG")}
+
+
 def test_call_variants_deletion():
     assembly_seq = "AGTA"
     wildtype_seq = "AGCTA"  # Deletion at position 3