harmonize_seqlengths: Forces seqlengths in GRanges to be that in header of bam_file
In jrboyd/seqqc:

Description Usage Arguments Value Examples

This is important to avoid errors while using the gr to query the bam_file.

1	harmonize_seqlengths(gr, bam_file)

`gr`	GRanges
`bam_file`	bam_file

GRanges with seqlengths matching bam_file

peak_files = dir(system.file("extdata", package = "seqqc"), pattern = "Peak$", full.names = TRUE)
peak_grs = seqsetvis::easyLoad_narrowPeak(peak_files)
query_gr = resize(seqsetvis::ssvOverlapIntervalSets(peak_grs), 6e2, fix = "center")

bam_files = dir(system.file("extdata", package = "seqqc"), pattern = "^M.+bam$", full.names = TRUE)

query_gr.harmonized = seqqc:::harmonize_seqlengths(query_gr, bam_files[1])
seqlengths(query_gr.harmonized)