off-by-one kmer coverage print
Brought to you by:
mariocaccamo,
zamin
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.