while () { our(@F) = split(/\s+/, $_, 0); push @{$r{join ' ' x 8, @F[0..3]};}, [@F[4, 6]]; sub END { foreach $k (keys %r) { my($x, $y); map {$x += $$_[0]; $y += $$_[1];} @{$r{$k};}; my @g = split(/\s+/,$k); print OUT "$g[0]\t@g[1]\t@g[2]\t@g[3]\t", $x / scalar(@{$r{$k};}), "\t$y\n"; } } } #### 1 111 C T 1 0 6 1 136 G A 1 0 6 1 136 G A 1 0 9 1 244 C CT 1 0 2 1 262 A G 1 0 2 1 268 A C 1 0 2 1 268 A C 1 0 4 1 270 C T 1 0 2 #### map {$x += $$_[0]; $y += $$_[1];} @{$r{$k};};