Skip to content

Commit

Permalink
Bug: was printing opposite strand to truth. whoops.
Browse files Browse the repository at this point in the history
  • Loading branch information
Torsten Seemann committed Sep 28, 2013
1 parent 664d62e commit 8c2cca2
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion barrnap
Original file line number Diff line number Diff line change
Expand Up @@ -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

Expand Down

0 comments on commit 8c2cca2

Please sign in to comment.