| harmonize_seqlengths {seqsetvis} | R Documentation |
ensures compatibility between seqlength of gr and bam_file based on header
harmonize_seqlengths(gr, bam_file)
gr |
GRanges, object to harmonize seqlengths for |
bam_file |
character, a path to a valid bam file |
gr with seqlengths matching bam_file
library(GenomicRanges)
gr = GRanges("chr1", IRanges(1, 100))
#seqlengths has not been set
seqlengths(gr)
bam = system.file("extdata/test.bam", package = "seqsetvis")
gr2 = harmonize_seqlengths(gr, bam)
#seqlengths now set
seqlengths(gr2)