Test case where parts assemble into an enzyme site

Author: veghp

tests/data/assemblies/fragment_1.gb

@@ -0,0 +1,13 @@
+LOCUS       fragment_1                62 bp    DNA     linear   UNK 01-JAN-1980
+DEFINITION  .
+ACCESSION   fragment_1
+VERSION     fragment_1
+KEYWORDS    .
+SOURCE      
+  ORGANISM  .
+            .
+FEATURES             Location/Qualifiers
+ORIGIN
+        1 cgtctcacta actaagcagc ccggctttcg ccctattgag tctaccgcgg cgtcttgaga
+       61 cg
+//

tests/data/assemblies/fragment_2.gb

@@ -0,0 +1,13 @@
+LOCUS       fragment_2                80 bp    DNA     linear   UNK 01-JAN-1980
+DEFINITION  .
+ACCESSION   fragment_2
+VERSION     fragment_2
+KEYWORDS    .
+SOURCE      
+  ORGANISM  .
+            .
+FEATURES             Location/Qualifiers
+ORIGIN
+        1 cgtctcagtc tcaaatcaat agatgtactg tagaatcctg atagaaaagg cagaagcccg
+       61 gtaatcgcgc taatgagacg
+//

tests/test_type2S_assembly.py

@@ -22,7 +22,7 @@
 #     )
 #     simulation = assembly.simulate(sequence_repository=repo)
 #     assert len(simulation.construct_records) == 1
-#     # We expect two errors: 
+#     # We expect two errors:
 #     assert len(simulation.errors) == 2
 
 #     assembly = dc.Type2sRestrictionAssembly(
@@ -48,3 +48,14 @@ def test_single_assembly_with_wrong_enzyme(tmpdir):
 #     )
 #     simulation = assembly.simulate(sequence_repository=repo)
 #     assert len(simulation.construct_records) == 4
+
+
+def test_assembling_into_new_enzyme_site():
+    assembly = dc.Type2sRestrictionAssembly(
+        parts=["fragment_1", "fragment_2"], enzyme="BsmBI"
+    )
+    simulation = assembly.simulate(sequence_repository=repo)
+    assert len(simulation.construct_records) == 1
+    assert len(simulation.errors) == 1
+    assert simulation.errors[0].message == "Assembly creates a new enzyme site"
+