[ create a new paste ] login | about

Project: bioinformatics
Link: http://bioinformatics.codepad.org/CQu9XCmJ    [ raw code | output | fork ]

proch - Perl, pasted on Apr 23:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
$quality = ';;;;;;;;;;;9;7;;.7;393333';
$avg = avgqual($quality);
print "$quality
$avg
";

sub avgqual {
	my $string = shift;
	my $sum;
	for ($i=0; $i<length($string); $i++) {
		$q = substr($string, $i, 1);

		$sum+=ord($q) - 33;
	}
	return $sum/length($string);
}


Output:
1
2
;;;;;;;;;;;9;7;;.7;393333
23.4


Create a new paste based on this one


Comments: