/*
* 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.statham2014mods.processed;
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 InspectCrossover {
public static void main(String[] args) throws Exception { //debug
File dir= Dirs.getFile("dir");
File d= new File(dir+"/other_datasets/statham2014_mcf7-mods/locate");
File d2= new File(dir+"/data/human_magnus/h2az/norm/peak/locate/ensGene_hg19_exon");
File tarf= new File(d+
"/GSM1383850_MCF7_NDR_intergenic.bed"
);
File reff= new File(d2+"/peak_norm_wig_mcf7_chip-h2az_time-0_hg19"+
"/locate_intergenic.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();
}
}