From a515a08c6187c6fbf384da794fbc9a0189d89c53 Mon Sep 17 00:00:00 2001 From: Auden Cote-L'Heureux <52716489+AudenCote@users.noreply.github.com> Date: Fri, 6 Oct 2023 10:41:29 -0400 Subject: [PATCH] Fixing same contig/OG number bug in 6_FilterPartials.py --- 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 d243472..f1248de 100644 --- a/PTL1/Transcriptomes/Scripts/6_FilterPartials.py +++ b/PTL1/Transcriptomes/Scripts/6_FilterPartials.py @@ -301,7 +301,7 @@ def filter_NTD_data(args, OGLenDB): good_NTD_seqs = []; good_AA_seqs = [] for rec in starting_NTD_seqs: og_number = re.split('OG.{1}_', rec)[-1][:6] - og_prefix = rec.split(og_number)[0][-4:] + 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: