From 1a6adddacbd13ec3d17b019a454aae15dfa637ef Mon Sep 17 00:00:00 2001 From: Vignan Thumma <126392049+vignan8@users.noreply.github.com> Date: Sat, 25 Mar 2023 12:10:34 +0530 Subject: [PATCH] Update opt.py --- opt.py | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/opt.py b/opt.py index c909e4d..42b8046 100755 --- a/opt.py +++ b/opt.py @@ -3,10 +3,14 @@ from vaccine.load import dat as vaccine from corona import corona +# Replace thymine (T) with uracil (U) in virus sequence virus = virus.replace("T", "U") + +# Replace pseudouridine (Ψ) with uracil (U) in vaccine sequence vaccine = vaccine.replace("Ψ", "U") """ +# Find the start position of the vaccine sequence in the virus sequence for i in range(len(virus)-len(vaccine)): mm = virus[i:i+len(vaccine)] mr = sum([c1 == c2 for c1,c2 in zip(mm, vaccine)]) / len(vaccine)