txSeqsFromFa: Get transcript sequence from a GrangesList and a faFile or...

Description Usage Arguments Details Value See Also

View source: R/ranges_helpers.R

Description

For each GRanges object, find the sequence of it from faFile or BSgenome.

Usage

1
txSeqsFromFa(grl, faFile, is.sorted = FALSE, keep.names = TRUE)

Arguments

grl

a GRangesList object

faFile

FaFile, BSgenome, fasta/index file path or an ORFik experiment. This file is usually used to find the transcript sequences from some GRangesList.

is.sorted

a speedup, if you know the grl ranges are sorted

keep.names

a logical, default (TRUE), if FALSE: return as character vector without names.

Details

A wrapper around extractTranscriptSeqs that works for ORFik experiment input. For debug of errors do: which(!(unique(seqnamesPerGroup(grl, FALSE)) This happens usually when the grl contains chromsomes that the fasta file does not have. A normal error is that mitocondrial chromosome is called MT vs chrM even though they have same seqlevelsStyle. The above line will give you which chromosome it is missing.

Value

a DNAStringSet of the transcript sequences

See Also

Other ExtendGenomicRanges: asTX(), coveragePerTiling(), extendLeaders(), extendTrailers(), reduceKeepAttr(), tile1(), windowPerGroup()


ORFik documentation built on March 27, 2021, 6 p.m.