The bioBakery help forum

Average Read Coverage per Species Output

I would like to explore some count modelling statistics but MetaPlhan only outputs proportions. Could an option be included to output the median read coverage of all of the marker genes for each species? A median read coverage statistic wouldn’t be susceptible to differences in gene length between marker genes, unlike the plain counts which are proportional to both gene copies and gene length. -t rel_ab_w_read_stats is not what I want because it is larger for larger genomes but I want to estimate bacterial abundance independent of marker genes lengths.