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

Updating header and Hook name check CheckSetup.py

Browse Source
This commit is contained in:
Auden Cote-L'Heureux 2024-01-26 11:46:06 -05:00 committed by GitHub
parent df54de3589
commit bf669310fd
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 12 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
PTL1/Genomes/Scripts/CheckSetup.py
View File

@ -1,3 +1,10 @@
# Last updated Nov 2023
# Author: Auden Cote-L'Heureux
# This script is run as the first step of the PhyloToL 6 Part 1 GENOMES pipeline,
# before any sequence data are actually processed. It checks to ensure that the input
# CDS files and databases are properly located and formatted.
import os, sys, re
from Bio import SeqIO
@ -21,6 +28,11 @@ def check_databases(params):
fasta = [file for file in os.listdir(params.databases + '/db_OG') if file.endswith('.fasta')]
dmnd = [file for file in os.listdir(params.databases + '/db_OG') if file.endswith('.dmnd')]
if len(fasta) == 1 and len(dmnd) == 1:
if fasta[0].split('.fasta')[0] != dmnd[0].split('.dmnd')[0]:
print('\nERROR: The file names (except for the extensions) of the OG reference (Hook) database .fasta and .dmnd databases must match!\n')
exit()
if len(fasta) == 0:
print('\nERROR: No Hook fasta file found in the Databases/db_OG folder\n')
exit()