Hey, quick question. I need the final number of re...
# taxprofiler
h
Hey, quick question. I need the final number of reads of the analysis ready fastqs, although I did not output the reads themselves, so I can't verify this. Is this the same as the "reads_unmapped" column in the multiqc_samtools_stats.txt file in the output?
j
Actuaylly no, that's not supported at the moment... although it would make sense to support it!
If you could make an issue, I _think_I could add it quite quickly
You can't use the bbduk or samtools bits as these both happen before run merging... so if you do run-merging it could be misleading
If you can send your command I can tell you roughly where that info may apply
If you are running host-removal, then yeah probably the
samtools stats unmapped
is the closest you wil lget
But it would be useful to have the actual numbers in the analysis ready fastqs regardless what you are running for preprocessing (e.g. wiht
seqkit
maybe)
h
Thanks for your reply, James! I have added an issue on this. For this run, I did not have runmerging activated, as I only had one run for each sample (which is usually the case for us).
👍 1