open( IN, "<".$ARGV[0] );
open( OUT, ">".$ARGV[1] );

$L = <IN>;
print OUT "Gene\tNonzero_Counts\tMean_Counts\n";

while( $L = <IN> ) {
    chomp( $L );
    @D = split( /\t/, $L );

    $sum = 0;
    $total = 0;
    $nonzero = 0;

    for( $i = 1; $i < scalar @D; $i++ ) {
	$sum += $D[$i];
	$total++;
	if( $D[$i] > 0 ) {
	    $nonzero++;
	}
    }

    print OUT $D[0]."\t".$nonzero."\t".( $sum / $total )."\n";
}

close( IN );
close( OUT );