import
random
def
generateDNASequence():
l
=
[
'C'
,
'A'
,
'G'
,
'T'
]
res
=
""
for
i
in
range
(
0
,
40
):
res
=
res
+
random.choice(l)
return
res
def
applyGammaRadiation(dna):
pos
=
random.randint(
1
,
100
)
cdna
=
''
l
=
[
'C'
,
'A'
,
'G'
,
'T'
]
if
(pos >
50
):
changepos
=
random.randint(
0
,
39
)
dl
=
[]
dl[:
0
]
=
dna
ch
=
""
+
dl[changepos]
l.remove(ch)
ms
=
random.choice(l)
cl
=
[]
cl[:
0
]
=
dna
cl[changepos]
=
ms
cdna
=
''.join([
str
(e)
for
e
in
cl])
else
:
cdna
=
dna
return
cdna
def
detectMutation(dna, cdna):
count
=
0
for
x, y
in
zip
(dna, cdna):
if
x
=
=
y:
count
=
count
+
1
else
:
break
return
count
dna
=
generateDNASequence()
print
(dna
+
" (Original DNA)"
)
cdna
=
applyGammaRadiation(dna)
print
(cdna
+
" (DNA after radiation)"
)
count
=
detectMutation(dna, cdna)
if
count
=
=
40
:
print
(
"No Mutation detected"
)
else
:
pos
=
"^"
print
(pos.rjust(count
+
1
))
print
(
"Mutation detected at pos = "
, (count
+
1
))