Beefy Boxes and Bandwidth Generously Provided by pair Networks
"be consistent"

Re: Genbank file parsing

by blazar (Canon)
on Jan 11, 2005 at 14:00 UTC ( #421273=note: print w/ replies, xml ) Need Help??

in reply to Genbank file parsing

Dear monks, My question is related to file parsing within while loops. I have a large file containing many individual files all in the same format.
I don't think so. I suppose you have a large file containing info related to many individual files.
This is how the file looks: How can I extract all the info between ORIGIN and // for each record??
# extract of file:  

FEATURES             Location/Qualifiers    
                     /note="blah blah"   
COUNT                200
        1 lots of nice info
       61 lots of nice info
      121 lots of nice info
If you can rely on this format, here's how I'd do it:
#!/usr/bin/perl -ln use strict; use warnings; if ($_ eq 'ORIGIN') { local $/='//'; print <>; }
#!/usr/bin/perl -ln use strict; use warnings; print if $_ eq 'ORIGIN' .. $_ eq '//' __END__
Of course these are intended to be as minimal examples: adapt the techniques shown here to your needs.
This is my code so far:
open (FILE, $ARGV[0]) or die "unable to open FILE\n";
Why not using <> in the first place? Also, you'd better:
  • use lexical FHs,
  • use the three args form of open(),
  • put relevant info in the error message (i.e. at least include $!).
Note: I skipped the rest

Comment on Re: Genbank file parsing
Select or Download Code

Log In?

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

How do I use this? | Other CB clients
Other Users?
Others avoiding work at the Monastery: (6)
As of 2015-10-06 21:44 GMT
Find Nodes?
    Voting Booth?

    Does Humor Belong in Programming?

    Results (161 votes), past polls