nf-core / nanoseq

Nanopore demultiplexing, QC and alignment pipeline
https://nf-co.re/nanoseq
MIT License
179 stars 80 forks source link

Different Chromosome counts between samtools and bambu #278

Open fbndr98 opened 1 month ago

fbndr98 commented 1 month ago

Description of the bug

Hello, I hope to get help or maybe better an explanation of the following issue. I have run the nanoseq pipeline using standard settings on a couple of samples in the meantime. I'm interested in the counts per chromosome but by chance I realized, that the counts differ between samtools (.idxstats in the samtools output or the raw index counts of the mutliqc report) and the bambu output counts_gene.txt if gene counts for each chromosome are summed. Not too much, but the numbers do differ. I'm a bit confused why this is the case and would be happy for suggestions which statistic is the more accurate or more reliable. Thanks!

Command used and terminal output

No response

Relevant files

No response

System information

No response