notebook.community

Edit and run 





In [354]:



    


%pylab inline
%config InlineBackend.figure_format = 'retina'
import py
import pandas as pd
import seaborn as sns
from collections import OrderedDict
import glob
import os
import re
from IPython.display import FileLink
import pyprind
from khmer import Countgraph
from khmer._oxli.partitioning import StreamingPartitioner, Component
from khmer._oxli.traversal import Traverser
from khmer._oxli.hashing import Kmer
from khmer._oxli.app import PartitioningApp
from khmer import reverse_hash
import json
import cython
import cythonmagic
import networkx as nx
from khmer import reverse_complement



    


















    






Populating the interactive namespace from numpy and matplotlib

















In [2]:



    


%load_ext cython



    


















    






/home/camille/miniconda/envs/dev.py3.khmer.assembly/lib/python3.5/site-packages/Cython/Distutils/old_build_ext.py:30: UserWarning: Cython.Distutils.old_build_ext does not properly handle dependencies and is deprecated.
  "Cython.Distutils.old_build_ext does not properly handle dependencies "

















In [3]:



    


sns.set_style('ticks')



    











In [4]:



    


sp = StreamingPartitioner.load('../data/SRR1300451.partitioner.2.json')



    


















    






('Loading', '../data/SRR1300451.partitioner.2.graph', 'as CountGraph')

















In [7]:



    


components = sp.component_dict
lump = None
for cid, comp in components.items():
    if lump is None or len(comp) > len(lump):
        lump = comp



    











In [8]:



    


rep_comp = None
for cid, comp in components.items():
    if len(comp) > 1000 and len(comp) < 2000:
        rep_comp = comp
        break



    











In [355]:



    


lump_counts = Component.tag_counts(lump, sp.graph)



    











In [356]:



    


sns.distplot(lump_counts, kde=False)



    


















    
Out[356]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f3be9660a58>










    
























In [357]:



    


lump_degrees = []
for tag in lump:
    lump_degrees.append(trav.degree(Kmer.create(tag, sp.graph.ksize())))



    











In [358]:



    


lump_kmers = pd.DataFrame({'tag': list(lump), 
                           'kmer': [reverse_hash(tag, sp.graph.ksize()) for tag in lump],
                           'count': lump_counts,
                           'degree': lump_degrees})



    











In [359]:



    


lump_kmers['n_A'] = lump_kmers.kmer.str.count('A')
lump_kmers['n_C'] = lump_kmers.kmer.str.count('C')
lump_kmers['n_G'] = lump_kmers.kmer.str.count('G')
lump_kmers['n_T'] = lump_kmers.kmer.str.count('T')



    











In [360]:



    


lump_kmers['GC'] = (lump_kmers['n_C'] + lump_kmers['n_G']) / sp.graph.ksize()



    











In [361]:



    


sns.distplot(lump_kmers.GC, kde=False, bins=sp.graph.ksize())



    


















    
Out[361]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f3c202479b0>










    
























In [362]:



    


lump_kmers.groupby('GC').mean()['degree'].plot()
ylabel('Mean $k$-mer degree')



    


















    
Out[362]:








<matplotlib.text.Text at 0x7f3c20168630>










    
























In [363]:



    


lump_kmers = pd.melt(lump_kmers, 
                     id_vars=['kmer', 'tag', 'count', 'degree', 'GC'], 
                     value_vars=['n_A', 'n_C', 'n_G', 'n_T'], 
                     value_name='nucl_count')



    











In [364]:



    


lump_kmers.head(10)



    


















    
Out[364]:










  
    

      
      kmer
      tag
      count
      degree
      GC
      variable
      nucl_count
    

  
  
    

      0
      AAAAAAAAAAAAAAAAAAAAAAAAA
      0
      255
      8
      0.00
      n_A
      25
    

    

      1
      AAAAAAAAAAAAAAAAAAAAAAAAT
      1
      200
      8
      0.00
      n_A
      24
    

    

      2
      AAAAAAAAAAAAAAAAAAAAAAAAC
      2
      255
      8
      0.04
      n_A
      24
    

    

      3
      AAAAAAAAAAAAAAAAAAAAAAAAG
      3
      226
      8
      0.04
      n_A
      24
    

    

      4
      AAAAAAAAAAAAAAAAAAAAAAATA
      4
      106
      8
      0.00
      n_A
      24
    

    

      5
      AAAAAAAAAAAAAAAAAAAAAAATT
      5
      71
      8
      0.00
      n_A
      23
    

    

      6
      AAAAAAAAAAAAAAAAAAAAAAATG
      7
      3
      8
      0.04
      n_A
      23
    

    

      7
      AAAAAAAAAAAAAAAAAAAAAAACA
      8
      38
      8
      0.04
      n_A
      24
    

    

      8
      AAAAAAAAAAAAAAAAAAAAAAACT
      9
      226
      8
      0.04
      n_A
      23
    

    

      9
      AAAAAAAAAAAAAAAAAAAAAAACC
      10
      22
      8
      0.08
      n_A
      23
    

  




















In [365]:



    


lump_nucl_content = lump_kmers.groupby(['variable', 'nucl_count']).mean()[['count', 'degree']]



    











In [366]:



    


lump_nucl_content.unstack(level=0).plot(y='count', xlim=[15,25])
ylabel('Mean $k$-mer count')
xlabel('Nucleotide Count')
title('Lump Tag Coverage by Nucleotide Content')
savefig('lump_tag_cov_by_nucl.png', dpi=600)



    


















    
























In [367]:



    


FileLink('images/lump_tag_cov_by_nucl.png')



    


















    
Out[367]:






Path (images/lump_tag_cov_by_nucl.png) doesn't exist. It may still be in the process of being generated, or you may have the incorrect path.
















In [368]:



    


lump_nucl_content.unstack(level=0).plot(y='degree')
ylabel('Mean Degree')
xlabel('Nucleotide Count')
title('Lump Tag Degrees by Nucleotide Content')
savefig('images/lump_tag_deg_by_nucl.png', dpi=600)



    


















    
























In [369]:



    


FileLink('lump_tag_deg_by_nucl.png')



    


















    
Out[369]:






lump_tag_deg_by_nucl.png

















In [370]:



    


lump_kmers.degree.value_counts(sort=False).plot(kind='bar')



    


















    
Out[370]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f3c20082080>










    
























In [301]:



    


def extract_local_compact_dbg(start, graph, max_breadth=None, verbose=True):

    from khmer._oxli.assembly import CompactingAssembler
    import networkx as nx
    from collections import deque
    from khmer import reverse_complement
    
    def kmer_str(kmer):
        return '{0}/{1}'.format(kmer.kmer, reverse_complement(kmer.kmer))

    K = graph.ksize()
    G = nx.Graph()
    q = deque()
    breadth = deque()
    seen = set()
    
    kmers = set()
    def add_kmers(path):
        for i in range(0, len(path)-K+1):
            kmers.add(hash(Kmer(path[i:i+K])))

    trav = Traverser(graph)
    asm = CompactingAssembler(graph)
    path = asm.assemble(Kmer(start))
    
    def add(src, dst, path):
        path_cov, _, _ = graph.get_median_count(path)
        if verbose:
            add_kmers(path)
            add_kmers(src.kmer)
            add_kmers(dst.kmer)
        G.add_edge(hash(src), hash(dst), path=path, Weight=len(path)-K+1, coverage=path_cov)
        G.add_node(hash(src), kmer=kmer_str(src), coverage=graph.get(src.kmer))
        G.add_node(hash(dst), kmer=kmer_str(dst), coverage=graph.get(dst.kmer))

    src = Kmer(path[:K])
    dst = Kmer(path[-K:])
    if hash(src) != hash(dst):
        add(src, dst, path)
        q.append(src)
        breadth.append(0)
    q.append(dst)
    breadth.append(0)

    n = 0
    while(len(q) > 0):
        src = q.popleft()
        cur_breadth = breadth.popleft()

        if n % 1000 == 0:
            if verbose:
                print(n, 'iterations, ', len(G.nodes()), 'nodes,', len(q), 'on queue, breadth at', cur_breadth)
            else:
                if n % 50000 == 0:
                    print(n, 'iterations, ', len(G.nodes()), 'nodes,', len(q), 'on queue, breadth at', cur_breadth)
        if max_breadth is not None and cur_breadth > max_breadth:
            if verbose:
                print('max breadth reached at', len(G.nodes()), 'nodes')
            break
        n += 1

        for neighbor in trav.neighbors(src):
            if hash(neighbor) in seen:
                continue
            if hash(neighbor) == hash(src):
                continue
            seen.add(hash(neighbor))

            path = asm.assemble(neighbor)
            #print(neighbor, path, len(path))
            if len(path) == K:
                add(src, neighbor, path)
                seen.add(hash(src))
                q.append(neighbor)
                breadth.append(cur_breadth + 1)
                continue

            left = Kmer(path[:K])
            right = Kmer(path[-K:])
            if hash(neighbor) == hash(left):
                dst = right
                seen.add(hash(Kmer(path[-K-1:-1])))
            elif hash(src) == hash(left):
                dst = right
            elif hash(neighbor) == hash(right):
                dst = left
                seen.add(hash(Kmer(path[1:K+1])))
            elif hash(src) == hash(right):
                dst = left
            else:
                print(src, trav.degree(src), '\n', 
                      neighbor, reverse_complement(neighbor.kmer), trav.degree(neighbor), '\n',
                      left, right, path, seen)
                raise RuntimeError()

            add(src, dst, path)
            q.append(dst)
            breadth.append(cur_breadth + 1)
    # normalize edge weights
    max_weight = max([G.get_edge_data(*e)['Weight'] for e in G.edges()])
    for e in G.edges():
        data = G.get_edge_data(*e)
        data['Weight'] = float(data['Weight']) / max_weight

    if verbose:
        print(n, 'iterations, ', len(G.nodes()), 'nodes,', len(q), 'on queue, breadth at', 
              cur_breadth, len(kmers), 'k-mers')
    return G



    











In [12]:



    


max(list(nx.degree(G).values()))



    


















    
Out[12]:








8

















In [96]:



    


G_lump.get_edge_data(0, 281474976710656)['Weight'] = 4



    











In [103]:



    


[G_lump.get_edge_data(*e)['Weight'] for e in G_lump.edges()]



    


















    
Out[103]:








[0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.07086614173228346,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.2047244094488189,
 0.007874015748031496,
 0.007874015748031496,
 0.07086614173228346,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.023622047244094488,
 0.023622047244094488,
 0.007874015748031496,
 0.12598425196850394,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.16535433070866143,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.11023622047244094,
 0.10236220472440945,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.2125984251968504,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.047244094488188976,
 0.007874015748031496,
 0.06299212598425197,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.023622047244094488,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.10236220472440945,
 0.08661417322834646,
 0.10236220472440945,
 0.07086614173228346,
 0.007874015748031496,
 0.023622047244094488,
 0.023622047244094488,
 0.031496062992125984,
 0.13385826771653545,
 0.023622047244094488,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.023622047244094488,
 0.023622047244094488,
 0.15748031496062992,
 0.1732283464566929,
 0.023622047244094488,
 0.03937007874015748,
 0.13385826771653545,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.07086614173228346,
 0.07874015748031496,
 0.007874015748031496,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.07874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.06299212598425197,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.09448818897637795,
 0.007874015748031496,
 0.007874015748031496,
 0.07874015748031496,
 0.09448818897637795,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.16535433070866143,
 0.015748031496062992,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.023622047244094488,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.2047244094488189,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.015748031496062992,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.08661417322834646,
 0.007874015748031496,
 0.007874015748031496,
 0.1732283464566929,
 0.05511811023622047,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.14173228346456693,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.047244094488188976,
 0.023622047244094488,
 0.023622047244094488,
 0.007874015748031496,
 0.1732283464566929,
 0.007874015748031496,
 0.12598425196850394,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.1889763779527559,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.06299212598425197,
 0.015748031496062992,
 0.023622047244094488,
 0.09448818897637795,
 0.14173228346456693,
 0.023622047244094488,
 0.10236220472440945,
 0.14960629921259844,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.1968503937007874,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.07086614173228346,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.031496062992125984,
 0.023622047244094488,
 0.023622047244094488,
 0.14960629921259844,
 0.1889763779527559,
 0.007874015748031496,
 0.007874015748031496,
 0.07874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.015748031496062992,
 0.2047244094488189,
 0.007874015748031496,
 0.007874015748031496,
 0.14960629921259844,
 0.007874015748031496,
 0.007874015748031496,
 0.13385826771653545,
 0.18110236220472442,
 0.07086614173228346,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.07874015748031496,
 0.10236220472440945,
 0.03937007874015748,
 0.007874015748031496,
 0.05511811023622047,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.06299212598425197,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.09448818897637795,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.07874015748031496,
 0.007874015748031496,
 0.16535433070866143,
 0.015748031496062992,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.11811023622047244,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.13385826771653545,
 0.1732283464566929,
 0.16535433070866143,
 0.023622047244094488,
 0.1732283464566929,
 0.023622047244094488,
 0.023622047244094488,
 0.007874015748031496,
 0.023622047244094488,
 0.031496062992125984,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.031496062992125984,
 0.03937007874015748,
 0.023622047244094488,
 0.007874015748031496,
 0.13385826771653545,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.2047244094488189,
 0.03937007874015748,
 0.031496062992125984,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.023622047244094488,
 0.023622047244094488,
 0.023622047244094488,
 0.015748031496062992,
 0.007874015748031496,
 0.007874015748031496,
 0.15748031496062992,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.015748031496062992,
 0.007874015748031496,
 0.1968503937007874,
 0.14960629921259844,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.07086614173228346,
 0.11023622047244094,
 0.007874015748031496,
 0.015748031496062992,
 0.12598425196850394,
 0.11023622047244094,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.1732283464566929,
 0.047244094488188976,
 0.14960629921259844,
 0.09448818897637795,
 0.047244094488188976,
 0.05511811023622047,
 0.007874015748031496,
 0.007874015748031496,
 0.14173228346456693,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.18110236220472442,
 0.18110236220472442,
 0.18110236220472442,
 0.06299212598425197,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.14173228346456693,
 0.023622047244094488,
 0.12598425196850394,
 0.2125984251968504,
 0.015748031496062992,
 0.07874015748031496,
 0.031496062992125984,
 0.07874015748031496,
 0.007874015748031496,
 0.015748031496062992,
 0.03937007874015748,
 0.007874015748031496,
 0.13385826771653545,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.16535433070866143,
 0.05511811023622047,
 0.07086614173228346,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.11023622047244094,
 0.047244094488188976,
 0.007874015748031496,
 0.06299212598425197,
 0.14960629921259844,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.16535433070866143,
 0.023622047244094488,
 0.12598425196850394,
 0.1889763779527559,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.18110236220472442,
 0.007874015748031496,
 0.007874015748031496,
 0.015748031496062992,
 0.15748031496062992,
 0.1889763779527559,
 0.07086614173228346,
 0.015748031496062992,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.05511811023622047,
 0.007874015748031496,
 0.03937007874015748,
 0.007874015748031496,
 0.09448818897637795,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.11023622047244094,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.03937007874015748,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.047244094488188976,
 0.007874015748031496,
 0.031496062992125984,
 0.031496062992125984,
 0.1889763779527559,
 0.007874015748031496,
 0.10236220472440945,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.08661417322834646,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.07086614173228346,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.13385826771653545,
 0.13385826771653545,
 0.07874015748031496,
 0.14173228346456693,
 0.2047244094488189,
 0.1889763779527559,
 0.007874015748031496,
 0.1968503937007874,
 0.13385826771653545,
 0.15748031496062992,
 0.11811023622047244,
 0.007874015748031496,
 0.007874015748031496,
 0.10236220472440945,
 0.007874015748031496,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.023622047244094488,
 0.14960629921259844,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.09448818897637795,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.13385826771653545,
 0.07874015748031496,
 0.007874015748031496,
 0.09448818897637795,
 0.007874015748031496,
 0.007874015748031496,
 0.05511811023622047,
 0.023622047244094488,
 0.12598425196850394,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.07086614173228346,
 0.023622047244094488,
 0.023622047244094488,
 0.007874015748031496,
 0.16535433070866143,
 0.031496062992125984,
 0.11023622047244094,
 0.03937007874015748,
 0.07874015748031496,
 0.11023622047244094,
 0.015748031496062992,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.03937007874015748,
 0.07086614173228346,
 0.007874015748031496,
 0.007874015748031496,
 0.05511811023622047,
 0.11023622047244094,
 0.03937007874015748,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.10236220472440945,
 0.11023622047244094,
 0.1732283464566929,
 0.14960629921259844,
 0.007874015748031496,
 0.023622047244094488,
 0.03937007874015748,
 0.06299212598425197,
 0.007874015748031496,
 0.2125984251968504,
 0.10236220472440945,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.031496062992125984,
 0.007874015748031496,
 0.07086614173228346,
 0.023622047244094488,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.03937007874015748,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.10236220472440945,
 0.047244094488188976,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.03937007874015748,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.14173228346456693,
 0.023622047244094488,
 0.023622047244094488,
 0.06299212598425197,
 0.007874015748031496,
 0.11023622047244094,
 0.007874015748031496,
 0.007874015748031496,
 0.047244094488188976,
 0.031496062992125984,
 0.05511811023622047,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.031496062992125984,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.03937007874015748,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.07086614173228346,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.031496062992125984,
 0.007874015748031496,
 0.28346456692913385,
 0.18110236220472442,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.023622047244094488,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.007874015748031496,
 0.06299212598425197,
 0.023622047244094488,
 0.023622047244094488,
 0.03937007874015748,
 0.08661417322834646,
 0.11023622047244094,
 0.031496062992125984,
 0.13385826771653545,
 0.2440944881889764,
 ...]

















In [139]:



    


G_lump = extract_local_compact_dbg('A' * 25, sp.graph, max_breadth=50)



    


















    






