Empty files were also published from the module samtools/fastq
Closed this issue · 0 comments
LilyAnderssonLee commented
Description of the bug
A bug was found in Taxprofiler
versions newer than 1.0.1. Three FASTQ files were published by the module samtools/fastq
after converting the unmapped human reads BAM file generated from minimap2
, and two of them are empty. All these three files were passed to the classifiers and it causes errors.
Files under the folder results/samtools/fastq
-rw-rw----+ 1 lili.li users 28 Jan 19 10:34 sample1.unmapped_1.fastq.gz
-rw-rw----+ 1 lili.li users 55M Jan 19 10:34 sample1.unmapped_other.fastq.gz
-rw-rw----+ 1 lili.li users 28 Jan 19 10:34 sample.unmapped_singleton.fastq.gz
Error message:
Command executed:
kraken2 \
--db k2_pluspf_20230605 \
--threads 12 \
--report sample1_k2_pluspf.kraken2.kraken2.report.txt \
--gzip-compressed \
--unclassified-out sample1_k2_pluspf.kraken2.unclassified.fastq \
--classified-out sample1_k2_pluspf.kraken2.classified.fastq \
--output /dev/null \
\
\
sample1.unmapped_1.fastq.gz
pigz -p 12 *.fastq
cat <<-END_VERSIONS > versions.yml
"NFCORE_TAXPROFILER:TAXPROFILER:PROFILING:KRAKEN2_KRAKEN2":
kraken2: $(echo $(kraken2 --version 2>&1) | sed 's/^.*Kraken version //; s/ .*$//')
pigz: $( pigz --version 2>&1 | sed 's/pigz //g' )
END_VERSIONS
Command exit status:
0
Command output:
(empty)
Command error:
Loading database information... done.
0 sequences (0.00 Mbp) processed in 0.040s (0.0 Kseq/m, 0.00 Mbp/m).
0 sequences classified (-nan%)
0 sequences unclassified (-nan%)
pigz: skipping: *.fastq does not exist
Command used and terminal output
nextflow run nf-core/taxprofiler -r 1.1.2 -profile hasta,singularity --input cases/20231227/samplesheet_20231227.csv --databases cases/20231227/databases_20231227.csv --outdir results/20231227 --save_preprocessed_reads --perform_longread_qc --longread_qc_skipqualityfilter --perform_longread_hostremoval --hostremoval_reference references/GCF_000001405.39_GRCh38.p13_genomic.fna --save_hostremoval_index --save_hostremoval_unmapped --run_kraken2 --kraken2_save_reads --kraken2_save_readclassification --run_centrifuge --centrifuge_save_reads --run_krona -params-file scripts/params.json -resume
Relevant files
No response
System information
Taxprofiler V.1.1.2