INNER CODE UNIT · Python

to_add

HadrienG/InSilicoSeq · iss/error_models/__init__.py:141

            to_add = self.read_length - len(mut_seq)
            if orientation == "forward":
                for i in range(to_add):
                    if read_end + i >= len(full_sequence):
                        nucl_to_add = "A"
                    else:
                        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):

View source record →