@nucs=("A","C","G","T"); $size=5000; for ($i=0; $i<$size; $i++) { $seqR .= $nucs[int(rand(4))]; } print "Seq($size): $seqR\n"; if (/ATG([ACGT][ACGT][ACGT]){3,5000}(TAA|TAG|TGA)/) { print "This seq. might contain a coding region\n" } else{ print "This sequence most liklely does not contin a coding region\n" }