#!/usr/bin/perl

# Use category=1 for all.

$dir = "/ddn/gs1/home/grimmsa/ruifang/HiC/fromBaylor/Neva_peaks_and_blocks";
@diets = ("LF", "HF", "HFLF");

foreach $diet (@diets) {
  open(IN, "$dir/$diet\_blocks.txt");
  open(TMP, ">$diet.domain.tmp");
  while(<IN>) {
    chomp $_; @ar = split/\t/, $_;
    next if ($ar[1] eq "x1");
    print TMP "$ar[0]\t$ar[1]\t$ar[2]\t1\n";
  }
  close(IN); close(TMP);
  system "sort -k1,1 -k2,2n -k3,3n $diet.domain.tmp > $diet.domain.bed";
  system "rm $diet.domain.tmp";
  system "/ddn/gs1/home/grimmsa/tools/tabix-0.2.6/bgzip $diet.domain.bed";
  system "/ddn/gs1/home/grimmsa/tools/tabix-0.2.6/tabix -p bed $diet.domain.bed.gz";
}