91c925b21b7c9c02bcc035c75d463162cd7e7d97 angie Mon Dec 5 08:08:23 2011 -0800 Oops, forgot to run 'make test' before committing, and overlookeda test program... diff --git src/lib/tests/vcfParseTest.c src/lib/tests/vcfParseTest.c index c414f40..c6a5fd1 100644 --- src/lib/tests/vcfParseTest.c +++ src/lib/tests/vcfParseTest.c @@ -15,31 +15,31 @@ "usage:\n" " vcfParseTest fileOrUrl.vcf.gz seqName start end\n" "\n" "fileOrUrl.vcf.gz needs to have been compressed by tabix, and index file\n" "fileOrUrl.vcf.gz.tbi must exist.\n" ); } static struct optionSpec options[] = { {NULL, 0}, }; void vcfParseTest(char *fileOrUrl, char *seqName, int start, int end) /* vcfParseTest - Parse VCF header and data lines in given position range.. */ { -struct vcfFile *vcff = vcfTabixFileMayOpen(fileOrUrl, seqName, start, end, 100); +struct vcfFile *vcff = vcfTabixFileMayOpen(fileOrUrl, seqName, start, end, 100, -1); if (vcff == NULL) errAbort("Failed to parse \"%s\" and/or its index file \"%s.tbi\"", fileOrUrl, fileOrUrl); int recCount = slCount(vcff->records); printf("Finished parsing \"%s\" items in %s:%d-%d, got %d data rows\n", fileOrUrl, seqName, start+1, end, recCount); if (recCount > 0) printf("First (up to) 100 rows in range:\n"); int i = 0; struct vcfRecord *rec = vcff->records; while (rec != NULL && i < 100) { printf("%s\t%d\t%d\t%s:%s/%s\t%s\n", rec->chrom, rec->chromStart, rec->chromEnd, rec->name, rec->alleles[0], rec->alleles[1], rec->qual); rec = rec->next;