#!/usr/bin/perl

$rgb = "127,0,255";

$refGenome = "hg38";
$idColumns = 5;

$path = `pwd`; $path =~ s/\s//g;
$path =~ s/\/data\/wade\///;

$ls = `ls -1 H3K27ac.*.bigBed`;
foreach $bbfile (split/\n/, $ls) {
  @ar = split/\./, $bbfile;
  @arX = splice(@ar, 0, $idColumns);
  $id = join(".", @arX);
  print "track type=bigBed name=\"$id peaks\" db=$refGenome visibility=dense color=$rgb bigDataUrl=https://snpinfo.niehs.nih.gov/ucscview/$path/$bbfile\n";
}


#  if ($id =~ /H3K27ac/) { $rgb = "127,0,255"; $vr = 10; }
#  elsif ($id =~ /H3K4me/) { $rgb = "255,128,0"; $vr = 20; }
#  elsif ($id =~ /H3K27me3/) { $rgb = "0,128,255"; $vr = 10; }
#  elsif ($id =~ /TET2/) { $rgb = "0,153,0"; $vr = 10; }
#  elsif ($id =~ /IgG/) { $rgb = "102,102,102"; $vr = 10; }

#H3K27ac.Abcam.1ug.1_100.CaCl2_30-rep1.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Abcam.1ug.1_100.CaCl2_30-rep2.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Invitrogen.1_50.1_100.CaCl2_15-rep1.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Invitrogen.1_50.1_100.CaCl2_15-rep2.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Invitrogen.1_50.1_100.CaCl2_30-rep1.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Invitrogen.1_50.1_100.CaCl2_30-rep2.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Invitrogen.1_50.1_200.CaCl2_30-rep1.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed
#H3K27ac.Invitrogen.1_50.1_200.CaCl2_30-rep2.homer-peakcalls.s500md1000-PE-fdre5.ExclBlklist.bigBed