Skip to content

Commit

Permalink
biodb: orth-ss fix strand bug
Browse files Browse the repository at this point in the history
  • Loading branch information
Weigang Qiu committed Oct 29, 2023
1 parent 4ebca0d commit 591b993
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion bin/biodb
Original file line number Diff line number Diff line change
Expand Up @@ -172,7 +172,7 @@ sub export_orth_bbss {
# my $id = $orf->{'locus'} . "|" . $orf->{'strain'};
my $id = "Bb_" . $orf->{'strain'};
my $seq = Bio::Seq->new(-id => $id, -seq => $orf->{'seq'});
$seq = $seq->revcom() if !$orf->{'strand'};
$seq = $seq->revcom() if !$ref->{'strand'}; # revcom if B31 ORF is on the opposite strand
$out->write_seq($seq);
}

Expand Down

0 comments on commit 591b993

Please sign in to comment.