#!/usr/bin/perl
#Naga_betrapally_6_3_16
use strict;
use warnings;
use Cwd;
use FileHandle;
my $name;                         														# Naming of variables, remembers the name of the file
my $cwd = getcwd();              														#Takes in the current directory of files
my $opendir = "$cwd";           														#Used to open the current directory and remembers the location
opendir (DIR, $opendir) or die $!;          											#Open directory command
my @directory = readdir DIR;     														#Grabbing all the fas files stored in the location
my $scalar = scalar(@directory);                										#Counting the number of files and printing it in the next line
#print "The number of sample files in this folder: $scalar \n";
# for (my $i=0; $i<scalar(@directory); $i++){
# $name = "$directory[$i]";
# print "The file read in: $name\n";
# my $digits = substr $name, 0, 4;                  									#Grabbing the first 4 characters to name the sample run for MEGA
# my $tmp = "$opendir/$directory[$i]";

# }

for(my$i=0;$i<@directory;$i++){
print $directory[$i];
&average($directory[$i]);
}

 #average(@directory);

sub average {
my $header;
my @array_lines;
my $n=0;
my $filename=@_; 																			#Takes in the sequence file from the user input with the year in the header
my @average_array;
my $sum =0;
my $avarage;
	 next if($filename =~ /^\.$/);
     next if($filename =~ /^\.\.$/);
open (FILE,$filename) or die "No input file provided or the input file does not exist in the path, '$_' $!";
print "$filename\n";
while (my $file = <FILE>){
chomp $file;
if ($file =~/^Name/){
$header = $file;
}
else{
$array_lines[$n] = $file;  
chop $array_lines[$n];
$n++;
}
}					
for(my $i=0;$i<@array_lines;$i++){
my @tmp = split(/\t/,($array_lines[$i]));
#print $tmp[3];
push @average_array, $tmp[3];
my $scalar = scalar(@average_array);
#print $scalar;
}	
foreach (@average_array){
	#print "$_\n";
	my $tmp = $_;
	$tmp =~ s/,//;
	#print "$tmp\n";
	$sum = $sum+$tmp;
	}	
#print "The sum of all the mapping is $sum\n";
$avarage = $sum/11;
#print "THe average of the mapping is $avarage\n";
close (FILE);
}