Perhaps not the best efficient way, but it does what you ask:
use warnings;
use Data::Dumper;
my $fseq = 'CCCCGCGC';
my @nsub = ('CCCCG', 'CCCGC', 'CGCGC');
my @results;
for (1 .. $#nsub)
{
my $current = $nsub[$_];
my $previous = $nsub[$_ -1];
# the "#" is used as a separator.
# put a different character (or sequence of characters)
# if you believe that it could also be in your strings
if ( "$previous#$current" =~ /(\w+)#\1/ )
{
my $found = $1;
printf "%d -> %s (%s) %d -> %s \n",
$_ -1,
$previous,
$found,
$_,
$current;
$current =~ s/^$found/"-" x length($found)/e;
$previous =~ s/$found$/"-" x length($found)/e;
push @results, [ $_ -1, $previous];
push @results, [ $_, $current];
}
else
{
printf "%d -> no overlap\n", $_
}
}
print Data::Dumper->Dump([ \@results], ['result']);
Result: (adjusted)
0 -> CCCCG (CCCG) 1 -> CCCGC
1 -> CCCGC (CGC) 2 -> CGCGC
$result = [
[ 0, 'C----' ],
[ 1, '----C' ],
[ 1, 'CC---' ],
[ 2, '---GC' ]
];
-
Are you posting in the right place? Check out Where do I post X? to know for sure.
-
Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
<code> <a> <b> <big>
<blockquote> <br /> <dd>
<dl> <dt> <em> <font>
<h1> <h2> <h3> <h4>
<h5> <h6> <hr /> <i>
<li> <nbsp> <ol> <p>
<small> <strike> <strong>
<sub> <sup> <table>
<td> <th> <tr> <tt>
<u> <ul>
-
Snippets of code should be wrapped in
<code> tags not
<pre> tags. In fact, <pre>
tags should generally be avoided. If they must
be used, extreme care should be
taken to ensure that their contents do not
have long lines (<70 chars), in order to prevent
horizontal scrolling (and possible janitor
intervention).
-
Want more info? How to link
or How to display code and escape characters
are good places to start.