diff --git a/mutalyzer/mapping.py b/mutalyzer/mapping.py index e0ae81adfa79088b804bf1add7d66520ce74ed93..78356ef9c83285c5b6c2b89637693d0c1607e7a6 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: