Convert a splice junction ID into a GRanges object
junc_to_gr(junc_id)
A GRanges
object
junc_to_gr("chr5:10-50:+")
#> GRanges object with 1 range and 0 metadata columns:
#> seqnames ranges strand
#> <Rle> <IRanges> <Rle>
#> [1] chr5 10-50 +
#> -------
#> seqinfo: 1 sequence from an unspecified genome; no seqlengths
junc_id <- c("chr5:10-50:+", "chr_special1:5-20:-")
junc_to_gr(junc_id)
#> GRanges object with 2 ranges and 0 metadata columns:
#> seqnames ranges strand
#> <Rle> <IRanges> <Rle>
#> [1] chr5 10-50 +
#> [2] chr_special1 5-20 -
#> -------
#> seqinfo: 2 sequences from an unspecified genome; no seqlengths