digest.py

from Bio.Seq import Seq
from Bio import SeqIO
from Bio.Restriction import *

rb=
rb = DpnII + HinfI + DdeI + MseI
dpnii_dig=[]
for record in SeqIO.parse("GCA_905146935.1_idScaPyra1.1_genomic.fna", "fasta"):
	my_seq = Seq(str(record.seq))
	dpnii_dig.append(DpnII.catalyse(my_seq))

test_1=list(sum(dpnii_dig, ()))


hinfi_dig=[]
for i in range(len(test_1)):
	newSeq= test_1[i]
	hinfi_dig.append(HinfI.catalyse(my_seq))
print(hinfi_dig)

test_2=list(sum(hinfi_dig, ()))

ddei_dig=[]
for i in range(len(test_1)):
	newSeq= test_1[i]
	hinfi_dig.append(HinfI.catalyse(my_seq))