Dear Perlmonks, I have not yet had time to delve into BioPerl (although I can't wait!).

Here's a brief introduction to Bioperl for you:

#!/usr/bin/perl -w

use strict;
use Bio::SeqIO;

my @files   = <*.seq>;
my $pattern = '[AN]{10}';

foreach my $file (@files) {
   my $seqio = Bio::SeqIO->new( -file => $file );
   while ( my $seqobj = $seqio->next_seq() ) {
      my $raw_seq = $seqobj->seq();
      if ( $raw_seq =~ /$pattern/o ) {
         print "Found a polyA tail in $file.\n";
      }
   }
}
[download]

This seems to work for me.

If you have any questions about how it works, just ask. I made a few assumptions:

• If there were a sequence file with 4 sequences, two of which had poly-A tails, the program would print:
  

  Found a polyA tail in myfile
  Found a polyA tail in myfile
  [download]

• I'm assuming that you're using some sort of format that Bioperl supports.

Also, is there community interest in an Intro to Bioperl Tutorial? The documentation at doc.bioperl.org seems to be very geared towards developers and other experienced Bioperl users. If there is interest I'd be happy to write it, alhough I might need some help (I'm a long time reader of PM, short time member).

---

---

Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!

Titles consisting of a single word are discouraged, and in most cases are disallowed outright.

Read Where should I post X? if you're not absolutely sure you're posting in the right place.

Please read these before you post! —

• How do I compose an effective node title?
• How do I post a question effectively?
• Markup in the Monastery

Posts may use any of the Perl Monks Approved HTML tags:

a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr

You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)

	For:		Use:
	&		&
	<		&lt;
	>		&gt;
	[		&#91;
	]		&#93;

Link using PerlMonks shortcuts! What shortcuts can I use for linking?

See Writeup Formatting Tips and other pages linked from there for more info.