src/hg/pslCDnaGenomeMatch/pslCDnaGenomeMatch.c 1.24
1.24 2009/08/18 23:24:05 baertsch
fix counter
Index: src/hg/pslCDnaGenomeMatch/pslCDnaGenomeMatch.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/pslCDnaGenomeMatch/pslCDnaGenomeMatch.c,v
retrieving revision 1.23
retrieving revision 1.24
diff -b -B -U 4 -r1.23 -r1.24
--- src/hg/pslCDnaGenomeMatch/pslCDnaGenomeMatch.c 18 Aug 2009 22:31:32 -0000 1.23
+++ src/hg/pslCDnaGenomeMatch/pslCDnaGenomeMatch.c 18 Aug 2009 23:24:05 -0000 1.24
@@ -788,8 +788,9 @@
name, score, nextBestScore, maxScore);
if (maxScore == score)
{
maxHits ++;
+ nextBestScore = score;
verbose(3,"%s score %d == maxScore %d next %d \n",
name, score, maxScore, nextBestScore);
}
else
@@ -1329,8 +1330,9 @@
if (psl->strand[0] == '+' && psl->strand[1] == '-')
pslRc(psl);
pslTabOut(psl, outFile);
outputCount++;
+ filterCount++;
return;
}
/* one loci for each place the mrna alignments to the genome */
for (align = alignList ; align != NULL ; align= align->next)
@@ -1445,9 +1447,9 @@
doOneMrna(lastName, subList);
alignFreeList(&subList);
//pslFreeList(&pslList);
verbose(1,"Wrote %d alignments out of %d \n", outputCount, aliCount);
-verbose(1,"Filtered %d out of %d mRNAs \n", filterCount, mrnaCount);
+verbose(1,"Kept %d out of %d mRNAs \n", filterCount, mrnaCount);
if (computeSS)
{
fprintf(outFile, "%8d %8d %8d %8d %8d %8d %8d %8d %8d %8d %8d %8d %8d %8d %8d %8d \n",
getHistogram('a','a'), getHistogram('a','c'), getHistogram('a','g'), getHistogram('a','t'),