https://bioinformaticsonline.com/snippets/view/27267/for-a-given-dna-sequence-find-its-rna-transcript-find-its-reverse-complement-and-check-if-the-reverse-complement-contains-a-start-codon? https://bioinformaticsonline.com/snippets/view/27267/for-a-given-dna-sequence-find-its-rna-transcript-find-its-reverse-complement-and-check-if-the-reverse-complement-contains-a-start-codon Tue, 10 May 2016 11:42:30 -0500 https://bioinformaticsonline.com/snippets/view/27267/for-a-given-dna-sequence-find-its-rna-transcript-find-its-reverse-complement-and-check-if-the-reverse-complement-contains-a-start-codon #!/usr/bin/perl use strict; use warnings; my $DNA = "GATTACACAT"; #transcribe DNA to RNA - T changes to U my $RNA = $DNA; $RNA =~ s/T/U/g; print "RNA sequence is $RNA\n"; #find the reverse complement of $DNA using substitution operator #first - reverse the sequence my $rcDNA = reverse($DNA); $rcDNA =~ s/T/A/g; $rcDNA =~ s/A/T/g; $rcDNA =~ s/G/C/g; $rcDNA =~ s/C/G/g; print "Reverse complement of $DNA is $rcDNA\n"; #did it work? #find the reverse complement of $DNA using translation operator #first - reverse the sequence $rcDNA = reverse($DNA); #then - complement the sequence $rcDNA =~ tr/ACGT/TGCA/; #then - print the reverse complement print "Reverse complement of $DNA is $rcDNA\n"; #look for a start codon in the reverse sequence if($rcDNA =~ /ATG/){ print "Start codon found\n"; } else{ print "Start codon not found\n"; }]]> Nishi Singh