DNA Reading + Digestion site

from Bio.Seq import Seq
from Bio.Alphabet import IUPAC


DNA0 = input("Plak de DNA sequentie")
seq = DNA0.upper()
enz0 = input("plak de enzym naam:")
enz = enz0.upper()

y= seq[0]
y_minus = 0

restriction_site = ""
restriction_site = 0
fragment = ""
enz = ""


A = DNA0.count("A")
T = DNA0.count("T")
C = DNA0.count("C")
G = DNA0.count("G")

my_seq = Seq(DNA0, IUPAC.unambiguous_dna)

def transcribe(sequence):
    return sequence.replace('T', 'U')

def digestion (seq,enz):
    if (enz == "XBAI"): recognition_site = "TCTAGA"
    elif (enz == "ECORI"): recognition_site = "GAATTC"
    else:
        print("Enzym niet herkend")
    fragment = (seq.split(recognition_site))
    print("De digestie van het DNA fragment"+ enz +", De resulterende DNA fragmenten na digestie")
    for No_fr, i in enumerate(fragment):
        print(No_fr,i, end= "\n")
    str_regDNA0 = "(_^_)".join(map(str,fragment))
    str_regDNA1 = str(str_regDNA0[0:-1])
    print("\nHet aantal DNA fragmenten na disgestie ziet er als volgt uit")
    return(str_regDNA1)

def translate_rna(s):
    codon2aa = {"AAA": "K", "AAC": "N", "AAG": "K", "AAU": "N",
                "ACA": "T", "ACC": "T", "ACG": "T", "ACU": "T",
                "AGA": "R", "AGC": "S", "AGG": "R", "AGU": "S",
                "AUA": "I", "AUC": "I", "AUG": "M", "AUU": "I",

                "CAA": "Q", "CAC": "H", "CAG": "Q", "CAU": "H",
                "CCA": "P", "CCC": "P", "CCG": "P", "CCU": "P",
                "CGA": "R", "CGC": "R", "CGG": "R", "CGU": "R",
                "CUA": "L", "CUC": "L", "CUG": "L", "CUU": "L",

                "GAA": "E", "GAC": "D", "GAG": "E", "GAU": "D",
                "GCA": "A", "GCC": "A", "GCG": "A", "GCU": "A",
                "GGA": "G", "GGC": "G", "GGG": "G", "GGU": "G",
                "GUA": "V", "GUC": "V", "GUG": "V", "GUU": "V",

                "UAA": "_", "UAC": "Y", "UAG": "_", "UAU": "T",
                "UCA": "S", "UCC": "S", "UCG": "S", "UCU": "S",
                "UGA": "_", "UGC": "C", "UGG": "W", "UGU": "C",
                "UUA": "L", "UUC": "F", "UUG": "L", "UUU": "F"}

    l = [codon2aa.get(s[n:n + 3], 'X') for n in range(0, len(s), 3)]
    return "".join(l)

RNA = transcribe(DNA0)

print("De volledige sequentie is :", (A + C + T + G), "bp lang.")
print("Het aantal A bp in de sequentie is :", A)
print("Het aantal T bp in de sequentie is :", T)
print("Het aantal C bp in de sequentie is :", C)
print("Het aantal G bp in de sequentie is :", G)
print("De GC ratio van de sequence is :", (G + C) / (G + C + T + A), '\n')

print("De complimentaire DNA sequentie is:")
print (my_seq.complement(), '\n')

print("De reversecomplimentaire DNA sequentie is:")
print(my_seq.reverse_complement(), '\n')

print("De RNA sequentie is:")
print(RNA, '\n')

print("De Eiwitsequentie is:")
protein = translate_rna(RNA)
print(protein, '\n')
print("Het aantal aminozuren van het eiwit is:")
print(len(protein) - 1, '\n\n\n')
print(digestion(DNA0, enz0),'\n\n\n')
print("Deze software is het eigendom van Loeke van der Linden, De restrictie software is eigendom van Kevin Harrijvan")
print("Deze Software kan daarom niet instaan voor de juistheid of (on)volledigheid van de gepubliceerde informatie.")
print("Loeke aanvaardt geen enkele verantwoordelijkheid of aansprakelijkheid voor eventuele schade, van welke aard dan ook, die op welke wijze dan ook voortvloeit uit het gebruik")