--- a/extract_genomic_dna.py	Wed Jan 20 07:37:31 2016 -0500
+++ b/extract_genomic_dna.py	Wed Mar 02 09:13:24 2016 -0500
@@ -17,6 +17,7 @@
 parser.add_argument('--reference_genome_source', dest='reference_genome_source', help="Source of reference genome file")
 parser.add_argument('--reference_genome', dest='reference_genome', help="Reference genome file")
 parser.add_argument('--output_format', dest='output_format', help="Output format")
+parser.add_argument('--description_field_delimiter', dest='description_field_delimiter', default=None, help="Fasta description field delimiter")
 parser.add_argument('--output', dest='output', help="Output dataset")
 args = parser.parse_args()
 
@@ -159,7 +160,8 @@
         if input_is_gff:
             start, end = egdu.convert_bed_coords_to_gff([start, end])
         fields = [args.genome, str(chrom), str(start), str(end), strand]
-        meta_data = "_".join(fields)
+        field_delimiter = extract_genomic_dna_utils.get_description_field_delimiter(args.description_field_delimiter)
+        meta_data = field_delimiter.join(fields)
         if name.strip():
             out.write(">%s %s\n" % (meta_data, name))
         else: