#!/usr/bin/perl
use strict;
use warnings;
use 5.014;
my @prot = qw/ D K /;
my ($acid_cnt, $base_cnt, $neutral_cnt);
while(@prot)
{
my $aa = shift (@prot);
if($aa =~/[DNEQ]/gi)
{
++$acid_cnt;
}
if($aa=~/[KRH]/gi)
{
++$base_cnt;
}
say "$aa pos: ", pos($aa) // 'pos reset';
if($aa=~/[DNEQKRH]/gi)
{
++$neutral_cnt;
}
}
####
C:\Old_Data\perlp>perl t7.pl
aa D pos: pos reset
aa K pos: 1
####
C:\Old_Data\perlp>perl t9.pl test.fas
Name: >DROME_HH_Q02937
Number of acidic amino acids:33
Number of basic amino acids:35
Number of neutral amino acids:136
Name: >DROME_HH_Q02938
Number of acidic amino acids:18
Number of basic amino acids:17
Number of neutral amino acids:69
Name: >DROTME_HH_Q02936
Number of acidic amino acids:14
Number of basic amino acids:18
Number of neutral amino acids:67