YuuMJ1997/EukPhylo
1
0
Fork 0
mirror of http://43.156.76.180:8026/YuuMJ/EukPhylo.git synced 2025-12-27 21:50:26 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fixing same contig/OG number bug in 6_FilterPartials.py

Browse Source
This commit is contained in:
Auden Cote-L'Heureux 2023-10-06 10:41:29 -04:00 committed by GitHub
parent 694ece4648
commit a515a08c61
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
PTL1/Transcriptomes/Scripts/6_FilterPartials.py
View File

@ -301,7 +301,7 @@ def filter_NTD_data(args, OGLenDB):
good_NTD_seqs = []; good_AA_seqs = [] good_NTD_seqs = []; good_AA_seqs = []
for rec in starting_NTD_seqs: for rec in starting_NTD_seqs:
og_number = re.split('OG.{1}_', rec)[-1][:6] 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 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 len(starting_AA_seqs[rec]) >= 0.33*OGLenDB[og] and '-'.join(rec.split('_')[1:]) not in short_from_translation: