From bacb094028126e44eaacbceb0576e136bbd77d10 Mon Sep 17 00:00:00 2001 From: Auden Cote-L'Heureux <52716489+AudenCote@users.noreply.github.com> Date: Mon, 27 Nov 2023 14:11:07 -0500 Subject: [PATCH] Cleaning up length filter --- PTL1/Transcriptomes/Scripts/6_FilterPartials.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PTL1/Transcriptomes/Scripts/6_FilterPartials.py b/PTL1/Transcriptomes/Scripts/6_FilterPartials.py index f1248de..808e407 100644 --- a/PTL1/Transcriptomes/Scripts/6_FilterPartials.py +++ b/PTL1/Transcriptomes/Scripts/6_FilterPartials.py @@ -304,7 +304,7 @@ def filter_NTD_data(args, OGLenDB): og_prefix = rec.split(og_number)[-2][-4:] og = og_prefix + og_number - if len(starting_AA_seqs[rec]) <= 1.5*OGLenDB[og] and len(starting_AA_seqs[rec]) >= 0.33*OGLenDB[og] and '-'.join(rec.split('_')[1:]) not in short_from_translation: + if len(starting_AA_seqs[rec]) <= 1.5*OGLenDB[og] and 3*len(starting_AA_seqs[rec]) >= OGLenDB[og] and '-'.join(rec.split('_')[1:]) not in short_from_translation: good_NTD_seqs.append((rec, starting_NTD_seqs[rec])) good_AA_seqs.append((rec, starting_AA_seqs[rec]))