extract_sequence_between_boundaries: Extract sequence between two samples boundaries

In adnaniazi/tailfindr: Poly(A) Tail Length Estimator for Oxford Nanopore RNA and DNA Reads

Extract sequence between two samples boundaries

Description

The functions takes two sample numbers – start and end – as inputs and finds their corresponding closest matching events in the event table. It then finds the FASTA sequence between these bounding events.

Usage

extract_sequence_between_boundaries(event_data, start, end)

Arguments

event_data	a data frame or tibble. The event table of a read. Must contain start and model_state columns
start	a numeric. Starting sample number
end	a numeric. Ending sample number

Value

A character string is returned containing FastQ sequencing

Examples

## Not run: 
event_data <- data.frame(start=seq(0, 150, 5),
                         model_state=sample(c("A", "T", "C", "G"), 1))
sequence <- extract_sequence_between_boundaries(event_data, 98, 131)

## End(Not run)

adnaniazi/tailfindr documentation built on March 23, 2024, 7:07 a.m.