[Bioperl-l] While running without stop
    Dave Messina 
    David.Messina at sbc.su.se
       
    Thu May 19 19:17:01 UTC 2011
    
    
  
I'm not familiar with these newfangled @-s and @+s, so I can't help with
that.
But the traditional way to do this is with pos():
while ( my $seq = $in->next_seq() ) {
    my $seqstring = $seq->seq;
    while ( $seqstring =~ m/(ATCGA)/g ) {
        my $end = pos($seqstring);
        my $start = $end - length($&) + 1;
        print "ATCGA commence à la position $start\nse termine juste avant
la position $end\n";
    }
}
Dave
    
    
More information about the Bioperl-l
mailing list