def AT(seq):
at = sum(seq.count(x) for x in ["a", "t", "A", "T"])
try:
c = round(at / len(seq), 2)
return c
except ZeroDivisionError:
return 0.0
def GC(seq):
gc = sum(seq.count(h) for h in ["g", "c", "G", "C"])
try:
y = round((gc / len(seq))*100, 2)
return y
except ZeroDivisionError:
return 0.0
dna_seq = "GTCTCTACCTTGACAGACCTCCAGCCGTACATGCGACAGTTCGTGGCTCACCTGCAGGAGACCAGCCCGCTGAGGGATGCCGTCGTCATCGAGCAG"
at_content = AT(dna_seq)
print(at_content)
gc_content = GC(dna_seq)
print(gc_content)
ZGVmIEFUKHNlcSk6CiAgICBhdCA9IHN1bShzZXEuY291bnQoeCkgZm9yIHggaW4gWyJhIiwgInQiLCAiQSIsICJUIl0pCiAgICB0cnk6CiAgICAgICAgYyA9IHJvdW5kKGF0IC8gbGVuKHNlcSksIDIpCiAgICAgICAgcmV0dXJuIGMKICAgIGV4Y2VwdCBaZXJvRGl2aXNpb25FcnJvcjoKICAgICAgICByZXR1cm4gMC4wCgpkZWYgR0Moc2VxKToKICAgIGdjID0gc3VtKHNlcS5jb3VudChoKSBmb3IgaCBpbiBbImciLCAiYyIsICJHIiwgIkMiXSkKICAgIHRyeToKICAgICAgICB5ID0gcm91bmQoKGdjIC8gbGVuKHNlcSkpKjEwMCwgMikKICAgICAgICByZXR1cm4geQogICAgZXhjZXB0IFplcm9EaXZpc2lvbkVycm9yOgogICAgICAgIHJldHVybiAwLjAKCmRuYV9zZXEgPSAiR1RDVENUQUNDVFRHQUNBR0FDQ1RDQ0FHQ0NHVEFDQVRHQ0dBQ0FHVFRDR1RHR0NUQ0FDQ1RHQ0FHR0FHQUNDQUdDQ0NHQ1RHQUdHR0FUR0NDR1RDR1RDQVRDR0FHQ0FHIgphdF9jb250ZW50ID0gQVQoZG5hX3NlcSkKcHJpbnQoYXRfY29udGVudCkKZ2NfY29udGVudCA9IEdDKGRuYV9zZXEpCnByaW50KGdjX2NvbnRlbnQp