Is your feature request related to a problem? Please describe.
parse.grl("chr1:36568608-36569851:+", seqlengths = seqlengths(hg19))
would give the output of
GRangesList object of length 1:
[[1]]
GRanges object with 1 range and 0 metadata columns:
seqnames ranges strand
<Rle> <IRanges> <Rle>
[1] chr1 1-249250621 +
-------
seqinfo: 84 sequences from an unspecified genome
due to the trailing strand info.
Describe the solution you'd like
Detect if there are 3 components after splitting the string by ":", and if the last element is length one of "+", "-", then parse that part as strand value.
Describe alternatives you've considered
Or throw an error if don't want to handle it.
Is your feature request related to a problem? Please describe.
would give the output of
due to the trailing strand info.
Describe the solution you'd like
Detect if there are 3 components after splitting the string by ":", and if the last element is length one of "+", "-", then parse that part as strand value.
Describe alternatives you've considered
Or throw an error if don't want to handle it.