Skip to content
Snippets Groups Projects
Commit e79c7b35 authored by Andrew Lyon's avatar Andrew Lyon
Browse files

Actually correctly handle completely ordered sequences in run_iupred

parent a9c15e8a
Branches
No related merge requests found
......@@ -321,6 +321,7 @@ def run_iupred(fasta, iupred_path="~/bin/IUPred/"):
diso_pos.append(int(line.split()[0]))
if not diso_pos:
regions.append([])
continue
diso_indices = np.array(diso_pos) - 1
is_non_continuous = np.diff(np.append(diso_indices[0]-999, diso_indices)) != 1
slice_start_indices = diso_indices[is_non_continuous]
......
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