From 1c012dc92814916749009f236c895a007e76fc1d Mon Sep 17 00:00:00 2001
From: Martijn Vermaat <martijn@vermaat.name>
Date: Tue, 23 Feb 2016 23:24:52 +0100
Subject: [PATCH] Only run import-reference on mitochondrial refs
---
mutalyzer/mapping.py | 3 +++
1 file changed, 3 insertions(+)
diff --git a/mutalyzer/mapping.py b/mutalyzer/mapping.py
index e0ae81ad..78356ef9 100644
--- a/mutalyzer/mapping.py
+++ b/mutalyzer/mapping.py
@@ -860,6 +860,9 @@ def import_from_reference(assembly, reference):
retriever = Retriever.GenBankRetriever(output)
record = retriever.loadrecord(reference)
+ if record.molType != 'm':
+ raise ValueError('Only mitochondial references are supported')
+
select_transcript = len(record.geneList) > 1
for gene in record.geneList:
--
GitLab