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

Updating headers in 6_FilterPartials.py

Browse Source
This commit is contained in:
Auden Cote-L'Heureux 2025-03-19 09:20:37 -04:00 committed by GitHub
parent d1b9a64e60
commit 8487f1d836
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -5,7 +5,8 @@
# First, all sequences shorter than 33% or longer than 150% the average length of sequences
# from the same OG in the Hook database are removed. Then, for each transcriptomic sample,
# all sequences within an OG are compared at the nucleotide level to the sequence with the
# highest “score” (defined as k-mer coverage multiplied by length). The script should be run
# highest “score” (defined as k-mer coverage multiplied by length) using BLAST, and sequences that
# are 98% identical to the master sequence are removed. The script should be run
# as part of the EukPhylo Part 1 pipeline using the script wrapper.py. It requires that the
# structure of the 'Output' folder be as output by script 5, and that the Databases/db_OG folder
# contains a .fasta file containing all amino acid sequences in the OG reference database (Hook)