Beefy Boxes and Bandwidth Generously Provided by pair Networks
XP is just a number
 
PerlMonks  

Re^3: Perl regular expression for amino acid sequence

by dws (Chancellor)
on Dec 01, 2004 at 21:13 UTC ( [id://411575]=note: print w/replies, xml ) Need Help??


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

I need to know what the found pattern was and I need to know if it's repeated which is why two regexp wont work I'm afraid.

Two regexes will work just fine. Use the first to do coarse filtering, and the second to filter.

while ($seq{$k} =~ /([QGYN]{3,6})/g) { next if $1 =~ m/QQQ|GGG|YYY|NNN/; print "\n$k"; print $1." begins at position ", (pos($seq{$k})-length($s)) , "\n"; }
This has the benefit of being blindingly obvious about what you're doing.

Oops: ikegami is correct. This is blindingly wrong.

Replies are listed 'Best First'.
Re^4: Perl regular expression for amino acid sequence
by ikegami (Patriarch) on Dec 01, 2004 at 21:25 UTC

    It has the benefit of being blindingly unobvious about what it's not doing. It doesn't detect 'GNN', 'NNGYGY' and 'NNGNN' givn the input 'xxxxxxxGNNNxxxxxxxNNNGYGYxxxxxxxGYGYNNNxxxxxxxNNNGNNNxxxxxxx'

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://411575]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others having a coffee break in the Monastery: (5)
As of 2024-04-25 06:58 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found