From 6cb294bd922c2cb888d59ab26bf988f283bb11e3 Mon Sep 17 00:00:00 2001
From: Magali Ruffier <mr6@ebi.ac.uk>
Date: Mon, 4 Dec 2017 11:55:39 +0000
Subject: [PATCH] ENSCORESW-2517: only compare on same strand
---
misc-scripts/xref_mapping/XrefParser/RefSeqCoordinateParser.pm | 1 +
1 file changed, 1 insertion(+)
diff --git a/misc-scripts/xref_mapping/XrefParser/RefSeqCoordinateParser.pm b/misc-scripts/xref_mapping/XrefParser/RefSeqCoordinateParser.pm
index 1573606f19..b723c8dba1 100644
--- a/misc-scripts/xref_mapping/XrefParser/RefSeqCoordinateParser.pm
+++ b/misc-scripts/xref_mapping/XrefParser/RefSeqCoordinateParser.pm
@@ -255,6 +255,7 @@ sub run_script {
# Create a range registry for all the exons of the ensembl transcript
foreach my $transcript(@$transcripts) {
+ if ($transcript->strand != $transcript_of->strand) { next; }
my $exons = $transcript->get_all_Exons();
my $rr2 = Bio::EnsEMBL::Mapper::RangeRegistry->new();
my $rr4 = Bio::EnsEMBL::Mapper::RangeRegistry->new();
--
GitLab