You can count the number of occurrences of a particular character in a string with tr:
use 5.010; + my $str1 = 'MALSSTAATTSSKLKLSNPPSLSHTFTASASASVSNSTSFR'; my $a_count = $str1 =~ tr/A/A/; # 6
length will give you the length of the entire string.
Now to actually pull out the uppercase sequence from your sample input, are you reading lines from a file? Something like this would probably work:
#!/usr/bin/env perl use 5.010; for (<>) { if (/^>/) { # Header } elsif (/^[A-Z]+$/) { # Protein my $a = tr/A/A/; say "A: $a, length: " . length; } }
Then simply run it with script.pl < protein.txt. Modify the say ... line to taste, or more likely, replace it with the rest of your logic. You can also choose to parse the header if needed, in the # Header section.
You could of course modify this to actually open the file in your script with open instead, if that is more desirable:
open my $fh, '<', $filename or die "Couldn't open $filename: $!"; for (<$fh>) {
use strict; use warnings; omitted for brevity.
|
---|
Replies are listed 'Best First'. | |
---|---|
Re^2: How to count the length of a sequence of alphabets and number of occurence of a particular alphabet in the sequence?
by davi54 (Sexton) on Oct 07, 2019 at 18:08 UTC | |
by rjt (Curate) on Oct 07, 2019 at 18:30 UTC | |
by stevieb (Canon) on Oct 07, 2019 at 18:15 UTC | |
by davi54 (Sexton) on Oct 11, 2019 at 15:05 UTC | |
by hippo (Bishop) on Oct 11, 2019 at 15:47 UTC | |
by davi54 (Sexton) on Oct 11, 2019 at 23:32 UTC | |
|
In Section
Seekers of Perl Wisdom