48cc0aa29826d6f6e99357a12eda728d76a7d132 markd Fri May 17 02:31:33 2024 -0700 fixed error in pslSpliceJunctions TSV header diff --git src/utils/pslSpliceJunctions/pslSpliceJunctions.c src/utils/pslSpliceJunctions/pslSpliceJunctions.c index 8e63e6e..c20db7a 100644 --- src/utils/pslSpliceJunctions/pslSpliceJunctions.c +++ src/utils/pslSpliceJunctions/pslSpliceJunctions.c @@ -113,31 +113,31 @@ { processPslGap(transPsl, iBlk, tGapIdx, genomeFh, juncsTsvFh); tGapIdx++; } } } static void pslSpliceJunctions(char *pslFile, char *genome2bitFile, char *junctionsTsv) /* pslSpliceJunctions - Extract splice junctions for a PSL file. */ { struct psl *transPsls = pslLoadAll(pslFile); slSort(&transPsls, pslCmpTarget); struct twoBitFile *genomeFh = twoBitOpen(genome2bitFile); FILE *juncsTsvFh = mustOpen(junctionsTsv, "w"); -fprintf(juncsTsvFh, "qName\tqStart\tqEnd\tqSize\ttName\ttStart\tEnd\tSize\tstrand\tiBlk\ttGapIdx\ttGapSize\tjunc5p\tjunc3p\tisCannon\n"); +fprintf(juncsTsvFh, "qName\tqStart\tqEnd\tqSize\ttName\ttStart\ttEnd\ttSize\tstrand\tiBlk\ttGapIdx\ttGapSize\tjunc5p\tjunc3p\tisCannon\n"); for (struct psl *transPsl = transPsls; transPsl != NULL; transPsl = transPsl->next) processPsl(transPsl, genomeFh, juncsTsvFh); carefulClose(&juncsTsvFh); twoBitClose(&genomeFh); } int main(int argc, char *argv[]) /* Process command line. */ { optionInit(&argc, argv, options); if (argc != 4) usage(); pslSpliceJunctions(argv[1], argv[2], argv[3]); return 0;