open (FH,"<prosite.dat"); while (<FH>) { if (/^PA (.+)\.$/) { $c++; print "$c\t$1\n"; $pattern=$1; $pattern =~ s/-//g; $pattern =~ s/\{/[^/g; $pattern =~ s/\}/]/g; $pattern =~ s/\(/\{/g; $pattern =~ s/\)/\}/g; $pattern =~ s/x/[ARNDCEQGHILKMFPSTWYV]/g; print "$c\t$pattern\n"; last if $c > 2; } }
You must be logged in to post a comment.