$seq_file=$ARGV[0];
open (IN,"$seq_file") || die "Cannot open $seq_file";
$fs=0;
$ss=0;
$tss=0;
while(<IN>){
 chomp;
 @arr=split(/\s+/,$_);
 $fs += $arr[2];
 $ss += $arr[3];
 $tss += $arr[8]; 
 }
close(IN);
open(OUT,">>settingsc") || die "Cannot open settingsc";
print OUT "FS = $fs\n";
print OUT "SS = $ss\n";
print OUT "Total_Silent_Sites = $tss\n";
close(OUT);

