Beefy Boxes and Bandwidth Generously Provided by pair Networks
go ahead... be a heretic
 
PerlMonks  

Re: Re: Re: Re: Re: BioInformatics - polyA tail search

by BrowserUk (Pope)
on Sep 02, 2003 at 19:14 UTC ( #288391=note: print w/replies, xml ) Need Help??


in reply to Re: Re: Re: Re: BioInformatics - polyA tail search
in thread BioInformatics - polyA tail search

Something like this will get you started.

perl -nle" print "$ARGV:($./$+[0]): $1" if m[([AN]{10,}]g;" file*

This will print lines like

filename:(10/50): ANNNANANAAN

where the first number is theline in the file and the second is the offset within the line.

For unix you need to swap "s for 's, and under Win32 you would need to add BEGIN{ @ARGV=map{ glob } @ARGV } to expand the wildcard filespec supplied on the comand line.

See perlrun for the switches used, and perlre for the regex.


Examine what is said, not who speaks.
"Efficiency is intelligent laziness." -David Dunham
"When I'm working on a problem, I never think about beauty. I think only how to solve the problem. But when I have finished, if the solution is not beautiful, I know it is wrong." -Richard Buckminster Fuller
If I understand your problem, I can solve it! Of course, the same can be said for you.

Log In?
Username:
Password:

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

How do I use this? | Other CB clients
Other Users?
Others studying the Monastery: (6)
As of 2020-05-31 22:41 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    If programming languages were movie genres, Perl would be:















    Results (177 votes). Check out past polls.

    Notices?