Sgb_to_gtdb_profile.py keep read count estimation from -t rel_ab_w_read_stats

Microbial community profiling MetaPhlAn
1

Dear MetaPhlAn developers,

I am running MetaPhlAn4 v. 4.0.6 with addition of the parameter -t rel_ab_w_read_stats to obtain the estimate of the number of reads coming from each clade. I have been testing the sgb_to_gtdb_profile.py script and found that it only preserves the relative abundances, not the read count estimation. Is there currently a way to preserve the read counts or is this something you would consider making an option with this script?

Thank you in advance and thank you for developing these very nice programs!

Best regards

2

Dear @sarmol
In the current version of the script this is not possible, but we will make it available in the next version of the code

3

Dear @aitor.blancomiguez
That sounds great, thank you.