sourmash
sourmash copied to clipboard
does `tax metagenome -F lineage_summary` use abundance-weighted output?
it looks like it does not, which is surprising to me.
Specifically, on one sample (JRH_1) I see in gather output:
the recovered matches hit 73.1% of the abundance-weighted query. the recovered matches hit 36.3% of the query k-mers (unweighted).
and then in the lineage_summary TSV file I see that 63.688% are unclassified.
So actually I'm not sure that number makes sense either. Huh.
(The kreport format looks fine! 26.9% unclassified!)