04b879361227c44e9602945464b83e5346b5dd4f
markd
  Tue Jul 19 22:38:07 2022 -0700
add option to output autoSql style header to bigBedNamedItems

diff --git src/utils/bigBedNamedItems/bigBedNamedItems.c src/utils/bigBedNamedItems/bigBedNamedItems.c
index 36827d1..9704bbb 100644
--- src/utils/bigBedNamedItems/bigBedNamedItems.c
+++ src/utils/bigBedNamedItems/bigBedNamedItems.c
@@ -1,72 +1,80 @@
 /* bigBedNamedItems - Extract item(s) of given name(s) from bigBed. */
 
 /* Copyright (C) 2013 The Regents of the University of California 
  * See kent/LICENSE or http://genome.ucsc.edu/license/ for licensing information. */
 #include "common.h"
 #include "linefile.h"
 #include "hash.h"
 #include "options.h"
 #include "localmem.h"
 #include "udc.h"
 #include "bPlusTree.h"
 #include "bigBed.h"
 #include "obscure.h"
 #include "zlibFace.h"
+#include "bigBedCmdSupport.h"
 
 char *field = "name";
+boolean header = FALSE;
 
 void usage()
 /* Explain usage and exit. */
 {
 errAbort(
   "bigBedNamedItems - Extract item of given name from bigBed\n"
   "usage:\n"
   "   bigBedNamedItems file.bb name output.bed\n"
   "options:\n"
   "   -nameFile - if set, treat name parameter as file full of space delimited names\n"
   "   -field=fieldName - use index on field name, default is \"name\"\n"
+  "   -header - output a autoSql-style header (starts with '#').\n"
   );
 }
 
 /* Command line validation table. */
 static struct optionSpec options[] = {
    {"nameFile", OPTION_BOOLEAN},
    {"field", OPTION_STRING},
+   {"header", OPTION_BOOLEAN},
    {NULL, 0},
 };
 
 void bigBedNamedItems(char *bigBedFile, char *name, char *outFile)
 /* bigBedNamedItems - Extract item(s) of given name(s) from bigBed. */
 {
 struct bbiFile *bbi = bigBedFileOpen(bigBedFile);
 FILE *f = mustOpen(outFile, "w");
+if (header)
+    bigBedCmdOutputHeader(bbi, f);
+
 struct lm *lm = lmInit(0);
 struct bigBedInterval *intervalList;
 int fieldIx;
 struct bptFile *bpt = bigBedOpenExtraIndex(bbi, field, &fieldIx);
 if (optionExists("nameFile"))
     {
     int wordCount = 0;
     char **words;
     char *buf;
     readAllWords(name, &words, &wordCount, &buf);
     intervalList = bigBedMultiNameQuery(bbi, bpt, fieldIx, words, wordCount, lm);
     }
 else
     {
     intervalList = bigBedNameQuery(bbi, bpt, fieldIx, name, lm);
     }
 bigBedIntervalListToBedFile(bbi, intervalList, f);
 carefulClose(&f);
 }
 
 int main(int argc, char *argv[])
 /* Process command line. */
 {
 optionInit(&argc, argv, options);
 if (argc != 4)
     usage();
 field = optionVal("field", field);
+header = optionExists("header");
 bigBedNamedItems(argv[1], argv[2], argv[3]);
 return 0;
 }