Re^4: Perl regular expression for amino acid sequence


We don't bite newbies here... much
	PerlMonks

Re^4: Perl regular expression for amino acid sequence

by Roy Johnson (Monsignor)

on Dec 01, 2004 at 20:53 UTC ( [id://411565]=note: print w/replies, xml )

Need Help??

in reply to Re^3: Perl regular expression for amino acid sequence
in thread Perl regular expression for amino acid sequence

The fix is something like:

my $cur;
while ($seq{$k} =~ /([QGYN]{3,6})/g) {
   $cur = $1;
   pos($seq{$k}) -= length($cur);
   $cur =~ s/(.)\1\1.*/$1$1/;
   if (length($cur) >= 3) {
      pos($seq{$k}) += length($cur);
   }
   else { ++pos($seq{$k}); next }
   print "\n$k";
   print $cur." begins at position ", (pos($seq{$k})-length($s)) , "\n
+";
}
[download]

Caution: Contents may have been coded under pressure.