#!/usr/bin/perl

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

open(FH, $infile);

my $total = 0;

my $total_first = 0;

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

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