rangeset

package main

import (
	"fmt"

	"github.com/brentp/go-athenaeum/rangeset"
)

func main() {

	chromLength := 1000
	genes, _ := rangeset.New(0, chromLength)
	genes.SetRange(100, 200)
	genes.SetRange(250, 350)
	genes.SetRange(900, 950)

	regions, _ := rangeset.New(0, chromLength)
	regions.SetRange(150, 250)
	regions.SetRange(500, 925)

	fmt.Println("genes          :", genes.Ranges())
	fmt.Println("regions        :", regions.Ranges())
	fmt.Println("genes ∩ regions:", rangeset.Intersection(genes, regions).Ranges())
	fmt.Println("genes ∪ regions:", rangeset.Union(genes, regions).Ranges())
	fmt.Println("genes - regions:", rangeset.Difference(genes, regions).Ranges())
	fmt.Println("regions - genes:", rangeset.Difference(regions, genes).Ranges())

}

Output:

genes          : [Range(100-200) Range(250-350) Range(900-950)]
regions        : [Range(150-250) Range(500-925)]
genes ∩ regions: [Range(150-200) Range(900-925)]
genes ∪ regions: [Range(100-350) Range(500-950)]
genes - regions: [Range(100-150) Range(250-350) Range(925-950)]
regions - genes: [Range(200-250) Range(500-900)]