read coverage calculation

[Maisie Nash]

Hi,

I was wondering specifically how the read coverage is caululated for each bin created by MaxBin2? I could not find this in the ReadMe!

Any information would be greatly apreciated.

Thanks!

[Yu-Wei Wu]

Dear Maisie,
I did not notice that there is a comment in the discussion until today. Sorry for my late reply. The reads coverage of each bin were calculated as follows. Assuming that there are n contigs c1 to cn, each with coverages a1 to an. So the coverage for c1 is a1, c2 is a2, etc. The bin abundance was estimated by

Bin_abund = sum(a_i * len(c_i)) / sum(len(c_i))

where i ranges from 1 to n, and len() indicates the contig length.

Hope that helps.