global lock
lock = l
l = multiprocessing.Lock()
# --------------------------------------------------------------------------- #
# This is the Wrapper script that runs the processing code found in #
# hisatgenotype_modules/hisatgenotype_typing_process #
# --------------------------------------------------------------------------- #
if __name__ == '__main__':
parser = ArgumentParser(
description="Extract variants from multiple sequence alignments")
# Add Arguments
arguments.args_databases(parser, location=True)
arguments.args_var_gaps(parser)
arguments.args_extract_vars(parser)
arguments.args_no_partial(parser)
arguments.args_common(parser)
args = parser.parse_args()
if args.locus_list == "":
locus_list = []
else:
locus_list = args.locus_list.split(',')
if args.inter_gap > args.intra_gap:
print("Error: --inter-gap (%d) must be smaller than --intra-gap (%d)" \
% (args.inter_gap, args.intra_gap),
file=sys.stderr)
sys.exit(1)
print("Error: indexing failed! "\
"Perhaps, you may have forgotten to build %s executables?" \
% aligner,
file=sys.stderr)
sys.exit(1)
# --------------------------------------------------------------------------- #
# Main function to build index and run script #
# --------------------------------------------------------------------------- #
if __name__ == '__main__':
parser = ArgumentParser(
description="Build genotype genome")
# Add Arguments
arguments.args_databases(parser)
arguments.args_var_gaps(parser)
arguments.args_set_aligner(parser,
False) # no missmatch option
arguments.args_build_genome(parser)
arguments.args_common(parser)
args = parser.parse_args()
if args.inter_gap > args.intra_gap:
print("Error: --inter-gap (%d) must be smaller than --intra-gap (%d)" \
% (args.inter_gap, args.intra_gap),
file=sys.stderr)
sys.exit(1)
if not args.base_fname:
args.base_fname = 'genotype_genome'
import sys
import os
from argparse import ArgumentParser, FileType
from hisatgenotype_modules.hisatgenotype_typing_process import extract_reads
import hisatgenotype_modules.hisatgenotype_args as arguments
# --------------------------------------------------------------------------- #
# This is the Wrapper script that runs the processing code found in #
# hisatgenotype_modules/hisatgenotype_typing_process #
# --------------------------------------------------------------------------- #
if __name__ == '__main__':
parser = ArgumentParser(description='Extract reads')
# Add Arguments
arguments.args_databases(parser, True,
True) # Add option to change genotype_genome name
arguments.args_input_output(parser)
arguments.args_aligner_inputs(parser)
arguments.args_set_aligner(parser, False) # No missmatch option
arguments.args_single_end(parser)
arguments.args_extract_reads(parser)
arguments.args_common(parser)
args = parser.parse_args()
if args.locus_list:
print(
"--locus-list option not implemented in this script yet, skipping")
if not args.graph_index:
print("--linear-index not implemented yet, skipping")