Bubble caller: don't report serial bubbles

mcveanlab / mccortex

De novo genome assembly and multisample variant calling

https://github.com/mcveanlab/mccortex/wiki

MIT License

113 stars 25 forks source link

Closed noporpoise closed 8 years ago

noporpoise commented 9 years ago

  _    _    _  
_/a\__/b\__/c\_
 \_/  \_/  \_/ 
  A    B    C

We currently report adjacent bubbles as one big bubble. For the given example we report abc|ABC, bc|BC, c|C, and the reverse calls.

Instead we should only report stand alone bubbles: a|A, b|B, c|C.

noporpoise commented 8 years ago

Fixed.