DaehwanKimLab / centrifuge

Classifier for metagenomic sequences
GNU General Public License v3.0
246 stars 73 forks source link

Abundances on different taxonomic levels conflict #143

Open GrantHov opened 6 years ago

GrantHov commented 6 years ago

I want to analyze abundances on genus level. However I see that genus levels abundances (and read counts) are less than abundances of species in the corresponding genus, which seems a bit counter intuitive to me.

For example, if I do "grep Homo A_Y_T2_24706_ACTTGAtr_report.tsv", I get Homo 9605 genus 0 12 0 0.0 Homo sapiens 9606 species 3238442024 738484 395792 1.26543e-06

Shouln't the software summarize abundances in hierarchical manner? Please let me know if I miss something. Cheers,