Skip to content

Commit

Permalink
Add a comment (and asserts) about a strange test to flag bubbles as c…
Browse files Browse the repository at this point in the history
…ontigs.
  • Loading branch information
brianwalenz committed Sep 15, 2020
1 parent 276f344 commit a5f397b
Showing 1 changed file with 9 additions and 0 deletions.
9 changes: 9 additions & 0 deletions src/bogart/AS_BAT_Instrumentation.C
Original file line number Diff line number Diff line change
Expand Up @@ -272,11 +272,20 @@ classifyTigsAsUnassembled(TigVector &tigs,
bool r4 = classifyRule4(utg, F, nCoverage, bCoverage, lowcovFraction, lowcovDepth);

// If either first or last read is flagged as a bubble, it's not unassembled.
//
// Sept 2020 - bri isn't sure if this is being clever about something or
// not. When bubbles are assigned - search for setBubble - we're both
// setting the tig _isBubble to true and labelling all reads as
// isBubble(). Could we then be merging a bubble into something else?
// The asserts were added.

if ((OG->isBubble(utg->firstRead()->ident) == true) ||
(OG->isBubble(utg-> lastRead()->ident) == true))
bb = true;

if (bb == true) assert(utg->_isBubble == true);
if (bb == false) assert(utg->_isBubble == false);

// If flagged, we're done, just move on.

if ((bb == false) && (rr == false) && (rs || r1 || r2 || r3 || r4))
Expand Down

0 comments on commit a5f397b

Please sign in to comment.