Fix output of b37 in header (fixes #1121)

There should be a better way to do this, but for now the solution
has been to regenerate the header for the file within each python
module that is returning GFF data.

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)
 

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