good chemistry is complicated, and a little bit messy -LW |
|
PerlMonks |
Re: simple but stuckby stephen (Priest) |
on Apr 25, 2002 at 15:27 UTC ( [id://161995]=note: print w/replies, xml ) | Need Help?? |
If I'm reading what you're trying to do correctly, you're trying to read sequences from a FastA format file, count the length of each sequence, and print out the ones that are more than 250 bases long. If possible, I'd suggest using the bioperl module. That does all of the work of reading fasta format, and so will make your work easier. If you can't install modules yourself, perhaps you can get your system administrator to install it for you. Here's an example of what I think you're trying to do, using the BioPerl modules.
stephen Update: Removed redundant length() call.
In Section
Seekers of Perl Wisdom
|
|