#!/usr/bin/perl

@A_WGS = ("15");
#@A_WGS = ("20", "25", "30");
@A_GERP = ("nogerp", "gerp0", "gerp1", "gerp2");
@A_COLORS = (
"102,102,0",
"153,153,30",
"204,0,0",
"255,0,0",
"255,0,204",
"255,204,204",
"255,153,0",
"255,204,0",
"255,255,0",
"204,255,0",
"0,255,0",
"53,128,0",
"0,0,204",
"102,153,255",
"153,204,255",
"0,255,255",
"204,255,255",
"153,0,204",
"204,51,255",
"204,153,255",
"102,102,102",
"153,153,153",
"204,204,204");

foreach $WGS (@A_WGS) {
	foreach $GERP (@A_GERP) {
		@A_FILES = <counts_matched.WGS$WGS.$GERP.maf1pct.*.bigwig>;

		for ($i = 0; $i < @A_FILES; $i++) {
			$FILE = $A_FILES[$i];
			($GRP) = $FILE =~ /maf1pct\.(.*)\.bigwig/;
			$COLOR = $A_COLORS[$i];

			print "track type=bigWig name=WGS$WGS.$GERP.maf1pct.$GRP bigDataUrl=http://pastinenlab.genome.mcgill.ca/hg19/immunoseq/$FILE visibility=2 maxHeightPixels=20 color=$COLOR\n";
		}
	}
}