diff --git a/mutalyzer/describe.py b/mutalyzer/describe.py
index 555894661475b0ed114568f06dbf4d72b50e6c90..c76ff4c4811ba3a2a9f4c92abd0dcc8766721b88 100644
--- a/mutalyzer/describe.py
+++ b/mutalyzer/describe.py
@@ -624,6 +624,8 @@ def describe_dna(s1, s2):
         #if
     #for
 
+    if not description:
+        return Allele([DNAVar()])
     return description
 #describe_dna
 
diff --git a/tests/test_describe.py b/tests/test_describe.py
index f1ccb6180a6c62aedb9a3b63dc7444febbf953a3..080f3716d35955df13c659cdd77b5e3ca93b0bc2 100644
--- a/tests/test_describe.py
+++ b/tests/test_describe.py
@@ -25,7 +25,7 @@ class TestDescribe(MutalyzerTest):
         """
         reference = "ACGTCGATTCGCTAGCTTCGGGGGATAGATAGAGATATAGAGAT"
 
-        result = describe.describe(reference, sample)
+        result = describe.describe_dna(reference, sample)
         assert result[0].type == expected[0]
         assert result[0].start == expected[1]
         assert result[0].end == expected[2]
@@ -39,41 +39,37 @@ class TestDescribe(MutalyzerTest):
         """
         Test 1.
         """
-        result = describe.describe(
+        result = describe.describe_dna(
             'ATGATGATCAGATACAGTGTGATACAGGTAGTTAGACAA',
             'ATGATTTGATCAGATACATGTGATACCGGTAGTTAGGACAA')
-        description = describe.allele_description(result)
-        assert description == '[5_6insTT;17del;26A>C;35dup]'
+        assert unicode(result) == '[5_6insTT;17del;26A>C;35dup]'
 
     def test2(self):
         """
         Test 2.
         """
-        result = describe.describe(
+        result = describe.describe_dna(
             'TAAGCACCAGGAGTCCATGAAGAAGATGGCTCCTGCCATGGAATCCCCTACTCTACTGTG',
             'TAAGCACCAGGAGTCCATGAAGAAGCTGGATCCTCCCATGGAATCCCCTACTCTACTGTG')
-        description = describe.allele_description(result)
-        assert description == '[26A>C;30C>A;35G>C]'
+        assert unicode(result) == '[26A>C;30C>A;35G>C]'
 
     def test3(self):
         """
         Test 3.
         """
-        result = describe.describe(
+        result = describe.describe_dna(
             'TAAGCACCAGGAGTCCATGAAGAAGATGGCTCCTGCCATGGAATCCCCTACTCTA',
             'TAAGCACCAGGAGTCCATGAAGAAGCCATGTCCTGCCATGGAATCCCCTACTCTA')
-        description = describe.allele_description(result)
-        assert description == '[26_29inv;30C>G]'
+        assert unicode(result) == '[26_29inv;30C>G]'
 
     def test4(self):
         """
         Test 4.
         """
-        result = describe.describe(
+        result = describe.describe_dna(
             'TAAGCACCAGGAGTCCATGAAGAAGATGGCTCCTGCCATGGAATCCCCTACTCTA',
             'TAAGCACCAGGAGTCCATGAAGAAGCCATGTCCTGCCATGAATCCCCTACTCTA')
-        description = describe.allele_description(result)
-        assert description == '[26_29inv;30C>G;41del]'
+        assert unicode(result) == '[26_29inv;30C>G;41del]'
 
     def test5(self):
         """