diff --git a/jcvi/formats/genbank.py b/jcvi/formats/genbank.py
index <HASH>..<HASH> 100644
--- a/jcvi/formats/genbank.py
+++ b/jcvi/formats/genbank.py
@@ -49,11 +49,6 @@ class MultiGenBank(BaseFile):
             for f in rf:
                 self.print_gffline(gff_fw, f, seqid)
                 nfeats += 1
-
-                for sf in f.sub_features:
-                    self.print_gffline(gff_fw, sf, seqid, parent=f)
-                    nfeats += 1
-
             nrecs += 1
 
         logging.debug(
Remove .sub_features extraction as no longer needed (thanks @chen1i6c<I>, #<I>) (#<I>)