From bb08f5a3572854615e9b2cda2548d7535bf074bb Mon Sep 17 00:00:00 2001 From: Jesse Luehrs Date: Mon, 19 Nov 2012 03:31:00 -0600 Subject: another solution --- DBRU.pl | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) create mode 100644 DBRU.pl diff --git a/DBRU.pl b/DBRU.pl new file mode 100644 index 0000000..476dbcb --- /dev/null +++ b/DBRU.pl @@ -0,0 +1,17 @@ +#!/usr/bin/env perl +use strict; +use warnings; +use 5.016; + +use List::MoreUtils 'uniq'; + +chomp(my @k1mers = <>); +@k1mers = uniq map { $_, reverse_complement($_) } @k1mers; +my $k = length($k1mers[0]) - 1; +say for map { '(' . substr($_, 0, $k) . ', ' . substr($_, 1, $k) . ')' } + @k1mers; + +sub reverse_complement { + my ($str) = @_; + return scalar reverse $str =~ tr/ACGT/TGCA/r; +} -- cgit v1.2.3