0 iterations,  0 nodes, 0 on queue, breadth at 0
1000 iterations,  2075 nodes, 1237 on queue, breadth at 5
2000 iterations,  3317 nodes, 1766 on queue, breadth at 6
3000 iterations,  4479 nodes, 2272 on queue, breadth at 7
4000 iterations,  5519 nodes, 2659 on queue, breadth at 7
5000 iterations,  6410 nodes, 2913 on queue, breadth at 8
6000 iterations,  7315 nodes, 3272 on queue, breadth at 8
7000 iterations,  8196 nodes, 3578 on queue, breadth at 8
8000 iterations,  9047 nodes, 3817 on queue, breadth at 9
9000 iterations,  9863 nodes, 4049 on queue, breadth at 9
10000 iterations,  10684 nodes, 4255 on queue, breadth at 9
11000 iterations,  11484 nodes, 4474 on queue, breadth at 9
12000 iterations,  12358 nodes, 4763 on queue, breadth at 10
13000 iterations,  13120 nodes, 4960 on queue, breadth at 10
14000 iterations,  13921 nodes, 5152 on queue, breadth at 10
15000 iterations,  14744 nodes, 5375 on queue, breadth at 10
16000 iterations,  15553 nodes, 5626 on queue, breadth at 10
17000 iterations,  16403 nodes, 5877 on queue, breadth at 11
18000 iterations,  17180 nodes, 6063 on queue, breadth at 11
19000 iterations,  18053 nodes, 6340 on queue, breadth at 11
20000 iterations,  18818 nodes, 6498 on queue, breadth at 11
21000 iterations,  19555 nodes, 6657 on queue, breadth at 11
22000 iterations,  20336 nodes, 6853 on queue, breadth at 12
23000 iterations,  21196 nodes, 7105 on queue, breadth at 12
24000 iterations,  21958 nodes, 7291 on queue, breadth at 12
25000 iterations,  22836 nodes, 7586 on queue, breadth at 12
26000 iterations,  23564 nodes, 7726 on queue, breadth at 12
27000 iterations,  24295 nodes, 7859 on queue, breadth at 12
28000 iterations,  24988 nodes, 7957 on queue, breadth at 12
29000 iterations,  25833 nodes, 8210 on queue, breadth at 13
30000 iterations,  26815 nodes, 8568 on queue, breadth at 13
31000 iterations,  27541 nodes, 8711 on queue, breadth at 13
32000 iterations,  28305 nodes, 8913 on queue, breadth at 13
33000 iterations,  29157 nodes, 9134 on queue, breadth at 13
34000 iterations,  29901 nodes, 9310 on queue, breadth at 13
35000 iterations,  30612 nodes, 9448 on queue, breadth at 13
36000 iterations,  31304 nodes, 9571 on queue, breadth at 13
37000 iterations,  32132 nodes, 9791 on queue, breadth at 14
38000 iterations,  33057 nodes, 10146 on queue, breadth at 14
39000 iterations,  34012 nodes, 10462 on queue, breadth at 14
40000 iterations,  34723 nodes, 10567 on queue, breadth at 14
41000 iterations,  35469 nodes, 10738 on queue, breadth at 14
42000 iterations,  36365 nodes, 10997 on queue, breadth at 14
43000 iterations,  37114 nodes, 11137 on queue, breadth at 14
44000 iterations,  37896 nodes, 11292 on queue, breadth at 14
45000 iterations,  38602 nodes, 11412 on queue, breadth at 14
46000 iterations,  39256 nodes, 11454 on queue, breadth at 14
47000 iterations,  40094 nodes, 11719 on queue, breadth at 15
48000 iterations,  41012 nodes, 12072 on queue, breadth at 15
49000 iterations,  42126 nodes, 12592 on queue, breadth at 15
50000 iterations,  42809 nodes, 12641 on queue, breadth at 15
51000 iterations,  43592 nodes, 12819 on queue, breadth at 15
52000 iterations,  44543 nodes, 13157 on queue, breadth at 15
53000 iterations,  45382 nodes, 13386 on queue, breadth at 15
54000 iterations,  46131 nodes, 13566 on queue, breadth at 15
55000 iterations,  46901 nodes, 13707 on queue, breadth at 15
56000 iterations,  47602 nodes, 13835 on queue, breadth at 15
57000 iterations,  48245 nodes, 13933 on queue, breadth at 15
58000 iterations,  49105 nodes, 14196 on queue, breadth at 16
59000 iterations,  49895 nodes, 14358 on queue, breadth at 16
60000 iterations,  50725 nodes, 14590 on queue, breadth at 16
61000 iterations,  51738 nodes, 14996 on queue, breadth at 16
62000 iterations,  52614 nodes, 15337 on queue, breadth at 16
63000 iterations,  53370 nodes, 15525 on queue, breadth at 16
64000 iterations,  54069 nodes, 15653 on queue, breadth at 16
65000 iterations,  55039 nodes, 16023 on queue, breadth at 16
66000 iterations,  55862 nodes, 16272 on queue, breadth at 16
67000 iterations,  56574 nodes, 16380 on queue, breadth at 16
68000 iterations,  57361 nodes, 16590 on queue, breadth at 16
69000 iterations,  58189 nodes, 16784 on queue, breadth at 16
70000 iterations,  58844 nodes, 16847 on queue, breadth at 16
71000 iterations,  59522 nodes, 16898 on queue, breadth at 16
72000 iterations,  60358 nodes, 17131 on queue, breadth at 17
73000 iterations,  61174 nodes, 17364 on queue, breadth at 17
74000 iterations,  62001 nodes, 17654 on queue, breadth at 17
75000 iterations,  62732 nodes, 17827 on queue, breadth at 17
76000 iterations,  63731 nodes, 18238 on queue, breadth at 17
77000 iterations,  64602 nodes, 18543 on queue, breadth at 17
78000 iterations,  65260 nodes, 18593 on queue, breadth at 17
79000 iterations,  65903 nodes, 18643 on queue, breadth at 17
80000 iterations,  66709 nodes, 18858 on queue, breadth at 17
81000 iterations,  67576 nodes, 19134 on queue, breadth at 17
82000 iterations,  68376 nodes, 19342 on queue, breadth at 17
83000 iterations,  69081 nodes, 19460 on queue, breadth at 17
84000 iterations,  69901 nodes, 19726 on queue, breadth at 17
85000 iterations,  70633 nodes, 19853 on queue, breadth at 17
86000 iterations,  71347 nodes, 19965 on queue, breadth at 17
87000 iterations,  72024 nodes, 20026 on queue, breadth at 17
88000 iterations,  72738 nodes, 20158 on queue, breadth at 17
89000 iterations,  73510 nodes, 20322 on queue, breadth at 18
90000 iterations,  74466 nodes, 20672 on queue, breadth at 18
91000 iterations,  75308 nodes, 20916 on queue, breadth at 18
92000 iterations,  75961 nodes, 20995 on queue, breadth at 18
93000 iterations,  76563 nodes, 20984 on queue, breadth at 18
94000 iterations,  77528 nodes, 21419 on queue, breadth at 18
95000 iterations,  78313 nodes, 21674 on queue, breadth at 18
96000 iterations,  79112 nodes, 21874 on queue, breadth at 18
97000 iterations,  79813 nodes, 21944 on queue, breadth at 18
98000 iterations,  80434 nodes, 21974 on queue, breadth at 18
99000 iterations,  81250 nodes, 22210 on queue, breadth at 18
100000 iterations,  82114 nodes, 22488 on queue, breadth at 18
101000 iterations,  82872 nodes, 22723 on queue, breadth at 18
102000 iterations,  83667 nodes, 22880 on queue, breadth at 18
103000 iterations,  84370 nodes, 22969 on queue, breadth at 18
104000 iterations,  85109 nodes, 23145 on queue, breadth at 18
105000 iterations,  85823 nodes, 23238 on queue, breadth at 18
106000 iterations,  86569 nodes, 23352 on queue, breadth at 18
107000 iterations,  87260 nodes, 23445 on queue, breadth at 18
108000 iterations,  87939 nodes, 23538 on queue, breadth at 18
109000 iterations,  88695 nodes, 23671 on queue, breadth at 19
110000 iterations,  89712 nodes, 24165 on queue, breadth at 19
111000 iterations,  90486 nodes, 24351 on queue, breadth at 19
112000 iterations,  91282 nodes, 24588 on queue, breadth at 19
113000 iterations,  91883 nodes, 24574 on queue, breadth at 19
114000 iterations,  92567 nodes, 24629 on queue, breadth at 19
115000 iterations,  93352 nodes, 24863 on queue, breadth at 19
116000 iterations,  94128 nodes, 25148 on queue, breadth at 19
117000 iterations,  94934 nodes, 25448 on queue, breadth at 19
118000 iterations,  95594 nodes, 25578 on queue, breadth at 19
119000 iterations,  96293 nodes, 25652 on queue, breadth at 19
120000 iterations,  96942 nodes, 25704 on queue, breadth at 19
121000 iterations,  97757 nodes, 25910 on queue, breadth at 19
122000 iterations,  98544 nodes, 26123 on queue, breadth at 19
123000 iterations,  99310 nodes, 26349 on queue, breadth at 19
124000 iterations,  99978 nodes, 26451 on queue, breadth at 19
125000 iterations,  100824 nodes, 26667 on queue, breadth at 19
126000 iterations,  101521 nodes, 26784 on queue, breadth at 19
127000 iterations,  102279 nodes, 26921 on queue, breadth at 19
128000 iterations,  102979 nodes, 27044 on queue, breadth at 19
129000 iterations,  103718 nodes, 27191 on queue, breadth at 19
130000 iterations,  104365 nodes, 27214 on queue, breadth at 19
131000 iterations,  105049 nodes, 27286 on queue, breadth at 19
132000 iterations,  105783 nodes, 27436 on queue, breadth at 19
133000 iterations,  106610 nodes, 27685 on queue, breadth at 20
134000 iterations,  107512 nodes, 28027 on queue, breadth at 20
135000 iterations,  108177 nodes, 28119 on queue, breadth at 20
136000 iterations,  108941 nodes, 28310 on queue, breadth at 20
137000 iterations,  109621 nodes, 28389 on queue, breadth at 20
138000 iterations,  110281 nodes, 28451 on queue, breadth at 20
139000 iterations,  110946 nodes, 28561 on queue, breadth at 20
140000 iterations,  111642 nodes, 28755 on queue, breadth at 20
141000 iterations,  112318 nodes, 28928 on queue, breadth at 20
142000 iterations,  112946 nodes, 29036 on queue, breadth at 20
143000 iterations,  113516 nodes, 29082 on queue, breadth at 20
144000 iterations,  114193 nodes, 29143 on queue, breadth at 20
145000 iterations,  114809 nodes, 29126 on queue, breadth at 20
146000 iterations,  115476 nodes, 29194 on queue, breadth at 20
147000 iterations,  116333 nodes, 29423 on queue, breadth at 20
148000 iterations,  117039 nodes, 29617 on queue, breadth at 20
149000 iterations,  117749 nodes, 29773 on queue, breadth at 20
150000 iterations,  118366 nodes, 29867 on queue, breadth at 20
151000 iterations,  119152 nodes, 30025 on queue, breadth at 20
152000 iterations,  119854 nodes, 30124 on queue, breadth at 20
153000 iterations,  120568 nodes, 30250 on queue, breadth at 20
154000 iterations,  121323 nodes, 30438 on queue, breadth at 20
155000 iterations,  122073 nodes, 30596 on queue, breadth at 20
156000 iterations,  122901 nodes, 30804 on queue, breadth at 20
157000 iterations,  123451 nodes, 30796 on queue, breadth at 20
158000 iterations,  124125 nodes, 30874 on queue, breadth at 20
159000 iterations,  124822 nodes, 31038 on queue, breadth at 20
160000 iterations,  125499 nodes, 31142 on queue, breadth at 20
161000 iterations,  126332 nodes, 31405 on queue, breadth at 21
162000 iterations,  127089 nodes, 31629 on queue, breadth at 21
163000 iterations,  127689 nodes, 31681 on queue, breadth at 21
164000 iterations,  128431 nodes, 31814 on queue, breadth at 21
165000 iterations,  129172 nodes, 32001 on queue, breadth at 21
166000 iterations,  129823 nodes, 32037 on queue, breadth at 21
167000 iterations,  130510 nodes, 32117 on queue, breadth at 21
168000 iterations,  131169 nodes, 32256 on queue, breadth at 21
169000 iterations,  131777 nodes, 32363 on queue, breadth at 21
170000 iterations,  132424 nodes, 32471 on queue, breadth at 21
171000 iterations,  132972 nodes, 32491 on queue, breadth at 21
172000 iterations,  133395 nodes, 32381 on queue, breadth at 21
173000 iterations,  134077 nodes, 32454 on queue, breadth at 21
174000 iterations,  134740 nodes, 32495 on queue, breadth at 21
175000 iterations,  135352 nodes, 32541 on queue, breadth at 21
176000 iterations,  136216 nodes, 32856 on queue, breadth at 21
177000 iterations,  136802 nodes, 32901 on queue, breadth at 21
178000 iterations,  137469 nodes, 33020 on queue, breadth at 21
179000 iterations,  138144 nodes, 33177 on queue, breadth at 21
180000 iterations,  138713 nodes, 33201 on queue, breadth at 21
181000 iterations,  139505 nodes, 33365 on queue, breadth at 21
182000 iterations,  140173 nodes, 33439 on queue, breadth at 21
183000 iterations,  140871 nodes, 33533 on queue, breadth at 21
184000 iterations,  141620 nodes, 33755 on queue, breadth at 21
185000 iterations,  142307 nodes, 33880 on queue, breadth at 21
186000 iterations,  143117 nodes, 34089 on queue, breadth at 21
187000 iterations,  143856 nodes, 34213 on queue, breadth at 21
188000 iterations,  144400 nodes, 34171 on queue, breadth at 21
189000 iterations,  145030 nodes, 34180 on queue, breadth at 21
190000 iterations,  145632 nodes, 34241 on queue, breadth at 21
191000 iterations,  146292 nodes, 34306 on queue, breadth at 21
192000 iterations,  147128 nodes, 34540 on queue, breadth at 22
193000 iterations,  147870 nodes, 34753 on queue, breadth at 22
194000 iterations,  148482 nodes, 34814 on queue, breadth at 22
195000 iterations,  149030 nodes, 34802 on queue, breadth at 22
196000 iterations,  149725 nodes, 34973 on queue, breadth at 22
197000 iterations,  150439 nodes, 35098 on queue, breadth at 22
198000 iterations,  151075 nodes, 35137 on queue, breadth at 22
199000 iterations,  151716 nodes, 35140 on queue, breadth at 22
200000 iterations,  152386 nodes, 35278 on queue, breadth at 22
201000 iterations,  152931 nodes, 35318 on queue, breadth at 22
202000 iterations,  153584 nodes, 35431 on queue, breadth at 22
203000 iterations,  154195 nodes, 35532 on queue, breadth at 22
204000 iterations,  154681 nodes, 35402 on queue, breadth at 22
205000 iterations,  155265 nodes, 35348 on queue, breadth at 22
206000 iterations,  156007 nodes, 35455 on queue, breadth at 22
207000 iterations,  156638 nodes, 35464 on queue, breadth at 22
208000 iterations,  157233 nodes, 35449 on queue, breadth at 22
209000 iterations,  158018 nodes, 35709 on queue, breadth at 22
210000 iterations,  158563 nodes, 35685 on queue, breadth at 22
211000 iterations,  159204 nodes, 35721 on queue, breadth at 22
212000 iterations,  159866 nodes, 35801 on queue, breadth at 22
213000 iterations,  160499 nodes, 35878 on queue, breadth at 22
214000 iterations,  161223 nodes, 36000 on queue, breadth at 22
215000 iterations,  161902 nodes, 36041 on queue, breadth at 22
216000 iterations,  162516 nodes, 36032 on queue, breadth at 22
217000 iterations,  163215 nodes, 36131 on queue, breadth at 22
218000 iterations,  163894 nodes, 36241 on queue, breadth at 22
219000 iterations,  164573 nodes, 36311 on queue, breadth at 22
220000 iterations,  165385 nodes, 36529 on queue, breadth at 22
221000 iterations,  166147 nodes, 36687 on queue, breadth at 22
222000 iterations,  166734 nodes, 36642 on queue, breadth at 22
223000 iterations,  167320 nodes, 36622 on queue, breadth at 22
224000 iterations,  167951 nodes, 36696 on queue, breadth at 22
225000 iterations,  168577 nodes, 36714 on queue, breadth at 22
226000 iterations,  169281 nodes, 36794 on queue, breadth at 23
227000 iterations,  170066 nodes, 37078 on queue, breadth at 23
228000 iterations,  170747 nodes, 37175 on queue, breadth at 23
229000 iterations,  171306 nodes, 37110 on queue, breadth at 23
230000 iterations,  171850 nodes, 37022 on queue, breadth at 23
231000 iterations,  172450 nodes, 37108 on queue, breadth at 23
232000 iterations,  173211 nodes, 37267 on queue, breadth at 23
233000 iterations,  173836 nodes, 37278 on queue, breadth at 23
234000 iterations,  174499 nodes, 37273 on queue, breadth at 23
235000 iterations,  175183 nodes, 37425 on queue, breadth at 23
236000 iterations,  175722 nodes, 37423 on queue, breadth at 23
237000 iterations,  176325 nodes, 37527 on queue, breadth at 23
238000 iterations,  176929 nodes, 37617 on queue, breadth at 23
239000 iterations,  177521 nodes, 37684 on queue, breadth at 23
240000 iterations,  178040 nodes, 37597 on queue, breadth at 23
241000 iterations,  178790 nodes, 37736 on queue, breadth at 23
242000 iterations,  179470 nodes, 37780 on queue, breadth at 23
243000 iterations,  180100 nodes, 37812 on queue, breadth at 23
244000 iterations,  180714 nodes, 37886 on queue, breadth at 23
245000 iterations,  181463 nodes, 38054 on queue, breadth at 23
246000 iterations,  182084 nodes, 38106 on queue, breadth at 23
247000 iterations,  182709 nodes, 38198 on queue, breadth at 23
248000 iterations,  183346 nodes, 38308 on queue, breadth at 23
249000 iterations,  183945 nodes, 38378 on queue, breadth at 23
250000 iterations,  184722 nodes, 38541 on queue, breadth at 23
251000 iterations,  185360 nodes, 38551 on queue, breadth at 23
252000 iterations,  186019 nodes, 38582 on queue, breadth at 23
253000 iterations,  186715 nodes, 38684 on queue, breadth at 23
254000 iterations,  187372 nodes, 38831 on queue, breadth at 23
255000 iterations,  188039 nodes, 38868 on queue, breadth at 23
256000 iterations,  188841 nodes, 39079 on queue, breadth at 23
257000 iterations,  189523 nodes, 39238 on queue, breadth at 23
258000 iterations,  190202 nodes, 39347 on queue, breadth at 23
259000 iterations,  190766 nodes, 39314 on queue, breadth at 23
260000 iterations,  191387 nodes, 39342 on queue, breadth at 23
261000 iterations,  191981 nodes, 39382 on queue, breadth at 23
262000 iterations,  192614 nodes, 39425 on queue, breadth at 23
263000 iterations,  193386 nodes, 39601 on queue, breadth at 24
264000 iterations,  194048 nodes, 39707 on queue, breadth at 24
265000 iterations,  194770 nodes, 39839 on queue, breadth at 24
266000 iterations,  195338 nodes, 39800 on queue, breadth at 24
267000 iterations,  195912 nodes, 39743 on queue, breadth at 24
268000 iterations,  196497 nodes, 39806 on queue, breadth at 24
269000 iterations,  197156 nodes, 39922 on queue, breadth at 24
270000 iterations,  197788 nodes, 39917 on queue, breadth at 24
271000 iterations,  198474 nodes, 40002 on queue, breadth at 24
272000 iterations,  199171 nodes, 40098 on queue, breadth at 24
273000 iterations,  199779 nodes, 40138 on queue, breadth at 24
274000 iterations,  200312 nodes, 40131 on queue, breadth at 24
275000 iterations,  200940 nodes, 40227 on queue, breadth at 24
276000 iterations,  201477 nodes, 40242 on queue, breadth at 24
277000 iterations,  202010 nodes, 40204 on queue, breadth at 24
278000 iterations,  202692 nodes, 40248 on queue, breadth at 24
279000 iterations,  203349 nodes, 40285 on queue, breadth at 24
280000 iterations,  204030 nodes, 40360 on queue, breadth at 24
281000 iterations,  204655 nodes, 40385 on queue, breadth at 24
282000 iterations,  205292 nodes, 40508 on queue, breadth at 24
283000 iterations,  206026 nodes, 40702 on queue, breadth at 24
284000 iterations,  206586 nodes, 40726 on queue, breadth at 24
285000 iterations,  207198 nodes, 40808 on queue, breadth at 24
286000 iterations,  207811 nodes, 40846 on queue, breadth at 24
287000 iterations,  208428 nodes, 40939 on queue, breadth at 24
288000 iterations,  209111 nodes, 41003 on queue, breadth at 24
289000 iterations,  209815 nodes, 41124 on queue, breadth at 24
290000 iterations,  210437 nodes, 41151 on queue, breadth at 24
291000 iterations,  211086 nodes, 41192 on queue, breadth at 24
292000 iterations,  211718 nodes, 41251 on queue, breadth at 24
293000 iterations,  212402 nodes, 41407 on queue, breadth at 24
294000 iterations,  213124 nodes, 41529 on queue, breadth at 24
295000 iterations,  213948 nodes, 41723 on queue, breadth at 24
296000 iterations,  214432 nodes, 41609 on queue, breadth at 24
297000 iterations,  215138 nodes, 41713 on queue, breadth at 24
298000 iterations,  215680 nodes, 41670 on queue, breadth at 24
299000 iterations,  216292 nodes, 41684 on queue, breadth at 24
300000 iterations,  216946 nodes, 41738 on queue, breadth at 24
301000 iterations,  217575 nodes, 41808 on queue, breadth at 24
302000 iterations,  218202 nodes, 41812 on queue, breadth at 25
303000 iterations,  219001 nodes, 42016 on queue, breadth at 25
304000 iterations,  219589 nodes, 42019 on queue, breadth at 25
305000 iterations,  220317 nodes, 42161 on queue, breadth at 25
306000 iterations,  220895 nodes, 42106 on queue, breadth at 25
307000 iterations,  221514 nodes, 42106 on queue, breadth at 25
308000 iterations,  222056 nodes, 42133 on queue, breadth at 25
309000 iterations,  222638 nodes, 42120 on queue, breadth at 25
310000 iterations,  223259 nodes, 42109 on queue, breadth at 25
311000 iterations,  223919 nodes, 42140 on queue, breadth at 25
312000 iterations,  224668 nodes, 42311 on queue, breadth at 25
313000 iterations,  225294 nodes, 42430 on queue, breadth at 25
314000 iterations,  225841 nodes, 42422 on queue, breadth at 25
315000 iterations,  226504 nodes, 42533 on queue, breadth at 25
316000 iterations,  227040 nodes, 42542 on queue, breadth at 25
317000 iterations,  227627 nodes, 42566 on queue, breadth at 25
318000 iterations,  228286 nodes, 42623 on queue, breadth at 25
319000 iterations,  228988 nodes, 42725 on queue, breadth at 25
320000 iterations,  229682 nodes, 42781 on queue, breadth at 25
321000 iterations,  230244 nodes, 42758 on queue, breadth at 25
322000 iterations,  230867 nodes, 42828 on queue, breadth at 25
323000 iterations,  231523 nodes, 42976 on queue, breadth at 25
324000 iterations,  232067 nodes, 42979 on queue, breadth at 25
325000 iterations,  232714 nodes, 43054 on queue, breadth at 25
326000 iterations,  233117 nodes, 42847 on queue, breadth at 25
327000 iterations,  233725 nodes, 42955 on queue, breadth at 25
328000 iterations,  234322 nodes, 43052 on queue, breadth at 25
329000 iterations,  235016 nodes, 43137 on queue, breadth at 25
330000 iterations,  235685 nodes, 43245 on queue, breadth at 25
331000 iterations,  236252 nodes, 43191 on queue, breadth at 25
332000 iterations,  236966 nodes, 43277 on queue, breadth at 25
333000 iterations,  237573 nodes, 43284 on queue, breadth at 25
334000 iterations,  238154 nodes, 43364 on queue, breadth at 25
335000 iterations,  238785 nodes, 43456 on queue, breadth at 25
336000 iterations,  239687 nodes, 43775 on queue, breadth at 25
337000 iterations,  240279 nodes, 43745 on queue, breadth at 25
338000 iterations,  240838 nodes, 43681 on queue, breadth at 25
339000 iterations,  241479 nodes, 43740 on queue, breadth at 25
340000 iterations,  242023 nodes, 43673 on queue, breadth at 25
341000 iterations,  242600 nodes, 43641 on queue, breadth at 25
342000 iterations,  243205 nodes, 43712 on queue, breadth at 25
343000 iterations,  243929 nodes, 43851 on queue, breadth at 25
344000 iterations,  244636 nodes, 43889 on queue, breadth at 26
345000 iterations,  245458 nodes, 44122 on queue, breadth at 26
346000 iterations,  246190 nodes, 44250 on queue, breadth at 26
347000 iterations,  246872 nodes, 44372 on queue, breadth at 26
348000 iterations,  247480 nodes, 44380 on queue, breadth at 26
349000 iterations,  248146 nodes, 44381 on queue, breadth at 26
350000 iterations,  248707 nodes, 44428 on queue, breadth at 26
351000 iterations,  249248 nodes, 44397 on queue, breadth at 26
352000 iterations,  249834 nodes, 44352 on queue, breadth at 26
353000 iterations,  250565 nodes, 44416 on queue, breadth at 26
354000 iterations,  251265 nodes, 44539 on queue, breadth at 26
355000 iterations,  251935 nodes, 44677 on queue, breadth at 26
356000 iterations,  252514 nodes, 44724 on queue, breadth at 26
357000 iterations,  253070 nodes, 44765 on queue, breadth at 26
358000 iterations,  253690 nodes, 44813 on queue, breadth at 26
359000 iterations,  254254 nodes, 44825 on queue, breadth at 26
360000 iterations,  254803 nodes, 44769 on queue, breadth at 26
361000 iterations,  255496 nodes, 44863 on queue, breadth at 26
362000 iterations,  256230 nodes, 44988 on queue, breadth at 26
363000 iterations,  256886 nodes, 45036 on queue, breadth at 26
364000 iterations,  257503 nodes, 45059 on queue, breadth at 26
365000 iterations,  258067 nodes, 45058 on queue, breadth at 26
366000 iterations,  258690 nodes, 45191 on queue, breadth at 26
367000 iterations,  259222 nodes, 45167 on queue, breadth at 26
368000 iterations,  259854 nodes, 45220 on queue, breadth at 26
369000 iterations,  260348 nodes, 45116 on queue, breadth at 26
370000 iterations,  260878 nodes, 45090 on queue, breadth at 26
371000 iterations,  261435 nodes, 45082 on queue, breadth at 26
372000 iterations,  262058 nodes, 45101 on queue, breadth at 26
373000 iterations,  262727 nodes, 45205 on queue, breadth at 26
374000 iterations,  263402 nodes, 45232 on queue, breadth at 26
375000 iterations,  264051 nodes, 45306 on queue, breadth at 26
376000 iterations,  264794 nodes, 45437 on queue, breadth at 26
377000 iterations,  265370 nodes, 45480 on queue, breadth at 26
378000 iterations,  265919 nodes, 45467 on queue, breadth at 26
379000 iterations,  266558 nodes, 45514 on queue, breadth at 26
380000 iterations,  267421 nodes, 45821 on queue, breadth at 26
381000 iterations,  267955 nodes, 45714 on queue, breadth at 26
382000 iterations,  268616 nodes, 45833 on queue, breadth at 26
383000 iterations,  269163 nodes, 45793 on queue, breadth at 26
384000 iterations,  269752 nodes, 45760 on queue, breadth at 26
385000 iterations,  270405 nodes, 45809 on queue, breadth at 26
386000 iterations,  270974 nodes, 45829 on queue, breadth at 26
387000 iterations,  271652 nodes, 45906 on queue, breadth at 26
388000 iterations,  272405 nodes, 46048 on queue, breadth at 27
389000 iterations,  273196 nodes, 46372 on queue, breadth at 27
390000 iterations,  273948 nodes, 46570 on queue, breadth at 27
391000 iterations,  274667 nodes, 46752 on queue, breadth at 27
392000 iterations,  275221 nodes, 46713 on queue, breadth at 27
393000 iterations,  275915 nodes, 46811 on queue, breadth at 27
394000 iterations,  276515 nodes, 46853 on queue, breadth at 27
395000 iterations,  277089 nodes, 46870 on queue, breadth at 27
396000 iterations,  277623 nodes, 46767 on queue, breadth at 27
397000 iterations,  278328 nodes, 46848 on queue, breadth at 27
398000 iterations,  279054 nodes, 46921 on queue, breadth at 27
399000 iterations,  279766 nodes, 47049 on queue, breadth at 27
400000 iterations,  280354 nodes, 47086 on queue, breadth at 27
401000 iterations,  280895 nodes, 47082 on queue, breadth at 27
402000 iterations,  281466 nodes, 47112 on queue, breadth at 27
403000 iterations,  282107 nodes, 47167 on queue, breadth at 27
404000 iterations,  282715 nodes, 47208 on queue, breadth at 27
405000 iterations,  283268 nodes, 47156 on queue, breadth at 27
406000 iterations,  283971 nodes, 47224 on queue, breadth at 27
407000 iterations,  284700 nodes, 47368 on queue, breadth at 27
408000 iterations,  285367 nodes, 47403 on queue, breadth at 27
409000 iterations,  285995 nodes, 47451 on queue, breadth at 27
410000 iterations,  286584 nodes, 47471 on queue, breadth at 27
411000 iterations,  287157 nodes, 47510 on queue, breadth at 27
412000 iterations,  287695 nodes, 47498 on queue, breadth at 27
413000 iterations,  288219 nodes, 47466 on queue, breadth at 27
414000 iterations,  288809 nodes, 47436 on queue, breadth at 27
415000 iterations,  289395 nodes, 47445 on queue, breadth at 27
416000 iterations,  289931 nodes, 47454 on queue, breadth at 27
417000 iterations,  290545 nodes, 47428 on queue, breadth at 27
418000 iterations,  291248 nodes, 47574 on queue, breadth at 27
419000 iterations,  291917 nodes, 47587 on queue, breadth at 27
420000 iterations,  292555 nodes, 47592 on queue, breadth at 27
421000 iterations,  293234 nodes, 47645 on queue, breadth at 27
422000 iterations,  293931 nodes, 47790 on queue, breadth at 27
423000 iterations,  294461 nodes, 47794 on queue, breadth at 27
424000 iterations,  295057 nodes, 47811 on queue, breadth at 27
425000 iterations,  295825 nodes, 47971 on queue, breadth at 27
426000 iterations,  296560 nodes, 48142 on queue, breadth at 27
427000 iterations,  297044 nodes, 48021 on queue, breadth at 27
428000 iterations,  297621 nodes, 48013 on queue, breadth at 27
429000 iterations,  298164 nodes, 47958 on queue, breadth at 27
430000 iterations,  298854 nodes, 48024 on queue, breadth at 27
431000 iterations,  299463 nodes, 48059 on queue, breadth at 27
432000 iterations,  300074 nodes, 48097 on queue, breadth at 27
433000 iterations,  300707 nodes, 48136 on queue, breadth at 27
434000 iterations,  301409 nodes, 48243 on queue, breadth at 28
435000 iterations,  302051 nodes, 48385 on queue, breadth at 28
436000 iterations,  302738 nodes, 48551 on queue, breadth at 28
437000 iterations,  303459 nodes, 48686 on queue, breadth at 28
438000 iterations,  304056 nodes, 48688 on queue, breadth at 28
439000 iterations,  304683 nodes, 48661 on queue, breadth at 28
440000 iterations,  305343 nodes, 48701 on queue, breadth at 28
441000 iterations,  305946 nodes, 48730 on queue, breadth at 28
442000 iterations,  306469 nodes, 48680 on queue, breadth at 28
443000 iterations,  307076 nodes, 48648 on queue, breadth at 28
444000 iterations,  307840 nodes, 48762 on queue, breadth at 28
445000 iterations,  308567 nodes, 48851 on queue, breadth at 28
446000 iterations,  309260 nodes, 49020 on queue, breadth at 28
447000 iterations,  309798 nodes, 49022 on queue, breadth at 28
448000 iterations,  310348 nodes, 49052 on queue, breadth at 28
449000 iterations,  310914 nodes, 49062 on queue, breadth at 28
450000 iterations,  311547 nodes, 49139 on queue, breadth at 28
451000 iterations,  312168 nodes, 49223 on queue, breadth at 28
452000 iterations,  312759 nodes, 49202 on queue, breadth at 28
453000 iterations,  313419 nodes, 49243 on queue, breadth at 28
454000 iterations,  314133 nodes, 49341 on queue, breadth at 28
455000 iterations,  314785 nodes, 49344 on queue, breadth at 28
456000 iterations,  315437 nodes, 49411 on queue, breadth at 28
457000 iterations,  316039 nodes, 49419 on queue, breadth at 28
458000 iterations,  316556 nodes, 49407 on queue, breadth at 28
459000 iterations,  317176 nodes, 49502 on queue, breadth at 28
460000 iterations,  317724 nodes, 49465 on queue, breadth at 28
461000 iterations,  318220 nodes, 49376 on queue, breadth at 28
462000 iterations,  318869 nodes, 49407 on queue, breadth at 28
463000 iterations,  319379 nodes, 49365 on queue, breadth at 28
464000 iterations,  319966 nodes, 49404 on queue, breadth at 28
465000 iterations,  320650 nodes, 49473 on queue, breadth at 28
466000 iterations,  321383 nodes, 49609 on queue, breadth at 28
467000 iterations,  322037 nodes, 49616 on queue, breadth at 28
468000 iterations,  322719 nodes, 49691 on queue, breadth at 28
469000 iterations,  323425 nodes, 49779 on queue, breadth at 28
470000 iterations,  324060 nodes, 49853 on queue, breadth at 28
471000 iterations,  324599 nodes, 49827 on queue, breadth at 28
472000 iterations,  325211 nodes, 49825 on queue, breadth at 28
473000 iterations,  325907 nodes, 49893 on queue, breadth at 28
474000 iterations,  326607 nodes, 50026 on queue, breadth at 28
475000 iterations,  327096 nodes, 49873 on queue, breadth at 28
476000 iterations,  327618 nodes, 49770 on queue, breadth at 28
477000 iterations,  328209 nodes, 49765 on queue, breadth at 28
478000 iterations,  328927 nodes, 49874 on queue, breadth at 28
479000 iterations,  329512 nodes, 49883 on queue, breadth at 28
480000 iterations,  330128 nodes, 49944 on queue, breadth at 28
481000 iterations,  330725 nodes, 49938 on queue, breadth at 28
482000 iterations,  331406 nodes, 50011 on queue, breadth at 29
483000 iterations,  331923 nodes, 49933 on queue, breadth at 29
484000 iterations,  332482 nodes, 49963 on queue, breadth at 29
485000 iterations,  333198 nodes, 50158 on queue, breadth at 29
486000 iterations,  333835 nodes, 50210 on queue, breadth at 29
487000 iterations,  334431 nodes, 50183 on queue, breadth at 29
488000 iterations,  334916 nodes, 50091 on queue, breadth at 29
489000 iterations,  335670 nodes, 50221 on queue, breadth at 29
490000 iterations,  336244 nodes, 50240 on queue, breadth at 29
491000 iterations,  336773 nodes, 50166 on queue, breadth at 29
492000 iterations,  337496 nodes, 50289 on queue, breadth at 29
493000 iterations,  338240 nodes, 50425 on queue, breadth at 29
494000 iterations,  338966 nodes, 50493 on queue, breadth at 29
495000 iterations,  339608 nodes, 50598 on queue, breadth at 29
496000 iterations,  340154 nodes, 50587 on queue, breadth at 29
497000 iterations,  340709 nodes, 50585 on queue, breadth at 29
498000 iterations,  341288 nodes, 50632 on queue, breadth at 29
499000 iterations,  341924 nodes, 50698 on queue, breadth at 29
500000 iterations,  342501 nodes, 50700 on queue, breadth at 29
501000 iterations,  343078 nodes, 50681 on queue, breadth at 29
502000 iterations,  343716 nodes, 50711 on queue, breadth at 29
503000 iterations,  344381 nodes, 50800 on queue, breadth at 29
504000 iterations,  345078 nodes, 50875 on queue, breadth at 29
505000 iterations,  345673 nodes, 50861 on queue, breadth at 29
506000 iterations,  346328 nodes, 50907 on queue, breadth at 29
507000 iterations,  346850 nodes, 50859 on queue, breadth at 29
508000 iterations,  347469 nodes, 50944 on queue, breadth at 29
509000 iterations,  348002 nodes, 50931 on queue, breadth at 29
510000 iterations,  348538 nodes, 50869 on queue, breadth at 29
511000 iterations,  349178 nodes, 50842 on queue, breadth at 29
512000 iterations,  349743 nodes, 50826 on queue, breadth at 29
513000 iterations,  350299 nodes, 50814 on queue, breadth at 29
514000 iterations,  350911 nodes, 50817 on queue, breadth at 29
515000 iterations,  351543 nodes, 50887 on queue, breadth at 29
516000 iterations,  352234 nodes, 50956 on queue, breadth at 29
517000 iterations,  352950 nodes, 51040 on queue, breadth at 29
518000 iterations,  353647 nodes, 51118 on queue, breadth at 29
519000 iterations,  354308 nodes, 51204 on queue, breadth at 29
520000 iterations,  354922 nodes, 51290 on queue, breadth at 29
521000 iterations,  355482 nodes, 51294 on queue, breadth at 29
522000 iterations,  356053 nodes, 51268 on queue, breadth at 29
523000 iterations,  356782 nodes, 51379 on queue, breadth at 29
524000 iterations,  357366 nodes, 51416 on queue, breadth at 29
525000 iterations,  357925 nodes, 51347 on queue, breadth at 29
526000 iterations,  358467 nodes, 51276 on queue, breadth at 29
527000 iterations,  359087 nodes, 51294 on queue, breadth at 29
528000 iterations,  359712 nodes, 51322 on queue, breadth at 29
529000 iterations,  360324 nodes, 51364 on queue, breadth at 29
530000 iterations,  360884 nodes, 51374 on queue, breadth at 29
531000 iterations,  361453 nodes, 51365 on queue, breadth at 29
532000 iterations,  362128 nodes, 51448 on queue, breadth at 30
533000 iterations,  362616 nodes, 51345 on queue, breadth at 30
534000 iterations,  363057 nodes, 51245 on queue, breadth at 30
535000 iterations,  363698 nodes, 51388 on queue, breadth at 30
536000 iterations,  364323 nodes, 51395 on queue, breadth at 30
537000 iterations,  364923 nodes, 51418 on queue, breadth at 30
538000 iterations,  365378 nodes, 51254 on queue, breadth at 30
539000 iterations,  366119 nodes, 51328 on queue, breadth at 30
540000 iterations,  366712 nodes, 51357 on queue, breadth at 30
541000 iterations,  367274 nodes, 51327 on queue, breadth at 30
542000 iterations,  368005 nodes, 51449 on queue, breadth at 30
543000 iterations,  368724 nodes, 51535 on queue, breadth at 30
544000 iterations,  369395 nodes, 51600 on queue, breadth at 30
545000 iterations,  370036 nodes, 51653 on queue, breadth at 30
546000 iterations,  370675 nodes, 51734 on queue, breadth at 30
547000 iterations,  371249 nodes, 51763 on queue, breadth at 30
548000 iterations,  371802 nodes, 51773 on queue, breadth at 30
549000 iterations,  372400 nodes, 51806 on queue, breadth at 30
550000 iterations,  372922 nodes, 51733 on queue, breadth at 30
551000 iterations,  373512 nodes, 51751 on queue, breadth at 30
552000 iterations,  374196 nodes, 51828 on queue, breadth at 30
553000 iterations,  374786 nodes, 51775 on queue, breadth at 30
554000 iterations,  375455 nodes, 51886 on queue, breadth at 30
555000 iterations,  376150 nodes, 51955 on queue, breadth at 30
556000 iterations,  376723 nodes, 51940 on queue, breadth at 30
557000 iterations,  377321 nodes, 51936 on queue, breadth at 30
558000 iterations,  377884 nodes, 51921 on queue, breadth at 30
559000 iterations,  378528 nodes, 52035 on queue, breadth at 30
560000 iterations,  379077 nodes, 52046 on queue, breadth at 30
561000 iterations,  379602 nodes, 51975 on queue, breadth at 30
562000 iterations,  380277 nodes, 52053 on queue, breadth at 30
563000 iterations,  380859 nodes, 52102 on queue, breadth at 30
564000 iterations,  381407 nodes, 52092 on queue, breadth at 30
565000 iterations,  382046 nodes, 52108 on queue, breadth at 30
566000 iterations,  382657 nodes, 52160 on queue, breadth at 30
567000 iterations,  383388 nodes, 52254 on queue, breadth at 30
568000 iterations,  384100 nodes, 52379 on queue, breadth at 30
569000 iterations,  384805 nodes, 52443 on queue, breadth at 30
570000 iterations,  385410 nodes, 52464 on queue, breadth at 30
571000 iterations,  386054 nodes, 52531 on queue, breadth at 30
572000 iterations,  386596 nodes, 52529 on queue, breadth at 30
573000 iterations,  387207 nodes, 52546 on queue, breadth at 30
574000 iterations,  387950 nodes, 52630 on queue, breadth at 30
575000 iterations,  388441 nodes, 52523 on queue, breadth at 30
576000 iterations,  388987 nodes, 52492 on queue, breadth at 30
577000 iterations,  389497 nodes, 52398 on queue, breadth at 30
578000 iterations,  390090 nodes, 52374 on queue, breadth at 30
579000 iterations,  390681 nodes, 52315 on queue, breadth at 30
580000 iterations,  391322 nodes, 52357 on queue, breadth at 30
581000 iterations,  391835 nodes, 52286 on queue, breadth at 30
582000 iterations,  392348 nodes, 52224 on queue, breadth at 30
583000 iterations,  392959 nodes, 52255 on queue, breadth at 31
584000 iterations,  393469 nodes, 52099 on queue, breadth at 31
585000 iterations,  393943 nodes, 52023 on queue, breadth at 31
586000 iterations,  394414 nodes, 51865 on queue, breadth at 31
587000 iterations,  395011 nodes, 51901 on queue, breadth at 31
588000 iterations,  395576 nodes, 51883 on queue, breadth at 31
589000 iterations,  396158 nodes, 51805 on queue, breadth at 31
590000 iterations,  396825 nodes, 51817 on queue, breadth at 31
591000 iterations,  397459 nodes, 51884 on queue, breadth at 31
592000 iterations,  398054 nodes, 51932 on queue, breadth at 31
593000 iterations,  398801 nodes, 52067 on queue, breadth at 31
594000 iterations,  399433 nodes, 52071 on queue, breadth at 31
595000 iterations,  400110 nodes, 52109 on queue, breadth at 31
596000 iterations,  400778 nodes, 52113 on queue, breadth at 31
597000 iterations,  401391 nodes, 52186 on queue, breadth at 31
598000 iterations,  401967 nodes, 52236 on queue, breadth at 31
599000 iterations,  402555 nodes, 52302 on queue, breadth at 31
600000 iterations,  403108 nodes, 52311 on queue, breadth at 31
601000 iterations,  403727 nodes, 52355 on queue, breadth at 31
602000 iterations,  404285 nodes, 52376 on queue, breadth at 31
603000 iterations,  404926 nodes, 52426 on queue, breadth at 31
604000 iterations,  405638 nodes, 52529 on queue, breadth at 31
605000 iterations,  406295 nodes, 52544 on queue, breadth at 31
606000 iterations,  406907 nodes, 52573 on queue, breadth at 31
607000 iterations,  407514 nodes, 52600 on queue, breadth at 31
608000 iterations,  408085 nodes, 52569 on queue, breadth at 31
609000 iterations,  408683 nodes, 52562 on queue, breadth at 31
610000 iterations,  409222 nodes, 52581 on queue, breadth at 31
611000 iterations,  409870 nodes, 52700 on queue, breadth at 31
612000 iterations,  410411 nodes, 52672 on queue, breadth at 31
613000 iterations,  410937 nodes, 52585 on queue, breadth at 31
614000 iterations,  411601 nodes, 52637 on queue, breadth at 31
615000 iterations,  412184 nodes, 52664 on queue, breadth at 31
616000 iterations,  412729 nodes, 52642 on queue, breadth at 31
617000 iterations,  413395 nodes, 52688 on queue, breadth at 31
618000 iterations,  413984 nodes, 52729 on queue, breadth at 31
619000 iterations,  414677 nodes, 52818 on queue, breadth at 31
620000 iterations,  415345 nodes, 52862 on queue, breadth at 31
621000 iterations,  416019 nodes, 52937 on queue, breadth at 31
622000 iterations,  416713 nodes, 53018 on queue, breadth at 31
623000 iterations,  417345 nodes, 53074 on queue, breadth at 31
624000 iterations,  417954 nodes, 53117 on queue, breadth at 31
625000 iterations,  418541 nodes, 53137 on queue, breadth at 31
626000 iterations,  419122 nodes, 53128 on queue, breadth at 31
627000 iterations,  419818 nodes, 53195 on queue, breadth at 31
628000 iterations,  420337 nodes, 53073 on queue, breadth at 31
629000 iterations,  420895 nodes, 52976 on queue, breadth at 31
630000 iterations,  421469 nodes, 52938 on queue, breadth at 31
631000 iterations,  422098 nodes, 52957 on queue, breadth at 31
632000 iterations,  422774 nodes, 53052 on queue, breadth at 31
633000 iterations,  423347 nodes, 53027 on queue, breadth at 31
634000 iterations,  423844 nodes, 52920 on queue, breadth at 31
635000 iterations,  424352 nodes, 52831 on queue, breadth at 31
636000 iterations,  424994 nodes, 52813 on queue, breadth at 32
637000 iterations,  425375 nodes, 52628 on queue, breadth at 32
638000 iterations,  425986 nodes, 52626 on queue, breadth at 32
639000 iterations,  426553 nodes, 52617 on queue, breadth at 32
640000 iterations,  427115 nodes, 52570 on queue, breadth at 32
641000 iterations,  427707 nodes, 52538 on queue, breadth at 32
642000 iterations,  428546 nodes, 52726 on queue, breadth at 32
643000 iterations,  429146 nodes, 52762 on queue, breadth at 32
644000 iterations,  429728 nodes, 52741 on queue, breadth at 32
645000 iterations,  430440 nodes, 52868 on queue, breadth at 32
646000 iterations,  431058 nodes, 52835 on queue, breadth at 32
647000 iterations,  431784 nodes, 52948 on queue, breadth at 32
648000 iterations,  432465 nodes, 52973 on queue, breadth at 32
649000 iterations,  433063 nodes, 53012 on queue, breadth at 32
650000 iterations,  433618 nodes, 53050 on queue, breadth at 32
651000 iterations,  434158 nodes, 53015 on queue, breadth at 32
652000 iterations,  434691 nodes, 53020 on queue, breadth at 32
653000 iterations,  435290 nodes, 53082 on queue, breadth at 32
654000 iterations,  435935 nodes, 53132 on queue, breadth at 32
655000 iterations,  436522 nodes, 53140 on queue, breadth at 32
656000 iterations,  437255 nodes, 53288 on queue, breadth at 32
657000 iterations,  437894 nodes, 53268 on queue, breadth at 32
658000 iterations,  438523 nodes, 53299 on queue, breadth at 32
659000 iterations,  439107 nodes, 53275 on queue, breadth at 32
660000 iterations,  439736 nodes, 53235 on queue, breadth at 32
661000 iterations,  440357 nodes, 53290 on queue, breadth at 32
662000 iterations,  440881 nodes, 53232 on queue, breadth at 32
663000 iterations,  441372 nodes, 53199 on queue, breadth at 32
664000 iterations,  441991 nodes, 53276 on queue, breadth at 32
665000 iterations,  442554 nodes, 53264 on queue, breadth at 32
666000 iterations,  443175 nodes, 53243 on queue, breadth at 32
667000 iterations,  443771 nodes, 53230 on queue, breadth at 32
668000 iterations,  444342 nodes, 53265 on queue, breadth at 32
669000 iterations,  444928 nodes, 53239 on queue, breadth at 32
670000 iterations,  445627 nodes, 53301 on queue, breadth at 32
671000 iterations,  446141 nodes, 53236 on queue, breadth at 32
672000 iterations,  446762 nodes, 53253 on queue, breadth at 32
673000 iterations,  447383 nodes, 53283 on queue, breadth at 32
674000 iterations,  448077 nodes, 53344 on queue, breadth at 32
675000 iterations,  448752 nodes, 53383 on queue, breadth at 32
676000 iterations,  449355 nodes, 53405 on queue, breadth at 32
677000 iterations,  449961 nodes, 53463 on queue, breadth at 32
678000 iterations,  450583 nodes, 53522 on queue, breadth at 32
679000 iterations,  451139 nodes, 53483 on queue, breadth at 32
680000 iterations,  451789 nodes, 53498 on queue, breadth at 32
681000 iterations,  452467 nodes, 53551 on queue, breadth at 32
682000 iterations,  453052 nodes, 53548 on queue, breadth at 32
683000 iterations,  453632 nodes, 53527 on queue, breadth at 32
684000 iterations,  454218 nodes, 53513 on queue, breadth at 32
685000 iterations,  454844 nodes, 53516 on queue, breadth at 32
686000 iterations,  455407 nodes, 53451 on queue, breadth at 32
687000 iterations,  455962 nodes, 53427 on queue, breadth at 32
688000 iterations,  456493 nodes, 53323 on queue, breadth at 32
689000 iterations,  457157 nodes, 53350 on queue, breadth at 33
690000 iterations,  457703 nodes, 53278 on queue, breadth at 33
691000 iterations,  458286 nodes, 53273 on queue, breadth at 33
692000 iterations,  458870 nodes, 53252 on queue, breadth at 33
693000 iterations,  459353 nodes, 53119 on queue, breadth at 33
694000 iterations,  460094 nodes, 53218 on queue, breadth at 33
695000 iterations,  460908 nodes, 53394 on queue, breadth at 33
696000 iterations,  461488 nodes, 53394 on queue, breadth at 33
697000 iterations,  462121 nodes, 53428 on queue, breadth at 33
698000 iterations,  462676 nodes, 53416 on queue, breadth at 33
699000 iterations,  463350 nodes, 53436 on queue, breadth at 33
700000 iterations,  463984 nodes, 53456 on queue, breadth at 33
701000 iterations,  464613 nodes, 53432 on queue, breadth at 33
702000 iterations,  465198 nodes, 53465 on queue, breadth at 33
703000 iterations,  465743 nodes, 53444 on queue, breadth at 33
704000 iterations,  466304 nodes, 53451 on queue, breadth at 33
705000 iterations,  466819 nodes, 53429 on queue, breadth at 33
706000 iterations,  467398 nodes, 53420 on queue, breadth at 33
707000 iterations,  468106 nodes, 53540 on queue, breadth at 33
708000 iterations,  468709 nodes, 53611 on queue, breadth at 33
709000 iterations,  469359 nodes, 53676 on queue, breadth at 33
710000 iterations,  470037 nodes, 53786 on queue, breadth at 33
711000 iterations,  470705 nodes, 53801 on queue, breadth at 33
712000 iterations,  471274 nodes, 53768 on queue, breadth at 33
713000 iterations,  471922 nodes, 53786 on queue, breadth at 33
714000 iterations,  472539 nodes, 53844 on queue, breadth at 33
715000 iterations,  473161 nodes, 53864 on queue, breadth at 33
716000 iterations,  473640 nodes, 53795 on queue, breadth at 33
717000 iterations,  474195 nodes, 53795 on queue, breadth at 33
718000 iterations,  474781 nodes, 53813 on queue, breadth at 33
719000 iterations,  475443 nodes, 53826 on queue, breadth at 33
720000 iterations,  476183 nodes, 53916 on queue, breadth at 33
721000 iterations,  476788 nodes, 53966 on queue, breadth at 33
722000 iterations,  477374 nodes, 54000 on queue, breadth at 33
723000 iterations,  478080 nodes, 54094 on queue, breadth at 33
724000 iterations,  478691 nodes, 54139 on queue, breadth at 33
725000 iterations,  479234 nodes, 54067 on queue, breadth at 33
726000 iterations,  479841 nodes, 54073 on queue, breadth at 33
727000 iterations,  480430 nodes, 54087 on queue, breadth at 33
728000 iterations,  481072 nodes, 54123 on queue, breadth at 33
729000 iterations,  481679 nodes, 54114 on queue, breadth at 33
730000 iterations,  482346 nodes, 54199 on queue, breadth at 33
731000 iterations,  482910 nodes, 54248 on queue, breadth at 33
732000 iterations,  483525 nodes, 54276 on queue, breadth at 33
733000 iterations,  484129 nodes, 54237 on queue, breadth at 33
734000 iterations,  484837 nodes, 54328 on queue, breadth at 33
735000 iterations,  485425 nodes, 54306 on queue, breadth at 33
736000 iterations,  485997 nodes, 54240 on queue, breadth at 33
737000 iterations,  486581 nodes, 54235 on queue, breadth at 33
738000 iterations,  487243 nodes, 54251 on queue, breadth at 33
739000 iterations,  487847 nodes, 54282 on queue, breadth at 33
740000 iterations,  488364 nodes, 54260 on queue, breadth at 33
741000 iterations,  488943 nodes, 54214 on queue, breadth at 33
742000 iterations,  489570 nodes, 54236 on queue, breadth at 34
743000 iterations,  490148 nodes, 54202 on queue, breadth at 34
744000 iterations,  490817 nodes, 54311 on queue, breadth at 34
745000 iterations,  491333 nodes, 54256 on queue, breadth at 34
746000 iterations,  491810 nodes, 54103 on queue, breadth at 34
747000 iterations,  492474 nodes, 54113 on queue, breadth at 34
748000 iterations,  493286 nodes, 54246 on queue, breadth at 34
749000 iterations,  493935 nodes, 54291 on queue, breadth at 34
750000 iterations,  494564 nodes, 54359 on queue, breadth at 34
751000 iterations,  495103 nodes, 54251 on queue, breadth at 34
752000 iterations,  495677 nodes, 54219 on queue, breadth at 34
753000 iterations,  496316 nodes, 54228 on queue, breadth at 34
754000 iterations,  496921 nodes, 54211 on queue, breadth at 34
755000 iterations,  497544 nodes, 54214 on queue, breadth at 34
756000 iterations,  498043 nodes, 54205 on queue, breadth at 34
757000 iterations,  498562 nodes, 54208 on queue, breadth at 34
758000 iterations,  499094 nodes, 54166 on queue, breadth at 34
759000 iterations,  499632 nodes, 54184 on queue, breadth at 34
760000 iterations,  500298 nodes, 54306 on queue, breadth at 34
761000 iterations,  500898 nodes, 54353 on queue, breadth at 34
762000 iterations,  501463 nodes, 54328 on queue, breadth at 34
763000 iterations,  502180 nodes, 54456 on queue, breadth at 34
764000 iterations,  502739 nodes, 54388 on queue, breadth at 34
765000 iterations,  503447 nodes, 54474 on queue, breadth at 34
766000 iterations,  504051 nodes, 54458 on queue, breadth at 34
767000 iterations,  504739 nodes, 54528 on queue, breadth at 34
768000 iterations,  505311 nodes, 54502 on queue, breadth at 34
769000 iterations,  505954 nodes, 54544 on queue, breadth at 34
770000 iterations,  506490 nodes, 54552 on queue, breadth at 34
771000 iterations,  507033 nodes, 54524 on queue, breadth at 34
772000 iterations,  507594 nodes, 54538 on queue, breadth at 34
773000 iterations,  508332 nodes, 54681 on queue, breadth at 34
774000 iterations,  509033 nodes, 54793 on queue, breadth at 34
775000 iterations,  509642 nodes, 54830 on queue, breadth at 34
776000 iterations,  510205 nodes, 54837 on queue, breadth at 34
777000 iterations,  510922 nodes, 54905 on queue, breadth at 34
778000 iterations,  511527 nodes, 54949 on queue, breadth at 34
779000 iterations,  512071 nodes, 54903 on queue, breadth at 34
780000 iterations,  512647 nodes, 54891 on queue, breadth at 34
781000 iterations,  513216 nodes, 54828 on queue, breadth at 34
782000 iterations,  513723 nodes, 54734 on queue, breadth at 34
783000 iterations,  514374 nodes, 54828 on queue, breadth at 34
784000 iterations,  515053 nodes, 54916 on queue, breadth at 34
785000 iterations,  515593 nodes, 54883 on queue, breadth at 34
786000 iterations,  516182 nodes, 54889 on queue, breadth at 34
787000 iterations,  516808 nodes, 54907 on queue, breadth at 34
788000 iterations,  517427 nodes, 54916 on queue, breadth at 34
789000 iterations,  518088 nodes, 54978 on queue, breadth at 34
790000 iterations,  518700 nodes, 54991 on queue, breadth at 34
791000 iterations,  519290 nodes, 54943 on queue, breadth at 34
792000 iterations,  519922 nodes, 54956 on queue, breadth at 34
793000 iterations,  520554 nodes, 54992 on queue, breadth at 34
794000 iterations,  521052 nodes, 54910 on queue, breadth at 34
795000 iterations,  521687 nodes, 54948 on queue, breadth at 34
796000 iterations,  522230 nodes, 54888 on queue, breadth at 35
797000 iterations,  522896 nodes, 54888 on queue, breadth at 35
798000 iterations,  523551 nodes, 54963 on queue, breadth at 35
799000 iterations,  524092 nodes, 54913 on queue, breadth at 35
800000 iterations,  524611 nodes, 54773 on queue, breadth at 35
801000 iterations,  525287 nodes, 54845 on queue, breadth at 35
802000 iterations,  526111 nodes, 55050 on queue, breadth at 35
803000 iterations,  526824 nodes, 55119 on queue, breadth at 35
804000 iterations,  527378 nodes, 55095 on queue, breadth at 35
805000 iterations,  528007 nodes, 55091 on queue, breadth at 35
806000 iterations,  528502 nodes, 54939 on queue, breadth at 35
807000 iterations,  529142 nodes, 54943 on queue, breadth at 35
808000 iterations,  529837 nodes, 54999 on queue, breadth at 35
809000 iterations,  530444 nodes, 54954 on queue, breadth at 35
810000 iterations,  530956 nodes, 54868 on queue, breadth at 35
811000 iterations,  531481 nodes, 54821 on queue, breadth at 35
812000 iterations,  532035 nodes, 54783 on queue, breadth at 35
813000 iterations,  532520 nodes, 54713 on queue, breadth at 35
814000 iterations,  533122 nodes, 54755 on queue, breadth at 35
815000 iterations,  533731 nodes, 54793 on queue, breadth at 35
816000 iterations,  534257 nodes, 54760 on queue, breadth at 35
817000 iterations,  534855 nodes, 54771 on queue, breadth at 35
818000 iterations,  535530 nodes, 54844 on queue, breadth at 35
819000 iterations,  536210 nodes, 54885 on queue, breadth at 35
820000 iterations,  536878 nodes, 54939 on queue, breadth at 35
821000 iterations,  537536 nodes, 54990 on queue, breadth at 35
822000 iterations,  538085 nodes, 54927 on queue, breadth at 35
823000 iterations,  538696 nodes, 54923 on queue, breadth at 35
824000 iterations,  539323 nodes, 54973 on queue, breadth at 35
825000 iterations,  539893 nodes, 55003 on queue, breadth at 35
826000 iterations,  540381 nodes, 54920 on queue, breadth at 35
827000 iterations,  540952 nodes, 54886 on queue, breadth at 35
828000 iterations,  541704 nodes, 55039 on queue, breadth at 35
829000 iterations,  542325 nodes, 55037 on queue, breadth at 35
830000 iterations,  542884 nodes, 55067 on queue, breadth at 35
831000 iterations,  543430 nodes, 55048 on queue, breadth at 35
832000 iterations,  544188 nodes, 55174 on queue, breadth at 35
833000 iterations,  544724 nodes, 55127 on queue, breadth at 35
834000 iterations,  545255 nodes, 55003 on queue, breadth at 35
835000 iterations,  545817 nodes, 54962 on queue, breadth at 35
836000 iterations,  546399 nodes, 54899 on queue, breadth at 35
837000 iterations,  546905 nodes, 54781 on queue, breadth at 35
838000 iterations,  547527 nodes, 54810 on queue, breadth at 35
839000 iterations,  548196 nodes, 54887 on queue, breadth at 35
840000 iterations,  548748 nodes, 54878 on queue, breadth at 35
841000 iterations,  549345 nodes, 54904 on queue, breadth at 35
842000 iterations,  549970 nodes, 54954 on queue, breadth at 35
843000 iterations,  550661 nodes, 54986 on queue, breadth at 35
844000 iterations,  551245 nodes, 54952 on queue, breadth at 35
845000 iterations,  551800 nodes, 54920 on queue, breadth at 35
846000 iterations,  552465 nodes, 54951 on queue, breadth at 35
847000 iterations,  553076 nodes, 54965 on queue, breadth at 35
848000 iterations,  553681 nodes, 54974 on queue, breadth at 35
849000 iterations,  554203 nodes, 54928 on queue, breadth at 35
850000 iterations,  554828 nodes, 54958 on queue, breadth at 35
851000 iterations,  555387 nodes, 54905 on queue, breadth at 36
852000 iterations,  556156 nodes, 55047 on queue, breadth at 36
853000 iterations,  556806 nodes, 55087 on queue, breadth at 36
854000 iterations,  557307 nodes, 55009 on queue, breadth at 36
855000 iterations,  557922 nodes, 55034 on queue, breadth at 36
856000 iterations,  558630 nodes, 55134 on queue, breadth at 36
857000 iterations,  559393 nodes, 55265 on queue, breadth at 36
858000 iterations,  560107 nodes, 55295 on queue, breadth at 36
859000 iterations,  560625 nodes, 55268 on queue, breadth at 36
860000 iterations,  561267 nodes, 55308 on queue, breadth at 36
861000 iterations,  561788 nodes, 55184 on queue, breadth at 36
862000 iterations,  562451 nodes, 55225 on queue, breadth at 36
863000 iterations,  563083 nodes, 55210 on queue, breadth at 36
864000 iterations,  563669 nodes, 55151 on queue, breadth at 36
865000 iterations,  564202 nodes, 55170 on queue, breadth at 36
866000 iterations,  564746 nodes, 55178 on queue, breadth at 36
867000 iterations,  565345 nodes, 55172 on queue, breadth at 36
868000 iterations,  565858 nodes, 55130 on queue, breadth at 36
869000 iterations,  566371 nodes, 55090 on queue, breadth at 36
870000 iterations,  566916 nodes, 55090 on queue, breadth at 36
871000 iterations,  567464 nodes, 55030 on queue, breadth at 36
872000 iterations,  568048 nodes, 55055 on queue, breadth at 36
873000 iterations,  568615 nodes, 54997 on queue, breadth at 36
874000 iterations,  569358 nodes, 55096 on queue, breadth at 36
875000 iterations,  569954 nodes, 55098 on queue, breadth at 36
876000 iterations,  570598 nodes, 55146 on queue, breadth at 36
877000 iterations,  571184 nodes, 55115 on queue, breadth at 36
878000 iterations,  571845 nodes, 55169 on queue, breadth at 36
879000 iterations,  572417 nodes, 55220 on queue, breadth at 36
880000 iterations,  572950 nodes, 55172 on queue, breadth at 36
881000 iterations,  573485 nodes, 55137 on queue, breadth at 36
882000 iterations,  574121 nodes, 55187 on queue, breadth at 36
883000 iterations,  574833 nodes, 55268 on queue, breadth at 36
884000 iterations,  575439 nodes, 55254 on queue, breadth at 36
885000 iterations,  575981 nodes, 55229 on queue, breadth at 36
886000 iterations,  576542 nodes, 55209 on queue, breadth at 36
887000 iterations,  577198 nodes, 55227 on queue, breadth at 36
888000 iterations,  577778 nodes, 55204 on queue, breadth at 36
889000 iterations,  578396 nodes, 55187 on queue, breadth at 36
890000 iterations,  578986 nodes, 55159 on queue, breadth at 36
891000 iterations,  579543 nodes, 55118 on queue, breadth at 36
892000 iterations,  580064 nodes, 55018 on queue, breadth at 36
893000 iterations,  580686 nodes, 55032 on queue, breadth at 36
894000 iterations,  581367 nodes, 55133 on queue, breadth at 36
895000 iterations,  581912 nodes, 55105 on queue, breadth at 36
896000 iterations,  582515 nodes, 55101 on queue, breadth at 36
897000 iterations,  583238 nodes, 55178 on queue, breadth at 36
898000 iterations,  583980 nodes, 55301 on queue, breadth at 36
899000 iterations,  584516 nodes, 55234 on queue, breadth at 36
900000 iterations,  585064 nodes, 55193 on queue, breadth at 36
901000 iterations,  585706 nodes, 55231 on queue, breadth at 36
902000 iterations,  586349 nodes, 55233 on queue, breadth at 36
903000 iterations,  586951 nodes, 55240 on queue, breadth at 36
904000 iterations,  587488 nodes, 55201 on queue, breadth at 36
905000 iterations,  588082 nodes, 55172 on queue, breadth at 36
906000 iterations,  588661 nodes, 55122 on queue, breadth at 37
907000 iterations,  589401 nodes, 55262 on queue, breadth at 37
908000 iterations,  590036 nodes, 55295 on queue, breadth at 37
909000 iterations,  590494 nodes, 55155 on queue, breadth at 37
910000 iterations,  591094 nodes, 55112 on queue, breadth at 37
911000 iterations,  591736 nodes, 55139 on queue, breadth at 37
912000 iterations,  592410 nodes, 55186 on queue, breadth at 37
913000 iterations,  593123 nodes, 55259 on queue, breadth at 37
914000 iterations,  593668 nodes, 55193 on queue, breadth at 37
915000 iterations,  594262 nodes, 55206 on queue, breadth at 37
916000 iterations,  594832 nodes, 55105 on queue, breadth at 37
917000 iterations,  595476 nodes, 55152 on queue, breadth at 37
918000 iterations,  596111 nodes, 55160 on queue, breadth at 37
919000 iterations,  596656 nodes, 55096 on queue, breadth at 37
920000 iterations,  597237 nodes, 55128 on queue, breadth at 37
921000 iterations,  597797 nodes, 55158 on queue, breadth at 37
922000 iterations,  598408 nodes, 55243 on queue, breadth at 37
923000 iterations,  598920 nodes, 55197 on queue, breadth at 37
924000 iterations,  599352 nodes, 55026 on queue, breadth at 37
925000 iterations,  599832 nodes, 54942 on queue, breadth at 37
926000 iterations,  600408 nodes, 54951 on queue, breadth at 37
927000 iterations,  600979 nodes, 54910 on queue, breadth at 37
928000 iterations,  601542 nodes, 54873 on queue, breadth at 37
929000 iterations,  602233 nodes, 54945 on queue, breadth at 37
930000 iterations,  602806 nodes, 54933 on queue, breadth at 37
931000 iterations,  603428 nodes, 54907 on queue, breadth at 37
932000 iterations,  604071 nodes, 54931 on queue, breadth at 37
933000 iterations,  604699 nodes, 54953 on queue, breadth at 37
934000 iterations,  605296 nodes, 54967 on queue, breadth at 37
935000 iterations,  605804 nodes, 54941 on queue, breadth at 37
936000 iterations,  606354 nodes, 54934 on queue, breadth at 37
937000 iterations,  606956 nodes, 54927 on queue, breadth at 37
938000 iterations,  607603 nodes, 54962 on queue, breadth at 37
939000 iterations,  608217 nodes, 54982 on queue, breadth at 37
940000 iterations,  608770 nodes, 54979 on queue, breadth at 37
941000 iterations,  609314 nodes, 54947 on queue, breadth at 37
942000 iterations,  610010 nodes, 54983 on queue, breadth at 37
943000 iterations,  610525 nodes, 54892 on queue, breadth at 37
944000 iterations,  611092 nodes, 54916 on queue, breadth at 37
945000 iterations,  611707 nodes, 54919 on queue, breadth at 37
946000 iterations,  612304 nodes, 54894 on queue, breadth at 37
947000 iterations,  612789 nodes, 54748 on queue, breadth at 37
948000 iterations,  613401 nodes, 54771 on queue, breadth at 37
949000 iterations,  614064 nodes, 54821 on queue, breadth at 37
950000 iterations,  614566 nodes, 54790 on queue, breadth at 37
951000 iterations,  615156 nodes, 54786 on queue, breadth at 37
952000 iterations,  615825 nodes, 54881 on queue, breadth at 37
953000 iterations,  616475 nodes, 54931 on queue, breadth at 37
954000 iterations,  617128 nodes, 54956 on queue, breadth at 37
955000 iterations,  617689 nodes, 54866 on queue, breadth at 37
956000 iterations,  618294 nodes, 54858 on queue, breadth at 37
957000 iterations,  618909 nodes, 54881 on queue, breadth at 37
958000 iterations,  619486 nodes, 54858 on queue, breadth at 37
959000 iterations,  620069 nodes, 54875 on queue, breadth at 37
960000 iterations,  620687 nodes, 54907 on queue, breadth at 37
961000 iterations,  621243 nodes, 54795 on queue, breadth at 38
962000 iterations,  621863 nodes, 54810 on queue, breadth at 38
963000 iterations,  622480 nodes, 54784 on queue, breadth at 38
964000 iterations,  622948 nodes, 54626 on queue, breadth at 38
965000 iterations,  623586 nodes, 54637 on queue, breadth at 38
966000 iterations,  624167 nodes, 54605 on queue, breadth at 38
967000 iterations,  624811 nodes, 54594 on queue, breadth at 38
968000 iterations,  625408 nodes, 54532 on queue, breadth at 38
969000 iterations,  626018 nodes, 54527 on queue, breadth at 38
970000 iterations,  626587 nodes, 54498 on queue, breadth at 38
971000 iterations,  627176 nodes, 54432 on queue, breadth at 38
972000 iterations,  627738 nodes, 54362 on queue, breadth at 38
973000 iterations,  628321 nodes, 54292 on queue, breadth at 38
974000 iterations,  628840 nodes, 54154 on queue, breadth at 38
975000 iterations,  629397 nodes, 54156 on queue, breadth at 38
976000 iterations,  629963 nodes, 54147 on queue, breadth at 38
977000 iterations,  630533 nodes, 54159 on queue, breadth at 38
978000 iterations,  631037 nodes, 54107 on queue, breadth at 38
979000 iterations,  631492 nodes, 53973 on queue, breadth at 38
980000 iterations,  631919 nodes, 53804 on queue, breadth at 38
981000 iterations,  632496 nodes, 53824 on queue, breadth at 38
982000 iterations,  633036 nodes, 53768 on queue, breadth at 38
983000 iterations,  633617 nodes, 53718 on queue, breadth at 38
984000 iterations,  634275 nodes, 53764 on queue, breadth at 38
985000 iterations,  634832 nodes, 53730 on queue, breadth at 38
986000 iterations,  635452 nodes, 53757 on queue, breadth at 38
987000 iterations,  636049 nodes, 53739 on queue, breadth at 38
988000 iterations,  636688 nodes, 53772 on queue, breadth at 38
989000 iterations,  637259 nodes, 53790 on queue, breadth at 38
990000 iterations,  637760 nodes, 53727 on queue, breadth at 38
991000 iterations,  638294 nodes, 53694 on queue, breadth at 38
992000 iterations,  638908 nodes, 53719 on queue, breadth at 38
993000 iterations,  639479 nodes, 53672 on queue, breadth at 38
994000 iterations,  640094 nodes, 53663 on queue, breadth at 38
995000 iterations,  640631 nodes, 53637 on queue, breadth at 38
996000 iterations,  641182 nodes, 53653 on queue, breadth at 38
997000 iterations,  641858 nodes, 53720 on queue, breadth at 38
998000 iterations,  642391 nodes, 53602 on queue, breadth at 38
999000 iterations,  642922 nodes, 53539 on queue, breadth at 38
1000000 iterations,  643556 nodes, 53524 on queue, breadth at 38
1001000 iterations,  644143 nodes, 53490 on queue, breadth at 38
1002000 iterations,  644732 nodes, 53412 on queue, breadth at 38
1003000 iterations,  645331 nodes, 53432 on queue, breadth at 38
1004000 iterations,  645951 nodes, 53493 on queue, breadth at 38
1005000 iterations,  646449 nodes, 53409 on queue, breadth at 38
1006000 iterations,  647099 nodes, 53442 on queue, breadth at 38
1007000 iterations,  647712 nodes, 53472 on queue, breadth at 38
1008000 iterations,  648295 nodes, 53444 on queue, breadth at 38
1009000 iterations,  648947 nodes, 53503 on queue, breadth at 38
1010000 iterations,  649527 nodes, 53444 on queue, breadth at 38
1011000 iterations,  650156 nodes, 53452 on queue, breadth at 38
1012000 iterations,  650702 nodes, 53396 on queue, breadth at 38
1013000 iterations,  651339 nodes, 53394 on queue, breadth at 38
1014000 iterations,  651918 nodes, 53421 on queue, breadth at 38
1015000 iterations,  652531 nodes, 53441 on queue, breadth at 38
1016000 iterations,  653065 nodes, 53368 on queue, breadth at 39
1017000 iterations,  653555 nodes, 53248 on queue, breadth at 39
1018000 iterations,  654244 nodes, 53259 on queue, breadth at 39
1019000 iterations,  654786 nodes, 53152 on queue, breadth at 39
1020000 iterations,  655386 nodes, 53140 on queue, breadth at 39
1021000 iterations,  655965 nodes, 53051 on queue, breadth at 39
1022000 iterations,  656568 nodes, 53007 on queue, breadth at 39
1023000 iterations,  657237 nodes, 53022 on queue, breadth at 39
1024000 iterations,  657761 nodes, 52958 on queue, breadth at 39
1025000 iterations,  658403 nodes, 52970 on queue, breadth at 39
1026000 iterations,  658964 nodes, 52917 on queue, breadth at 39
1027000 iterations,  659534 nodes, 52848 on queue, breadth at 39
1028000 iterations,  660077 nodes, 52756 on queue, breadth at 39
1029000 iterations,  660640 nodes, 52716 on queue, breadth at 39
1030000 iterations,  661238 nodes, 52775 on queue, breadth at 39
1031000 iterations,  661824 nodes, 52782 on queue, breadth at 39
1032000 iterations,  662362 nodes, 52758 on queue, breadth at 39
1033000 iterations,  662827 nodes, 52611 on queue, breadth at 39
1034000 iterations,  663346 nodes, 52553 on queue, breadth at 39
1035000 iterations,  663874 nodes, 52506 on queue, breadth at 39
1036000 iterations,  664432 nodes, 52428 on queue, breadth at 39
1037000 iterations,  665073 nodes, 52405 on queue, breadth at 39
1038000 iterations,  665661 nodes, 52385 on queue, breadth at 39
1039000 iterations,  666230 nodes, 52353 on queue, breadth at 39
1040000 iterations,  666909 nodes, 52390 on queue, breadth at 39
1041000 iterations,  667469 nodes, 52359 on queue, breadth at 39
1042000 iterations,  668046 nodes, 52343 on queue, breadth at 39
1043000 iterations,  668589 nodes, 52325 on queue, breadth at 39
1044000 iterations,  669096 nodes, 52256 on queue, breadth at 39
1045000 iterations,  669625 nodes, 52217 on queue, breadth at 39
1046000 iterations,  670239 nodes, 52250 on queue, breadth at 39
1047000 iterations,  670762 nodes, 52159 on queue, breadth at 39
1048000 iterations,  671339 nodes, 52155 on queue, breadth at 39
1049000 iterations,  671826 nodes, 52097 on queue, breadth at 39
1050000 iterations,  672372 nodes, 52069 on queue, breadth at 39
1051000 iterations,  672918 nodes, 51984 on queue, breadth at 39
1052000 iterations,  673465 nodes, 51960 on queue, breadth at 39
1053000 iterations,  674092 nodes, 51956 on queue, breadth at 39
1054000 iterations,  674781 nodes, 52028 on queue, breadth at 39
1055000 iterations,  675377 nodes, 52037 on queue, breadth at 39
1056000 iterations,  676013 nodes, 52070 on queue, breadth at 39
1057000 iterations,  676564 nodes, 52001 on queue, breadth at 39
1058000 iterations,  677132 nodes, 52000 on queue, breadth at 39
1059000 iterations,  677765 nodes, 52039 on queue, breadth at 39
1060000 iterations,  678346 nodes, 52015 on queue, breadth at 39
1061000 iterations,  678981 nodes, 51984 on queue, breadth at 39
1062000 iterations,  679553 nodes, 51947 on queue, breadth at 39
1063000 iterations,  680158 nodes, 51983 on queue, breadth at 39
1064000 iterations,  680795 nodes, 51975 on queue, breadth at 39
1065000 iterations,  681340 nodes, 51876 on queue, breadth at 39
1066000 iterations,  681936 nodes, 51859 on queue, breadth at 39
1067000 iterations,  682554 nodes, 51913 on queue, breadth at 39
1068000 iterations,  683185 nodes, 51949 on queue, breadth at 39
1069000 iterations,  683663 nodes, 51790 on queue, breadth at 40
1070000 iterations,  684214 nodes, 51650 on queue, breadth at 40
1071000 iterations,  684940 nodes, 51723 on queue, breadth at 40
1072000 iterations,  685447 nodes, 51608 on queue, breadth at 40
1073000 iterations,  686103 nodes, 51645 on queue, breadth at 40
1074000 iterations,  686678 nodes, 51583 on queue, breadth at 40
1075000 iterations,  687306 nodes, 51573 on queue, breadth at 40
1076000 iterations,  688035 nodes, 51630 on queue, breadth at 40
1077000 iterations,  688592 nodes, 51606 on queue, breadth at 40
1078000 iterations,  689224 nodes, 51616 on queue, breadth at 40
1079000 iterations,  689764 nodes, 51518 on queue, breadth at 40
1080000 iterations,  690344 nodes, 51470 on queue, breadth at 40
1081000 iterations,  690886 nodes, 51371 on queue, breadth at 40
1082000 iterations,  691450 nodes, 51354 on queue, breadth at 40
1083000 iterations,  692045 nodes, 51379 on queue, breadth at 40
1084000 iterations,  692596 nodes, 51356 on queue, breadth at 40
1085000 iterations,  693155 nodes, 51365 on queue, breadth at 40
1086000 iterations,  693685 nodes, 51292 on queue, breadth at 40
1087000 iterations,  694273 nodes, 51289 on queue, breadth at 40
1088000 iterations,  694832 nodes, 51266 on queue, breadth at 40
1089000 iterations,  695476 nodes, 51296 on queue, breadth at 40
1090000 iterations,  696064 nodes, 51275 on queue, breadth at 40
1091000 iterations,  696631 nodes, 51223 on queue, breadth at 40
1092000 iterations,  697211 nodes, 51183 on queue, breadth at 40
1093000 iterations,  697795 nodes, 51188 on queue, breadth at 40
1094000 iterations,  698334 nodes, 51144 on queue, breadth at 40
1095000 iterations,  698821 nodes, 51077 on queue, breadth at 40
1096000 iterations,  699317 nodes, 50985 on queue, breadth at 40
1097000 iterations,  699853 nodes, 50974 on queue, breadth at 40
1098000 iterations,  700409 nodes, 50936 on queue, breadth at 40
1099000 iterations,  700920 nodes, 50804 on queue, breadth at 40
1100000 iterations,  701500 nodes, 50760 on queue, breadth at 40
1101000 iterations,  701996 nodes, 50680 on queue, breadth at 40
1102000 iterations,  702542 nodes, 50644 on queue, breadth at 40
1103000 iterations,  703112 nodes, 50568 on queue, breadth at 40
1104000 iterations,  703708 nodes, 50563 on queue, breadth at 40
1105000 iterations,  704356 nodes, 50584 on queue, breadth at 40
1106000 iterations,  704995 nodes, 50629 on queue, breadth at 40
1107000 iterations,  705579 nodes, 50575 on queue, breadth at 40
1108000 iterations,  706186 nodes, 50577 on queue, breadth at 40
1109000 iterations,  706706 nodes, 50514 on queue, breadth at 40
1110000 iterations,  707251 nodes, 50498 on queue, breadth at 40
1111000 iterations,  707874 nodes, 50522 on queue, breadth at 40
1112000 iterations,  708492 nodes, 50517 on queue, breadth at 40
1113000 iterations,  709116 nodes, 50490 on queue, breadth at 40
1114000 iterations,  709678 nodes, 50433 on queue, breadth at 40
1115000 iterations,  710254 nodes, 50368 on queue, breadth at 40
1116000 iterations,  710883 nodes, 50341 on queue, breadth at 40
1117000 iterations,  711485 nodes, 50315 on queue, breadth at 40
1118000 iterations,  712106 nodes, 50331 on queue, breadth at 40
1119000 iterations,  712695 nodes, 50349 on queue, breadth at 40
1120000 iterations,  713321 nodes, 50407 on queue, breadth at 40
1121000 iterations,  713922 nodes, 50355 on queue, breadth at 41
1122000 iterations,  714490 nodes, 50280 on queue, breadth at 41
1123000 iterations,  715117 nodes, 50258 on queue, breadth at 41
1124000 iterations,  715640 nodes, 50164 on queue, breadth at 41
1125000 iterations,  716400 nodes, 50296 on queue, breadth at 41
1126000 iterations,  716919 nodes, 50140 on queue, breadth at 41
1127000 iterations,  717501 nodes, 50063 on queue, breadth at 41
1128000 iterations,  718170 nodes, 50125 on queue, breadth at 41
1129000 iterations,  718786 nodes, 50166 on queue, breadth at 41
1130000 iterations,  719370 nodes, 50058 on queue, breadth at 41
1131000 iterations,  719947 nodes, 50003 on queue, breadth at 41
1132000 iterations,  720539 nodes, 49915 on queue, breadth at 41
1133000 iterations,  721166 nodes, 49957 on queue, breadth at 41
1134000 iterations,  721788 nodes, 50025 on queue, breadth at 41
1135000 iterations,  722317 nodes, 49983 on queue, breadth at 41
1136000 iterations,  722856 nodes, 50001 on queue, breadth at 41
1137000 iterations,  723417 nodes, 49974 on queue, breadth at 41
1138000 iterations,  724066 nodes, 50074 on queue, breadth at 41
1139000 iterations,  724640 nodes, 50074 on queue, breadth at 41
1140000 iterations,  725251 nodes, 50035 on queue, breadth at 41
1141000 iterations,  725813 nodes, 49991 on queue, breadth at 41
1142000 iterations,  726442 nodes, 50009 on queue, breadth at 41
1143000 iterations,  727037 nodes, 49979 on queue, breadth at 41
1144000 iterations,  727570 nodes, 49858 on queue, breadth at 41
1145000 iterations,  728098 nodes, 49787 on queue, breadth at 41
1146000 iterations,  728596 nodes, 49664 on queue, breadth at 41
1147000 iterations,  729123 nodes, 49637 on queue, breadth at 41
1148000 iterations,  729617 nodes, 49561 on queue, breadth at 41
1149000 iterations,  730219 nodes, 49551 on queue, breadth at 41
1150000 iterations,  730769 nodes, 49438 on queue, breadth at 41
1151000 iterations,  731311 nodes, 49375 on queue, breadth at 41
1152000 iterations,  731833 nodes, 49322 on queue, breadth at 41
1153000 iterations,  732363 nodes, 49287 on queue, breadth at 41
1154000 iterations,  732971 nodes, 49262 on queue, breadth at 41
1155000 iterations,  733580 nodes, 49308 on queue, breadth at 41
1156000 iterations,  734201 nodes, 49317 on queue, breadth at 41
1157000 iterations,  734868 nodes, 49389 on queue, breadth at 41
1158000 iterations,  735475 nodes, 49351 on queue, breadth at 41
1159000 iterations,  736025 nodes, 49289 on queue, breadth at 41
1160000 iterations,  736553 nodes, 49229 on queue, breadth at 41
1161000 iterations,  737088 nodes, 49195 on queue, breadth at 41
1162000 iterations,  737793 nodes, 49252 on queue, breadth at 41
1163000 iterations,  738387 nodes, 49222 on queue, breadth at 41
1164000 iterations,  738994 nodes, 49165 on queue, breadth at 41
1165000 iterations,  739628 nodes, 49182 on queue, breadth at 41
1166000 iterations,  740266 nodes, 49212 on queue, breadth at 41
1167000 iterations,  740896 nodes, 49227 on queue, breadth at 41
1168000 iterations,  741524 nodes, 49230 on queue, breadth at 41
1169000 iterations,  742140 nodes, 49241 on queue, breadth at 41
1170000 iterations,  742717 nodes, 49236 on queue, breadth at 41
1171000 iterations,  743372 nodes, 49275 on queue, breadth at 42
1172000 iterations,  744043 nodes, 49297 on queue, breadth at 42
1173000 iterations,  744598 nodes, 49251 on queue, breadth at 42
1174000 iterations,  745161 nodes, 49198 on queue, breadth at 42
1175000 iterations,  745853 nodes, 49236 on queue, breadth at 42
1176000 iterations,  746528 nodes, 49255 on queue, breadth at 42
1177000 iterations,  747157 nodes, 49224 on queue, breadth at 42
1178000 iterations,  747820 nodes, 49241 on queue, breadth at 42
1179000 iterations,  748373 nodes, 49224 on queue, breadth at 42
1180000 iterations,  748980 nodes, 49228 on queue, breadth at 42
1181000 iterations,  749544 nodes, 49177 on queue, breadth at 42
1182000 iterations,  750158 nodes, 49152 on queue, breadth at 42
1183000 iterations,  750786 nodes, 49187 on queue, breadth at 42
1184000 iterations,  751377 nodes, 49200 on queue, breadth at 42
1185000 iterations,  751937 nodes, 49156 on queue, breadth at 42
1186000 iterations,  752500 nodes, 49130 on queue, breadth at 42
1187000 iterations,  753058 nodes, 49120 on queue, breadth at 42
1188000 iterations,  753644 nodes, 49135 on queue, breadth at 42
1189000 iterations,  754229 nodes, 49157 on queue, breadth at 42
1190000 iterations,  754857 nodes, 49185 on queue, breadth at 42
1191000 iterations,  755480 nodes, 49194 on queue, breadth at 42
1192000 iterations,  756137 nodes, 49245 on queue, breadth at 42
1193000 iterations,  756762 nodes, 49216 on queue, breadth at 42
1194000 iterations,  757308 nodes, 49134 on queue, breadth at 42
1195000 iterations,  757876 nodes, 49101 on queue, breadth at 42
1196000 iterations,  758389 nodes, 49050 on queue, breadth at 42
1197000 iterations,  758902 nodes, 48989 on queue, breadth at 42
1198000 iterations,  759464 nodes, 48966 on queue, breadth at 42
1199000 iterations,  760105 nodes, 48984 on queue, breadth at 42
1200000 iterations,  760651 nodes, 48899 on queue, breadth at 42
1201000 iterations,  761196 nodes, 48867 on queue, breadth at 42
1202000 iterations,  761754 nodes, 48817 on queue, breadth at 42
1203000 iterations,  762365 nodes, 48769 on queue, breadth at 42
1204000 iterations,  762896 nodes, 48698 on queue, breadth at 42
1205000 iterations,  763541 nodes, 48718 on queue, breadth at 42
1206000 iterations,  764128 nodes, 48679 on queue, breadth at 42
1207000 iterations,  764763 nodes, 48691 on queue, breadth at 42
1208000 iterations,  765361 nodes, 48692 on queue, breadth at 42
1209000 iterations,  765895 nodes, 48612 on queue, breadth at 42
1210000 iterations,  766474 nodes, 48642 on queue, breadth at 42
1211000 iterations,  767108 nodes, 48676 on queue, breadth at 42
1212000 iterations,  767683 nodes, 48612 on queue, breadth at 42
1213000 iterations,  768304 nodes, 48583 on queue, breadth at 42
1214000 iterations,  768963 nodes, 48614 on queue, breadth at 42
1215000 iterations,  769531 nodes, 48572 on queue, breadth at 42
1216000 iterations,  770153 nodes, 48578 on queue, breadth at 42
1217000 iterations,  770789 nodes, 48568 on queue, breadth at 42
1218000 iterations,  771388 nodes, 48570 on queue, breadth at 42
1219000 iterations,  771951 nodes, 48549 on queue, breadth at 42
1220000 iterations,  772668 nodes, 48647 on queue, breadth at 42
1221000 iterations,  773370 nodes, 48680 on queue, breadth at 43
1222000 iterations,  773939 nodes, 48608 on queue, breadth at 43
1223000 iterations,  774525 nodes, 48563 on queue, breadth at 43
1224000 iterations,  775178 nodes, 48584 on queue, breadth at 43
1225000 iterations,  775899 nodes, 48626 on queue, breadth at 43
1226000 iterations,  776599 nodes, 48661 on queue, breadth at 43
1227000 iterations,  777273 nodes, 48678 on queue, breadth at 43
1228000 iterations,  777777 nodes, 48571 on queue, breadth at 43
1229000 iterations,  778347 nodes, 48536 on queue, breadth at 43
1230000 iterations,  778964 nodes, 48520 on queue, breadth at 43
1231000 iterations,  779555 nodes, 48472 on queue, breadth at 43
1232000 iterations,  780187 nodes, 48507 on queue, breadth at 43
1233000 iterations,  780775 nodes, 48519 on queue, breadth at 43
1234000 iterations,  781345 nodes, 48513 on queue, breadth at 43
1235000 iterations,  781943 nodes, 48560 on queue, breadth at 43
1236000 iterations,  782532 nodes, 48575 on queue, breadth at 43
1237000 iterations,  783125 nodes, 48611 on queue, breadth at 43
1238000 iterations,  783677 nodes, 48575 on queue, breadth at 43
1239000 iterations,  784304 nodes, 48577 on queue, breadth at 43
1240000 iterations,  784942 nodes, 48578 on queue, breadth at 43
1241000 iterations,  785587 nodes, 48616 on queue, breadth at 43
1242000 iterations,  786188 nodes, 48609 on queue, breadth at 43
1243000 iterations,  786764 nodes, 48562 on queue, breadth at 43
1244000 iterations,  787405 nodes, 48586 on queue, breadth at 43
1245000 iterations,  787906 nodes, 48505 on queue, breadth at 43
1246000 iterations,  788429 nodes, 48493 on queue, breadth at 43
1247000 iterations,  789019 nodes, 48497 on queue, breadth at 43
1248000 iterations,  789634 nodes, 48496 on queue, breadth at 43
1249000 iterations,  790138 nodes, 48390 on queue, breadth at 43
1250000 iterations,  790664 nodes, 48370 on queue, breadth at 43
1251000 iterations,  791244 nodes, 48365 on queue, breadth at 43
1252000 iterations,  791825 nodes, 48298 on queue, breadth at 43
1253000 iterations,  792432 nodes, 48329 on queue, breadth at 43
1254000 iterations,  793083 nodes, 48338 on queue, breadth at 43
1255000 iterations,  793717 nodes, 48371 on queue, breadth at 43
1256000 iterations,  794312 nodes, 48370 on queue, breadth at 43
1257000 iterations,  794847 nodes, 48309 on queue, breadth at 43
1258000 iterations,  795409 nodes, 48283 on queue, breadth at 43
1259000 iterations,  795986 nodes, 48270 on queue, breadth at 43
1260000 iterations,  796587 nodes, 48265 on queue, breadth at 43
1261000 iterations,  797197 nodes, 48235 on queue, breadth at 43
1262000 iterations,  797720 nodes, 48108 on queue, breadth at 43
1263000 iterations,  798355 nodes, 48137 on queue, breadth at 43
1264000 iterations,  798924 nodes, 48105 on queue, breadth at 43
1265000 iterations,  799531 nodes, 48094 on queue, breadth at 43
1266000 iterations,  800137 nodes, 48068 on queue, breadth at 43
1267000 iterations,  800717 nodes, 48064 on queue, breadth at 43
1268000 iterations,  801343 nodes, 48083 on queue, breadth at 43
1269000 iterations,  801951 nodes, 48093 on queue, breadth at 44
1270000 iterations,  802714 nodes, 48237 on queue, breadth at 44
1271000 iterations,  803207 nodes, 48081 on queue, breadth at 44
1272000 iterations,  803762 nodes, 48007 on queue, breadth at 44
1273000 iterations,  804423 nodes, 48034 on queue, breadth at 44
1274000 iterations,  805168 nodes, 48147 on queue, breadth at 44
1275000 iterations,  805805 nodes, 48135 on queue, breadth at 44
1276000 iterations,  806431 nodes, 48131 on queue, breadth at 44
1277000 iterations,  806999 nodes, 48092 on queue, breadth at 44
1278000 iterations,  807592 nodes, 48039 on queue, breadth at 44
1279000 iterations,  808220 nodes, 48039 on queue, breadth at 44
1280000 iterations,  808853 nodes, 48043 on queue, breadth at 44
1281000 iterations,  809425 nodes, 48029 on queue, breadth at 44
1282000 iterations,  809990 nodes, 48044 on queue, breadth at 44
1283000 iterations,  810576 nodes, 48052 on queue, breadth at 44
1284000 iterations,  811164 nodes, 48070 on queue, breadth at 44
1285000 iterations,  811713 nodes, 48054 on queue, breadth at 44
1286000 iterations,  812308 nodes, 48064 on queue, breadth at 44
1287000 iterations,  812936 nodes, 48100 on queue, breadth at 44
1288000 iterations,  813572 nodes, 48104 on queue, breadth at 44
1289000 iterations,  814241 nodes, 48156 on queue, breadth at 44
1290000 iterations,  814872 nodes, 48173 on queue, breadth at 44
1291000 iterations,  815509 nodes, 48171 on queue, breadth at 44
1292000 iterations,  816093 nodes, 48128 on queue, breadth at 44
1293000 iterations,  816713 nodes, 48149 on queue, breadth at 44
1294000 iterations,  817233 nodes, 48118 on queue, breadth at 44
1295000 iterations,  817779 nodes, 48088 on queue, breadth at 44
1296000 iterations,  818403 nodes, 48108 on queue, breadth at 44
1297000 iterations,  819022 nodes, 48088 on queue, breadth at 44
1298000 iterations,  819512 nodes, 48001 on queue, breadth at 44
1299000 iterations,  820003 nodes, 47927 on queue, breadth at 44
1300000 iterations,  820624 nodes, 47926 on queue, breadth at 44
1301000 iterations,  821178 nodes, 47907 on queue, breadth at 44
1302000 iterations,  821916 nodes, 47994 on queue, breadth at 44
1303000 iterations,  822537 nodes, 48029 on queue, breadth at 44
1304000 iterations,  823142 nodes, 48045 on queue, breadth at 44
1305000 iterations,  823698 nodes, 47993 on queue, breadth at 44
1306000 iterations,  824236 nodes, 47951 on queue, breadth at 44
1307000 iterations,  824771 nodes, 47944 on queue, breadth at 44
1308000 iterations,  825369 nodes, 47923 on queue, breadth at 44
1309000 iterations,  825922 nodes, 47879 on queue, breadth at 44
1310000 iterations,  826474 nodes, 47780 on queue, breadth at 44
1311000 iterations,  827079 nodes, 47744 on queue, breadth at 44
1312000 iterations,  827644 nodes, 47656 on queue, breadth at 44
1313000 iterations,  828259 nodes, 47656 on queue, breadth at 44
1314000 iterations,  828862 nodes, 47645 on queue, breadth at 44
1315000 iterations,  829408 nodes, 47577 on queue, breadth at 44
1316000 iterations,  830026 nodes, 47612 on queue, breadth at 44
1317000 iterations,  830655 nodes, 47629 on queue, breadth at 45
1318000 iterations,  831352 nodes, 47692 on queue, breadth at 45
1319000 iterations,  831892 nodes, 47608 on queue, breadth at 45
1320000 iterations,  832487 nodes, 47568 on queue, breadth at 45
1321000 iterations,  833120 nodes, 47590 on queue, breadth at 45
1322000 iterations,  833847 nodes, 47682 on queue, breadth at 45
1323000 iterations,  834453 nodes, 47639 on queue, breadth at 45
1324000 iterations,  835101 nodes, 47662 on queue, breadth at 45
1325000 iterations,  835654 nodes, 47644 on queue, breadth at 45
1326000 iterations,  836331 nodes, 47687 on queue, breadth at 45
1327000 iterations,  836928 nodes, 47670 on queue, breadth at 45
1328000 iterations,  837571 nodes, 47653 on queue, breadth at 45
1329000 iterations,  838174 nodes, 47681 on queue, breadth at 45
1330000 iterations,  838731 nodes, 47664 on queue, breadth at 45
1331000 iterations,  839331 nodes, 47684 on queue, breadth at 45
1332000 iterations,  839868 nodes, 47676 on queue, breadth at 45
1333000 iterations,  840391 nodes, 47633 on queue, breadth at 45
1334000 iterations,  840985 nodes, 47658 on queue, breadth at 45
1335000 iterations,  841624 nodes, 47707 on queue, breadth at 45
1336000 iterations,  842243 nodes, 47677 on queue, breadth at 45
1337000 iterations,  842834 nodes, 47642 on queue, breadth at 45
1338000 iterations,  843442 nodes, 47644 on queue, breadth at 45
1339000 iterations,  844035 nodes, 47628 on queue, breadth at 45
1340000 iterations,  844667 nodes, 47642 on queue, breadth at 45
1341000 iterations,  845335 nodes, 47713 on queue, breadth at 45
1342000 iterations,  845853 nodes, 47666 on queue, breadth at 45
1343000 iterations,  846411 nodes, 47669 on queue, breadth at 45
1344000 iterations,  847040 nodes, 47717 on queue, breadth at 45
1345000 iterations,  847700 nodes, 47738 on queue, breadth at 45
1346000 iterations,  848213 nodes, 47631 on queue, breadth at 45
1347000 iterations,  848758 nodes, 47581 on queue, breadth at 45
1348000 iterations,  849389 nodes, 47555 on queue, breadth at 45
1349000 iterations,  849969 nodes, 47555 on queue, breadth at 45
1350000 iterations,  850645 nodes, 47603 on queue, breadth at 45
1351000 iterations,  851232 nodes, 47581 on queue, breadth at 45
1352000 iterations,  851825 nodes, 47529 on queue, breadth at 45
1353000 iterations,  852414 nodes, 47511 on queue, breadth at 45
1354000 iterations,  852950 nodes, 47456 on queue, breadth at 45
1355000 iterations,  853479 nodes, 47409 on queue, breadth at 45
1356000 iterations,  854064 nodes, 47373 on queue, breadth at 45
1357000 iterations,  854623 nodes, 47320 on queue, breadth at 45
1358000 iterations,  855136 nodes, 47154 on queue, breadth at 45
1359000 iterations,  855794 nodes, 47203 on queue, breadth at 45
1360000 iterations,  856437 nodes, 47195 on queue, breadth at 45
1361000 iterations,  857042 nodes, 47165 on queue, breadth at 45
1362000 iterations,  857613 nodes, 47110 on queue, breadth at 45
1363000 iterations,  858210 nodes, 47105 on queue, breadth at 45
1364000 iterations,  858785 nodes, 47083 on queue, breadth at 45
1365000 iterations,  859471 nodes, 47150 on queue, breadth at 46
1366000 iterations,  860097 nodes, 47098 on queue, breadth at 46
1367000 iterations,  860756 nodes, 47072 on queue, breadth at 46
1368000 iterations,  861328 nodes, 47033 on queue, breadth at 46
1369000 iterations,  862010 nodes, 47083 on queue, breadth at 46
1370000 iterations,  862705 nodes, 47132 on queue, breadth at 46
1371000 iterations,  863370 nodes, 47139 on queue, breadth at 46
1372000 iterations,  863953 nodes, 47105 on queue, breadth at 46
1373000 iterations,  864502 nodes, 47087 on queue, breadth at 46
1374000 iterations,  865170 nodes, 47110 on queue, breadth at 46
1375000 iterations,  865762 nodes, 47069 on queue, breadth at 46
1376000 iterations,  866464 nodes, 47162 on queue, breadth at 46
1377000 iterations,  867039 nodes, 47148 on queue, breadth at 46
1378000 iterations,  867626 nodes, 47161 on queue, breadth at 46
1379000 iterations,  868183 nodes, 47139 on queue, breadth at 46
1380000 iterations,  868704 nodes, 47097 on queue, breadth at 46
1381000 iterations,  869250 nodes, 47041 on queue, breadth at 46
1382000 iterations,  869791 nodes, 47037 on queue, breadth at 46
1383000 iterations,  870403 nodes, 47036 on queue, breadth at 46
1384000 iterations,  871010 nodes, 47010 on queue, breadth at 46
1385000 iterations,  871672 nodes, 47054 on queue, breadth at 46
1386000 iterations,  872268 nodes, 47059 on queue, breadth at 46
1387000 iterations,  872850 nodes, 47038 on queue, breadth at 46
1388000 iterations,  873450 nodes, 47053 on queue, breadth at 46
1389000 iterations,  874050 nodes, 47061 on queue, breadth at 46
1390000 iterations,  874620 nodes, 47091 on queue, breadth at 46
1391000 iterations,  875193 nodes, 47094 on queue, breadth at 46
1392000 iterations,  875809 nodes, 47102 on queue, breadth at 46
1393000 iterations,  876415 nodes, 47082 on queue, breadth at 46
1394000 iterations,  876957 nodes, 47068 on queue, breadth at 46
1395000 iterations,  877588 nodes, 47082 on queue, breadth at 46
1396000 iterations,  878199 nodes, 47091 on queue, breadth at 46
1397000 iterations,  878839 nodes, 47123 on queue, breadth at 46
1398000 iterations,  879458 nodes, 47136 on queue, breadth at 46
1399000 iterations,  880062 nodes, 47100 on queue, breadth at 46
1400000 iterations,  880617 nodes, 47057 on queue, breadth at 46
1401000 iterations,  881190 nodes, 47020 on queue, breadth at 46
1402000 iterations,  881739 nodes, 46994 on queue, breadth at 46
1403000 iterations,  882322 nodes, 46972 on queue, breadth at 46
1404000 iterations,  882885 nodes, 46918 on queue, breadth at 46
1405000 iterations,  883485 nodes, 46842 on queue, breadth at 46
1406000 iterations,  884120 nodes, 46877 on queue, breadth at 46
1407000 iterations,  884739 nodes, 46874 on queue, breadth at 46
1408000 iterations,  885347 nodes, 46879 on queue, breadth at 46
1409000 iterations,  885962 nodes, 46859 on queue, breadth at 46
1410000 iterations,  886567 nodes, 46902 on queue, breadth at 46
1411000 iterations,  887122 nodes, 46856 on queue, breadth at 46
1412000 iterations,  887729 nodes, 46871 on queue, breadth at 47
1413000 iterations,  888403 nodes, 46934 on queue, breadth at 47
1414000 iterations,  889139 nodes, 47031 on queue, breadth at 47
1415000 iterations,  889779 nodes, 47057 on queue, breadth at 47
1416000 iterations,  890456 nodes, 47090 on queue, breadth at 47
1417000 iterations,  891096 nodes, 47119 on queue, breadth at 47
1418000 iterations,  891756 nodes, 47117 on queue, breadth at 47
1419000 iterations,  892385 nodes, 47122 on queue, breadth at 47
1420000 iterations,  892937 nodes, 47094 on queue, breadth at 47
1421000 iterations,  893542 nodes, 47076 on queue, breadth at 47
1422000 iterations,  894192 nodes, 47063 on queue, breadth at 47
1423000 iterations,  894870 nodes, 47117 on queue, breadth at 47
1424000 iterations,  895452 nodes, 47087 on queue, breadth at 47
1425000 iterations,  896033 nodes, 47087 on queue, breadth at 47
1426000 iterations,  896623 nodes, 47111 on queue, breadth at 47
1427000 iterations,  897140 nodes, 47084 on queue, breadth at 47
1428000 iterations,  897730 nodes, 47077 on queue, breadth at 47
1429000 iterations,  898236 nodes, 47018 on queue, breadth at 47
1430000 iterations,  898834 nodes, 47003 on queue, breadth at 47
1431000 iterations,  899446 nodes, 46987 on queue, breadth at 47
1432000 iterations,  900117 nodes, 47021 on queue, breadth at 47
1433000 iterations,  900691 nodes, 47029 on queue, breadth at 47
1434000 iterations,  901260 nodes, 46956 on queue, breadth at 47
1435000 iterations,  901840 nodes, 46926 on queue, breadth at 47
1436000 iterations,  902393 nodes, 46901 on queue, breadth at 47
1437000 iterations,  902939 nodes, 46897 on queue, breadth at 47
1438000 iterations,  903514 nodes, 46892 on queue, breadth at 47
1439000 iterations,  904112 nodes, 46868 on queue, breadth at 47
1440000 iterations,  904726 nodes, 46839 on queue, breadth at 47
1441000 iterations,  905271 nodes, 46832 on queue, breadth at 47
1442000 iterations,  905896 nodes, 46821 on queue, breadth at 47
1443000 iterations,  906492 nodes, 46816 on queue, breadth at 47
1444000 iterations,  907061 nodes, 46784 on queue, breadth at 47
1445000 iterations,  907665 nodes, 46763 on queue, breadth at 47
1446000 iterations,  908254 nodes, 46736 on queue, breadth at 47
1447000 iterations,  908799 nodes, 46683 on queue, breadth at 47
1448000 iterations,  909389 nodes, 46671 on queue, breadth at 47
1449000 iterations,  909951 nodes, 46662 on queue, breadth at 47
1450000 iterations,  910531 nodes, 46661 on queue, breadth at 47
1451000 iterations,  911104 nodes, 46624 on queue, breadth at 47
1452000 iterations,  911739 nodes, 46611 on queue, breadth at 47
1453000 iterations,  912339 nodes, 46591 on queue, breadth at 47
1454000 iterations,  912910 nodes, 46518 on queue, breadth at 47
1455000 iterations,  913491 nodes, 46493 on queue, breadth at 47
1456000 iterations,  914059 nodes, 46469 on queue, breadth at 47
1457000 iterations,  914641 nodes, 46492 on queue, breadth at 47
1458000 iterations,  915237 nodes, 46454 on queue, breadth at 47
1459000 iterations,  915837 nodes, 46417 on queue, breadth at 48
1460000 iterations,  916501 nodes, 46453 on queue, breadth at 48
1461000 iterations,  917198 nodes, 46521 on queue, breadth at 48
1462000 iterations,  917862 nodes, 46612 on queue, breadth at 48
1463000 iterations,  918518 nodes, 46644 on queue, breadth at 48
1464000 iterations,  919107 nodes, 46617 on queue, breadth at 48
1465000 iterations,  919727 nodes, 46596 on queue, breadth at 48
1466000 iterations,  920412 nodes, 46635 on queue, breadth at 48
1467000 iterations,  920996 nodes, 46600 on queue, breadth at 48
1468000 iterations,  921618 nodes, 46570 on queue, breadth at 48
1469000 iterations,  922261 nodes, 46587 on queue, breadth at 48
1470000 iterations,  922938 nodes, 46631 on queue, breadth at 48
1471000 iterations,  923536 nodes, 46644 on queue, breadth at 48
1472000 iterations,  924122 nodes, 46638 on queue, breadth at 48
1473000 iterations,  924697 nodes, 46632 on queue, breadth at 48
1474000 iterations,  925240 nodes, 46565 on queue, breadth at 48
1475000 iterations,  925775 nodes, 46553 on queue, breadth at 48
1476000 iterations,  926275 nodes, 46500 on queue, breadth at 48
1477000 iterations,  926854 nodes, 46486 on queue, breadth at 48
1478000 iterations,  927467 nodes, 46460 on queue, breadth at 48
1479000 iterations,  928131 nodes, 46454 on queue, breadth at 48
1480000 iterations,  928687 nodes, 46406 on queue, breadth at 48
1481000 iterations,  929312 nodes, 46386 on queue, breadth at 48
1482000 iterations,  929918 nodes, 46395 on queue, breadth at 48
1483000 iterations,  930456 nodes, 46310 on queue, breadth at 48
1484000 iterations,  930995 nodes, 46330 on queue, breadth at 48
1485000 iterations,  931563 nodes, 46307 on queue, breadth at 48
1486000 iterations,  932174 nodes, 46320 on queue, breadth at 48
1487000 iterations,  932813 nodes, 46350 on queue, breadth at 48
1488000 iterations,  933379 nodes, 46356 on queue, breadth at 48
1489000 iterations,  934015 nodes, 46414 on queue, breadth at 48
1490000 iterations,  934624 nodes, 46416 on queue, breadth at 48
1491000 iterations,  935227 nodes, 46396 on queue, breadth at 48
1492000 iterations,  935830 nodes, 46394 on queue, breadth at 48
1493000 iterations,  936424 nodes, 46396 on queue, breadth at 48
1494000 iterations,  937028 nodes, 46386 on queue, breadth at 48
1495000 iterations,  937588 nodes, 46358 on queue, breadth at 48
1496000 iterations,  938105 nodes, 46310 on queue, breadth at 48
1497000 iterations,  938708 nodes, 46299 on queue, breadth at 48
1498000 iterations,  939350 nodes, 46313 on queue, breadth at 48
1499000 iterations,  939923 nodes, 46256 on queue, breadth at 48
1500000 iterations,  940520 nodes, 46241 on queue, breadth at 48
1501000 iterations,  941149 nodes, 46236 on queue, breadth at 48
1502000 iterations,  941718 nodes, 46135 on queue, breadth at 48
1503000 iterations,  942253 nodes, 46069 on queue, breadth at 48
1504000 iterations,  942840 nodes, 46072 on queue, breadth at 48
1505000 iterations,  943473 nodes, 46079 on queue, breadth at 48
1506000 iterations,  944154 nodes, 46147 on queue, breadth at 49
1507000 iterations,  944792 nodes, 46168 on queue, breadth at 49
1508000 iterations,  945461 nodes, 46226 on queue, breadth at 49
1509000 iterations,  946059 nodes, 46230 on queue, breadth at 49
1510000 iterations,  946670 nodes, 46205 on queue, breadth at 49
1511000 iterations,  947191 nodes, 46076 on queue, breadth at 49
1512000 iterations,  947810 nodes, 46031 on queue, breadth at 49
1513000 iterations,  948441 nodes, 46039 on queue, breadth at 49
1514000 iterations,  949119 nodes, 46144 on queue, breadth at 49
1515000 iterations,  949756 nodes, 46176 on queue, breadth at 49
1516000 iterations,  950404 nodes, 46181 on queue, breadth at 49
1517000 iterations,  950986 nodes, 46173 on queue, breadth at 49
1518000 iterations,  951579 nodes, 46199 on queue, breadth at 49
1519000 iterations,  952128 nodes, 46170 on queue, breadth at 49
1520000 iterations,  952713 nodes, 46187 on queue, breadth at 49
1521000 iterations,  953291 nodes, 46192 on queue, breadth at 49
1522000 iterations,  953846 nodes, 46138 on queue, breadth at 49
1523000 iterations,  954361 nodes, 46053 on queue, breadth at 49
1524000 iterations,  955033 nodes, 46073 on queue, breadth at 49
1525000 iterations,  955651 nodes, 46052 on queue, breadth at 49
1526000 iterations,  956213 nodes, 46008 on queue, breadth at 49
1527000 iterations,  956821 nodes, 45995 on queue, breadth at 49
1528000 iterations,  957440 nodes, 45995 on queue, breadth at 49
1529000 iterations,  958047 nodes, 45980 on queue, breadth at 49
1530000 iterations,  958603 nodes, 45980 on queue, breadth at 49
1531000 iterations,  959170 nodes, 45963 on queue, breadth at 49
1532000 iterations,  959784 nodes, 45988 on queue, breadth at 49
1533000 iterations,  960462 nodes, 46030 on queue, breadth at 49
1534000 iterations,  961020 nodes, 46026 on queue, breadth at 49
1535000 iterations,  961564 nodes, 45997 on queue, breadth at 49
1536000 iterations,  962239 nodes, 46047 on queue, breadth at 49
1537000 iterations,  962788 nodes, 46021 on queue, breadth at 49
1538000 iterations,  963455 nodes, 46072 on queue, breadth at 49
1539000 iterations,  964044 nodes, 46052 on queue, breadth at 49
1540000 iterations,  964660 nodes, 46028 on queue, breadth at 49
1541000 iterations,  965234 nodes, 46020 on queue, breadth at 49
1542000 iterations,  965749 nodes, 45945 on queue, breadth at 49
1543000 iterations,  966332 nodes, 45916 on queue, breadth at 49
1544000 iterations,  966953 nodes, 45954 on queue, breadth at 49
1545000 iterations,  967548 nodes, 45898 on queue, breadth at 49
1546000 iterations,  968131 nodes, 45858 on queue, breadth at 49
1547000 iterations,  968807 nodes, 45909 on queue, breadth at 49
1548000 iterations,  969396 nodes, 45899 on queue, breadth at 49
1549000 iterations,  969969 nodes, 45857 on queue, breadth at 49
1550000 iterations,  970586 nodes, 45878 on queue, breadth at 49
1551000 iterations,  971183 nodes, 45864 on queue, breadth at 49
1552000 iterations,  971780 nodes, 45869 on queue, breadth at 50
1553000 iterations,  972380 nodes, 45841 on queue, breadth at 50
1554000 iterations,  973021 nodes, 45836 on queue, breadth at 50
1555000 iterations,  973612 nodes, 45793 on queue, breadth at 50
1556000 iterations,  974226 nodes, 45740 on queue, breadth at 50
1557000 iterations,  974750 nodes, 45610 on queue, breadth at 50
1558000 iterations,  975370 nodes, 45566 on queue, breadth at 50
1559000 iterations,  975986 nodes, 45591 on queue, breadth at 50
1560000 iterations,  976607 nodes, 45666 on queue, breadth at 50
1561000 iterations,  977226 nodes, 45675 on queue, breadth at 50
1562000 iterations,  977894 nodes, 45708 on queue, breadth at 50
1563000 iterations,  978431 nodes, 45607 on queue, breadth at 50
1564000 iterations,  979009 nodes, 45574 on queue, breadth at 50
1565000 iterations,  979551 nodes, 45540 on queue, breadth at 50
1566000 iterations,  980134 nodes, 45518 on queue, breadth at 50
1567000 iterations,  980668 nodes, 45531 on queue, breadth at 50
1568000 iterations,  981273 nodes, 45514 on queue, breadth at 50
1569000 iterations,  981847 nodes, 45495 on queue, breadth at 50
1570000 iterations,  982543 nodes, 45554 on queue, breadth at 50
1571000 iterations,  983172 nodes, 45559 on queue, breadth at 50
1572000 iterations,  983755 nodes, 45523 on queue, breadth at 50
1573000 iterations,  984351 nodes, 45512 on queue, breadth at 50
1574000 iterations,  985013 nodes, 45553 on queue, breadth at 50
1575000 iterations,  985628 nodes, 45569 on queue, breadth at 50
1576000 iterations,  986200 nodes, 45597 on queue, breadth at 50
1577000 iterations,  986759 nodes, 45604 on queue, breadth at 50
1578000 iterations,  987366 nodes, 45588 on queue, breadth at 50
1579000 iterations,  988039 nodes, 45620 on queue, breadth at 50
1580000 iterations,  988608 nodes, 45600 on queue, breadth at 50
1581000 iterations,  989143 nodes, 45541 on queue, breadth at 50
1582000 iterations,  989781 nodes, 45553 on queue, breadth at 50
1583000 iterations,  990301 nodes, 45496 on queue, breadth at 50
1584000 iterations,  990947 nodes, 45523 on queue, breadth at 50
1585000 iterations,  991551 nodes, 45496 on queue, breadth at 50
1586000 iterations,  992202 nodes, 45523 on queue, breadth at 50
1587000 iterations,  992788 nodes, 45518 on queue, breadth at 50
1588000 iterations,  993371 nodes, 45543 on queue, breadth at 50
1589000 iterations,  993981 nodes, 45555 on queue, breadth at 50
1590000 iterations,  994587 nodes, 45517 on queue, breadth at 50
1591000 iterations,  995168 nodes, 45454 on queue, breadth at 50
1592000 iterations,  995785 nodes, 45444 on queue, breadth at 50
1593000 iterations,  996451 nodes, 45475 on queue, breadth at 50
1594000 iterations,  996987 nodes, 45397 on queue, breadth at 50
1595000 iterations,  997569 nodes, 45385 on queue, breadth at 50
1596000 iterations,  998138 nodes, 45423 on queue, breadth at 50
max breadth reached at 998654 nodes
1596946 iterations,  998654 nodes, 45349 on queue, breadth at 51 8250658 k-mers

















