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.