sufforest / SolidBin

A genome binning method for contig binning, based on semi-supervised spectral clustering method.
GNU General Public License v2.0
13 stars 3 forks source link

Key Error when seed contig is 1000 bp long #6

Closed Vini2 closed 4 years ago

Vini2 commented 4 years ago

When building the composition profiles, by default SolidBin keeps contigs longer than 1000 bp. However, when running test_getmarker.pl to get the seeds, the minimum contig length is set to 1000 bp. I came across a case where I got a seed contig which is 1000 bp long but its composition profile was not created. This resulted in a Key Error at line 304. So I think the minimum contig length should be set to 1001 in line 141 of SolidBin.py. Setting the minimum contig length to 1001 solved the Key Error. Please correct me if I'm wrong.

Thank you!

ziyewang commented 4 years ago

I think that you are correct. Thanks very much.

Best, ziye