From 8c2cca22f9c68e72e7050dbc85d47564a86e8071 Mon Sep 17 00:00:00 2001 From: Torsten Seemann Date: Sat, 28 Sep 2013 12:37:59 +1000 Subject: [PATCH] Bug: was printing opposite strand to truth. whoops. --- barrnap | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/barrnap b/barrnap index 4c4b87e..4aa558c 100755 --- a/barrnap +++ b/barrnap @@ -71,7 +71,7 @@ foreach (@hits) { my @x = split ' ', $_; err("bad line in nhmmer output - @x") unless defined $x[6] and $x[6] =~ m/^\d+$/; - my($begin,$end,$strand) = $x[6] < $x[7] ? ($x[6],$x[7],'-') : ($x[7],$x[6],'+'); + my($begin,$end,$strand) = $x[6] < $x[7] ? ($x[6],$x[7],'+') : ($x[7],$x[6],'-'); my($seqid, $gene, $prod) = ($x[0], $x[2], $x[2]); $prod =~ s/_r/ ribosomal /; # convert the short ID to an English string