In [280]:



    


G_lump_full = extract_local_compact_dbg('A' * 25, sp.graph, verbose=False)



    


















    






0 iterations,  0 nodes, 0 on queue, breadth at 0
50000 iterations,  42809 nodes, 12641 on queue, breadth at 15
100000 iterations,  82114 nodes, 22488 on queue, breadth at 18
150000 iterations,  118366 nodes, 29867 on queue, breadth at 20
200000 iterations,  152386 nodes, 35278 on queue, breadth at 22
250000 iterations,  184722 nodes, 38541 on queue, breadth at 23
300000 iterations,  216946 nodes, 41738 on queue, breadth at 24
350000 iterations,  248707 nodes, 44428 on queue, breadth at 26
400000 iterations,  280354 nodes, 47086 on queue, breadth at 27
450000 iterations,  311547 nodes, 49139 on queue, breadth at 28
500000 iterations,  342501 nodes, 50700 on queue, breadth at 29
550000 iterations,  372922 nodes, 51733 on queue, breadth at 30
600000 iterations,  403108 nodes, 52311 on queue, breadth at 31
650000 iterations,  433618 nodes, 53050 on queue, breadth at 32
700000 iterations,  463984 nodes, 53456 on queue, breadth at 33
750000 iterations,  494564 nodes, 54359 on queue, breadth at 34
800000 iterations,  524611 nodes, 54773 on queue, breadth at 35
850000 iterations,  554828 nodes, 54958 on queue, breadth at 35
900000 iterations,  585064 nodes, 55193 on queue, breadth at 36
950000 iterations,  614566 nodes, 54790 on queue, breadth at 37
1000000 iterations,  643556 nodes, 53524 on queue, breadth at 38
1050000 iterations,  672372 nodes, 52069 on queue, breadth at 39
1100000 iterations,  701500 nodes, 50760 on queue, breadth at 40
1150000 iterations,  730769 nodes, 49438 on queue, breadth at 41
1200000 iterations,  760651 nodes, 48899 on queue, breadth at 42
1250000 iterations,  790664 nodes, 48370 on queue, breadth at 43
1300000 iterations,  820624 nodes, 47926 on queue, breadth at 44
1350000 iterations,  850645 nodes, 47603 on queue, breadth at 45
1400000 iterations,  880617 nodes, 47057 on queue, breadth at 46
1450000 iterations,  910531 nodes, 46661 on queue, breadth at 47
1500000 iterations,  940520 nodes, 46241 on queue, breadth at 48
1550000 iterations,  970586 nodes, 45878 on queue, breadth at 49
1600000 iterations,  1000626 nodes, 45348 on queue, breadth at 51
1650000 iterations,  1030354 nodes, 44614 on queue, breadth at 52
1700000 iterations,  1060201 nodes, 44022 on queue, breadth at 53
1750000 iterations,  1089699 nodes, 43543 on queue, breadth at 54
1800000 iterations,  1118872 nodes, 42375 on queue, breadth at 55
1850000 iterations,  1148292 nodes, 41268 on queue, breadth at 56
1900000 iterations,  1177217 nodes, 40098 on queue, breadth at 57
1950000 iterations,  1206432 nodes, 38531 on queue, breadth at 59
2000000 iterations,  1235686 nodes, 37739 on queue, breadth at 60
2050000 iterations,  1264742 nodes, 36654 on queue, breadth at 61
2100000 iterations,  1293983 nodes, 35484 on queue, breadth at 63
2150000 iterations,  1322943 nodes, 34334 on queue, breadth at 64
2200000 iterations,  1351804 nodes, 33118 on queue, breadth at 66
2250000 iterations,  1380383 nodes, 31994 on queue, breadth at 67
2300000 iterations,  1409466 nodes, 31139 on queue, breadth at 69
2350000 iterations,  1438334 nodes, 30281 on queue, breadth at 70
2400000 iterations,  1467507 nodes, 29591 on queue, breadth at 72
2450000 iterations,  1496835 nodes, 29003 on queue, breadth at 74
2500000 iterations,  1526456 nodes, 28992 on queue, breadth at 75
2550000 iterations,  1555916 nodes, 28484 on queue, breadth at 77
2600000 iterations,  1584903 nodes, 27624 on queue, breadth at 79
2650000 iterations,  1613876 nodes, 26653 on queue, breadth at 81
2700000 iterations,  1642911 nodes, 26065 on queue, breadth at 83
2750000 iterations,  1671512 nodes, 25088 on queue, breadth at 84
2800000 iterations,  1700318 nodes, 24124 on queue, breadth at 87
2850000 iterations,  1729411 nodes, 23453 on queue, breadth at 89
2900000 iterations,  1758261 nodes, 22843 on queue, breadth at 91
2950000 iterations,  1786999 nodes, 22047 on queue, breadth at 93
3000000 iterations,  1815948 nodes, 21362 on queue, breadth at 95
3050000 iterations,  1844745 nodes, 20582 on queue, breadth at 98
3100000 iterations,  1873273 nodes, 19810 on queue, breadth at 100
3150000 iterations,  1902028 nodes, 19352 on queue, breadth at 103
3200000 iterations,  1930524 nodes, 18593 on queue, breadth at 105
3250000 iterations,  1959083 nodes, 17959 on queue, breadth at 108
3300000 iterations,  1987762 nodes, 17430 on queue, breadth at 111
3350000 iterations,  2016942 nodes, 17274 on queue, breadth at 114
3400000 iterations,  2046092 nodes, 16957 on queue, breadth at 117
3450000 iterations,  2074787 nodes, 16445 on queue, breadth at 120
3500000 iterations,  2103478 nodes, 15650 on queue, breadth at 123
3550000 iterations,  2132185 nodes, 15092 on queue, breadth at 126
3600000 iterations,  2161428 nodes, 15188 on queue, breadth at 129
3650000 iterations,  2190340 nodes, 14896 on queue, breadth at 132
3700000 iterations,  2219362 nodes, 14723 on queue, breadth at 136
3750000 iterations,  2247899 nodes, 13945 on queue, breadth at 139
3800000 iterations,  2276601 nodes, 13347 on queue, breadth at 143
3850000 iterations,  2305187 nodes, 12865 on queue, breadth at 147
3900000 iterations,  2333835 nodes, 12411 on queue, breadth at 151
3950000 iterations,  2361786 nodes, 11172 on queue, breadth at 155
4000000 iterations,  2389960 nodes, 10456 on queue, breadth at 160
4050000 iterations,  2418139 nodes, 9669 on queue, breadth at 165
4100000 iterations,  2446455 nodes, 9303 on queue, breadth at 170
4150000 iterations,  2474639 nodes, 8695 on queue, breadth at 175
4200000 iterations,  2502681 nodes, 7973 on queue, breadth at 181
4250000 iterations,  2531134 nodes, 7661 on queue, breadth at 188
4300000 iterations,  2559334 nodes, 7050 on queue, breadth at 194
4350000 iterations,  2587431 nodes, 6601 on queue, breadth at 202
4400000 iterations,  2615483 nodes, 6063 on queue, breadth at 210
4450000 iterations,  2643758 nodes, 5862 on queue, breadth at 218
4500000 iterations,  2671929 nodes, 5454 on queue, breadth at 227
4550000 iterations,  2700175 nodes, 4893 on queue, breadth at 236
4600000 iterations,  2728494 nodes, 4443 on queue, breadth at 247
4650000 iterations,  2757000 nodes, 4160 on queue, breadth at 258
4700000 iterations,  2785007 nodes, 3533 on queue, breadth at 271
4750000 iterations,  2813042 nodes, 3006 on queue, breadth at 286
4800000 iterations,  2841280 nodes, 2641 on queue, breadth at 304
4850000 iterations,  2869508 nodes, 2279 on queue, breadth at 325
4900000 iterations,  2897574 nodes, 1708 on queue, breadth at 350
4950000 iterations,  2925624 nodes, 1267 on queue, breadth at 382
5000000 iterations,  2953765 nodes, 851 on queue, breadth at 428
5050000 iterations,  2982088 nodes, 536 on queue, breadth at 500
5100000 iterations,  3010900 nodes, 249 on queue, breadth at 641

















