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:
