os.makedirs(os.path.dirname(args.output_prefix)) else: logger.error("The output prefix must be defined") end(1) ### Check that query files exist if args.query: QueryFile = args.query AliQueryFile = args.query if not os.path.isfile(QueryFile): logger.error(QueryFile+" (-q) is not a file.") end(1) elif not os.stat(QueryFile).st_size: logger.error(QueryFile+" (-q) is empty.") end(1) elif ApytramNeeds.count_sequences(QueryFile) !=1: logger.warning("%s (-q) contains more than one query. They are %s sequences." %(QueryFile,ApytramNeeds.count_sequences(QueryFile))) # If there are multiple probes, align them for the future coverage counter # Use Mafft start_mafft_time = time.time() MafftProcess = Aligner.Mafft(QueryFile) MafftProcess.QuietOption = True MafftProcess.AutoOption = True (MafftResult, err) = MafftProcess.get_output() AliQueryFile = "%s/References.ali.fasta" %TmpDirName ApytramNeeds.write_in_file(MafftResult,AliQueryFile) logger.debug("mafft --- %s seconds ---" % (time.time() - start_mafft_time)) # If the -pep option is used, the -q option must be precised if args.query_pep: