Skip to content

Commit 4fd8027

Browse files
Bug fix in extract_all_unmapped
1 parent f062efc commit 4fd8027

File tree

1 file changed

+3
-2
lines changed

1 file changed

+3
-2
lines changed

kb_python/extract.py

+3-2
Original file line numberDiff line numberDiff line change
@@ -351,8 +351,9 @@ def extract(
351351

352352
if extract_all_unmapped:
353353
# Save unmapped reads in a separate fastq file
354-
unmapped_fastq = os.path.join(out_dir, "all_unmapped")
355-
extract_matching_reads_by_header(extract_out_folder, fastq, unmapped_fastq)
354+
unmapped_fastq = os.path.join(out_dir, "all_unmapped/1.fastq.gz")
355+
mapped_fastq = os.path.join(extract_out_folder, "1.fastq.gz")
356+
extract_matching_reads_by_header(mapped_fastq, fastq, unmapped_fastq)
356357

357358
else:
358359
if not mm:

0 commit comments

Comments
 (0)