#!usr/bin/perl 
use strict;
use warnings;

if (! @ARGV) {
    print <<HELP;
Usage:

> basecount.pl <bases file>

HELP
    exit;
}

open my $dnaIn, '<', $ARGV[0] or die "Can't open bases file $ARGV[0]: $!\n";

my %counts;
my @baseList = qw(A T G C);

while (defined (my $line = <$dnaIn>)) {
    chomp $line;
    ++$counts{$_} for grep {/\S/} split '', $line;
}

my $bases;
my $errors;

$bases += $_ for @counts{@baseList};
$errors += $_ for map {$counts{$_}} grep {! /[ATGC]/} keys %counts;

print "Total bases: $bases\n";
print join (', ', map {"$_: $counts{$_}"} @baseList), "\n";
print "Errors: $errors\n" if $errors;