diff --git a/src/ecopcr.c b/src/ecopcr.c index 8892e1b..5f17e20 100644 --- a/src/ecopcr.c +++ b/src/ecopcr.c @@ -138,7 +138,7 @@ void printRepeat(ecoseq_t *seq, int32_t error1; int32_t error2; - int32_t ldeta,rdelta; + int32_t ldelta,rdelta; char *amplifia = NULL; int32_t amplength; @@ -204,8 +204,8 @@ void printRepeat(ecoseq_t *seq, superkingdom_name = "###"; } - ldeta=(pos1 <= delta)?pos1:delta; - rdelta(pos2+delta>=seqlength)?seqlength-pos2-1:delta; + ldelta=(pos1 <= delta)?pos1:delta; + rdelta=(pos2+delta>=seqlength)?seqlength-pos2-1:delta; amplifia = getSubSequence(seq->SQ,pos1-ldelta,pos2+rdelta); amplength= strlen(amplifia)-rdelta-ldelta; @@ -257,7 +257,7 @@ void printRepeat(ecoseq_t *seq, amplifia[amplength - o2->patlen - o1->patlen]=0; else { - delta=ldelta+rdelta+amplength-1; + delta=ldelta+rdelta+amplength; for (i=0;i