if args.start is not None: start = args.start if args.indexes is not None: createIndexes = args.indexes if args.queries is not None: runQueries = args.queries # Open results file and print header resultsFileName = 'results-pgsql' if resultsFileName != "": resultsFileName += '-' + tag resultsFileName += '.txt' resultsFile = open(resultsFileName, 'w') resultsFile.write(scriptVersion + '\n') result = Result() resultsFile.write(result.toHeader() + '\n') # Data files snpFilePath = 'snpData-chr{0}.txt' lociFilePath = 'lociData-chr{0}.txt' # Chromosome list chromosomes = ["21"] # dev list # If not in dev mode, iterate through all chromosomes if dev is False: chromosomes = ["1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","X","Y","MT"] # complete list if start != "1": # Allow restart from anywhere in chromosome list, sequentially as ordered above startList = [] hitMin = False for cur in chromosomes: