Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
def identify(self, options):
"""Identify marker genes in genomes.
Parameters
----------
options : argparse.Namespace
The CLI arguments input by the user.
"""
if options.genome_dir:
check_dir_exists(options.genome_dir)
if options.batchfile:
check_file_exists(options.batchfile)
make_sure_path_exists(options.out_dir)
genomes, tln_tables = self._genomes_to_process(
options.genome_dir, options.batchfile, options.extension)
self.genomes_to_process = genomes
markers = Markers(options.cpus)
markers.identify(genomes,
tln_tables,
options.out_dir,
options.prefix,
options.force)
def classify(self, options):
"""Determine taxonomic classification of genomes.
Parameters
----------
options : argparse.Namespace
The CLI arguments input by the user.
"""
# See ticket #255... perhaps an upstream version/OS issue?
if not hasattr(options, 'pplacer_cpus'):
options.pplacer_cpus = None
check_dir_exists(options.align_dir)
make_sure_path_exists(options.out_dir)
if options.scratch_dir:
make_sure_path_exists(options.scratch_dir)
genomes, _ = self._genomes_to_process(
options.genome_dir, options.batchfile, options.extension)
classify = Classify(options.cpus, options.pplacer_cpus)
classify.run(genomes,
options.align_dir,
options.out_dir,
options.prefix,
options.scratch_dir,
options.recalculate_red,
options.debug,
options.split_tree)
def align(self, options):
"""Create MSA from marker genes.
Parameters
----------
options : argparse.Namespace
The CLI arguments input by the user.
"""
check_dir_exists(options.identify_dir)
make_sure_path_exists(options.out_dir)
if not hasattr(options, 'outgroup_taxon'):
options.outgroup_taxon = None
markers = Markers(options.cpus, options.debug)
markers.align(options.identify_dir,
options.skip_gtdb_refs,
options.taxa_filter,
options.min_perc_aa,
options.custom_msa_filters,
options.skip_trimming,
options.rnd_seed,
options.cols_per_gene,
options.min_consensus,
options.max_consensus,