Hi, I have a list of composite regions (meaning each region is comprised of one or more non-overlapping/disjoint continuous intervals). Biologically, each number represents the genome-relative coordinates of a single nucleotide position of a transcript spanning one or more exon junctions. The numbers look something like this (for ease of reading, i've bolded the nucleotides at the exon junctions):

121328495,121328496,

121328497,121330673,121330674,121330675,121330676,121330677,121330678,121330991"

As you can see, there are gaps/skips in the region.

My goal is to use permutation test for the significance p-value compared to random of the overlap of these composite regions against a set of standard reference genomic ranges. So far, I've been thoroughly enjoying the regioneR package (by the way, it is SO useful. thank God it exists) and the permutation test in regioneR does exactly what I want to do, except it only seems to be able to do permutations of *continuous intervals*.

I was wondering if there's a way to still use regioneR to solve the problem? I can imagine it would be trivial if there were a way to "lock" the distances between **some** but **not all** intervals together - mutually locked intervals would represent composite regions, which could then be run through standard permtest.

However, if this is not possible, then is there any other way to accomplish this goal?

Regards, Angel