--- a/xmfa2gff3.py	Wed Aug 17 14:46:55 2016 -0400
+++ b/xmfa2gff3.py	Fri Sep 07 11:30:11 2018 -0400
@@ -1,12 +1,17 @@
 #!/usr/bin/env python
+import argparse
+import logging
 import sys
+
+from BCBio import GFF
 from Bio import SeqIO
 from Bio.Seq import Seq
+from Bio.SeqFeature import (
+    FeatureLocation,
+    SeqFeature
+)
 from Bio.SeqRecord import SeqRecord
-from Bio.SeqFeature import SeqFeature, FeatureLocation
-import argparse
-from BCBio import GFF
-import logging
+
 logging.basicConfig(level=logging.INFO)
 log = logging.getLogger(__name__)
 
@@ -131,6 +136,10 @@
                 # Ignore 0% identity sequences
                 if pid == 0:
                     continue
+
+                # Support for Biopython 1.68 and above, which removed sub_features
+                if not hasattr(other['feature'], "sub_features"):
+                    other['feature'].sub_features = []
                 other['feature'].sub_features.append(
                     SeqFeature(
                         FeatureLocation(real_start, real_end),