Method extendLeftSE

src/rnabloom/util/GraphUtils.java:6112–6204

(ArrayList<Kmer> kmers, 
                                            BloomFilterDeBruijnGraph graph, 
                                            int maxTipLen,
                                            float minKmerCov)

    }
        
    private static ArrayDeque<Kmer> extendLeftSE(ArrayList<Kmer> kmers, 
                                            BloomFilterDeBruijnGraph graph, 
                                            int maxTipLen,
                                            float minKmerCov) {
        // `kmers` is reversed
        
        final int k = graph.getK();
        final int numHash = graph.getMaxNumHash();
        final int readPairedKmersDist = graph.getReadPairedKmerDistance();
        final int numKmers = kmers.size();
        int maxExtensionLength = readPairedKmersDist - 2; // -1 for candidate k-mer; -1 for partner kmer on current sequence
        
        ArrayDeque<Kmer> candidates = kmers.get(numKmers-1).getPredecessors(k, numHash, graph);
        
        if (candidates.size() == 1) {
            Kmer c = candidates.peek();
            ArrayDeque<Kmer> e = naiveExtendLeftNoBackChecks(c, graph, maxTipLen, maxExtensionLength, minKmerCov);
            e.addFirst(c);
            return e;
        }
        
        final float pathMinCov = getMinimumKmerCoverage(kmers, Math.max(numKmers - readPairedKmersDist, 0), numKmers);
        float bestScore = 0;
        float bestCov = 0;
        ArrayDeque<Kmer> bestExtension = null;
        int[] result = new int[2];
        
        for (Kmer candidate : candidates) {
            ArrayDeque<Kmer> e = naiveExtendLeftNoBackChecks(candidate, graph, maxTipLen, maxExtensionLength, minKmerCov);
            e.addFirst(candidate);
            
            countKmerPairsReversedSE(e, kmers, 0, graph, result);
            
            int lastPartneredKmerIndex = result[1];
            
            if (lastPartneredKmerIndex >= 0 && result[0] > 0) {
                float cov = getMedianKmerCoverage(e);
                float score = Math.min(pathMinCov, cov) * result[0] / (lastPartneredKmerIndex+1);
                if (score > bestScore || (score == bestScore && cov > bestCov)) {
                    bestScore = score;
                    bestCov = cov;
                    bestExtension = e;
                    
                    Iterator<Kmer> itr = e.descendingIterator();
                    for (int i=e.size()-1; i>lastPartneredKmerIndex; --i) {
                        itr.next();
                        itr.remove();
                    }
                }
            }
            else {
                int gap = e.size();
                
                if (gap >= readPairedKmersDist-1 && result[0] == 0) {
                    continue;
                }

                // not enough supporting paired k-mers in first extension