In [104]:



    


nx.write_gml(G_lump, 'component-subsample-{0}.gml'.format(lump.component_id), )



    











In [281]:



    


lump_centrality = pd.DataFrame({'centrality':nx.centrality.eigenvector_centrality(G_lump_full, weight='coverage')})



    











In [128]:



    


print('\n'.join(['>{0}\n{1}'.format(n, G_lump_full.node[n]['kmer'][:25]) for n,d in nx.degree(G_lump_full).items() if d == 9]))



    


















    






>499604117652391
CAGGCTTTGACCAAGTCTGAGATCA
>460032190954739
TCCACTCTTCTTGTTAGCGGAGGAG
>368241324220894
GCACATGAGTAAACACAGGCGCTAA
>53808337628903
AAGAAGGAAAGCCGTCTCCGCGCTG
>697733847628470
GACGGGCGGTGTGTACAAAGGGCAG
>289730208293299
CTCGATGAAAGTGTATGTTGCATTA
>51754744265853
ACCATCCGAAGTCAATGTATCCGTT
>752856896119578
CCCGACGCAACTGGTGAACAGATCC
>216924397683620
AGATTTACCCCTTCCCCGAAGCCTA
>450553535362826
GGTTCACAACGACGGGGATCAGAGA
>183088559218282
GGGAGCACTGTGTTGGTATGGAGGT
>693087573631078
GAGATTGCCCAGGATTTCAAGACAG
>123542523473980
TCCTTTCCAGCTTCCATCAATTCAT
>529432943248808
TGCATCATAAGTCCAAGGATTCCCA
>108787783533728
TTGGTATGATCTATGTATCCGTGAT
>623860080632718
CAGTGTCTTCGGCTGGCGCAGCAGC
>658070847750927
CTTTCCAAGAAGGATGCAAAGAAGG
>11403879321502
GCAGCCTTACGTTCGTCCAAGGTTT
>411901806036348
TTGTCCTGGTTGAGTTAGATTTGGA
>489266180559647
TCGGAGGGAACCAGCTACTAGATGG
>295157793043504
TAAGATGATGATGTCCAAGTAAGAA
>148469477140894
GCAGAGGCATCAGCCTTGTTCATGT

















