more useful options | |
PerlMonks |
Re: BioInformatics - polyA tail searchby halley (Prior) |
on Sep 02, 2003 at 18:54 UTC ( [id://288387]=note: print w/replies, xml ) | Need Help?? |
Does a single file have a single sequence? You said these polyA tails are at (or near?) the end.
If the files have just one sequence and are small enough to fit in memory, slurp into a scalar, strip any irrelevant characters, and the regular expression qr/[AN]{10,}/ would find the tails. If the files are too large, or have more than one sequence, you'll have to work a little harder to fit in memory. Reading a couple lines at a time, searching the tailing line(s) for a matching tail. --
In Section
Seekers of Perl Wisdom
|
|