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.