/*
* To change this license header, choose License Headers in Project Properties.
* To change this template file, choose Tools | Templates
* and open the template in the editor.
*/
package data.human.test.enhancer;
import fork.lib.base.file.management.Dirs;
import fork.lib.bio.anno.genomic.BedReader;
import fork.lib.bio.anno.genomic.region.GenomicRegionsBuilder;
import fork.lib.bio.anno.genomic.region.GenomicRegionsComparator;
import fork.lib.bio.anno.genomic.region.GenomicRegionsCrossover;
import java.io.File;
/**
*
* @author man-mqbpjmg4
*/
public class Crossover {
public static void main(String[] args) throws Exception { //debug
File dir= Dirs.getFile("dir");
File tarf= new File(dir+"/data\\human_magnus\\h2az\\norm\\peak\\locate\\ensGene_hg19_exon\\peak_norm_wig_mcf7_chip-h2az_time-0_hg19/"+
"locate_exon.bed");
File reff= new File(dir+"/anno/bed/hg19/enhancers/human_sorted.bed");
GenomicRegionsBuilder gbe= new BedReader(tarf).getGenomicRegionsBuilder();
GenomicRegionsBuilder gbc= new BedReader(reff).getGenomicRegionsBuilder();
System.out.println(gbe.size());
System.out.println(gbc.size());
GenomicRegionsCrossover gc= new GenomicRegionsComparator(gbe, gbc).getCrossover();
gc.print();
}
}