Beefy Boxes and Bandwidth Generously Provided by pair Networks
Just another Perl shrine

Re^2: simple perl script to trim a FASTA file

by erix (Vicar)
on Apr 09, 2010 at 15:52 UTC ( #833851=note: print w/replies, xml ) Need Help??

in reply to Re: simple perl script to trim a FASTA file
in thread simple perl script to trim a FASTA file

You should indeed not edit the sequence files.

The previous answers (by almut and BrowserUK) just use __DATA__ because that makes it easy to show as compact example. And to use the examples they'll need a little editing.

For instance, in almut's script, replace the line:

while( <DATA> ) {
while( <> ) {

and invoke the script as:

perl sequences.fa > sequences_cleaned.fa

Replacing <DATA> with <> means it will not read the __DATA__ section of the program-file itself, but instead read from STDIN (aka standard input).


Log In?

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

How do I use this? | Other CB clients
Other Users?
Others pondering the Monastery: (4)
As of 2016-10-23 16:28 GMT
Find Nodes?
    Voting Booth?
    How many different varieties (color, size, etc) of socks do you have in your sock drawer?

    Results (301 votes). Check out past polls.