Commit 13b93f7b authored by Tizian Schulz's avatar Tizian Schulz
Browse files

Minor changes on comparison workflow

parent 0453dbdd
......@@ -238,7 +238,7 @@ rule search:
rule makeUdb:
input:
"%s/merged_{panID}_part{nr}.fasta" %config['dataDir']
"%s/merged_{panID}_part{nr}.fa" %config['dataDir']
output:
temp("UBLAST/dbs/{panID}_part{nr}.udb")
shell:
......@@ -249,7 +249,7 @@ rule makeChunks:
glist = "%s/chosen_{panID}.txt" %config['dataDir'],
seqs = expand("{assembly}", assembly=[i + "tmp" for i in glob("%s/*.fasta" %config['dataDir'])])
output:
temp("%s/merged_{panID}_part{nr}.fasta" %config['dataDir'])
temp("%s/merged_{panID}_part{nr}.fa" %config['dataDir'])
shell:
"python3 scripts/writeChunks.py {input.glist}"
......
......@@ -15,7 +15,7 @@ for n in colorNamesLine.split(' ')[:-1]:
if colorCounter > 0:
outfile.close()
outfileName = '/'.join(sys.argv[1].split('/')[:-1]) + "/merged_" + sys.argv[1].split(".txt")[0].split("chosen_")[1] + "_part" + str(int(colorCounter / CHUNK_SIZE) + 1) + ".fasta"
outfileName = '/'.join(sys.argv[1].split('/')[:-1]) + "/merged_" + sys.argv[1].split(".txt")[0].split("chosen_")[1] + "_part" + str(int(colorCounter / CHUNK_SIZE) + 1) + ".fa"
outfile = open(outfileName, 'w')
fastaFile = open(n.replace("fasta", "fastatmp").replace('\n',''), 'r')
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment