Beefy Boxes and Bandwidth Generously Provided by pair Networks
Keep It Simple, Stupid

Re: File::Wrap and Fasta format

by robsv (Curate)
on Mar 25, 2003 at 21:41 UTC ( #245796=note: print w/replies, xml ) Need Help??

in reply to File::Wrap and Fasta format

Here's another shorter regex (although it's not doing any substitution):
... print $fh ">$description\n"; print $fh "$_\n" foreach ($seq =~ /.{1,50}/g); ...

The fastest way I've found of dumping a FASTA sequence (with very large sequences, anyways), is with substr:
... print $fh ">$description\n"; print $fh substr($seq,0,50,'') . "\n" while ($seq); ...
The thing to watch here is that it's destructive - $seq will be empty when it finishes.

- robsv

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://245796]
and all is quiet...

How do I use this? | Other CB clients
Other Users?
Others making s'mores by the fire in the courtyard of the Monastery: (5)
As of 2016-12-03 13:53 GMT
Find Nodes?
    Voting Booth?
    On a regular basis, I'm most likely to spy upon:

    Results (56 votes). Check out past polls.