a44421a79fb36cc2036fe116b97ea3bc9590cd0c braney Fri Dec 2 09:34:39 2011 -0800 removed rcsid (#295) diff --git src/utils/bigBedInfo/bigBedInfo.c src/utils/bigBedInfo/bigBedInfo.c index 3d0fc66..4007b48 100644 --- src/utils/bigBedInfo/bigBedInfo.c +++ src/utils/bigBedInfo/bigBedInfo.c @@ -1,101 +1,100 @@ /* bigBedInfo - Show information about a bigBed file.. */ #include "common.h" #include "linefile.h" #include "hash.h" #include "options.h" #include "udc.h" #include "bbiFile.h" #include "bigBed.h" #include "obscure.h" #include "hmmstats.h" -static char const rcsid[] = "$Id: bigBedInfo.c,v 1.2 2009/11/13 19:02:38 kent Exp $"; void usage() /* Explain usage and exit. */ { errAbort( "bigBedInfo - Show information about a bigBed file.\n" "usage:\n" " bigBedInfo file.bb\n" "options:\n" " -udcDir=/dir/to/cache - place to put cache for remote bigBed/bigWigs\n" " -chroms - list all chromosomes and their sizes\n" " -zooms - list all zoom levels and theier sizes\n" " -as - get autoSql spec\n" ); } static struct optionSpec options[] = { {"udcDir", OPTION_STRING}, {"chroms", OPTION_BOOLEAN}, {"zooms", OPTION_BOOLEAN}, {"as", OPTION_BOOLEAN}, {NULL, 0}, }; void printLabelAndLongNumber(char *label, long long l) /* Print label: 1,234,567 format number */ { printf("%s: ", label); printLongWithCommas(stdout, l); printf("\n"); } void bigBedInfo(char *fileName) /* bigBedInfo - Show information about a bigBed file.. */ { struct bbiFile *bbi = bigBedFileOpen(fileName); printf("version: %d\n", bbi->version); printf("isCompressed: %s\n", (bbi->uncompressBufSize > 0 ? "yes" : "no")); printf("isSwapped: %d\n", bbi->isSwapped); printLabelAndLongNumber("itemCount", bigBedItemCount(bbi)); printLabelAndLongNumber("primaryDataSize", bbi->unzoomedIndexOffset - bbi->unzoomedDataOffset); if (bbi->levelList != NULL) { long long indexEnd = bbi->levelList->dataOffset; printLabelAndLongNumber("primaryIndexSize", indexEnd - bbi->unzoomedIndexOffset); } printf("zoomLevels: %d\n", bbi->zoomLevels); if (optionExists("zooms")) { struct bbiZoomLevel *zoom; for (zoom = bbi->levelList; zoom != NULL; zoom = zoom->next) printf("\t%d\t%d\n", zoom->reductionLevel, (int)(zoom->indexOffset - zoom->dataOffset)); } struct bbiChromInfo *chrom, *chromList = bbiChromList(bbi); printf("chromCount: %d\n", slCount(chromList)); if (optionExists("chroms")) for (chrom=chromList; chrom != NULL; chrom = chrom->next) printf("\t%s %d %d\n", chrom->name, chrom->id, chrom->size); if (optionExists("as")) { char *asText = bigBedAutoSqlText(bbi); if (asText == NULL) printf("as: n/a\n"); else { printf("as:\n"); printf("%s", asText); } } struct bbiSummaryElement sum = bbiTotalSummary(bbi); printLabelAndLongNumber("basesCovered", sum.validCount); printf("meanDepth (of bases covered): %f\n", sum.sumData/sum.validCount); printf("minDepth: %f\n", sum.minVal); printf("maxDepth: %f\n", sum.maxVal); printf("std of depth: %f\n", calcStdFromSums(sum.sumData, sum.sumSquares, sum.validCount)); } int main(int argc, char *argv[]) /* Process command line. */ { optionInit(&argc, argv, options); if (argc != 2) usage(); udcSetDefaultDir(optionVal("udcDir", udcDefaultDir())); bigBedInfo(argv[1]); if (verboseLevel() > 1) printVmPeak(); return 0; }