In [112]:



    


sns.distplot(list(nx.degree(G_lump_full).values()), kde=False)



    


















    
Out[112]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f5528188828>










    
























In [111]:



    


sns.distplot([len(G_lump_full.get_edge_data(*e)['path']) for e in G_lump_full.edges()], kde=False)



    


















    
Out[111]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f552aff1a20>










    
























In [9]:



    


sp.graph.calc_connected_graph_size('A' * 25)



    


















    
Out[9]:








25294103

















In [89]:



    


G_rep = extract_local_compact_dbg(reverse_hash(list(rep_comp)[0], 25), sp.graph)



    


















    






0 iterations,  2 nodes, 1 on queue, breadth at 0
474 iterations,  243 nodes, 0 on queue, breadth at 97 3849 k-mers

















In [49]:



    


sp.graph.calc_connected_graph_size(reverse_hash(list(rep_comp)[0], 25))



    


















    
Out[49]:








3849

















In [84]:



    


[(G_rep.node[i], i, G_rep.degree(i), sp.graph.kmer_degree(reverse_hash(i, 25))) for i in G_rep.nodes() if sp.graph.kmer_degree(reverse_hash(i, 25)) != G_rep.degree(i)]



    


















    
Out[84]:








[({'coverage': 13,
   'kmer': 'AAGAATTAGATGTAAAAAGCCATGA/TCATGGCTTTTTACATCTAATTCTT'},
  53133526972956,
  4,
  3),
 ({'coverage': 15,
   'kmer': 'TTTCTTTTAGGTGTCTAAGTAACTG/CAGTTACTTAGACACCTAAAAGAAA'},
  379696439374887,
  4,
  3),
 ({'coverage': 4,
   'kmer': 'ACACAGAATATAAAACTAGAATAGC/GCTATTCTAGTTTTATATTCTGTGT'},
  150362779103310,
  5,
  4),
 ({'coverage': 1,
   'kmer': 'GCCGCCTTTCTGTAGTATTGAATTT/AAATTCAATACTACAGAAAGGCGGC'},
  6052222078910,
  2,
  3),
 ({'coverage': 9,
   'kmer': 'CGCTACGTGAAAGAGAGTATGGGAT/ATCCCATACTCTCTTTCACGTAGCG'},
  116627258879291,
  2,
  3),
 ({'coverage': 6,
   'kmer': 'GATGGAGCGCTCCAAGGCGCAAGAG/CTCTTGCGCCTTGGAGCGCTCCATC'},
  675025376962182,
  4,
  3),
 ({'coverage': 7,
   'kmer': 'AAAGGCGGCTGTCCGAAGGTATCAT/ATGATACCTTCGGACAGCCGCCTTT'},
  17310709511265,
  2,
  3),
 ({'coverage': 1,
   'kmer': 'GAGAGGATGATAGATTTATTGAATG/CATTCAATAAATCTATCATCCTCTC'},
  587156494359142,
  2,
  3),
 ({'coverage': 4,
   'kmer': 'CCGAGAACATGAAATATGAATGATT/AATCATTCATATTTCATGTTCTCGG'},
  26767400228463,
  4,
  3),
 ({'coverage': 12,
   'kmer': 'TATTCAAACCGGAGTATCAAACGGC/GCCGTTTGATACTCCGGTTTGAATA'},
  305675990302910,
  4,
  3),
 ({'coverage': 30,
   'kmer': 'CCCGTCGAAGCCTCTTCGTCGTCAT/ATGACGACGAAGAGGCTTCGACGGG'},
  126216421600447,
  4,
  3),
 ({'coverage': 21,
   'kmer': 'CCGCCTTTCTGTAGTATTGAATTCC/GGAATTCAATACTACAGAAAGGCGG'},
  767830443257946,
  4,
  3),
 ({'coverage': 12,
   'kmer': 'TATTGCTAAAAGATAAAGAAAGCGA/TCGCTTTCTTTATCTTTTAGCAATA'},
  307382425010412,
  4,
  3),
 ({'coverage': 16,
   'kmer': 'GATGAACTCTTCGTATTCAAACCGG/CCGGTTTGAATACGAAGAGTTCATC'},
  771153102845318,
  4,
  3),
 ({'coverage': 14,
   'kmer': 'ATTAAGATGATATATCATCTGACTC/GAGTCAGATGATATATCATCTTAAT'},
  88815915525926,
  2,
  3),
 ({'coverage': 18,
   'kmer': 'AGCCGCCTTTCTGTAGTATTGAATT/AATTCAATACTACAGAAAGGCGGCT'},
  24208888315641,
  2,
  3),
 ({'coverage': 5,
   'kmer': 'CCGTCGATTAATCCATCATCGTTTT/AAAACGATGATGGATTAATCGACGG'},
  3054243158831,
  4,
  3),
 ({'coverage': 10,
   'kmer': 'TACGTTGGCGCTACGTGAAAGAGAG/CTCTCTTTCACGTAGCGCCAACGTA'},
  331498003039027,
  4,
  3),
 ({'coverage': 31,
   'kmer': 'GAGATGCAGATGACGACGAAGAGGC/GCCTCTTCGTCGTCATCTGCATCTC'},
  899276438946622,
  4,
  3),
 ({'coverage': 1,
   'kmer': 'GCAACATATGAACATCGAATGACGT/ACGTCATTCGATGTTCATATGTTGC'},
  200208949987166,
  2,
  3),
 ({'coverage': 10,
   'kmer': 'GAAACTTACACGTCTGTTGAGGCCA/TGGCCTCAACAGACGTGTAAGTTTC'},
  556767561335638,
  2,
  3),
 ({'coverage': 9,
   'kmer': 'CGGAAGATTCATGATTTGATTTCTG/CAGAAATCAAATCATGAATCTTCCG'},
  615829693274475,
  4,
  3),
 ({'coverage': 1,
   'kmer': 'ATTCAATAAATCTATCATCCTCTCC/GGAGAGGATGATAGATTTATTGAAT'},
  96826163751322,
  2,
  3),
 ({'coverage': 25,
   'kmer': 'CTCTTTCCTTTCCCGTCGAAGCCTC/GAGGCTTCGACGGGAAAGGAAAGAG'},
  674453121319846,
  4,
  3),
 ({'coverage': 1,
   'kmer': 'CTGGCCTCAACAGACGTGTAAGTTT/AAACTTACACGTCTGTTGAGGCCAG'},
  10207903469475,
  2,
  3),
 ({'coverage': 6,
   'kmer': 'CAACATATGAACATCGAATGACGTA/TACGTCATTCGATGTTCATATGTTG'},
  331527214207447,
  2,
  3),
 ({'coverage': 3,
   'kmer': 'AAAAAAGAGTATGCTGCAGATGGTA/TACCATCTGCAGCATACTCTTTTTT'},
  220244881908,
  4,
  3),
 ({'coverage': 1,
   'kmer': 'TCCCATACTCTCTTTCACGTAGCGT/ACGCTACGTGAAAGAGAGTATGGGA'},
  203606159315452,
  2,
  3),
 ({'coverage': 6,
   'kmer': 'AGTGAGATGGCTTCGAAGGGTGGGG/CCCCACCCTTCGAAGCCATCTCACT'},
  242751110184447,
  4,
  3)]

















