U _@sVdZdZdddddgZdddZd Zdd d ZddZddZedkrReedS)III )bedgff r TcCs|tks|tkr"|rdStdgS|dkr<|r2dStdgS|dkrV|rLdStdgS|t|t?}|t?}tdg}tD]L}|r||kr||S|tt||||d|tL}|tL}q||S)aS Uses the definition of a "genomic bin" described in Fig 7 of http://genome.cshlp.org/content/12/6/996.abstract. Parameters ---------- one : boolean If `one=True` (default), then only return the smallest bin that completely contains these coordinates (useful for assigning a single bin). If `one=False`, then return the set of *all* bins that overlap these coordinates (useful for looking for features that could intersect) fmt : 'gff' | 'bed' This specifies 1-based start coords (gff) or 0-based start coords (bed) rr) MAX_CHROM_SIZEset COORD_OFFSETS FIRST_SHIFTOFFSETSupdatelistrange NEXT_SHIFT)startstopfmtonebinsoffsetr,lib/python3.8/site-packages/gffutils/bins.pyr:s.       rc CsttD]\}}|}zt|d}Wntk r<|}YnXdt|t}|}|dd}}|dkrv|dd}}|dkr|dd}}d||f}d|}td jftqd S) zP Useful for debugging: how large is each bin, and what are the bin IDs? riZKbZMbZGbz(%s %s)z%s bpzTlevel: {i:1}; bins {binstart:<4} to {binstop:<4}; size: {actual_size:<12} {size:<6}N) enumerater IndexErrorrrprintformatlocals)irZbinstartZbinstopZbin_sizeZ actual_sizesuffixsizerrrprint_bin_sizess$  r&cCstdddddksttdddddks,ttdddddtdd d d dgksRttdddddttksnttd d dddksttd tdd tdddtdksttdddddksttd td tddtddksttdddddtdd d d dgksttdddddtd d dddtdgksDttd dddd ks\ttd tdd tddtdksttd tdd ttdddtdksttd tdd ttddtd ksttdddddksttddddtdgksttdddddks.ttddddtdgksLttdddddksdttddddtdgksttdddddksttdddddkstdS)Nrrr )rrr F)rrrrririiJiT)riiiiH"r )rAssertionErrorr rr differencerrrrrtestst&            r*__main__N)r T) rrrrr rr&r*__name__rrrr"s  IL