Add transcript strand annotation to SBS variants (defined in genes only)

annotate_transcript_strand(musica, genome_build, build_table = TRUE)

Arguments

musica

A musica object.

genome_build

Which genome build to use: hg19, hg38, or a custom TxDb object

build_table

Automatically build a table from this annotation

Value

None

Examples

data(musica)
annotate_transcript_strand(musica, 19)
#>   403 genes were dropped because they have exons located on both strands
#>   of the same reference sequence or on more than one reference sequence,
#>   so cannot be represented by a single genomic range.
#>   Use 'single.strand.genes.only=FALSE' to get all the genes in a
#>   GRangesList object, or use suppressMessages() to suppress this message.