In [43]:



    


nx.write_gml(G_rep, 'component-{0}.gml'.format(rep_comp.component_id), )



    











In [28]:



    


sns.distplot([len(G_rep.get_edge_data(*e)['path']) for e in G_rep.edges()], kde=False)



    


















    
Out[28]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f552b1f7c88>










    
























In [29]:



    


lyt = nx.layout.fruchterman_reingold_layout(G_rep, scale=3.0)
nx.draw_networkx_edges(G_rep, lyt)
nx.draw_networkx_nodes(G_rep, lyt)



    


















    
Out[29]:








<matplotlib.collections.PathCollection at 0x7f552b064a90>










    
























In [ ]:



    


%%cython



    











In [299]:



    


%%cython

from khmer._oxli.streaming import SplitPairedReader, FastxParser, Sequence
from khmer._oxli.partitioning import StreamingPartitioner
from khmer import reverse_hash

def kmers(sequence, K, which=[0,1,2,-27,-26,-25]):
    if which is None:
        which = range(0,len(sequence)-K+1)
    for i in which:
        yield sequence[i:i+K]

def read_degrees_by_position(left, right, part, limit=1000000):
    
    reader = SplitPairedReader(FastxParser(left), FastxParser(right))
    K = part.graph.ksize()

    data = {}
    for n, paired, first, second in reader:
        if n % 50000 == 0:
            print(n, '...', sep='')
        if limit and n > limit:
            break
        base = first.name
        fname = base + '.L'
        rname = base + '.R'
        data[fname] = []
        data[rname] = []
        for kmer in kmers(first.sequence, K):
            data[fname].append(part.graph.kmer_degree(kmer))
        for kmer in kmers(second.sequence, K):
            data[rname].append(part.graph.kmer_degree(kmer))
    
    return data

