675293a87172d440c9c78b1c1bdc54610bc35a77 braney Tue May 7 14:14:26 2024 -0700 changes to HPRC rearrangements track in response to SAB refs #33484 diff --git src/hg/utils/chainArrangeCollect/chainArrangeCollect.c src/hg/utils/chainArrangeCollect/chainArrangeCollect.c index d68e035..76edb77 100644 --- src/hg/utils/chainArrangeCollect/chainArrangeCollect.c +++ src/hg/utils/chainArrangeCollect/chainArrangeCollect.c @@ -32,31 +32,31 @@ static int count = 0; struct slName *names = hashSlNameFromHash(nameHash); int sizeQuery = bed->score; bed->score = slCount(names); struct dyString *dy = newDyString(100); for(; names; names = names->next) { dyStringAppend(dy, names->name); if (names->next) dyStringAppend(dy, ","); } bed->name = dy->string; // we're actually not outputting chainArrange structure because the label is coming // from an external program currently -fprintf(f, "%s %d %d %s%d.1 %d + %d %d 0 %s %d\n", bed->chrom, bed->chromStart, bed->chromEnd, prefix, count++, bed->score, bed->chromStart, bed->chromEnd, bed->name, sizeQuery); +fprintf(f, "%s\t%d\t%d\t%s%d.1\t%d\t+\t%d\t%d\t0\t%s\t%d\n", bed->chrom, bed->chromStart, bed->chromEnd, prefix, count++, bed->score, bed->chromStart, bed->chromEnd, bed->name, sizeQuery); } void chainArrangeCollect(char *prefix, char *inFile, char *outFile) /* chainArrangeCollect - collect overlapping beds into a single bed. */ { struct bed *allBeds = bedLoadAll(inFile); FILE *f = mustOpen(outFile, "w"); struct bed *bed, *prevBed = allBeds; prevBed->score = 1; struct hash *nameHash = newHash(0); hashStore(nameHash, prevBed->name); if (exact) { for(bed = prevBed->next; bed; bed = bed->next)