Menu

#7 off-by-one kmer coverage print

v1.0 (example)
open
nobody
None
5
2014-11-08
2014-11-08
No

Hi Zam,

When using Cortex during the hackathon, I ran into a bug when printing coverage information per unitig.

Line 10535 of CORTEX_release_v1.0.5.21/src/cortex_var/core/dB_graph_population.c:

for (i=0; i<len; i++)

should be

for (i=0; i<=len; i++)

as I believe that otherwise, it misses the coverage of the last kmer for each unitig.

Sample output (supernodes.fa) of the problem:

node_0 length:18 INFO:KMER=11
ATGGTTACCAAATTTTCC
Covg in Colour 0:
2 2 2 2 2 2 2

That sequence has length 18 so it contains exactly 18-11+1=8 11-mers, however only 7 coverages were printed.

Discussion


Log in to post a comment.

Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.