def tag_hexamer_decomposition(comp, K=25):
    
    decomp = {}
    for tag in comp:
        decomp[tag] = [hexamer for hexamer in kmers(reverse_hash(tag, K), 6, which=None)]

    return decomp



    











In [13]:



    


read_degrees = read_degrees_by_position('/store/biodb/MMETSP/data/Asterionellopsis_glacialis/SRR1300451/SRR1300451_1.fastq',
                                        '/store/biodb/MMETSP/data/Asterionellopsis_glacialis/SRR1300451/SRR1300451_2.fastq',
                                        sp)



    


















    






0...
50000...
100000...
150000...
200000...
250000...
300000...
350000...
400000...
450000...
500000...
550000...
600000...
650000...
700000...
750000...
800000...
850000...
900000...
950000...
1000000...

















In [14]:



    


read_degrees_df = pd.DataFrame(read_degrees)



    











In [28]:



    


read_degrees_melted = pd.melt(read_degrees_df.reset_index(), id_vars=['index'], value_name='degree',
                  value_vars=list(range(6)), var_name='position')



    











In [32]:



    


ttest_ind(read_degrees_melted.query('position in [0,1,2]').degree,
          read_degrees_melted.query('position in [3,4,5]').degree)



    


















    
Out[32]:








Ttest_indResult(statistic=11.34514920371045, pvalue=7.8447385924038063e-30)

















