#!/usr/bin/perl -w

my $Num_Bact_Virus_Chimera=0;
my $Num_Bact_Bact_Chimera=0;
my $Num_Virus_Virus_Chimera=0;
my $outfile='MC38_ChimeraTable.txt';
my $MG_Num=1;

open(OUT, ">$outfile") || die "Can't open outputfile $!\n";

print OUT join("\t",qw (MG_Num Num_Bact_Virus_Chimeras Num_Bact_Bact_Chimera
Num_Virus_Virus_Chimera), "\n") ;

while ($MG_Num <= 72) {

@files=</outputMC38/*.MG$MG_Num.*>;

foreach $file (@files) {

if ($file=qr/HybridViralBactContigsList.txt/){

$Num_Bact_Virus_Chimera = qx/wc -l $file/;
$Num_Bact_Virus_Chimera=$Num_Bact_Virus_Chimera-1;

}

if ($file=qr/HybridOnlyBactContigsList.txt/){

$Num_Bact_Bact_Chimera= qx/wc -l $file/;
$Num_Bact_Bact_Chimera=$Num_Bact_Bact_Chimera-1;
}

if ($file=qr/HybridOnlyVirusContigsList.txt/){

$Num_Virus_Virus_Chimera=qx/wc -l $file/;
$Num_Virus_Virus_Chimera=$Num_Virus_Virus_Chimera-1;
}

print OUT
join("\t",$MG_Num,$Num_Bact_Virus_Chimera,$Num_Bact_Bact_Chimera,$Num_Virus_Virus_Chimera,"\n");


}

}
$MG_Num=$MG_Num+1;