bac2b49a9cf03181c21ba22ab9ce4696f5735160 max Fri Jan 9 05:42:08 2026 -0800 in the SHANK2 locus, the Japan ToMMV VCF files has more than 1024 infos in a repeat, so making the buffer bigger, variant frequencies track diff --git src/lib/vcf.c src/lib/vcf.c index e4d30ca1a05..db148b4225f 100644 --- src/lib/vcf.c +++ src/lib/vcf.c @@ -531,31 +531,31 @@ vcff->headerString = dyStringCannibalize(&dyHeader); return vcff; } struct vcfFile *vcfFileFromHeader(char *name, char *headerString, int maxErr) /* Parse the VCF header string into a vcfFile object with no rows. * name is for error reporting. * If maxErr is non-negative then continue to parse until maxErr+1 errors have been found. * A maxErr less than zero does not stop and reports all errors. * Set maxErr to VCF_IGNORE_ERRS for silence. */ { struct lineFile *lf = lineFileOnString(name, TRUE, cloneString(headerString)); return vcfFileHeaderFromLineFile(lf, maxErr); } -#define VCF_MAX_INFO (4*1024) +#define VCF_MAX_INFO (4*4096) static void parseRefAndAlt(struct vcfFile *vcff, struct vcfRecord *record, char *ref, char *alt) /* Make an array of alleles, ref first, from the REF and comma-sep'd ALT columns. * Use the length of the reference sequence to set record->chromEnd. * Note: this trashes the alt argument, since this is expected to be its last use. */ { char *altAlleles[VCF_MAX_INFO]; int altCount = chopCommas(alt, altAlleles); record->alleleCount = 1 + altCount; record->alleles = vcfFileAlloc(vcff, record->alleleCount * sizeof(record->alleles[0])); record->alleles[0] = vcfFilePooledStr(vcff, ref); int i; for (i = 0; i < altCount; i++) record->alleles[1+i] = vcfFilePooledStr(vcff, altAlleles[i]); int refLen = strlen(ref);