sub Passref {
my ( $define, $seq, $count, $outfh ) = @_;
chomp($define);
my $len = length($$seq);
my %count = ( Miss => 0+( $$seq =~ m/^YCWLPHQIMTNKSRVADEGF/ig ) );
# YCWLPHQIMTNKSRVADEGF doesn't appear in test file
my %ratio = ( Miss => 0+( $count{Miss} / $len ) );
for my $letter (qw[ E D A V R S K N T M I Q H P L W C Y F ]) {
$count{$letter} = ( $$seq =~ m/$letter/g );
$ratio{$letter} = $count{$letter} / $len;
}
print {$outfh} join "\t", $define, $len,
@ratio{qw[ E D A V R S K N T M I Q H P L W C Y F Miss ]}, "\n";
}
|