#!/usr/bin/perl

my $infile = shift;
if(!( -f $infile)) { die "$infile not a valid file" }

open(FH, $infile);

my $total = 0;
my $pos = 0;

my $total_first = 0;
my $pos_first = 0;

while(<FH>) {
	if(/^\s*#/) { next }
	my($agi) = /(AT[1-5MC]G\d{5}\.\d+)/i;
	my ($num) = $agi =~ /\.(\d+)/;	
	$total++;
	if($num==1) { $total_first++ }
	if(/\s+Y\t+AT/) { 
		$pos++;
		if($num==1) { $pos_first++ }
	}
}

print "Total Gene Models: $total\n";
print "Total SigPositive: $pos\n";
print "Distinct Locii: $total_first\n";
print "Distinct SigPositive: $pos_first\n";