INNER CODE UNIT · Python
nucl_to_add
HadrienG/InSilicoSeq · iss/error_models/__init__.py:147
nucl_to_add = str(full_sequence[read_end + i])
mut_seq.append(nucl_to_add)
elif orientation == "reverse":
for i in range(to_add):
if read_start - 1 - i < 0:
nucl_to_add = "A"
else:
nucl_to_add = util.rev_comp(full_sequence[read_start - 1 - i])
mut_seq.append(nucl_to_add)
return Seq(mut_seq)
def introduce_indels(self, record, orientation, full_seq, bounds):
"""Introduce insertions or deletions in a sequence
Introduce insertion and deletion errors according to the probabilities
present in the indel choices list