Commit c0043541 authored by Yixiong Sun's avatar Yixiong Sun
Browse files

Fixed aiming for in sequences longer than 200 nucleotides. The positions are...

Fixed aiming for in sequences longer than 200 nucleotides. The positions are not always sorted so slicing can result in incorrect sequences
parent 47f94660
......@@ -397,8 +397,13 @@ def run_validation(module_number, BNs, module_sequences, target_sequences,test_i
new_seq += s
elif pdb_len > 200:
# NEED to sort aiming_for or else we will have incorrect pos
aiming_for = sorted(aiming_for)
first = aiming_for[0]
last = aiming_for[-1]
offset = max(0, first - 80)
if last - first < 300:
new_seq = seq[max(0, first - 80):last + 80]
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment