open(promotor, "H:\\bioinfo\\epd.txt"); while (<promotor>) { $c++; if (/^LOCUS(.+)(EP[0-9]+)/) { while ($seq =~ /(A[GC]GAC[ATCG]{2}TC)/g) {$positie=pos($seq)-length($1)+1; $tel++; print "$positie\t$1\t$tel\n$seq\n"};print ">$2\n";$seq=""} if (/^ /){ s/ //g; tr/actg/ACTG/; s/[0-9]//g; s/\n//g; $seq=$seq.$_ } #print; #last if $c > 50; }
You must be logged in to post a comment.