examples/Code/03-complement-file.pl

#!/usr/bin/perl -w
# Calculating the reverse complement of a strand of DNA

# read lines from file or STDIN
while ( $DNA = <> ) {

        # remove line ending
        chomp( $DNA );

        # Make a new (reverse) copy of the DNA
        $revcom = reverse $DNA;

        # Translate A->T, C->G, G->C, T->A
        $revcom =~ tr/ACGT/TGCA/;

        # Print the reverse complement DNA onto the screen
        print "$revcom\n";
}
1	dpavlin	45	#!/usr/bin/perl -w
2			# Calculating the reverse complement of a strand of DNA
3
4			# read lines from file or STDIN
5			while ( $DNA = <> ) {
6
7			# remove line ending
8			chomp( $DNA );
9
10			# Make a new (reverse) copy of the DNA
11			$revcom = reverse $DNA;
12
13			# Translate A->T, C->G, G->C, T->A
14			$revcom =~ tr/ACGT/TGCA/;
15
16			# Print the reverse complement DNA onto the screen
17			print "$revcom\n";
18			}