Re^2: motif finding


Clear questions and runnable code get the best and fastest answer
	PerlMonks

Re^2: motif finding

by educated_foo (Vicar)

on Jan 31, 2012 at 13:55 UTC ( [id://950975]=note: print w/replies, xml )

Need Help??

in reply to Re: motif finding
in thread motif finding

Or, even more Perl-ish, with a bit less extra work (e.g. only one //g loop):

use Term::ANSIColor;
open(READ,"<dna.txt");
$m = 'AGGGGG';
$_ = do { local $/; <READ> };   # read whole file
s/\s+//g;                       # remove blanks
s{$m}{                          # search the string
    push @c, 1 - length($m) + pos;               # remember position
    color('bold green').$m.color('reset');       # remember to reset!
}eg;
print "$_\n";                   # print transformed string
print "NUMBER OF SITES THE MOTIF ($m) IS PRESENT: ".@c."\n";
print "AND THE POSITION IN THE STRING IS:", join(',', @c), "\n\n";
[download]