diff --git a/src/libecoprimer/ecoprimer.h b/src/libecoprimer/ecoprimer.h index 0a9989a..0110e2c 100644 --- a/src/libecoprimer/ecoprimer.h +++ b/src/libecoprimer/ecoprimer.h @@ -179,6 +179,7 @@ typedef struct { uint32_t mind; //< minimum distance between primers uint32_t maxd; //< maximum distance between primers uint32_t sumd; //< distance sum + uint32_t amplifiacount; float yule; float quorumin; float quorumout; diff --git a/src/libecoprimer/pairs.c b/src/libecoprimer/pairs.c index 64d8905..08e3ac4 100644 --- a/src/libecoprimer/pairs.c +++ b/src/libecoprimer/pairs.c @@ -274,6 +274,7 @@ static void buildPrimerPairsForOneSeq(uint32_t seqid, current.maxd=DMAX; current.mind=DMAX; current.sumd=0; + current.amplifiacount=0; current.inexample=0; current.outexample=0; current.curseqid = 0; @@ -343,6 +344,7 @@ static void buildPrimerPairsForOneSeq(uint32_t seqid, { //pcurrent->inexample++; pcurrent->sumd+=distance; + pcurrent->amplifiacount++; if ((pcurrent->maxd==DMAX) || (distance > pcurrent->maxd)) pcurrent->maxd = distance;