in reply to motif finding
A more "Perl-ish" way of doing the same thing...
The eliminates counting characters one at a time, as in the $i loop in the original, in favor of using pattern matching on the entire character string. I have found that eliminating loop counters wherever possible greatly reduces the number of bugs in my code.use strict; use warnings; use Term::ANSIColor; use autodie; #Program to find motif site in a given protein sequence using files my $motif = "AGGGGG"; open( my $read, "<dna.txt" ); my @e = <$read>; $_ = join( " ", @e ); s/\s+//g; my @c; push @c, pos( ) - length( $motif ) + 1 while /$motif/g; s/$motif/color( 'bold green' ) . $motif . color( 'black' )/eg; print $_, "\n"; print "Number of sites the motif (AGGGGG) is present: ", scalar @c, "\ +n"; print "And the positions in the string are: ", join( ',', @c ), "\n\n" +;
|
---|
Replies are listed 'Best First'. | |
---|---|
Re^2: motif finding
by educated_foo (Vicar) on Jan 31, 2012 at 13:55 UTC | |
by Anonymous Monk on Oct 05, 2013 at 21:07 UTC | |
Re^2: motif finding
by RichardK (Parson) on Jan 31, 2012 at 14:19 UTC | |
by devi (Initiate) on Feb 03, 2012 at 06:08 UTC | |
by RichardK (Parson) on Feb 03, 2012 at 12:47 UTC |
In Section
Seekers of Perl Wisdom