In [36]:



    


sns.distplot(read_degrees_melted.query('position in [0,1,2]').degree, bins=8, kde=False)
sns.distplot(read_degrees_melted.query('position in [3,4,5]').degree, bins=8, kde=False)



    


















    
Out[36]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f45e6e19898>










    
























In [35]:



    


g = sns.FacetGrid(read_degrees_melted, col_wrap=3, col='position')
g.map(sns.distplot, 'degree', bins=8, kde=False)



    


















    
Out[35]:








<seaborn.axisgrid.FacetGrid at 0x7f45e709aa20>










    
























In [45]:



    


lump_tags = set(lump)



    











In [40]:



    


from khmer import Nodegraph



    











In [41]:



    


lump_filter = Nodegraph(25, 1e8, 4)
for tag in lump:
    lump_filter.count(reverse_hash(tag, 25))



    











In [233]:



    


def tag_read_positions(left, right, comp, K=25, limit=1000000):
    
    reader = SplitPairedReader(FastxParser(left), FastxParser(right))
    tags = set(comp)
    
    data = []
    for n, paired, first, second in reader:
        if n % 50000 == 0:
            print(n, '...', sep='')
        if limit and n > limit:
            break
        for i, kmer in enumerate(kmers(first.sequence, K, which=None)):
            tag = hash(Kmer(kmer))
            if tag in tags:
                data.append((tag,i,'L'))
        for i, kmer in enumerate(kmers(second.sequence, K, which=None)):
            tag = hash(Kmer(kmer))
            if tag in tags:
                data.append((tag, i + len(first.sequence), 'R'))
    return data



    











In [234]:



    


lump_tag_positions = tag_read_positions('/store/biodb/MMETSP/data/Asterionellopsis_glacialis/SRR1300451/SRR1300451_1.fastq',
                                        '/store/biodb/MMETSP/data/Asterionellopsis_glacialis/SRR1300451/SRR1300451_2.fastq',
                                        lump, limit=500000)



    


















    






0...
50000...
100000...
150000...
200000...
250000...
300000...
350000...
400000...
450000...
500000...

















In [238]:



    


lump_tag_positions = pd.DataFrame(lump_tag_positions, columns=['tag', 'position', 'pair'])



    











In [300]:



    


lump_hexamers = tag_hexamer_decomposition(lump)



    











In [284]:



    


lump_centrality.reset_index(inplace=True)
lump_centrality.rename(columns={'index': 'tag'}, inplace=True)



    











In [285]:



    


lump_centrality.head()



    


















    
Out[285]:










  
    

      
      tag
      centrality
    

  
  
    

      0
      0
      0.000173
    

    

      1
      1
      0.000078
    

    

      2
      2
      0.000114
    

    

      3
      3
      0.000091
    

    

      4
      4
      0.000011
    

  




















In [282]:



    


sns.distplot(log(lump_centrality['centrality']))



    


















    






/home/camille/miniconda/envs/dev.py3.khmer.assembly/lib/python3.5/site-packages/statsmodels/nonparametric/kdetools.py:20: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  y = X[:m/2+1] + np.r_[0,X[m/2+1:],0]*1j










    
Out[282]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f4224279780>










    
























In [259]:



    


lump_tag_positions



    


















    
Out[259]:










  
    

      
      tag
      position
      pair
    

  
  
    

      0
      848893025205450
      0
      L
    

    

      1
      17872380293928
      1
      L
    

    

      2
      71489521175712
      2
      L
    

    

      3
      104894234772756
      3
      L
    

    

      4
      17932431968773
      4
      L
    

    

      5
      71729727875094
      5
      L
    

    

      6
      286918911500378
      6
      L
    

    

      7
      21775739158888
      7
      L
    

    

      8
      87102956635555
      8
      L
    

    

      9
      348411826542222
      9
      L
    

    

      10
      267747399326264
      10
      L
    

    

      11
      415548276381594
      11
      L
    

    

      12
      666837022516710
      12
      L
    

    

      13
      247334954241548
      13
      L
    

    

      14
      674996011506270
      14
      L
    

    

      15
      579659547336910
      15
      L
    

    

      16
      66838375662392
      16
      L
    

    

      17
      267353502649570
      17
      L
    

    

      18
      526004237991322
      18
      L
    

    

      19
      412976036208486
      19
      L
    

    

      20
      222125566933122
      20
      L
    

    

      21
      236917234796022
      21
      L
    

    

      22
      176309350402086
      22
      L
    

    

      23
      507388536418399
      23
      L
    

    

      24
      126847134104599
      24
      L
    

    

      25
      24799357307271
      25
      L
    

    

      26
      194165269200002
      50
      R
    

    

      27
      668263103068491
      64
      R
    

    

      28
      421252598588718
      65
      R
    

    

      29
      542122190527912
      66
      R
    

    

      ...
      ...
      ...
      ...
    

    

      10271927
      406398407886866
      22
      L
    

    

      10271928
      307050254840192
      23
      L
    

    

      10271929
      102301112518147
      24
      L
    

    

      10271930
      409204450072590
      25
      L
    

    

      10271931
      298665135828455
      50
      R
    

    

      10271932
      74666283957113
      51
      R
    

    

      10271933
      18666570989278
      52
      R
    

    

      10271934
      4666642747319
      53
      R
    

    

      10271935
      204784904606806
      54
      R
    

    

      10271936
      819139618427226
      55
      R
    

    

      10271937
      1024758660023658
      56
      R
    

    

      10271938
      721334919566762
      57
      R
    

    

      10271939
      561854999061843
      58
      R
    

    

      10271940
      282359644641952
      59
      R
    

    

      10271941
      3538671725184
      60
      R
    

    

      10271942
      14154686900737
      61
      R
    

    

      10271943
      56618747602948
      62
      R
    

    

      10271944
      76414988839065
      63
      R
    

    

      10271945
      300578723920422
      64
      R
    

    

      10271946
      245900126060816
      65
      R
    

    

      10271947
      370629891133346
      66
      R
    

    

      10271948
      374132449493992
      67
      R
    

    

      10271949
      656483065794810
      68
      R
    

    

      10271950
      727070719870014
      69
      R
    

    

      10271951
      726049860370494
      70
      R
    

    

      10271952
      652399627796730
      71
      R
    

    

      10271953
      357798697501672
      72
      R
    

    

      10271954
      139179561843716
      73
      R
    

    

      10271955
      95279625813638
      74
      R
    

    

      10271956
      381118503254552
      75
      R
    

  



10271957 rows × 3 columns



















In [286]:



    


lump_stats = pd.merge(lump_centrality, lump_tag_positions, on='tag')



    











In [287]:



    


lump_stats['centrality'] = np.log(lump_stats['centrality'])



    











In [288]:



    


lump_stats['degree'] = lump_stats.tag.apply(lambda tag: sp.graph.kmer_degree(reverse_hash(tag, 25)))



    











In [303]:



    


from khmer import reverse_complement



    











In [304]:



    


lump_stats['kmer_f'] = lump_stats.tag.apply(lambda tag: reverse_hash(tag, 25))
lump_stats['kmer_r'] = lump_stats.kmer_f.apply(lambda kmer_f: reverse_complement(kmer_f))



    











In [337]:



    


freqs = {}
for i in range(100):
    freqs[i] = {'A': 0, 'C': 0, 'G': 0, 'T': 0}



    











In [338]:



    


def base_pos_count(row):
    for i, (b1,b2) in enumerate(zip(row.kmer_f, row.kmer_r)):
        freqs[i+row.position][b1] += 1
        freqs[i+row.position][b2] += 1



    











In [ ]:



    


_ = lump_stats.apply(base_pos_count, axis=1)



    











In [343]:



    


freqs_df = pd.DataFrame(freqs).T



    











In [349]:



    


freqs_df = freqs_df.divide(freqs_df.sum(axis=1), axis=0)



    











In [350]:



    


freqs_df.plot()



    


















    
Out[350]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f3c17cc9c18>










    
























In [305]:



    


lump_stats.head()



    


















    
Out[305]:










  
    

      
      tag
      centrality
      position
      pair
      degree
      kmer_f
      kmer_r
    

  
  
    

      0
      0
      -8.664243
      73
      R
      8
      AAAAAAAAAAAAAAAAAAAAAAAAA
      TTTTTTTTTTTTTTTTTTTTTTTTT
    

    

      1
      0
      -8.664243
      74
      R
      8
      AAAAAAAAAAAAAAAAAAAAAAAAA
      TTTTTTTTTTTTTTTTTTTTTTTTT
    

    

      2
      0
      -8.664243
      75
      R
      8
      AAAAAAAAAAAAAAAAAAAAAAAAA
      TTTTTTTTTTTTTTTTTTTTTTTTT
    

    

      3
      0
      -8.664243
      0
      L
      8
      AAAAAAAAAAAAAAAAAAAAAAAAA
      TTTTTTTTTTTTTTTTTTTTTTTTT
    

    

      4
      0
      -8.664243
      50
      R
      8
      AAAAAAAAAAAAAAAAAAAAAAAAA
      TTTTTTTTTTTTTTTTTTTTTTTTT
    

  




















In [276]:



    


figsize(12,8)
sns.distplot(lump_stats.query('pair == "L"').sample(n=500000)['position'], kde=False, bins=25)
title('distribution of start positions of lump $k$-mers in reads', )



    


















    
Out[276]:








<matplotlib.text.Text at 0x7f41d5f105f8>










    
























In [272]:



    


sns.jointplot(x='position', y='degree', data=lump_stats.sample(n=1000000).query('pair == "L"'))



    


















    
Out[272]:








<seaborn.axisgrid.JointGrid at 0x7f41f6cb3470>










    
























In [290]:



    


sns.violinplot(x='position', y='centrality', data=lump_stats.sample(n=1000000).query('pair == "L"'))



    


















    
Out[290]:








<matplotlib.axes._subplots.AxesSubplot at 0x7f421b313828>










    
























In [291]:



    


sns.jointplot('centrality', 'position', lump_stats.sample(n=1000000).query('pair == "L"'), kind='hexbin', 
              marginal_kws={'bins': 25})



    


















    
Out[291]:








<seaborn.axisgrid.JointGrid at 0x7f41aae6b208>










    
























In [223]:



    


lump_stats = pd.concat([lump_centrality.centrality, lump_avg_pos], axis=1)
lump_stats.rename(columns={0:'avg_pos'}, inplace=True)
lump_stats.dropna(inplace=True)



    











In [227]:



    


sns.jointplot('centrality', 'avg_pos', data=lump_stats, kind='kde')



    


















    






/home/camille/miniconda/envs/dev.py3.khmer.assembly/lib/python3.5/site-packages/statsmodels/nonparametric/kdetools.py:20: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  y = X[:m/2+1] + np.r_[0,X[m/2+1:],0]*1j










    
Out[227]:








<seaborn.axisgrid.JointGrid at 0x7f421d5611d0>










    
























In [228]:



    


sns.jointplot('centrality', 'avg_pos', data=lump_stats)



    


















    
Out[228]:








<seaborn.axisgrid.JointGrid at 0x7f421d561588>










    
























In [ ]:

Content source: camillescott/boink

Similar notebooks:

notebook.community | gallery | about