Skip to contents

Function to update the sequence library by extracting subsequences around the seed match.

Usage

.update_sequence_at(
  sequences,
  seed,
  mismatch_allowed = mismatch_,
  limit_ = limit
)

Arguments

sequences

A DNAStringSet object containing the sequences of the transcripts.

seed

A character string representing the seed sequence.

mismatch_allowed

An integer representing the maximum number of mismatches allowed in the seed match.

limit_

An integer representing the limit for sequence update.

Value

A DNAStringSet object containing the updated sequences and a data frame with the start sites of the seed matches.