UNIX for Dummies Questions & Answers

#!/usr/bin/perl
print "Enter minimum coverage: ";
chomp($min_cov=<STDIN>);
open I, "$ARGV[0]";
while (<I>){
  $pos=((split "[\t ]+",$_)[4]);
  $cov=((split "=",(split ";",(split "[\t ]+",$_)[8])[2])[1]);
  if ($cov >= $min_cov){
    $s{$pos}+=$cov;
    $n{$pos}++;
    $s_tot+=$cov;
    $n_tot++;
  }
}
END{
  for $i (keys %s){print "Mean_Coverage_for_position _$i = " . ($s{$i}/$n{$i}) . "\n"};
  print "Mean_Coverage = " . ($s_tot/$n_tot) . "\n"
}