c9b7830e130a43e817566cedfb1be819cc2164c5 galt Thu Feb 19 17:42:22 2026 -0800 fix oops do not want duplicate usage line. refs #28109 diff --git src/hg/utils/bigChainToChain/bigChainToChain.c src/hg/utils/bigChainToChain/bigChainToChain.c index 12764304f2d..caea7e87a4c 100644 --- src/hg/utils/bigChainToChain/bigChainToChain.c +++ src/hg/utils/bigChainToChain/bigChainToChain.c @@ -1,157 +1,154 @@ /* bigChainToChain - convert bigChain files back into a chain file. */ #include "common.h" #include "linefile.h" #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, chromHash, clBed, f, &processChromChunk); } else if (clPos != NULL) { genericBigToNonBigFromPos(bbi, chromHash, clPos, f, &processChromChunk); } 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 || 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; }