Commit b0d4e925 authored by Michele Clamp's avatar Michele Clamp
Browse files

Fixed strand bug in GenscanPeptide

parent 37bd46b6
......@@ -408,7 +408,7 @@ sub _make_homol {
my $newh = new Bio::EnsEMBL::SeqFeature(-start => $g1,
-end => $g2,
-strand => 1,
-strand => $pep_exon->strand,
);
$newh->primary_tag($pep_homol->primary_tag);
......@@ -483,7 +483,7 @@ sub _make_dna_homol {
my $newh = new Bio::EnsEMBL::SeqFeature(-start => $g1,
-end => $g2,
-strand => 1);
-strand => $pep_exon->strand * $pep_homol2->strand);
$newh->primary_tag($pep_homol2->primary_tag);
$newh->source_tag ($pep_homol2->source_tag);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment