Method layoutBackbones

src/rnabloom/olc/Layout.java:4060–4151

(String outFastaPath)

    }
    
    public boolean layoutBackbones(String outFastaPath) throws IOException {
        HashSet<String> containedSet = populateGraphFromOverlaps();
        
        if (!containedSet.isEmpty()) {
            printMessage("contained reads: " + NumberFormat.getInstance().format(containedSet.size()));
        }
        
        Set<String> vertexSet = graph.vertexSet();
        if (!vertexSet.isEmpty()) {
            int numDovetailReads = vertexSet.size();
            if (!stranded) {
                numDovetailReads = numDovetailReads/2;
            }
            printMessage("dovetail reads:  " + NumberFormat.getInstance().format(numDovetailReads));
        }
        
        Set<OverlapEdge> edgeSet = graph.edgeSet();
        
        if (edgeSet.size() > 1) {
            printMessage("G: |V|=" + NumberFormat.getInstance().format(vertexSet.size()) + " |E|=" + NumberFormat.getInstance().format(edgeSet.size()));

//            ArrayDeque<String> redundantNodeNames = removeRedundantNodes();
//            for (String n : redundantNodeNames) {
//                containedSet.add(getVertexName(n));
//            }
//            
//            printGraph();
//            numEdges = graph.edgeSet().size();
//            System.out.println("G: |V|=" + NumberFormat.getInstance().format(graph.vertexSet().size()) + " |E|=" + NumberFormat.getInstance().format(numEdges));

            resolveJunctions();
            printMessage("G: |V|=" + NumberFormat.getInstance().format(vertexSet.size()) + " |E|=" + NumberFormat.getInstance().format(edgeSet.size()));
        }
        
        // extract read sequences
        int numVertexReads = vertexSet.size();
        if (!stranded) {
            numVertexReads = numVertexReads/2;
        }
        
        HashMap<String, BitSequence> dovetailReadSeqs = new HashMap<>(numVertexReads);
        FastaReader fr = new FastaReader(seqFastaPath);
        FastaWriter fw = new FastaWriter(outFastaPath, false);
        long seqID = 0;
        long originalNumSeq = 0;
        while (fr.hasNext()) {
            ++originalNumSeq;
            String[] nameSeq = fr.nextWithName();
            String name = nameSeq[0];
            if (vertexSet.contains(name + '+')) {
                dovetailReadSeqs.put(name, new BitSequence(nameSeq[1]));
            }
            else if (!containedSet.contains(name)) {
                // this sequence either contains shorter sequences or it has no overlaps with others
                fw.write(Long.toString(++seqID), nameSeq[1]);
            }
        }
        fr.close();