R/annotate_variants.R
annotate_transcript_strand.Rd
Add transcript strand annotation to SBS variants (defined in genes only)
annotate_transcript_strand(musica, genome_build, build_table = TRUE)
musica | A |
---|---|
genome_build | Which genome build to use: hg19, hg38, or a custom TxDb object |
build_table | Automatically build a table from this annotation |
None
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.