#!/usr/bin/perl -w use strict; chomp(my @sequences = ); my (%sequences, $i); push @{$sequences{$_}}, $i++ for @sequences; for (keys %sequences) { my $complement = reverse; $complement =~ tr/GCAT/CGTA/; print "$_ complements sequence nr @{$sequences{$complement}}\n" if $sequences{$complement}; } __DATA__ GATTACA GATTACA TTAGCAA CATCATC TTGCTAA TGTAATC GATGATG