ó
£žbc @€ sä d d l m Z d d l Z d d l Z d d l Z d d l m Z d d l m Z m Z m
Z
m Z y d d l m Z Wn! e
k
r d d l m Z n Xd d l m Z m Z d d l m Z d d l m Z m Z d d l m Z d d
l m Z m Z m Z d d l m Z m Z m Z d d l! m" Z" d
Z# d Z$ d Z% d Z& d Z' d „ Z( d „ Z) d „ Z* d „ Z+ d „ Z, d „ Z- d „ Z. d „ Z/ d „ Z0 d „ Z1 d „ Z2 d „ Z3 d „ Z4 d „ Z5 d „ Z6 d „ Z7 d S(! iÿÿÿÿ( t with_statementN( t defaultdict( t joint existst dirnamet realpath( t OrderedDict( t qutilst qconfig( t get_aux_out_fpaths( t create_minimap_output_dirt parse_cs_tag( t get_chr_lengths_from_fastafile( t get_assembliest check_misassembled_blockst Alignment( t get_path_to_programt is_non_empty_filet relpath( t COVERAGE_FACTORs
circos.pngg{®Gáz¤?g¸ …ëQ¸®?iPà c C€ s( t | j ƒ ƒ } d } t | d ƒ } t | d ƒ h } x^ | j ƒ D]P \ } } | j d j d d | | d t | ƒ d g ƒ d ƒ t | | ƒ } qF WWd QXt | d
ƒ } t | d ƒ [} | j d ƒ | j d ƒ t j r | d
k r | j d ƒ nE | d k r | j d ƒ n) | d k r8| j d ƒ n
| j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d! ƒ | j d" ƒ | j d# ƒ Wd QX| | | f S($ Ni s reference.karyotype.txtt ws t chrt -t 0t lgreys
s
ideogram.confs
s
i s
default = 0r
i s default = 0.005r
id s default = 0.001r
s default = 0.0005r
s break = 0.005r
s
s thickness = 30p
s stroke_thickness = 2
s stroke_color = black
s fill = yes
s radius = 0.85r
s show_label = no
s label_font = default
s- label_radius = dims(ideogram,radius) + 0.05r
s label_size = 36
s label_parallel = yes
s band_stroke_thickness = 2
s show_bands = yes
s fill_bands = yes
s (
t lent keysR t opent itemst writet strt maxR t
prokaryote( t chr_lengthst
output_dirt num_chromosomest max_lent karyotype_fpatht out_ft namet seq_lent ideogram_fpath( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt create_ideogram$ sF 5
c C€ so t | d ƒ } t | d ƒ K} | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d
ƒ | j d ƒ | j d ƒ | j d
ƒ | j d ƒ | j d ƒ | d d( k rñ d | d } d } n d | } d } | j d t | ƒ d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d ƒ | j d! ƒ | j d" ƒ | j d ƒ | j d ƒ | j d ƒ | j d# ƒ | j d ƒ | j d$ ƒ | j d! ƒ | j d% ƒ | j d& | ƒ | j d ƒ | j d ƒ | j d' ƒ Wd QX| S() Ns
ticks.confR s show_ticks = yes
s show_tick_labels = yes
s show_grid = no
s
s skip_first_label = yes
s skip_last_label = no
s% radius = dims(ideogram,radius_outer)
s tick_separation = 2p
s min_label_distance_to_edge = 0p
s label_separation = 5p
s label_offset = 5p
s label_size = 12p
s thickness = 3p
i
i g ð?t Mbpt kbps label_multiplier = s
s
s
spacing = 1u
s color = dgrey
s size = 12p
s show_label = no
s format = %s
s
s
spacing = 5u
s color = black
s size = 18p
s show_label = yes
s label_size = 24p
s spacing = 10u
s size = 24p
s label_size = 32p
s suffix = " %s"
s i@B ( R R R R ( t chrom_unitsR" t ticks_fpathR& t label_multipliert suffix( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt create_ticks_confL s^
c C€ s t | d ƒ } d d d g } t | d ƒ æ } t ƒ } d d l m } xG | j j ƒ D]6 \ } } | | k r~ g | | s t |i i iÿÿÿÿi t *t +R s R s
( i i ( R Rk t fpathR R t detailed_contigs_reports_dirnameR
R R t show_snpsRQ R R RV RS RU R R R‘ R R9 R R ( R{ R• R| t root_dirR" t assembly_labelt aligner_dirpatht coords_basenamet _t coords_filtered_fpatht mismatches_fpatht mismatch_density_by_chromt coords_fileRc t s1R= t cigart ref_post opt n_basesR& t density_listRK RL R? t density( ( R| R• s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt create_mismatches_plot
sF
#
5D
?c € s, g } d } | s | | f Sx| D]û} t | | j d ƒ } t | j ƒ d k rZ q# n d } t ‡ ‡ f d † ƒ } t | d ƒ j}
xØ | j D]Í } | j r¾ | j | j k r¾ | j n | j } | | j k rã | j | n d } | sõ q” n xi t
| j ˆ t | j
ˆ d t | | ƒ ƒ ƒ D]3 }
|
t | | ƒ k r*| | |
c d 7 }
t |
j ƒ d k rj | j |
j | f ƒ | d 7} qj qj W| rÈ | j d | f ƒ n t | d ƒ X } | j t | j
ƒ ƒ d d d j g | D] \ } } d | | f ^ qƒ ƒ Wd QXt | d
ƒ }
t |
d ƒ % } | j d t | | ƒ d ƒ Wd QX|
| f S(
Ns
labels.txti R{ i t coverageR s
0 0 null t ,s
track%d=%ss
label.confsR z = 10
type = text
label_size = 30p
label_font = bold
label_parallel = yes
file = sÍ
r0 = eval(sprintf("%fr+5p", conf(conf(., track_idx)_pos)))
r1 = eval(sprintf("%fr+500p", conf(conf(., track_idx)_pos)))
condition = 1
value = eval(var(conf(., track_idx)))
( R R9 R8 R R Rº R¹ R R RP R R ( R! Rl R¾ t coverage_fpathR" t labels_txt_fpatht track_labelst plot_idxR? R{ RÀ R& Rv t labels_conf_fpath( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt
create_labels^ s& !
W c C€ sf | d d k r d } nI | d d k r2 d } n0 | d
k rG d } n | d k r\ d } n d
} | S( Ni i
i i N i i' iˆ i iè id i áõi áõi áõi@B ( ( R| R• ( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt set_window_size s c
C€ sþ t | d ƒ } t | d ƒ Ú } | j d t j r9 d n d | | f ƒ t j rc | j d ƒ n | j d ƒ x8 t | ƒ D]* \ } } | j d | d | j f ƒ q} W| j d
ƒ | rÎ | j d ƒ n | rô | j d | ré d
n d ƒ n Wd QX| S( Ns
legend.txtR sn 1) The outer circle represents reference sequence%s with GC (%%) heatmap [from %d%% (white) to %d%% (black)].
t sRB s= Color bars help to distinguish between different references.
s
2) Assembly tracks:
s assembly%d - %s
i sj Assembly tracks are combined with mismatches visualization: higher columns indicate larger mismatch rate.
sK
3) User-provided genes. A darker color indicates higher density of genes.
s:
%d) The inner circle represents read coverage histogram.
i i ( R R R R t is_combined_refR9 Rv (
Rl RŠ R‹ R¾ RË R" t legend_fpathR& R? R{ ( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt
create_legend s
"
&c + C€ sk t | d ƒ } t | ƒ s+ t j | ƒ n t | ƒ } t | | ƒ \ }
} } |
dB k rd dC }
n |
dD k ry dE }
n d }
t |
| ƒ } t | j ƒ ƒ } t | ƒ } t
| | ƒ \ } } g | D] } t | | | ƒ ^ qÈ } | sð d St
| | ƒ \ } } } } t | | | | ƒ \ } } g | D] } t | | | | | ƒ ^ q-} t | | | | ƒ \ } } t t | | | | g ƒ } t | | | | | ƒ \ } } t | d ƒ } d }! d }" t g t | ƒ }# | rýt |# d
<|# t g t | ƒ 7}# n | r t |# d
<|# j t ƒ n t |# d
ƒ |$ j d? ƒ |$ j d@ ƒ |$ j d- ƒ |$ j dA ƒ Wd QXt$ | | | | | | ƒ }* | |* f S(F Nt datai
i i i iè s circos.confgffffffî?i iÿÿÿÿR s+ <>
s <>
s karyotype = %s
s chromosomes_units = %d
s" chromosomes_display_default = yes
s track_width = s
s track%d_pos = %f
s
s dir = %s
s
file = %s
s
png = yes
s svg = no
s radius = 1500p
s angle_offset = -90
s auto_alpha_colors = yes
s auto_alpha_steps = 5
s background = white
s
s
s
s r0 = 1r - 50p
s r1 = 1r - 30p
s
s
t etcs housekeeping.confs max_points_per_track* = %d
s max_ideograms* = %d
s
s layers_overflow = collapse
s
s track_idx = track%d
s
s type = tile
s thickness = 50p
s stroke_thickness = 0
s layers = 1
s$ r0 = eval(sprintf("%.3fr",conf(tracks _pos) - conf(track_width)))
s$ r1 = eval(sprintf("%.3fr",conf(tracks _pos)))
s type = histogram
s thickness = 1
s fill_color = vlyellow
i s type = heatmap
s color = ylorbr-9
s fill_color = vlblue
s color = greys-6
s scale_log_base = 1.5
s r0 = 1r - 29p
s
r1 = 1r - 1p
s
i@B i † i † i' (% R R t ost makedirsR R* R1 R“ R RÑ Rr R RQ R RÆ R¸ R R t
MAX_POINTSRÐ t TRACK_INTERVALR t BIG_TRACK_INTERVALR8 R R R R t TRACK_WIDTHR’ t circos_png_fnameR RÓ RA R R9 RÕ (+ t ref_fpathRm Rn R" R… R¾ R” t loggerR† R! R$ R% R) R- R. R| R• Rl t
contig_pointsR{ t alignments_fpathsRŠ R‹ t gc_pointsR¿ t gene_pointst mismatches_fpathsR— t
cov_pointsR‰ t labels_fpathRÍ R} t radiusRÎ t track_intervalsR& R? R: Rv t alignments_confRÁ t circos_legend_fpath( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt create_conf¢ sþ %+
c C€ s& t | ƒ s t j | ƒ n t | | | | | | | | ƒ \ } } t d ƒ }
|
sv | j d | d | d ƒ d S|
d | g } t | d ƒ } t | d ƒ }
t | t ƒ } t j
| d t | d ƒ d
t |
d ƒ ƒ} | d k rt | ƒ r| | f S| j d | d
|
d ƒ d Sd S( Nt circossØ Circos is not installed!
If you want to create Circos plots, install Circos as described at http://circos.ca/tutorials/lessons/configuration/distribution_and_installation and run the following command:
circos -conf s
The plot legend is saved to s
s -confs
circos.logs
circos.errt stdoutR t stderri s& Circos diagram was not created. See s and s for details( NN( NN(
R RØ RÙ Rì R t warningRQ R RÞ R t call_subprocessR R ( Rß Rm Rn R… R¾ R” R" Rà R} Rë t circos_exect cmdlinet log_fpatht err_fpatht circos_png_fpatht return_code( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt do- s" ' -
(8 t
__future__R RØ t reRƒ t collectionsR t os.pathR R R R R t ImportErrort% quast_libs.site_packages.ordered_dictR6 R R t! quast_libs.ca_utils.align_contigsR t quast_libs.ca_utils.miscR
R t quast_libs.fastaparserR t quast_libs.icarus_utilsR
R R t quast_libs.qutilsR R R t quast_libs.reads_analyzerR RÞ RÝ RÛ RÜ RÚ R* R1 RA Ri Rr R R R R¸ RÆ RÈ RÐ RÑ RÕ Rì Rø ( ( ( s? /home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/circos.pyt sJ "
( 3 ( ( ! # ‹