Returns DNA sequences

Source: R/compute.R

Returns DNA sequences for given intervals

gseq.extract(intervals = NULL)

Arguments

intervals

intervals for which DNA sequence is returned

Value

An array of character strings representing DNA sequence.

Details

This function returns an array of sequence strings for each interval from 'intervals'. If intervals contain an additional 'strand' column and its value is '-1', the reverse-complementary sequence is returned.

See also

gextract

Examples


gdb.init_examples()
intervs <- gintervals(c(1, 2), 10000, 10020)
gseq.extract(intervs)
#> [1] "ccaggggccagcactgctcg" "GGGGCTCCCACCCGCCGTCC"