BioJulia / GenomeGraphs.jl

A modern genomics framework for julia
https://biojulia.dev
Other
67 stars 8 forks source link

Implement draft of perfect circle resolution for the dbg constructor #16

Closed TransGirlCodes closed 4 years ago

TransGirlCodes commented 4 years ago

Types of changes

This PR implements the following changes: (Please tick any or all of the following that are applicable)

:clipboard: Additional detail

Fixes #13

The DBG building process currently may fail to incorporate some kmers into a DBG if they exist in a perfect circle. Perfect circles don't have "ends" and so the kmers involved in them don't get incorporated into a DBG during construction at present

This commit implements a draft improvement that fixes this problem. After constructing unitigs and consuming kmers in the usual way. If there are any kmers left, they are left due to this circle problem, so it starts by picking one of the kmers involved in a circle and just keeps going forward until that same kmer is found again. The circle then is turned into a unitig - broken at an arbitrary point.

:ballot_box_with_check: Checklist