From 7a02f64351913d2692c6231fc0f454e6037246dc Mon Sep 17 00:00:00 2001
From: Remi Planel <rplanel@pasteur.fr>
Date: Fri, 3 Apr 2020 17:48:09 +0200
Subject: [PATCH] Fix bug shift one nt for guide
---
crisprbact/predict.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/crisprbact/predict.py b/crisprbact/predict.py
index 1b3821e..f0e09b8 100644
--- a/crisprbact/predict.py
+++ b/crisprbact/predict.py
@@ -35,13 +35,13 @@ def find_targets(seq):
for target in matching_targets:
matching_target = target.group(1)
start, end = target.span(1)
- start_min = 13
+ start_min = 14
if start >= start_min:
guide_start = start - start_min
- guide_end = start + 7
+ guide_end = end - 16 - 3
guide = seq_revcomp[guide_start:guide_end]
assert len(guide) == 20
- pos_seq_start = L - guide_start - 20
+ pos_seq_start = L - guide_end
pos_seq_stop = L - guide_start
pos_seq_pam = pos_seq_start - 3
yield dict(
--
GitLab