Skip to content

Commit 128e0cb

Browse files
Change ali_len to amino length for length adjustment of non-framshift hits
The length adjustment in p7_pli_postDomainDef_nonFrameshift was incorrectly using an ali_len in nucleotide positions. This has been fixed by dividing the length by 3.
1 parent 088a2a6 commit 128e0cb

File tree

1 file changed

+3
-3
lines changed

1 file changed

+3
-3
lines changed

src/p7_pipeline.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -2220,7 +2220,7 @@ p7_pli_postDomainDef_Frameshift(P7_PIPELINE *pli, P7_FS_PROFILE *gm_fs, P7_BG *b
22202220
ali_len = dom->jali - dom->iali + 1;
22212221
bitscore = dom->envsc;
22222222

2223-
if (ali_len < 8)
2223+
if (ali_len < 12)
22242224
{// anything less than this is a funny byproduct of the Forward score passing a very low threshold, but no reliable alignment existing that supports it
22252225
p7_alidisplay_Destroy(dom->ad);
22262226
p7_trace_fs_Destroy(dom->tr);
@@ -2379,8 +2379,8 @@ p7_pli_postDomainDef_nonFrameshift(P7_PIPELINE *pli, P7_OPROFILE *om, P7_BG *bg,
23792379
dom = pli->ddef->dcl + d;
23802380

23812381
env_len = dom->jenv - dom->ienv + 1;
2382-
ali_len = dom->jali - dom->iali + 1;
2383-
if (ali_len < 8)
2382+
ali_len = (dom->jali - dom->iali + 1) / 3;
2383+
if (ali_len < 4)
23842384
{ // anything less than this is a funny byproduct of the Forward score passing a very low threshold, but no reliable alignment existing that supports it
23852385
p7_alidisplay_Destroy(dom->ad);
23862386
p7_trace_fs_Destroy(dom->tr);

0 commit comments

Comments
 (0)