summaryrefslogtreecommitdiffstats
path: root/DBRU.pl
diff options
context:
space:
mode:
Diffstat (limited to 'DBRU.pl')
-rw-r--r--DBRU.pl17
1 files changed, 17 insertions, 0 deletions
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;
+}