diff --git a/server/call_missing.py b/server/call_missing.py
index a685b36..915cb55 100644
--- a/server/call_missing.py
+++ b/server/call_missing.py
@@ -167,7 +167,14 @@ def report_uncovered(gff_input, transcript_filename, genetests_filename,
     # Note: Start is 1-based, not 0-based as is in transcript files
     examined_regions = {}
 
+    header_done = False
     for record in gff_data:
+        if not header_done:
+            yield "##gff-version " + gff_data.data[0]
+            yield "##genome-build " + gff_data.data[1]
+            yield "# Produced by: call_missing.py"
+            header_done = True
+
         if f_out:
             yield str(record)
 
diff --git a/server/get_metadata.py b/server/get_metadata.py
index 91adbfe..94838bd 100755
--- a/server/get_metadata.py
+++ b/server/get_metadata.py
@@ -113,7 +113,16 @@ def genome_metadata(gff_input, genome_stats_file, progresstracker):
     progresstracker.metadata['chromosomes'] = list()
 
     # Progress through GFF input.
+    header_done = False
     for record in gff_data:
+        # Have to do this after calling the first record to
+        # get the iterator to read through the header data
+        if not header_done:
+            yield "##gff-version " + gff_data.data[0]
+            yield "##genome-build " + gff_data.data[1]
+            yield "# Produced by: get_metadata.py"
+            header_done = True
+
         # Record number of positions called.
         dist = (record.end - (record.start - 1))
         called_num += dist