Beefy Boxes and Bandwidth Generously Provided by pair Networks
good chemistry is complicated,
and a little bit messy -LW

Re: BioInformatics - polyA tail search

by johndageek (Hermit)
on Sep 03, 2003 at 22:38 UTC ( #288756=note: print w/replies, xml ) Need Help??

in reply to BioInformatics - polyA tail search

This may help. Not polished or thoroughly tested, but should get you to a starting point.


#!/usr/bin/perl -w use strict; ## set minimum length for polytail to search for my $polytail_min_len = 10; my $polytail_curr_len = $polytail_min_len; my $workline; my $polya_str_base; ## read data record for (<DATA>) { ## remove exteranious characters s/[\r\n]//g; # add new data to end of working data string $workline = $workline . $_; while (length($workline) > $polytail_curr_len) { $polya_str_base = substr($workline,0,$polytail_curr_len); ## remove desired characters from string $polya_str_base =~ s/[AN]//g; ## no characters left = all characters were in desired character set if (length($polya_str_base) == 0) { ## add a character from data set to string to test (ok - bump subscrip +t) $polytail_curr_len++; } else { ## a polytail of at least minimum length was found if ($polytail_curr_len > $polytail_min_len) { print substr($workline,0,$polytail_curr_len-1) . "\n"; ## trim characters of found polytail from working string $workline = substr($workline,$polytail_curr_len); ## reset length of string to test $polytail_curr_len = $polytail_min_len; } ## trim lead character from working field $workline = substr($workline,1); } } } __DATA__ 123456789abcdefghijklmnopqrstuvwxyz ACGGAAAAAAAAAAATCGGATCTGAATGTCTAGAGGGGTTCTCTCCCTTGGTGTGAGTCTAGCCCTGAAA +GTTGCANANAN NANANAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA +AAAAAAAAAAAA NTCCCACTCTANCGCGAAGCCAAATTTGTCGAGAGTACTCTGGGGGGAAGAGATCAGAATTGTGCAGACT +AATCCGTAACTGC CAAGTACTATTGGCCCTGTTCCAACCATCTAACCTCCTTATGATAACCATGCCACTAAATGGGTTCCTGG +ATCTGCACCT CATTCGCTCGCCTTATGGCCTCGGCTCTCTGCGTATCCACCCTCCTCGTCACCGCCATGCCCTTCGACCT +TCAGCGGGGG

Log In?

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

How do I use this? | Other CB clients
Other Users?
Others chanting in the Monastery: (6)
As of 2020-05-29 21:05 GMT
Find Nodes?
    Voting Booth?
    If programming languages were movie genres, Perl would be:

    Results (171 votes). Check out past polls.