From b0d4e925e5dc3be62c2e55384758e1659835451c Mon Sep 17 00:00:00 2001
From: Michele Clamp <michele@sanger.ac.uk>
Date: Wed, 9 Feb 2000 15:15:53 +0000
Subject: [PATCH] Fixed strand bug in GenscanPeptide
---
modules/Bio/EnsEMBL/Analysis/GenscanPeptide.pm | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/modules/Bio/EnsEMBL/Analysis/GenscanPeptide.pm b/modules/Bio/EnsEMBL/Analysis/GenscanPeptide.pm
index 6f83f97b8b..f5c59564d8 100755
--- a/modules/Bio/EnsEMBL/Analysis/GenscanPeptide.pm
+++ b/modules/Bio/EnsEMBL/Analysis/GenscanPeptide.pm
@@ -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);
--
GitLab