from opencloning.assembly2 import Assembly, assembly2str, assemble
from pydna.dseqrecord import Dseqrecord
template = Dseqrecord('cccgaggggaatcgaa')
insert = Dseqrecord('ggggaatcAcccgag')
asm = Assembly((template, insert), limit=5, use_all_fragments=True)
for a in asm.get_insertion_assemblies():
if a[0][0] == 1:
print(assembly2str(a))
prod = assemble(asm.fragments, a, is_insertion=True)
print(prod.seq)
print()
prints
('1[5:13]:2[0:8]', '2[9:15]:1[0:6]')
cccgaggggaatcAcccgaggggaatcgaa
This is the assembly:
1 cccgaggggaatcgaa
2 Acccgagggggaatc
cccgaggggaatcgaa 1[5:13]
||||||||
ggggaatcAcccgag 2[0:8]
ggggaatcAcccgag 2[9:15]
||||||
cccgaggggaatcgaa 1[0:6]
prints
This is the assembly: