https://bioinformaticsonline.com/snippets/view/44222/raku-script-to-find-palindrome-in-genomes? https://bioinformaticsonline.com/snippets/view/44222/raku-script-to-find-palindrome-in-genomes Tue, 07 Mar 2023 14:15:17 -0600 https://bioinformaticsonline.com/snippets/view/44222/raku-script-to-find-palindrome-in-genomes sub is-palindrome(Str $str) returns Bool { $str.=uc; # convert to uppercase $str.=subst:g/\s+//; # remove any spaces return $str eq $str.flip; } sub find-palindromes(Str $dna, Int $min-length, Int $max-length) { for $min-length..$max-length -> $length { for 0..^$dna.chars - $length -> $pos { my $substring = $dna.substr($pos, $length); if is-palindrome($substring) { say "Palindrome found at position $pos: $substring"; } } } } # Example usage my $dna = "GGATCCATGGCCTAGG"; # example DNA sequence find-palindromes($dna, 3, 8); # find palindromes with length between 3 and 8]]> BioStar