def shmir_from_fasta(siRNA, offtarget, regexp, original_frames, prefix):
siRNA2 = reverse_complement(siRNA)
frames = adjusted_frames(siRNA, siRNA2, 0, 0, deepcopy(original_frames)) # we do not have shifts here
shmirs = [frame.template() for frame in frames]
with allow_join_result():
foldings = group(fold.s(shmir, prefix=prefix).set(queue="subtasks") for shmir in shmirs).apply_async().get()
results = []
iter_frames = izip(frames, original_frames, foldings)
for frame, original_frame, folding in iter_frames:
score = score_from_transcript(frame, original_frame, folding["ss"], offtarget, regexp)
if validate_transcript_by_score(score):
results.append({"score": score, "frame": frame, "folding": folding, "found_sequence": siRNA})
return results
def test_validate_transcript_by_score_true(self):
score = {'structure': 70, 'all': 101}
self.assertTrue(validators.validate_transcript_by_score(score))
def test_validate_transcript_by_score_false(self):
score = {'structure': 70, 'all': 80}
self.assertFalse(validators.validate_transcript_by_score(score))
