Permalink
Browse files

Require at least two occurrences in the base sequence when making com…

…parative calls
  • Loading branch information...
1 parent 20519aa commit ef79e9e603d08cc78414d3e62cafaf241f624458 @jts committed Jul 8, 2012
Showing with 2 additions and 1 deletion.
  1. +2 −1 src/GraphDiff/GraphCompare.cpp
@@ -146,14 +146,15 @@ GraphCompareResult GraphCompare::process(const SequenceWorkItem& item)
count += rc_interval.size();
bool both_strands = interval.size() > 0 && rc_interval.size() > 0;
+ size_t min_base_coverage = m_parameters.bReferenceMode ? 1 : 2;
if(count >= m_parameters.minDiscoveryCount && count < m_parameters.maxDiscoveryCount && !variantAttempted && both_strands)
{
// Check if this k-mer is present in the other base index
size_t base_count = BWTAlgorithms::countSequenceOccurrences(kmer, m_parameters.baseIndex);
// k-mer present in the base read set, skip it
- if(base_count > 0)
+ if(base_count >= min_base_coverage)
continue;
if(m_parameters.verbose > 0)

0 comments on commit ef79e9e

Please sign in to comment.