#!/usr/bin/perl -w 

@files = qw(
by_rapano_ls021
by_rapano_ls024
by_rapa10_ls022
by_rapa10_ls025
by_rapa60_ls023
by_rapa60_ls026
prp4_ifhi_ls027
prp4_iflo_ls028
);

$tailm="_sacCer3.ih1.xpe.tophat.strandm.cov.trk.bb";
$tailp="_sacCer3.ih1.xpe.tophat.strandp.cov.trk.bb";

use DBI; 

$db_user = "donohue"; 
$db_passwd = ""; 
$database = "DBI:mysql:one_sacCer3"; 
$mdb = DBI->connect($database, $db_user, $db_passwd) || die ("couldn't connect to database"); 

open(IF,"ares_sc3_isoforms.bed");
#track type=bed name='ares saccer3 isoforms'  description='ares saccer3 isoforms'  genome=sacCer3 color=0,0,0  priority=034  visibility=full 
$a=<IF>;
#chrI	87285	87752	SNC1	0	+	87285	87752	0	2	102,252,	0,215,
while($a=<IF>){
@b=split(/\t/,$a);

#chrI	87285	87752	YAL030W	0	+	87285	87752	0	2	102,252,	0,215,
$chrom=$b[0];
$min=$b[1];
$max=$b[2];
$ann=$b[3];
$strand=$b[5];

$pos = "$chrom:$min-$max";
$duphash{$pos}++;
if($duphash{$pos}==1){
print "$ann\t$pos";
for($j=0;$j<8;$j++){

if($strand eq '-'){$tfile = "coverage/".$files[$j].$tailm;} 
else{$tfile = "coverage/".$files[$j].$tailp;} 

system("bigWigSummary $tfile $chrom $min $max 10 -type=mean > tmp.stats");
open(ST,"tmp.stats");
$mean=0;
while($a=<ST>){
chomp($a);
@b=split(/\t/,$a);
for($k=0;$k<10;$k++){if($b[$k] ne 'n/a'){$mean += $b[$k]/10;}}
$mean = abs($mean);
}
print "\t$mean";
close(ST);

}



print "\n";
}
}