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")