obitools/obitools3
2
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases 66 Wiki Activity

ecoPCR: fixed amplicon length computation bug

Browse Source
This commit is contained in:
Celine Mercier
2019-02-06 18:11:20 +01:00
parent 8f18907566
commit 6ca6d27ecb
1 changed files with 3 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/obi_ecopcr.c
View File

@ -1054,8 +1054,7 @@ int obi_ecopcr(const char* i_dms_name,
if (posj > posi) if (posj > posi)
length = posj - posi - o1->patlen - o2->patlen; length = posj - posi - o1->patlen - o2->patlen;
if (posj < posi) if (posj < posi)
//length = posj + apatseq->seqlen - posi - o1->patlen - o2->patlen; // TODO hein???? length = posj + apatseq->seqlen - posi - o1->patlen - o2->patlen;
length = posi - posj - o1->patlen - o2->patlen;
if ((length>0) && // For when primers touch or overlap if ((length>0) && // For when primers touch or overlap
(!min_len || (length >= min_len)) && (!min_len || (length >= min_len)) &&
(!max_len || (length <= max_len))) (!max_len || (length <= max_len)))
@ -1140,11 +1139,9 @@ int obi_ecopcr(const char* i_dms_name,
length = 0; length = 0;
if (posj > posi) if (posj > posi)
//length = posj - posi + 1 - o2->patlen - o1->patlen; /* - o1->patlen : deleted by <EC> (prior to the OBITools3) */ TODO ???? length = posj - posi + 1 - o2->patlen - o1->patlen; /* - o1->patlen : deleted by <EC> (prior to the OBITools3) */
length = posj - posi - o2->patlen - o1->patlen;
if (posj < posi) if (posj < posi)
//length = posj + apatseq->seqlen - posi - o1->patlen - o2->patlen; TODO ???? length = posj + apatseq->seqlen - posi - o1->patlen - o2->patlen;
length = posi - posj - o2->patlen - o1->patlen;
if ((length>0) && // For when primers touch or overlap if ((length>0) && // For when primers touch or overlap
(!min_len || (length >= min_len)) && (!min_len || (length >= min_len)) &&
(!max_len || (length <= max_len))) (!max_len || (length <= max_len)))