d9741ac7fe00ba4e2cc38bb8baab8c8b2ae6eec6
galt
  Fri Feb 13 18:23:42 2026 -0800
Updated the following seven utilities bigBedToBed, bigMafToMafj, bigWigToWig, bigWigToBedGraph, bigGenePredToGenePred, bigPslToPsl, bigChainToChain, added support for multiple -range parameters that uses and validates Marks desired bed and position values, in a new callback function in the library bigBedCmdSupport.c. Also added the ability to quickly check allowed chromosomes so users see and error when specifying incorrect chromsome names, and it works with all the filter options. These utils also have a new -skipChromCheck option that users can use if the want to skip the check. These have no tests. I will make tests for these soon. refs #28109, #36925

diff --git src/hg/utils/bigChainToChain/bigChainToChain.c src/hg/utils/bigChainToChain/bigChainToChain.c
index 928e33bf80c..12764304f2d 100644
--- src/hg/utils/bigChainToChain/bigChainToChain.c
+++ src/hg/utils/bigChainToChain/bigChainToChain.c
@@ -4,128 +4,154 @@
 #include "obscure.h"
 #include "hash.h"
 #include "chain.h"
 #include "bigBed.h"
 #include "bigChain.h"
 #include "chainNetDbLoad.h"
 #include "options.h"
 #include "basicBed.h"
 #include "bigBedCmdSupport.h"
 
 char *clChrom = NULL;
 int clStart = -1;
 int clEnd = -1;
 char *clBed = NULL;
 char *clPos = NULL;
+struct slName *clRange = NULL;
+struct hash *chromHash = NULL;
+boolean skipChromCheck = FALSE;
 
 char *inChain = NULL;
 char *inLinks = NULL;
 
 void usage()
 /* Explain usage and exit. */
 {
 errAbort(
   "bigChainToChain - convert bigChain files back into a chain file\n"
   "usage:\n"
   "   bigChainToChain bigChain.bb bigLinks.bb output.chain\n"
   "options:\n"
   "   -chrom=chr1 - if set restrict output to given chromosome\n"
   "   -start=N - if set, restrict output to only that over start\n"
   "   -end=N - if set, restict output to only that under end\n"
+  "   -range=\"chrom start end\" - if set, restrict output to only that within range from start to end. \n"
+  "          This range start is a half-open 0-based coordinate like used in BED files. \n"
+  "   -range=chrom:start-end - if set, restrict output to only that within range from start to end. \n"
+  "          This range start is a 1-based start position. \n"
+  "    Do not use range with chrom, start, and/or end options. \n"
+  "   -range may be specified multiple times for multiple ranges. \n"
   "   -bed=in.bed - restrict output to all regions in a BED file\n"
   "   -positions=in.pos - restrict output to all regions in a position file with 1-based start\n"
   "   -udcDir=/dir/to/cache - place to put cache for remote bigBed/bigWigs\n"
+  "   -skipChromCheck - skip checking chrom name.\n"
   );
 }
 
 /* Command line validation table. */
 static struct optionSpec options[] = {
    {"chrom", OPTION_STRING},
    {"start", OPTION_INT},
    {"end", OPTION_INT},
    {"bed", OPTION_STRING},
+   {"range", OPTION_STRING|OPTION_MULTI},
    {"positions", OPTION_STRING},
    {"udcDir", OPTION_STRING},
+   {"skipChromCheck", OPTION_BOOLEAN},
    {NULL, 0},
 };
 
 static void processChromChunk(struct bbiFile *bbi, char *chrom,
                               int start, int end, char *bedName, FILE *f)
 /* Output one chunk.  Only blocks where start is in the range will be written
  * to avoid outputting a block multiple tines.  */
 {
     char *chromName = chrom;
     struct chain  *chains = chainLoadIdRangeHub(NULL, inChain, inLinks, chromName, start, end, -1);
     chainWriteAll(chains, f);
 }
 
 
 void bigChainToChain(char *out)
 /* bigChainArrange - output a set of rearrangement breakpoints. */
 {
 struct bbiFile *bbi = bigBedFileOpen(inChain);
 FILE *f = mustOpen(out, "w");
+struct bbiChromInfo *chrom, *chromList = bbiChromList(bbi);
+if (!skipChromCheck)
+    chromHash = makeChromHash(chromList);
 
 if (clBed != NULL)
     {
-    genericBigToNonBigFromBed(bbi, clBed, f, &processChromChunk);
-    return;
+    genericBigToNonBigFromBed(bbi, chromHash, clBed, f, &processChromChunk);
     }
-if (clPos != NULL)
+else if (clPos != NULL)
     {
-    genericBigToNonBigFromPos(bbi, clPos, f, &processChromChunk);
-    return;
+    genericBigToNonBigFromPos(bbi, chromHash, clPos, f, &processChromChunk);
     }
-
-struct bbiChromInfo *chrom, *chromList = bbiChromList(bbi);
-
+else if (clRange != NULL)
+    {
+    genericBigToNonBigFromRange(bbi, chromHash, f, clRange, &processChromChunk);
+    }
+else
+    {
+    boolean chromFound;
     for (chrom = chromList; chrom != NULL; chrom = chrom->next)
 	{
 	if (clChrom != NULL && !sameString(clChrom, chrom->name))
 	    continue;
 
+	chromFound = TRUE;
 	int start = 0, end = chrom->size;
 	if (clStart >= 0)
 	    start = clStart;
 	if (clEnd >= 0)
 	    {
 	    end = clEnd;
 	    if (end > chrom->size)
 		end = chrom->size;
 	    }
 	if (start > end)
 	    errAbort("invalid range, start=%d > end=%d", start, end);
 
 	processChromChunk(bbi, chrom->name, start, end, NULL, f);
 	}
+    if (clChrom && !chromFound && !skipChromCheck)
+	errAbort("specified chrom %s not found in bigChain", clChrom);
+    }
+bbiChromInfoFreeList(&chromList);
+carefulClose(&f);
+bbiFileClose(&bbi);
 }
 
 int main(int argc, char *argv[])
 /* Process command line. */
 {
 optionInit(&argc, argv, options);
 if (argc != 4)
     usage();
 clChrom = optionVal("chrom", clChrom);
 clStart = optionInt("start", clStart);
 clEnd = optionInt("end", clEnd);
+clRange =optionMultiVal("range", clRange);
 clBed = optionVal("bed", clBed);
 clPos = optionVal("positions", clPos);
+skipChromCheck = optionExists("skipChromCheck");
 udcSetDefaultDir(optionVal("udcDir", udcDefaultDir()));
 
 if (argc != 3)
     usage();
 
-if ((clBed || clPos) && (clChrom || (clStart >= 0) || (clEnd >= 0)))
-    errAbort("-bed or -positions can not be used with -chrom -start or -end options");
-if (clBed && clPos)
-    errAbort("-bed and -positions can not be used together");
+if ((clBed || clPos || clRange) && (clChrom || (clStart >= 0) || (clEnd >= 0)))
+    errAbort("-bed or -positions or -range can not be used with -chrom -start or -end options");
+if ((clBed && clPos) || (clBed && clRange) || (clPos && clRange))
+    errAbort("-bed, -positions, and -range can not be used together");
 
 inChain = argv[1];
 inLinks = argv[2];
  
 bigChainToChain(argv[3]);
 
 if (verboseLevel() > 1)
     printVmPeak();
 return 0;
 }