From 11560560a2bcc0ffc521df53706306e7ede1f485 Mon Sep 17 00:00:00 2001 From: hjulienn <hanna.julienne@pasteur.fr> Date: Wed, 23 Mar 2022 13:12:04 +0100 Subject: [PATCH] results from simulation --- README.md | 6 ++-- reports/figures/meta_analysis_perf.png | Bin 30542 -> 26121 bytes .../Meta_analysis_imputation_performance.csv | 12 +++---- .../tables/random_masking_real_data_error.csv | 34 ++++++++++++++++++ src/visualization/Draw_Imputation_quality.R | 24 +++++++++---- src/visualization/Draw_loci_real_data.R | 2 +- .../Imputation_strategy_real_data.R | 3 -- 7 files changed, 61 insertions(+), 20 deletions(-) diff --git a/README.md b/README.md index f3b863c..5baf3d5 100644 --- a/README.md +++ b/README.md @@ -45,6 +45,7 @@ Project Organization │ │ │  ├── features │  │  └── add_hg_37_pos.py : perform liftover to hg37 (ad GnomAD is in hg37) + | |  |__ retrieve_filled_out.py : Retrieve european studies with more than 80% of filled out SNPS │ │ │  ├── models │ │ │ @@ -55,9 +56,7 @@ Project Organization impute Zscores using simulated data : mask 50 SNPs on 200, reimpute then save results assuming different sample size (based on hgcovid consortium or always 100) |__ Impution_real_data.py - - - + |__ Imputation_test_real_data.py : mask 10% of SNPs around significant loci and impute them back │ │ │  └── visualization <- Scripts to create exploratory and results oriented visualizations │  └── Draw_LD.R : Draw the LD matrix used to simulate Data @@ -65,6 +64,7 @@ Project Organization | (intrinsic variability due to sample size) | |__ Draw_signal_variability.R : draw imputed signal | ( variability due to sample size + imputation error) + |_ Imputation_strategy_real_data.R : Compare meta_analysis based on original Z-scores vs Imputed Z-scores │ └── tox.ini <- tox file with settings for running tox; see tox.readthedocs.io diff --git a/reports/figures/meta_analysis_perf.png b/reports/figures/meta_analysis_perf.png index 18fd66c3b67072e0284be3e3c840287cfa7af730..ac168eb5f20df2ccca6e05e8d92f948752795dfe 100644 GIT binary patch literal 26121 zcmcG#XH=9+ur`Xi0R;(4k_;j_NtP@ypyW6pk_15{D+~fdQgSws1O#NrSwKLtsAK_$ zoKX=Na+(1d!rUI*y7xKvp7Wiv)^~rnmdx8--PK)PRbBl&k$O5e$VeDS2nYzs)Neuz z2nf!gz(2OLXTUGv?e8eS9|U?@MylY?Q>RXyK7IPknKOU=^%o%_ArTP~F){JkvuDqp zJ4Zr7LP|<{{``3|GBR><ataCxN=nKL7cN}9c#(>Vikh05hK7cgmX?l=?$V`8^z`%$ z3=E8nj7&^S%*@OzEG(?7tZZy-?Ck89FJI>1;J9+-3MVHg7Z(>dH#ZLt4=*n-A0Hn- zKfi#0fS{nDkdTnDu<+HZS4BibL`6l##KgqK#U&&pBqb%Kq@=E0yCy9yEh8f%D=RA} zCnqm2ub`lysHk}T`gJ8GC1qu06%`c-1fr^{s-~uP<Hn7fH*cz|t7~XzXliO|X=!O| zYwPIf=<4d~>FMd~>l+vt7#bQH85u#LP-A0b6BCnLw{Dr5n%=&B+sw?&+}zy4!ot$h z(#p!}&Ye5f*48#QHh1sdwY9aiv$MN*@1DKA{r&s*9UL4S9UYyVoSdDVU0hsTU0vPW z+#Wo5@bKY7cXxLW4-Xg&=IQDA=+PrDFE4LzZyz5YI2`Wl>+9#|=kM<y5D*X;82I?{ z<0ns^1O){J2M33QgoK8MhJ}TNhlfW*L_|hLMny$MM@PrR#KgwN#>K@wefsp-vuDqr zKY#J!MSOgGLPA1fVq#KKQgU+g%a<=xQc_Y=Q`6GYUcGvio}QkOk&&61`TF(itgI{q z0+F4aos*N3o12@LmzSTPUr<m`SXfw8R8(ACTvAd}T3T9GR#skK{^re_ii(QL%F3#$ zs<&_7R##Wo)YR10*1mi9uCA`GzP`Spp`o#{v8k!4xw*NerRDwm_pPn1A3l8e`0?YX zPoLV_+S=ROJ32Z#J3G6&y1Ki&dwP04fBxLt+uPUI*Wcg&<;xc&5;-t1FgQ3kG&D3k zJUlWoGCDdsHa0dsK0YxqF*!LoH8nLoJv}osGdnx`_3PJf-@bkS{(Ww4?#GWG^Yily z3k!>jizpOoX=!PBd3j}JWp#CRZEbCReSKqNV{>zJYinzJdwXYR2aQH!Fqqxl-MziN zpFe-@@9$%=*n@+E!^6X)qod>FV;m0m{f&h`0RbC<Iz-9H`_)=HaU=Ece#_IsmMb{) zJsX6Gqj>6lL}U~H(ZKUBs+7QblvI^2T=*39xS&9v`Z+)CCz67*HzS!TFF<l`woEnG zuDdKMuDufT^y;hg3ZL)W+xNR;@733^G49ppVgKk6NfyE5dlb@R%mf6sED_|_&p$>H zd?NTi{O4K@MA&ce*x7KjY)z*sk9bXM+#4CA^@yL;!8vy!W3br%8{Jw1S;1Q=+X^!^ zZx3;91~d68g)(<$TXJ4xt@#}MTeIm1#6_jMy@IN#A~97qRzb4q^$VR07q?EE-sdY_ zNay>`7=fzvd6`bz78F(as?2>RX2sY2yYpxI7%Cb|I$j6(mY6oUFMp1l)R}dtnS8KS zdZxxH)Iw(bzH`VhIsK3@EkepVR~(tIGt9E#(;U2u>t)2$sX4EGZ8iBkHn5MCl>Pi> zL50e=zu;-H`2ypV+?x%Zx%9RM2zrbB6jE`?bl{*V4YhEgowa|0ul`PvgwvPslIJyv zwr?lDgqAVe2@qsgT&TH!I{hLAGOf(Oqu*h?G-2RK8mdD6S^bgOy~aY82Jg%1l6N(4 z$l$&$)Tp_ASWdw_Ow`QZMckga(^%!59HB77&V~)+FY;AHT7TMHiAgn{l0l57FT8?& z^x^P7g<6w5pT9rdcKYfEmdqPg(cvg<lft^eWuaLm{stGsQ(@;FSJoYE%&qLmvrcN_ z)|k)qW$yGpcAG1|c-wCYUGQ&IQ9)bwYz!PT4Y&~*V!vx)zkA<RbLlwK_`I3`?s8vK z@lxG}hKyGwa_p>`TxefREK1vm*JaQjQxMj=VNkzNIeIYV)Ld8*l?ur><1nE(86^Z& z<GwUBI!NEoxFImZiFoQ25uwvT!E`lk1jA78d?T|VERFO|7w>d~SNUjuP3$%QFI4xI zCeAO{$CA_gPTPWYGnV}xr?7RQ$$>97Qb*W-;nVu`;!uX(xKp9XXhbQqJLYkC?{lXR z3D(p%a5ev_`2pS)WQ%>qR$q-*WY{X5MHHF-(G$5J`lUu^{I?VH*1EE8qLy<7OcJs& zU#O30(vSLh&qS;%Hye~aaxJ^`$g@)2_hqkmGrZ9}z%CR;@R*Abafb2?n+$=-YhKmz zLc>T69#tqWVWDc_U&EY;ZX>h2i@EPgmZ6BtX2xX0d?|bhyV3=M#qNnb)ppQjb7O5@ zc%?CCk>1i>9%zZVpay3_+0`PgE}O)H<q|7nnr#tG;?ktII*)X5O^JpPVkimv#KTkx z{W&KHbvmN&F>2sCj0m!L9mYs0|H6pUf&fWRX~BU1z!c%CqK)p}rcB`a@YdZ>FYh9> zggS>d_o9*U`@9(VsGX3ZLzaQjh;Z>mrw`Q4GFn=u>R3t3s02i1fa!cyK~4-R>pg?T zV~KhKCz#r9Ia1{mlJGwba0YSX1qp(bP38<j^9&d{{S4COI{K!866_Sj{yNM7|G|T; zo-Bezoi9S@F)HRJq@Pj#4sR&xLhi}@BOqwWkZ0AD@`zt^*N?C)@zGfi{_|ZKO=N+V zVuM2p+3E=pVwcBZR>NUaI_$rHcr*G(+yCuk&PNE5;b%z>sXc|{P$vZ1rJq9PK0%$2 z=zD^aWP@M%H6~lVIR1k)<unoG=@l9X9e%}=i4bCJGTb!?69vr8CkGOLG>nO8&ZUx$ zqWsQ6^k+xtFEnw8cz0n$fA?2YgGkA9_m~T-{9Y_u2>MjmqKVWy%P`dQUwhjbfwJ@3 zO*LC*{_Bq<)1NyNeb4W<n|zsi7+?IcMFKw?u-3L&C{Hrj45QE8rNhu~vkYKQ@XJ-H z2aDw%De25j%V!wGc(NC0ggMOTy;0-j%K~OVYg&}5ENQoT#i7$Yujq1Ofkw4mE_H5A zNq1=^cJs`eJUM?x;anO%eaD`nTbvp$rZjH~!iErTEI$4=$ROkj*QNd_Qc<-4nEt+_ zKA+O%BDHcWmj%&6e;$d;=Ekos-z7*=eWo6D*08X(Lshk20$Z$G+NpiZ39WZ$(BiSS z*D1+qdee7H89A_V!-I~mvNHIMknUAQQ@AmB_==><K6(M6!RrAGjDu9w<ka6<{Wlvz zjF^1<d;R|Vv+@6zwHQH!H*UzdLOTh6+f>$N>^0?^CZXnkc99RUH?UEX_l)%S@ljDl zOQ$Y!P`fMrwg;ZI7~itlql<g0bAJxho$`kN(h~nKDGowp;ygQ#b_g(&*IJxKw|H;j zT~f(oZF)>uk9ppJa&AG1g>i2D#8g5)VdAZR?Qg{&Ly{76+~#s`^75@DLlvSmx+7Rf zN|1!NvPpj^j2J0*>Xjf}z>dqD1qNw?7FyUJ#XE}JV6dycs;uQc0n+%5|NHN91|1<q z;KgJ|X_aA1Y!Q8(mPk`d$Y>X`#!Y8mjg#*8847!gGFg%p=nUnn`*T!tRJpnM?81Lr zz%$CG<K|b48Y3wgnXg~`v)an&8{)0yDAy<j&40-YLAB3MO)_T8MmDWiDxDnpabE1v z$f$>C!9~S#bz&vFr?t&TNvP<4k-+UOUb~~lj@RiXIw;v>cG|P3K~aR-f3J76%zhAa zR6+y405Lfmgo;_c;NKwqn!6#`h8y17cAI|pmM(3T;9eW0#Yy)YXxZLV>2u^!DX0cF z@kbT?RT&EX5F81D|1Txrcq0_XnZ8K93;xatpXR(ruTYUeXO^*t_|y5cOEgA6^vw`6 zwJ@?pvZ*<<MeoqN!75<U-Vm)GA;kO_l75^+MmwN>RG7ySVS-A<^#`&*lF7lu1TM<6 zf!_wD&eaLup*{Hx_`WW3q#?n{cAPkG0ld%upSwmi8Kje+RVZrU*9w|)VC9s<B!4!N z|CU&(ip~jmFjcG|##Rrb#5-V}nuueB5cXy!2wiw10e;mC(qcj9AJq?hG+qJ9^msu; zOMdkFC;>wA?|q0yp-;JYn+25H8=+mKW@IW-7vj>)nU71^>-)aH>MLfeX9mXeVHT-b zwDc<Uk7IQ5%XWKAqXJX1d!oSyltt$qRB0Oj80a<ZHHWhE6rUaJ)_!7O6lw76c7oB$ zem$zEn2%(P6azn67JJr(DAD8PtF)6yeI`OusG(7M&Lh5(OH0d6nOz@RVURFU$!8Gd zI8|As%>n;b@B=}ZJ*s^*BY|64&wwE>G3|72!R1`GYM#OoMo**{rv~X5sapDqC=i>& z=FgZ}P>?hn4wqL*)P)*3YbVg3dBdlr<?s+I$t&6g(uqWBK2=`9Vn14F`fc>Gnb#MZ zx^qskEF{kLD~-d=JC1%g76k@Xt#~s@($-(axYkQ)9fu4Y-Q8qn_MLX&5h(7CK34}( z;siPupf>gxWFP&@!-_?*_+7X2R)LWR+`ya0`NPKcC#Xt%Z24EX`>&i#Q=#Fo(R{C8 z>k*P2A9jF)9_PSn&8>SU*I9=Cou}BwqBuv6T=rzG!v6VXQUGaDdJ?KYA5#vvirL4V za1a(@@!i^HEA05ZB3=S97}c>p8E>)>K!aeM%Ew)3)aInY=R8YnV44wX=;4LYR0?X4 zb&&tE>zL?cjTOoAfgUYnrol}@5V}w7N>mW*YO$F+64zZyQphGld}2y!d8&aXbce}# ziz-I{NYIM5ynZQP@V%wiuTWp9hVHJA_S%|!qHx#YxLkKj7U+pjc}y>3kG;wU8k^(g z?N)ZGc=5{x(jE?3NbPPzbybGsRitn}{se)%B4(AcDSiplySjgAXa8w5&+oRr5FAyV zv0o6>Bad}I`bqnhzjG%bWK?*zS>O%}Wb`NMi6^DSDPYm$$fe`jouGS4u-o#g2P)oF zkYsLA<Vt3PU7~W)SY*cOh=BCE4g@VilPU!z<wSb(yQ>5^MiEV@o^c6~7KlHOEE9fx zylqWsLFtRkl#<7BOHR*6YN452&Q0^N(%wMMroHLh!;D2|xJ4b)oJd#X%*hVY5TQP@ zb`h9el;YS$3{D@xkc#1`ZOFS(3)S&!QO6uE0yv}_5hOqNw#&qb63l0Fhk?TkqRVjO z!+!4*R06HMWCKY{Ha}$|%4e9;;+ebY5~6l1?G$iN!f@?Ne8tKV9{PHv7ZdeJUc_RT zd0aLZsCjdWgzgxEv&W!J$m?G~lDIx>#3NYz!Ki0hRapmf^cToXVzj|Q5z^A~GzHnK zMmH3ZpkQQIp#jd~cgtZVG>qwIT}3VIopRBoD3&qYyQmEBHTh9~FeXpxaGvjYa&^0| z6T>2JPSm$NUpA3xn69~UXJTJ+W#FJZb+0<Q_N??iWaV*QF%7T|rSp)$is`_S2D2aV zjUhs0Zr5WJ8KYM;ck-rF#`8fVqSFvM@j&D#!;Q*@+j3bgx4wFxLud}75hh-2wx@k7 zQ%nQ4f2@tIfc%oFBOF_<U{q6g`|}qom=IHhkhWp;F@dzIl~C@;;Bzgs`(d0baKHoQ z$jo6Y*zQ>%=M$pE%_rSWz9lkph5%&Ru(%zSUZ+u)^tsTZJ)P;%$Y)D^Y(NlIw8G>N zmSop2e16}s2xPHrL<rv*Lq^$TBmXs)zbGwIxo9m~diFzilEgRZ_$=>hA0KIjx`Pxw zA`2CbsSeqjc>x;b9mgHFzF&DHX@~!kj_L}Jtr~vh4Lr8#0|QbpH{i-?+ee!Zn(!z7 zSHQenUl(W}j6XG8XAm3j=2)y`*}+NkF!1A9k(8i-gH!8}YKfL}89&T!w?DbL=oC?S z=CEUBXMp=p&7*57ZGJE!3H>$XL7;~74busb%856j0*6&+$BgjmJVFfqK=dA$tH4Zw zFh+^~zzJLY?*Ry6kVrVdkxsrxTN{~DrYE-!e4-y<dH3I${ubrYih{ImK4}pC0;ML7 za7*8I7KGG(+t<T>d!RW8A}{_FlOT6yzNl@{Q%;8+s4RtBKZ|q-I#Kn)6O?XRIJ^m5 z8g?m$uZI0VlU`G6KXq~~xJ)@MBaJ&1A@nA&fNsoAW4&I=)07>*K+k=UrFUz)!YM88 z;PmF&5+Q+^m1QT_wTL^huH-)kx0t8Hj#NsAigq}9&_0EKkJZ-qj-tU)<)>JLYj9uz z9P6z!c^5a-b@3?CNa<Jag~e(+UDT*1Fsl`nVXHRDiSUK8V=8v;ob$gFT!6T3dY?f` zkA%Y`17#I&`%2<shkkYLe2F4Gx362h10q9(!VpXdJFEOl{cK*y>`-a1tuB>6x)0X; z7rr|-D*D(a0vH*%I*Es(1goE5D~nar72p<LXF^f0?9FOLY6hhaRmulFy&CBjEGqTU z59O2fQe|D;J@w8AQ(8Q@8EIMkba>KSup3Qx<C;}4DrODtxsB$Xd!)r#b+)dE-9pbe z%2M3zTrZ>MqVRd#$VnfaNneWD5;4hXE6?V$NP1Gta)~IHV&77Av3Vf&-b8^jb-E*> z9DmUX3v`MQy7m?_>+9K37dOsorCcxIbfdiv!6tLJrGjNq*K0-BJ~y`87)d91sj08A z-NNVcU!1wtx9fQIqWE4oHTxcJ4&Q!en_$2fKX6*oP9V4bvmGDr&DnNtKIkxofV>rB z?cGMP%cjQ0dKD7-_s+jBx+;>}wlta3d8O@ztNhywiJh#Y-(2OJtBv234-h3%7v&i2 z@MWq(^XdFmT|eB;HZci(U2<8AUy4vepVPe;Gry?>jXW7I>a--Ze6(D_Bc+=p5n*6- z<m4f>eSPEA&7?#cBS-C*b8k<Q4NgOlYhbeA_e1MdeFgkbWr9n=-BCsc5(C8(Xv7BV zpjLDC#sKkp-5886XM*W%+=AbvEOFX8?aJOuR>vUAf{hpa%5jIzQ?N{h`zR(EuhlD& ze72>i9o)Sb^9Z4v65cVSR`Fu~3)f-6(z6^nW*N=#=|>R(y6@gQQ(HW7ra6LoSpT>V zQz~lpDAq4a!#9i#nteI9^u5akE`sbQr6WP+KBt#ciM89BJAES2QRewZx_97!Jvn6b zTADXKlMppRlS9lq&q4Eh1S9{}PAZu>VG;*_z1cSIK5jRh6r_-y=QDm#4tyi!07;=- zFt?9NTF)T4NraBW{MMtMekZBsfR64r8z8o-D|#1|Uruj%M5qJ{v)_zCs<~{y^dHjs z^RZK)he<H`^KQ2g6!99sV?R8fQS>Bd&Saf-&oMeHm>-B{3bvmWR7px}(ad3AK)6}2 zY^AXrR!Ai7W{{DC(;nCAQER-VvrfC{dqW7aI7Y1dF?Xckh((&qv-YkK)}0HRIsU17 zN!PXlTxx)L{=*S9dnGjW7;(SQq*tFLZ?=t4{ln$%h8oA%7h1ojuvndHHUHA0h2vqd zmjne;+IdsF@o4qqZn|bieoRyy3Z~zW^j&S3Hnco9KpY_y>#->pxG!I0XYOl~Y<MnZ zI36~(Oajt0NpGUbE4aXk7BfS&HgojK+tV)j(U3Osn|a{e=UkG*bUc47uSuDHmn?%s z<i6FI!PM3#k;3fng7p<k9*i#bG^<{?{p*q-%5$}SHMJ8wyK$fli<9tn3_<mWkCjt9 zA*t@AH}{k5KJx>vR6DA*>b*I6lV1><OfXqSc+7488#p%&Y{9kY=AE<OzarrKpp27B z9h^U9PN^UP{J}usN6S&xJ&XK6vCJ3QN25aiMlq_ma>JcB`4!imfGIM_700nDHav*5 z0#||al6&-w2V%$4@#I4+*qo<P6)2(pz=BFY_IkNOV%_(n+7W$7;+5+fnS-pEIWS#v zhguGQYKXl6)XbRDg4UerMHW|m^NRPyQnxIM36}EkSRg@+DMQ%69RJ-by~RQGSYc@@ zAQJe#^ol%*B{Ow<vpLDNaRFKp+UW@w2ltre${<wb*vc!WOtD>~=D!SpMFmUgUOV1R z<6i7<31<Gm@`(DAA$lvMQEf@vCeT+X6tyiLusLOYV7B6AL4Dybq=cIPq`<JSiybLF zgC=rzkQCi{Q@_M4=C}~lbp|42lF1LW$1ac}B^DYU@u#aPuRV~8$VF8}q*<1OtcgL| zJL+dT3UfDUT@ze=c-NLMRjt8YCt-_>)_3pnVjn$4g_MMVc<8vf>cAcqxZAfaPR!qC zfgWJ=$P;{c{?=YP%^;;Faw&C#O3DjFTLua5u<<CXR0W;|#7cOZ`F+;%{+kO?iKyz} zn?1NCD*h%XbTTmWSOz(|v8Wbaf4$GDqs6dSr8b4S3VZCv2u9*2gH)-Fj!(U7WFXvB zx5ke-$82!)ey}P%YPjOR>mA+vO*Gp$5Negxp$R*3Px-MW_9}|i$1neQC-y|?%7pvH z+xhYHX1&X}3fH4e5#PZZZMV^jC0SSt$$RB{b#=EHY&fyFMzhtBm>)-mpcZFWu^zf3 zyYnsNS>#<o>{c1O?>ptPTalIz5}u;^w|+@|*NroJS8yx(n9iiP4BAP<jckFX2fA*W zAr%vW9f~5e%dooSbaeK}eCa?AY|C@LuJ2diq@$M(nvSjbe>IAGvaa}c|J!#P0j!*T zYk=V^j6DO3Ei!&xWW*BZHt;CvG@|FSi1&y;BJ=*t0~X<9Z{@so`k_FliJS4>Q8SzE z*&8ii`bF-Wjeh5fW(q{^Oeq$vy?epryVw>YP=n2X!FQlqE<mB=aN_Un*cp>BoH>RL z`=|8}o?bDeq3e0n7h!OYkN=LLTGNVXa^c&n+QF@#K5l9W^AkHlg)6wUK!N2`!byRp zQ`j*x{~n#cg=JahxO?F+(9$6}F)(4EvfEx?ogS(3{F?-p4qi8##;?JOrF2mf7)5zn z_MdBYi6^v)oq2QTP~&K;nA^AsRXWb25zbI8c7UJXC@QP-4Qd6$7pO}ekOC{Z6}4WM zDTQ=XIER#0MxVRwK>JehHN__Fnudnr6|ASv=abd9zW*~hd%aMYYfiHff_r1wP=}ox z%wV5f_J|OoMAQyD?oh0(OUTPX&2p>9w#R%+0{IHr;Aij276trpA}O@tgjZYUd7|S* zW<ANq%CPFEDAPG0ie|#PQZ`JthrSuv58n=4TFwM_^o%Hpy{3`tOU1iIR!y{Y;kJK8 zx=wxeMQaN6jG9Ta$=rwtD2L``dXFs36wiQ+aDXgh{&r%A-ho%%hC;E7_xCGUolK<0 zx{U`rXlq#jMb&<4oXknPuzsQl@6~1z>A2aUjb2c7QbHUY7O8dqbRb+X2kFpJTn4Ul z)%UPlt#jU0mAzhIiKxFc1Ho?<BBp~l$ned5sGlOji12tt{CTK>qI|@sC;voE5|X~^ z<ii^iEpjAC!c&Mn^riJnZz<E~6hwtO`0M-6U8l528Va^wqv4iuQm`G+Ir5C?+obc) z@Om3i=TLS0X5=M^36Vv)6s~XN13#M8XDOT<3jeXmNZ3r;MxhUdhY-Ym3=&viA?5$6 z@3L)P|E~Ru0rik6M@isyL%_R1nMhrwOIpeX%ZlLi;;-)Mm&!pf^t;Ee+gmP@b{!T* zdBj{Q$BkTeQFu|+@=QJNG-UA#T506zKq_}rL$Fp77NZL;+O<1Aqyt&cXmQ=9tJ%$0 z<L=zA=x_!_8D-er@P`eLTj8}kz3ejF$goQ-8T-dt)_xDCP^0)0IIVM|X?;zeGqB2f zE(0q!0zz#FN+Fq(cGW{o38rJi%bvWxCC~l%(5B=V=iG9#8(dCYO}E@=4{hCg?WIA= z_Hmc_h_y9i?|~}lfhfA?rSc`Sixq|2_friA&a`!>vY9Vay1r%i%ZQo#BGi&yhQm}Y z%jJfL2{zfH^QPOk&RVlaJ~J*nkj0mZ<-NPSuiAW?PZpEx_&&{UoLSO4Tw{jcqK-X? zBUfU@(EEdHfmc|Jy|Y0TO21m7mPT6c7P{XLoD>qo%4X(38NKIo53o&xRD?)+^*N`U zLe=@Oa=-8*DH4Ro#ixJKA9E|K(F7VT%$_)*{EXqz9q5j_8(#eE<%K$E&ry4kPgWd* zE;U%~(<bv0!(Cs2>LBP7TBRAXzIF1Pb!-U|-UUmfR(Wmj)y%=^_7#d)9BwW)2n9db z8~tL0$cM{D2x)%v4}p!~+;tqMw=!hc85T<SB8X1LQRK#AVhs*!B;POkv&-a4Fjr69 zM<~+j)f!5*(E2eyj1cN^ZOZG{$KL*O&3{~&1R?gKeapNdNM3PvapMXkf8lw^5B}5w zY)P@zbbx>D?kx8GahI3y0wS=K4%3ogxoEPAgSm)J1C_;I)UJfq%e~*WUOVucrz1r2 z`kM)VA){(3TAs6#J+c%uA>P$(_7$OP(o-xrdOtJg=QPdDLxQM!oP&#d6#KE?Y3uRq zggG(7DQ$f%a;tcP`4q-U{bwPK!cqK&<&-L}elAoK1wF0~;03)H8^&e1?pMCUvg5wT zHoM$=AFIq>-@4GZh&73OMs;?-A%F;%^!HpBmoSZc2LZz@+xPBJL-LKR=CsC=1@b6} z$V>)gd1w)aw`p$Gd~6?}pHP{EwpxUvD*dQ$tt&pl>|6AnMpRW@6`C>Hj(&d;dl=+; z8mPto#?eg=UV8kpC~$i^x)de6%VOxf=1F0rq9*{!UtQ_bUzg3CJQo#;8cDQ`+pnp; zkCi`@G0f7Cvz+O-_ni-`T$-?E-}nmGuHUSJzH-#k5fsD@e}8U>%dN>-p_!SQ6zNYs z_U2m0JNpfP7ot`5Va5IRHYn^7Vz|ZMaT?K|;C{3f{Q6FW<&O57#bmA}gk))kMupsW zL3ms6ToN(u#?SdIEca~FyQl!)JB9rLWtM5UfP;}I{IG~)#wNDB-&g?CHw)9O8j2a9 zWB%9g4sXw&45W+xQ5?ODGE?#lxc1rXlQFs?r?+UL!NXNu#Q5V>lF$f(ZUrjC|L(C8 zG3zz_FaoslUQ<7oH{F<x(UpT|_ZOLO5RL55QEpr?L0XNHif3HP*CP@Ug%piNS;GAr zBHNF+tv-_pe675S-`U-$45_3(t4qhFDp8jH@XbpSLhKSUpVB@-(Feg!+&tPig$f#r z*<k5;JLNb_f5$sVGGR+#9l-+vO&x4!2n?%+{!nKB=%xrHR`cdQ>-74SGDP|CvIivU zG)p#a&knm6S&D-D*Z-AqUA1%eZ6NapG%&<*e(bNRY+e*;YD7!^r5px}MXCV@-H5#| z6ZMR8@IHD45u;-4m?kh-JWb?N{*3=+!BGR7ZefS^P-|Y44K=7Xl@y9D)fQ<~s7s@1 z)>`ONxrSRqDR!gK-~x(&-on;X8gZVvXQaC(Q9ALGZ!mhl>b+=Y6V&M<7u4di5kH;X zx0wcBgJJ4Q-n;!eD)Dv<nw@RD5Mgu0?cFB!nLN?^rnxzKyx_qt(Toc7jrfkX4_4C0 zqy{6eJMi_#=*se=rd{oq<{GR;|2)f`T>HdlksAY!-uWWXtmG;?z4RlY1lH^3Qq%q$ zaG|l}V`eRzh9w=v#<iAKOwYMlU|x0`d9vV{&EV_4uH|Xf(G8Y?kynjMY3LdIuCG?2 zuAj8fwjWZ*UTW<4Qay3|%2k6{G5Nj^eQ-`CIoBsr=*DXlR%u+NWJ!0nOlG2`-O$#( zO%r|NEWNYqlVLBmAKq?jaYGyD$J8%FomDAjTt^+Kn!(v|I!4f;+5JgA7fSvR?qFUJ ze6WI>X7<ZVWv#f_WF(Bd>`|XrJxkDJ*Rj;H!W7`oau>NVT!iW9n<2?i?<`z{2QdFF z#%y&$wSSW_w@LlrebrD<L?TKS%_o0X=fJH2rz;8l3@OKjS)@087oR=dm-m8|3_&a! z^Ngigb!4wN;ILDs&>Re?gi~50D8-h%9H3~z9CUrP8z}tFSs(Cm$lMA3dgIdh#ur+5 za6cMGS5<M|S<PQY4rp(2ThHuw?>Uw$ebRv4#{78e4{I>a9wqfBM+%XXXeu>AGK-SC zQUiv(&J)cRvKQc0UO*%-b0^jM<{SnFV1~<M2VzqX=IwSl+q7W8-)_eSVAJz%F(ThZ zmSr+sFKHe<ctt*>ge|6k)EdsD{cx_)*;L#{upY!9BSC4uD7l1G{h_~@qRka-NBV?Z zt|F18z0EU857UN`{_dQpWK>!;U+mSqN5aAwONl+_#ADu1;%3N=FQbFazLDC;Ti4Si zyhrZ$OLqN;dWx!s=PFF>Y-`Ojr+hjjbzpkPQ*}cOya}a06?A;<8~%tXD{*OW9m~w< zog8L#WEkO4gV-Z|G<7F=rjtodVsc(V+YcfobR~Y9`Uy?Ka;@Efn-^8;%ShKReJy?_ z>P%GkRx&m8JZR*9il&9YG9Sd(kvw(ZWO789;JgmQu+*woe#}=_p(ni<O@wB+p7R+a z(lRCJjXpMGD{i*A`LT}zE0r)oCd+wzo#zjZk`0y$MzyBE9TF!V<OdSEr556FoZ1Ys z(=909anTebk&8`v<(`T&cY}%M4ia1kS=N*B-Z!-ng$=Fn`ts`W#QOI&HuxQ`ISn2! z&eJ-l5j$qfUU9lbhxN4F$V~C5+=q9cp)@i4u;&TM5BKn)tJe_C(BvM!<lEslYBo0K zTA7|tgbWs2YD-mh!mY~6m1YL8%XC<QM?9P1A!-YhVO#uJheD`36NKWa3T&=;UM216 z!z9cb5i+scx-h160?@fh{pMl=tK#hFj#!B#$P0z*1Dd9UG;zx<q8mmBx61<NOXVlf za|U=%xkH_Pocv6Zu`$tNE-nalDy4`v5aRH))Ii<7m>#nD*J?}}%B&@Rx0Fez-Yr|) zdeyVsyt+r~+8eyx<sk+U!i(vvr5oqHs3Cj(M(}g&Di|sw`Ijar>rpt!6$y}&KYfnZ z999v1=MEdWA`cd;QvwJrm%vcRDH+IO=+75cWV%dS6a&rg9f^_BmPj}DAD$B<Ua>m& zOkLRO1$2-r4Rs{GT5TQ7X7b-7@!b}dB^pu>Ke@$_=h{dSLl(Uufo8(S?$29LM<r_D zBA$<`eMj0btVn`5urIprVsSJ<e30AACey)zy$i2c_z{Z{yxf4e07+Jj$f)&Rgs#`? zoR7e(?x;AnI>TerXzfIVlx6|Xojn+*7KY>J>cmL9Z;0}3qiWon_r+bx=&xD|#Rw`l zg~bw+RTm)mfZ2Zyqd+%S<=x?U2H)0~B3yIp6g&Jup1ogte!~UG;?FhDJz4Rci7E_k zoHG_EW{#e%Cy{-}rb-Q==tfs~xS*4LSqF~Y=E3_IQ!*-g1OC(?Q~qJc>d4cihER!M zROMCV8nKVL+{<T=csDvpkp5yc+70jJuk)r4X5q-aC?Sj2R5P;@orQ9D8uM_+cW<-H z;4iw6{;EVYT{Ud2dno%1;=p!dr@KY|V^vk)NrzEg#f_Uwy$X}-S7t}Z5j{2$mi7(7 z)2#6DcV8qSd-t}b@9)qZ?fnF=qH%`kp0`C7`a77}LrCBxzBMA!(g!E(Ez(aSHNsxc zZJ6R^YkA=s1gr=zwmn*%w52IUM%HDGt)7~(p@zZA%(jKVw)50kq<jiVBganSd+`4L zaKlqq&x<}#ty3{jGhpqS0u@3t%Zk}L<K0`+S$)~Tj!V~Jx_Msm8MWo-K((1^qu*&= zCd1@;&Ja{HLMZrQslmg_3B=}e%IWwe^U+Pe<*Kl+K<rjjKwj|+d&|8%ZW2VejZ^=l zL3GD%jr+0k!o@kt>GZDZ`G9*QJ-nbMPcyxEalF4_@;X`%#eL9haBC;`-S*60MiLw# zv$~Bo%(HOpy<YgKMqBP23O6D4U^kRtjEVE7E_tDkqN@^=T#k5o;+@<K{sRCiw9iQI zII<}zZ6u#x*C9&-r$NQ1uP^iu94c~Wzb($Xu2H-^c=u{IGt}sN0H+boC<tM0yf47V z%K^VSQ~aGG(f16%H>5bqT-fIK6&Wlqt593KHj2WOsgUmJ*B46EeXcTSW00qoJy?ur z7&-d7bujt#xd3=qO<~aSiqm#i!Qy8gb(~b3bc&cs3H?mYzSwSgTJ88c?D-P3qq7lR z(~mF6$7!P{v}bz0x$Qm;LS^yLgv$Apb~9dwMR%t*=)5KCaZ(NgjVS#yk5=}*gk`3@ z8^^%`3rE!Q+)plHl`|E9Z@YI)EJlBkG6+U7v1709O_(oA;&>uh!Wf-HKX+kU<8u)v zI;-<TugothTtaxfLp7UxY`LEq2<-pk2Ku6-nc7ixI9+AHrfoJ#>f*o$c`jk|WQ90@ z*)fMz_`YPiC%=xnf+{*oDStk2!0~5sJ;oOmmfyTNB`HlULgydpH9Tm+M+I4YfNYW? z?;-;p#(_6*wN?7lyO?mW9wD?x4vUWu<=LlSEM+@`*BZKMp9dWgFUmCyAxdbYq2uH5 zKt-I@Z5k87-n*bwrgR_a<sc!aTk?d~WFhzILjt6<7FuSJnZ`aBh;D84ddV~Q@empz zbQZ_#zXM(S>4jEI=&}bUs)<GfO3NvVYUi`VDZ2)lLVw)o0@ZA9<eG0pevaGvDTIdr zl6wc9u{)JL$ExMyU_gRAvp-awnVr%57!|aMcuJd?FiC0;G*3c(uc71CB5%G!OpGA@ z(zdj3;YS2?p|db=eDHYUD?pQJqU~ys(!8KhN%1~lBYLfToV+6p74!wcD7B$*^ymmj zLpe>G!vhagZobR-T6!s)4NkDyo!VlCDcTuShAG`bdNJm{QZ47oF_DMV5*9`tN3HnI zqtN=4@)1wNDc#5Vg|~)?5u_iYnk`vllO3tj6it}rDAThdT_?W<Og*mZdP7+UY8T+W z7e08TJMnC4=%p5AdUoV7E(FD_4DOkx;OVRc70O?})aK=J9cK2d1ZB1U1=2EEGu$j^ z%LYj%M35>+9baUpTXQt$iYPM-<YSk{WNvFF#2sFHYY03#SY&Q^WB>{o8{6|2(vxLy zjfKy6u|I?&6b~X$x`ZrlmTdn~ee$vWnpcfisV?h>$0y%pkhKw2gd<WyhO4Txdn3oS zld`%?x}CP53=4WXBE2ScNG2vbwI0$<FNPdf7^DE#44#iQ&|y+8!P^z>_o5n5%uR?M z%tb`cX(VAGWW$Lu-b!4?a7uOKM?xbryKtA?FIgnW+VPIuIgOuSy{V1Sxx4Nx$7f4l zG}RHDXlM}893QK(7LBcJw+Av|xQpNF_wqWIPH+%PMcRAcMheVFuB6xq!{?`)ekSS} zSTG8mjSwP*(B(|6G#Bn!ho_H!pUdl__C=OgxYQN@tpdmTQz7+L)y3J60RQcLzF+eP zkeU>lrrihk;KYm5MjtPFAUiPT$`0%Ax|W1W?a$zmZfsTHKEo`fdk=sSKm4Z`G0>f! zrpfwo+0O&4x}XVwe`ZNdOCU|~IF?v$&uaF95^?}!u{HomzZ1}Bf`F}Idw)GyN{^9g z&PC{R%-;{^3I)^R(>qz!)gJH}{Z|4eT}HRFh7=QDciWr*RNp5-{oyAMWv)S6IXOIk zg@%HMvM!Tzd<Gg&JxWu)t%?{@^3i)jhFJ!uhLaKvj;NDqqD7@{bq7zRmq#}@sEev+ z(dI}J3^&^7Cq6-xGwoe5%jPbTfb+ar5_+H4z?Pl1&JiQLxF9Y$?b7L7zrh;3fufvN zqXm~nK6bGZhX_}W64W3uS-Q$%T>2G#LOhH#NGnR=`Ll<O1OSp5aN`7%X{!gDa~zi) z2XCZ*Bl|>ou$L{+6l4GX1s=Ju-)3&~PaAId%o|!}F6OVg?7vZ}oe3z<gZc{{7~4O< zRe;2_V!CH>$0Nj>>^HzL=zN6bfl2^5)_IHX4}er#7=Bdya$zyH*5WexpO##7O{=ng z{t91?PY{kLP*ZI|tk9&LRr#C#$G@SbuM1UK?Wz1r{fYjOCx~79yzNCtbM4+AgzuLj zL>1{A?z4Z-KYFD4*_+#Z+w1be4Y(#Bc>bfGU!DD_Khpl!@YeqV?{f;_c5UwpcE*Cw zjps4I^5B8J0+dHn!*pMk=b$rsc(@oqywIDnVGGfxHFS-9R8C+-c%14J)Z%*cV<AP& zd9ss609f{%5K)fp>F+=8#Tm`$z2<^HANZ<&iRd-}U}@-C=o_UtCJ^aXnH!U7bSLqu z0vL#Gp@vbl!7yn$_ttjPI;UGpwzv%e!MLvBaU&5yr@>296B)Uq2?0S$P~`EHN58zs z@cRVqQv~6BaoT&w#Km+3k7GVS$_K_^+e&LB1Z;*^k%8!J)}x>0{j5V}1O%dlyY+uV zy5F;7&;07}o$_eo*Z*vy6o$WI#mL?5<*1t-0&5DLMjz#y4hBEUb8fsJ=CojM;z?tj zQ5UX|QaUu{0mld&BK&Jl#Z{?|$-2>hTWifyy;p2M(g7_i5K&dM32FDF+GCDk$+78? zD5G0Xv(p}u%yC|Y)-@KNbqYh+*BE%@BRlfwj0U=@*2u<X(>O*4gH<9v8~2M&n&L&Y ze7%!HZu33Sf;cvNju<`O?DyW3BPkhkX~^xVgp)~}=jDEv{0;F*mv=u$PEA%12&OxW zOMY++Kfi3B<h|=$@AP_~sQvVpkDsPgGCBc<Ucs|y*Smj*`7O?Erl2ELn~0*>rah@| zw07LOc_|BGJ1n+5$CtdnD@C!|FnE?qS`_1;jr6i#UW2~;sowu$W@8#6_4Mko_lPl) z(6H1rSiv^-?HXJ5beIrqjpy#EOR=$URH;zZRi8rTQ5@wj-Sqk-3xAC*vG*!zT7!2w z(@^5a`&YKLc4`_;UEd>N6ET9!ONZeUO%iZ0`U)6bXHB*GZR@YWt2P&_BN@vI){h|J z7e1DYIo^J_5@G3;>S0W_7>g!+T|0U8)q)ot>HD+3ek}gI26egtqa-_IKTp5ul^=>M z0n0mV(Hw02XcEF|HMTIVa8F|f(zwK~3WoY<J}=t2=5LpEfy>=?c~BHdVLn@Tw<9j4 z`*P@Umakxo=0Wg{7J99vXtAZ}xmfyYW{OR<k<+N^EWx8YKbrZMrn=PbayV~M%M9Ko z*x@3tn0Z~3PIE8Kw}8qt_pSEca5e)`i{?^lkL89(dCI_j(MX+IwX);9aO&q!-3V9y zWF3F`^sUC<qfgFy47?IDvFD4XZV%6!BPM^Ty#gn*sLAfP_}EMM=10Q%bOYmP@0jFw zw36=2HV3hS6-gP}RnPd?cC}k;5wh4@q-+mAGY-LnC$Wt#A>#`y%MfE;Z5Vt>jQi1s zLjG<@x20z6t*WkP#l&*XaprT{6r7ifGxS-#>R{b)^hSz($hf=hSuU&tmb{K7rTc}n zr0uHi9S;t%M90^<%`YPNdL>f};JiJF$coPp#T3hzx|6n7_2;}*?OQy$+1^|$dsLBC z$e)-OXL7zlrf_V~`O4t%-266aGz87>)u*s+r6T(BEWv}#H*21_<?P^5GJ@+|)`-!3 zXF`H&Io1f=!=VYa7L(Y4>JOTUWw%RVg6kT&q=)>uMTT`jwFY?;AsR-p2E&>OBgBq9 z8tg?`w>b^#26TvWDWHaV7ZaVM?Ajo>g^a*MOqtsHWJmh#zup-}6L>jV+qS0P7MH|0 z{IGWa8|(U<46%_GuuxH3*qDWfpPyjqucI}|Sw}i$LFy^{+rOo0xGntkPXlqIr100` z7XZ{VNHnp)29Nbjd|e4ZBts2C9w{26iFj^B;w`GY&0%YsZjFK3ErT5DL`)Arhj7|K zC<sSTr%!;>VwHga*ok+n|3{9;e~V*>T!v5gJYCQI!>@QAG4qtNk_vGn;e@sU08TyP z{avsK7X1^XR16QK9I2W~hU5TH3E-M9bmbkyhrqOz(GdXs2K+aqYP4zU9seI915tow z{0S4{zX%n7An*T>QvTZspMWR-BWDDU-umCrG%9GCzBRR}?>YS$8_+MdGoFAd@4(SG z`)^i<=w)o;r?1e_o_{r)M!yn(H8B~W(>Xd~%3s&bahZf~P~u^l_|OGG?}ja=Psb`p zS&m*`Y<dpr*Iu}+n}0$Z1RcB9^&aEIFeMy53k$UAjWVzIDH5?TVGQbCzFIm(%>_kV z)KJ1xCAt9ZLf}c*XS*T+gUF&ztR1vHFXq%ORxU2nvKNJKl~2eP#gIb5;?p%H%ZAa- zK0h##s-$Y=X?NuKc}bM4tzR0uZbO5j5)5dlAxRp3B@9gg)gZSkFuD#n97c+5c5R5y zWJ%j+|DAp<jiAaEsxM1fmm*c3oX{@_$6)rh-%g3CYb4882cB=vRslSs^g(k7ZZl7$ zmTs`_F1NbTZ`uVQll*35{153q|J$wm|6thsU$(2%d{iiDZPAaQ-xLA83;y!#wi(Jo zzX?3&X^_NDh+6%B_#3}DK@Ld>w=a-DCH8+P2m^iZ5UP{Ve`p(D7bX5*q^JMGDhxr~ z&Z>B77u*vfae^Yc6N82jMnR6G;|Y*ciK0-tVhGgksIuse;kOkNljqH`$&g$;wv|sm z(!k7^(9g6h!V*d!RJn~eFrxxUBF_o+3`+fmcaFlze?c^nrU=401Tyo4+rydwukN3A zp2|=^6$ttsaw6zOtT*NrN52)q2m_$y92~zhbAq(>^Gn)$su24q4I#f@SvGhxev9g0 zK*SUMi_KDrP;`dtEYtk63Y}3Ici}Vt=2pa@$|(XH5#ir_HW(8BA;MflwA{aTZ0NrF zlZ5Fi$KOc*bi|%W?L>L*HBlWrKLWhNIze!+4f{?X=fv`2DOT9)m4DL)EVH7Iaxp&) zodRV|$&fSvmV`&8s-WeiJV%@QsxuUNAXcaT{H~0y*VHO+smySzl2a_AJ>gg286BZ0 zs}*Bh!x7_Q=y2eOm4g}N&k)8Du6|0m%%jA@2~znO-al*vniJB*=$7(YuZ!@kla|s6 zRpl}2)heYt$ZLPG7>G|Uq$i=`f9d|nA0jDYnie=(M%jOJBt{|v*x7BX`d5BS2ONpn z5>`Ur^8Y}R(EXyyx+VtB=O2~?#S%hzFXb1%<`-u7lhrxtL4edg0r~?v04>sC*k(3; zozg<#_j*&D{w)E3egCFl{p-+&Q2S5D*6Nhxow!H<xGdZ8TahGXuV=1%3YAVkPg5XW zY>m(_9A)gSXKu*^TK}vKLk-oxvR8oAvhP*2J788wK?_ih80X}7Us@<^he7Smc6Wco z6BV;ALJ{E3$~g>gemJILr+Uh_wZO}tNYpgpk10=eEEmw}k@8h%q%m-P3^7i~1N@_F zSr-9!E$W0LsSsVH{%;C~Q9|hJx8iun<v9?S6JH|`pYtaWPY3<&B2VvvXIKlocxu36 zqDrjpS+*<<P<hKnYgZ-%$W2^BYO0YM2?j+2%J?gBVt*1;v>G*ByZ)NrooWOW_QqhN zJ}oz9F+by6Z>--P)?InvO9(Erwy?n_>uv6Yb^MKXY{m7A?=+s<w&xqm-0-rGSymlE z@;loIMHK(OGyTsE1<;N@$eV4zKgkQQ*R%C@6yfq#s~>K7pBw^BVmoDY14iZB&A1g% z0lRrrp<%{I=r?h4#4KcZ&LHS6#y9mPgPxd<RGX3VU|Iau*TeLPAhXty>kM~6i3NK% zR%}J;{qH_hY3mj-4`n(qe4yZQe9OvlSoX!+i6u}C2Of?&;T+P*DivyOu>n{bqyv5> z=}%i!o5GY&I2j-MoleL`zs5K@y@09#s!G49AO97>@_*?5UkT7Q-O+0>b7ZRrg)dN+ zGZnzjhdgVgv|zzgVZrC=YdeLfPACyw;Y2M(&j50vW(M~0QK(U`KH3048*LF3qqJW% z0WD%00b%`iZr`X^6RirTiEGPiki%0ZvI6kmzFW+#S8f-k#-qrBV`#8)5kk7Qh!w9b zW?@z$q+rZJ+r~59;a3nV8}EPs!4vL;R1Ien^#=ex=3>V|Z-Q@uN_&3J=K^k+5cp2h zj|<7{)(>7j5^W`;M)2d=yXo<zh{2!~k0Q(sFf9mT)NRTVyt{1o6v6xzfVXNvbuCjS znBa+=*LjPKcyCHBXInWjS}d+ZZ$eE0wRc5rB%H(mt2eL8r*8+?n*%+LImRulvxP5; z#I3cj#MYL6;Qx&UICYDWD~}H!R$hRFRmI+P7{rtRl(4+eZ|K#GqcFai2x>@JFR=+q zzMD%bP})KW6jE8jl}6}i#STTB&@dFVv?|Tn^<5l7HH;$cpqLq%QU1g>zCyrSK4oT% z5k-WXtx4r;@A(&ZJMeqF$|!yYxJpi=NWqa_!3=|&o+!6^leb*9Av?LNvoH!(kBdi~ zE`G-|4EY>11_0t1HWR|MMOysWia$r@EhMnYqosG9wtt0Hes=5P3E^P0<j`EM?%9Le z1%4&hqMF?bFO9H9W1rtJGGI0U!CNt}L9z}T?GuN`JAzt`kLF~(Z@Pq9i$@n}E?{M- z*Okl+EMpXRb26dtXRR<<6aBQ|2`JpK;<BVkv@|way+aevQK0pYkrql3AFZg?rw_PM zd!e!9(c-gN$NBot@QL@!UR<#sH&E@nnz6-g9N6eSYXDGoutCmNwc9<%_?>4gAOrQS z&C2>I0XnMmO00$I9Ggt~aQz(dI<0K=)7jpApi5q)YG6B4(rCrvfv{voAdY6(=pj{8 zy6Jo9>c08TzQ_cvJ`-fG(uMnI;ds1H&<mI$Rdll^V?wf!DYs>T|Au1D3UPIJ;dd^q zc<JLgy2(fKG>Q}3s}CDux_(Al{`3h{4)e7F2r6-m7g}3NTi|{-egNwk&wi=DK1Zjc zJU;iPXlDRM2E(hN4}f!IctI9fdPJ282Zg9#8BQexw!$60J7WW=5~6I_#nGvERb2*c zdg!om6NTl2vWHaDx3dq-4FSv>5MTnc1=FaC)A*q0Ra?~fzw)%upM3A%1XLJ2w<Gyy zU}+Fa(m<$q!wGW`I~(&ZQs>|j!k0WC97O`Yb1#wFYZ~T8WunxZ0u;B>sAQ2#SKXZ{ zT0B7<%Uc38(%3;ONz<aAkvdN+$U}B06_LTeSR<HGR*xAe?IgZ1gqtGn<%%ERuYf*8 zisd4J)ph-U%w7MT#`U#QHz5?IAF4wKY6Ss@?1cC3u?d;AP-ypzWVSlthX5Butci}= z4LR(Vk-{7q{i0d`PZZ3F{Rr8F6#ZgQT!;0wUi5#VjZV<~9bG-1L1qJL1L(D;w@HPL z<gNuN!w65(8uB>QNUwqpyRA<BKb<f5=kQ`3GUxN}1B}uy>bmqLL^=KN2kN5&N~Lig z?7Iw2o?d3D@r~~qr%vd6;`T_Xt)>V00VxVRjH<$^wIK_!xjS?xv<XmJT@b*EPm+7h zYDUX+sv!%xu=rbyBR+nELf`#9gz(t|s(PiGsOI~vC{L;wZaYqP9zJ}nxmiO`J<EKj zcXT(vkn}~2freeHVp2uIG61y}T<6trHx)$r9F1~v5^8APDqU6`fqc01IMLcBDte$+ z{kgi(@?-tI@5Sdkd{=OSHyMqXx++Dx?De5bs#ed_-i-Y6c_*x71I{9&vVGMF_gl9u z_~>&Jx^pTYYY(qriu=v;wgHJ<w9c+tqj*;Rf*?OXkG1y8#mJ$s$U+S~Lpeu_bNQY! z*V5CpSlQO{z@%p@1H`J8B?3Pyi!_C?$0EHyWU;jV)`yH~>!-VVTZQ1K=xhelcTy~$ zxukL)UR=kpBPt&b`eR>bd7!2$;;$4vt41gWysfy#R)5`blj+j2MN@zwyUZeV=dYw` z#bH(}elq&RSOovGQiP)YI{y*?gBk?qN1MKX$|keuSjrMUD^qtD#k@`gj#)~pM|Fij zKnl3>r?lRkwkcrUttXwhNNf@k4N&8N3DJS}AkQ2AwkW{}WTYFv*E9iP>s=Qwmhc+_ z*qmr^{(91gd<`Q&wzErkH!sh?&4U$AA+?uTq}zG1=m!sh=RybPH*!0y$v99VJe1F8 zB!--SR2_uk+8gK0^bx7!vb6w$;#%Y1#k{?}==w0f=KRDvkWPG*_D;GW`lV_Ut#p_J zXhN<7Vd?=Lce9f;KY_E)6ud}Y`a>GfBQCWRNCQ$gXpT1Mb+0ZMv|{?|M}}ZTFF5Tq zxjz9)c6_D(>X|y=!44;33eGnc9=>K0Eqc_J)i=sk-)BFyW5US+h*EK_JN3pMi?MR? ztQ$zT;6qjukTF~%))g)_?mO^V^fxMe`a|lKK6p~3VK*82Qa!a(<#fdmOyyPZVmf2S zoB6U!@R_X;Emmdt`a?V;#ADM)hZU2xwiLdmGa9zFB5jM{=L2YXcSQ|cw`!p7-~5Nd zqu{zr42@s`pF3`rh_6Tu$a2phMj6A&BB=};G)%w|jknnps7?yr6F3vmr!L{$C%)oU z>ZQ}Dnm>)HWNMeo{ALlK^Q&2-t5A5W!Y^5e9|d0bCq@IZ$mN6n#dtiE(81t*hK7ey z!ajJ}RV!*08!gD3{Wb|y%Iw|og`PAZ%Mgy^ZgJho0ttSfy0kYJi_Wp!<!@pfBM6fj z*`Ja%976W>Rz3sZ;I~T5U3zzqiMt8`ZZDaNIqypFk7W2nD1~PFnDEFIWD<ktGqmmF z6)pU2R{RU`+a%`A!kMEwN!5`P7nP5SZb0KS4wQdNrJ*K)6~wa2(QvqjrQa8N%wY)b z5@B~YOO{@JbUh2#>@K(!)vRBxIc;FQzsbE<D8wHj6y4M(HF-S5Wk>{?`fm=6FMIjn z+a2dvCYK8OMAlzw0n9&B23c=%bakMQ20o{(Y26u7t-)5SU!XEXDk4|5-}*mu$)>c( zXeY}^6+H4j3S5;R9iRa%{Whsf>CXJL|1mZ+9hZwjcNsIYRQl?u>4UJBP1J8ZIZ<F| zIRT#0?*!ZDxg5@<Nb+R>eCPMgscWPJh@OFhR<SYTV+0!PhBa9RgV_;0?`c#qD)1t7 zK=-dP03nqF<cGy^B1i1x$O-+W+Xds4tiYFtXkm|nX}6g~#{>aj`4Exrw`DhUILGJh zr^B^VgR@kC*jNPEzuFV`V`m>uv=9SYgC07l;*3knGLyV#1x5xySR6K)uBAU)+LtYj z#xjk*x%=S^yW?4;w625QfuqxhdhGZ#y0@Y{Gy>qZBHI{wU3AT-Tu*S<$+^MVyEjb2 z!TFC_wgS!T#@NW{Y59>}PN6C0txjk`N$Xz?b6PPbkM6y@8qs&_E8m^YH@6}c^o6x} zg^`MiqO4|SCvfL^o-`CgzTt;9C0MklAOL{yW5XQWoN%+HbPGn7h+;<g1%%64!74ym zAVKsXI~3F?Mse$5(ZMLLhE<HL`oXA(?HxHt;8Row_0Vj#{9$InVDH;+O|_1@6(ZzF zyK{=(%d`Gw{d!0LPito$4^{j3@zf)g>|3%|2%!iel-rPOcZ6&sLW{BsGt6j33MrKR zo-mf|Wf{gOx5z$OMz+b8L5yX_GMbs^8oIlG&+GZ;`Q!O-UdK7-I_GlEnd>^=_xm#< zc3V=Nw+Xb_iD478p7p0Aizum3&o_*v@9Pe1RyD;`&orVG?J`a3g!zOFO}fc@!s~e@ zz~(Dl|0AtNxWWx3AXEIL?RxAN-<RgL&L@VokN%G&ijQT=)CW~k$5j_T&W9Msc$D<C zx7~B>cOdQ<6Q2>14~~TNQ5|nG7M9Z~!E08swi}+)(bg-!UbaFW6G)$YZ^rtHB)Q6_ zDSr2+@WJp6`f&tu=xlgr#Ipw9Ipq6yVaoezJC!i2&U7z~F%RT1gh+b3byLa=4S4N0 z>Fo8^&{!`O*gd->O&Z~qp2k*m8h4itxK<2#GUwFAAKB*zbwHRMmvG4kYEx@Ts}#(L zLC2tHIn}DlF`{~=X_1Ewbe(Q!p7!I}WhF4kb)>A=t<fTIMQxtuX|*J>V{>h+qgJNJ z*;zn;SL8^XHS<LjCTpI(srvJdlfVB3V@G_EWW)5)i#+L`Ma2D9vVQ&-9`j<|4$07x z#)oe8k_TG=qqltp*pwxZsxL0wP}<@VF5u3*@}xwnoAG(<&29FxvOM<pga(mNB{OJo zKe|+tYEwtj5Ag;C`W&&dpQhg-eO0%d1i`c}#>bg&af4YYk}wg_t(lH{BSpHXb4My< zZFvpGwWgf{Y?(x#rIGU^9pFYzxch9eAPbS@0B9e{-SCVesCNYTssa*))$Hb4>il)* z*_1?}#lJ)0Zk1*sRD1XYQNv5ojDjtm$zM7^1Bvr0B;{xZM2eCB7x)(VuUKeSAMMnv zJ0Tyt&PBhA1MnOlHp7Y1iP%2Jh`|9cC;hJH^H-g=2ugi7BrCVLPD~Dd;D%%c7hc;3 zAHry~EnnTi`E~PkvzbH-rDY;>S`@AI3V>owCX@iW@Sabv`;-xL3fW=9fxrTZWqiFx zK0&L+I?^omgqLCgK(>4^ljq%XZS`W2>zS6BtPq!Tt*)j~;gx(w*FzyV7Xx%fIL2)0 zdHX)?xCB!#vFNUG6dBVQk*YsyVQ)~=g=1MFm<;ja#s_C*<vB?VwKjF|Zq0<!346b| z+82MiZ?UlSdh=tFxE!pu+G8=$W?}pEEXJ^IUWs48wZmLg-t3J#9^F+I1EqdzO?uft z*1D>3>&GSF`o%QdiFU<10`MN#8(N`eN<SLwN`92;1yPdE68V6(poyo}+2{pRsWbc& zWBnzkI>18hJQ8t(8GNu_L;^gt4rcyk=PysheHVb2^BkQ@f34y#t@B58m&)s~4@cF_ zeEfu#R<y*#S*<-`MB_+2>^qIsv?wknuwB<#3r|!Yd&CVM@)|ltgpTw{acbS$FO?qP zxg$uIHRxn6-SlI!@)jPnW#;sM{rl{_)(=Q5*qxg|goccE!+h?xAMpxB(C*qRZQ$_L zc8mF1aJvt*=g&sT@*y6%=<;Ta56$a$q7nL>qNu=wc>9cL#jkm8Q>5C0?WdUc>8D4O zO`sa@(n4O1jkA0WXG*bz>5Lf1$OeZM-IA8I?QF~sy=V{DdT2d*7Z4;L_jYq4av8a6 zz6Gv8+3>^H(YWmlEq#{&*qCnQ<l8x@kIMACw|yiE#j|9s3vjA8^DRreIElxi*-7T! z#2F<V_{W%r7Tq1Ytf`@2#`<ULjDE*y&WQ+ugGeT=Z;AMYon!eQ^+m5%4SE>5j#+b3 zgFqNr7N)Q2yAbml<^v6d_X6j|xME{=t!a(0Ngz*T|HZO`;_HLO%E{&S^VIooIDXbf z<-!p-=wMpIr20kVneX&Q(D()Y(xgJ$>!71g<w^r_!iLrvQCT=!^2{kBE<z%^yEKe4 z)F{Mdl=5kYT+mD&?&5q3zGZT@KAPM|t^qoR9Y1SXY?wioj?-08vqO^sL+h?WgJ$v3 zu4(dxdhi38C5xgUeRL`PV6>$am*SA96+!~V5+mM|G^p?~-7hw`=s%wJQe#XWOh{uO zYR{1is(S2w$WuBtSr_&`JRd|0h@5FD*}v@`I0Vt}UTuMUV}DGKTgg-R<K}%21NGR~ z(HpWNdE}3ai;}g*IBCSU&-Tm0*tJ^Gx*K#M$9s(N;?X5P79#5mLBc6-L$IC$Ja(^c zcQEVWbt1bT$mkS@`gP1P8s>BEN{MnwdiVOm|MaVekPC!n#Zu9PLM$5g0UjyA+QUyH zs}5L|qF;N=LrjrlJ;asSbL%3tNfVp$uvc?9!bj_y%m7L(r5B!1(ks<4oj+=(I$Ejk z0l@H-X4Fhvnz&crm4`xl<YRZ>5tgO9x(zs~lI6nmQ^n>I4^)ozs`Xh%gN}h^czr3% zs__nqH)3sxMa_Vj*8`*}h1Z&LU=@x(rT>7^>eghylcLmx!bxP;7^f92?q8QH7tYWx zCd}pGLff>ASiRg1-9detaKVT0>ri7AE=hEOrAPq36W;ZvX+5Fwd&hw^#Q_vf^a}3t z_!h7=eh5YEM0VLgHZ6?Sznw3{zbW|Q!S{*BS;K<r<78CMy{hm{$s$P}ev;{)tkFL# zaVT&WY5laZKCkVC-sPNM?z+~)zc6wUVswtNH}?y}DJ6e~1-1QA{f=Wh?@qWv;c!Ew zj9RhjW>kP5iDCQbiXD_y)#xjPe#rYsC@oQoNJEE}r)zPa`TA4QK?H+z6B{Z}`KOkV zD-ywcCC7|I#M*VF(La$~-y`pF+=j)>=fRVpO2qbK;+7XmmkNNzfGD;tF{tD2)d{Xc zDBV-oImOW+bYP3+%7~5TtCJ{sDbi}F63!nJNaj#$(EuV1hMKBmztLZ+>Qa4)Wh5Xi zi)^Fd-OIv$XNsL0Alv=wj-TFl6xKB6IHxz2l%SVb>#hU~mK>e4!3`uH_+%HPgB{JL z2aCXz1|b!!tNjHD-X;5Y_p10;XMZ93O+&;?$jc6a#DWPWjB!l&fs75;&&O+6D?TH6 z6RfkGl=;+9E<7>5Gjjfd*xS<$3kYgmj(wd!HL<{5nReWl3wb#7QOquB7%x{u$BV*x zl3T3&f4NXvvkpBeA)a04$I-a@y^{69nb6nZwx@?@9YGv<&57bzE-0PSaqN`*BKf^_ z>d9ZErs%EK+wRbkcj#nG;s?=CZBZ$AWlVPly%q7<nRhfQL?;t&(m+n~Kv;Beo|~(_ zLy~@>WQUR`WAYcwxa>EDmYf%ePA`4Fpj}^H5OTt6#UqdVS5!w^nXa>tcqm2n*{#f# zqgnp6p2aAy`KY*Ajt+LHUvK=C7(^Ph@LK$ufS&K>9-3!7H>p?Yn02!)MfIIdEr&ZC z$4NvlWEq+ViqLQ3_-_4_%^1B3W;&{J%6Gw$Cy^L+nx{)>=`M?5du0giIbfKM-n_gl z_-OHbkVyMJ(M++^me0tLfp7ckWgi=q^U>uDq1-y>uV#-mVoKo8{(f$8Iuk{}(V~em z`>ahop0M19>-f%#j7F!C59;Kvp6{Z8egH#pWZlpRoTxO$x@b8ZG!=)e*S@X%Y`8$4 zlKvvs_kJRjrF@t}O5unhO~-N9uf)qw<{i#aBF*8F;c_>BijOyayU4)ugItqXsKAES z>aN`Q8$881EHBDS5rde;cJ{4*+C6p6NjG{i%M^F;?)Ax*Kvk2WTYy_HfcW~^<Ilk6 zbJ_F^dc4N6t<0#~o<Y#$k&+jr#A`17{dBSypnL&<-&6EH9)rBB3sEKUgPnDmR?afg zhbgg`@AbDwr6Tr&tc66D7X=ENpNQQp26vOz%A80)8!E+HybXhx=A0CmtWS{85?U&0 zRCCauZW!IQUywe!im=!#JL}-n>r1$P7}X2;L(YsUfKrDZ%bC|9+o5ONF}$M&Hlx!+ zQf<R}RPS-{l0RJ<|Ln!a=hoRwckVAUa#M)Xt8|<Fw~|r4srO&Y!yw2W33mJgHKB#t zm6@^Km^zR)(tz@f2<|3~DW=vJX@S9H4HBC6O>sVp%G3`qIZw$pk3nJ%<deVNW)%^{ z4e&GG8cvAv_$lqm#gi>Zj8d8vV8vVJMaTi4x&AIV@^Z9f!cBUqGkf@xu3N()MWAF6 z_|u5cjXHzG%%F^M#LrG(j+qW#8d?obOR?kW>iu=9J^tk^kbWa3bFxxRKJG<77=Rj5 zl)Fn7CR7$WB9BKecJ`G`J+d0DzCd4PW3tF}0lGvbtVA~~)7u-2&c{)zYd_f}g14g1 zU*a82E5M;6*Uk6Cz+h;3RMh@k$h9=|9WJm_XUf1L2AYUE0YLYHM)i~AeXNs8A=$e2 zABqowL#llx0BKu8TEt$^LKhw|TK5>5l}st{ImwCFw!~AFGD!A~(7@M|_G1b>nL(M? zR#~GFBW91?z?zmofJ0lHA;9c4t9CQi`3=3$ApAibP0DyTN7d{WEnBt%Hfc&T|0qK8 z6CW2P&7$CqG)A}gebw$UQ*ZrtD^ZZ$0@{!o!1PR8b~6{S&MIJaUToW4-6<Pi&Qm9) zG3vb&<P{nnXO@z~34qgdR5xC612b=!AWr};EapTi333$C8at#gqM<E=(d(TPcE7B` z31;WoB^mzSnH%?oGp2ays``27gB0|A%#4q%GqUn$%Vku#+3^;OGRm&GAwzqo%fM@> zj7fS}@ptP9RPGF>i`V4<Ho9`tZNCAI$|yT3ff9S$*S@AN6`%6Vk#}V%2dphL8cfS! z1SZTJPUUrZ?0Bsy+CIp3<3hC0QVz6g@(MkU8+RKRI7&4AsF_a^j<9V}AbjlkNM0!% zGjxD9tR*`py@)ax$ty-^DYg{66uB@G#EuFKk9!*E97yh}d&d;6O=HOgy72+H6VA)t z&@QQZCF<T&BJ@32k4-UeADX7>x6jt1+pTo<E})R-i8aH4si%RT<sf;k0NQAvi`&os zo{G*kFTt9Yv$EqZ{^`J$rZ(pyGPuMy;fA^lq7#!CNphD_z#fntc&HE_{OVN&e@kx* zlX?g4fPCL!P?-NWlX@E5Q%?eZ<*$AST6x<0JiwQCRRe<veZF9%K~Gfk&Q)`J*5zM< z*TQgkHM*QFH|A-+Cx&?<7W0!bji%}UnW5im*#SYf!zm)w%k-^ey#0wJ0A+&<KyTd~ zJHV0xnyMDU;x2}dKuUe&4}``W%CQWD#RO-Z{huW)+`pOaQipnjb=ZTNSpEA1a@`x- zWFS!I@J5U@R@v-9Ny84!-AS37eMNmB&j}%}3S4T8w6(ki2KRC>I<>2FO;9%uth0CM z$yb?q4y&~oD6@_NkF^Gm1Vf4fG9<!aeZ=zIE_{<;4Yc{qd)%v`rjI~T_$p63>@5u) zd5|E?X8*{Hi)9M9hTKFf$}6y%r0CwG$1N(nse=+1LWZX(e>+jmMu?BS8gVcoB9#g> zLqRA>tj|{Gx`8TR0KVb!C0;${%@;m!XpyY$H&B!O7}K!wy@kdf;9ha+wX{xDtdFi1 zHO%qd*vv6jdah6%j#}(0RZ0l_3}ZY8h0$UX7Dvr!V&?Xj_J2~{ym<>?<9&zlR>CNQ zB+UaTj@6%Ts6>IyA~4Z)q$qA>{Gq9-0_=v9c*0;|8loMnacRY-C6gN=qyj3i3N&F8 zy76h&8emq_F`26u;j^z2p#DnEb*DFz<v~rApT74x-RA}%?7L1Gl}^f2k={4*J9zbW z5l;QhvkGc435Jx>5`&k+Vk!Gs^Bwfk$4q_|&`uBZDPEYpsM$q92J=i1<@?@L-;Jy| z9Ql}3C2fmuzH|RE>m?G~iL$xSk-YvcqB4NXLH3zyQRJ!}#*#Y{tOmiaHd8z7STCE; z-wxY2xq^MO1ryd`9Tj$b^hK*1PVWP}C^XuK(^?>cr!u@TyZV<okiEQa!oz7)JxKq) zc(HB$@1)0ZK!w^aH5eVmAW9b@{mI`Aex42A@d!ZEhw-{M5Fz4TF^n9e|5!+MjVN|S z_>LSA*oXJL$0|2~SPEBMnJw~4+j`YHcUZg_=QeE*v;Z#%(!vHVle^AXi881}tkI1@ z$L<uA0X9ER?9{#FPtY-RY~4Zs0i)$%N9jr(O6k{8;gh$A|GOUoVY(@LPq+zshmlJd z&+yNcie2sQ-f)K%S$a>YZR-&L0#swH4#Ax%1i^oosKA=Drg*}I{$V84PQb?d@4Ow} zWK5VPnQv<mprUdxb}Y-crg~*WO7Gym%Y24Mp?Qf4-%;`NwGY35G+$Yt`v+Ozbb<q3 zD=lumx~lENVHaLI3M$O(#&3p#)hTaDS#zldWhh2@IvWx61<4}Bf+a;)hE;y`2;SGG zv>o?cz;9c)6Zg!j2^YB$uB+IDJ->L_0R%oj)NmX7l&6_LsxsML`x>Rl*QY2X>#2At zs%yIt@=t7~a8!2wallM$DxPnCv;(e9ToWbq7j-)4p?2pU&-Jv;7VNL_^klf(@%5FI zs?Zp6pO1;;@(WJ05=m@-<Nlf)qxxeUw=tZ6`EMoV!8u)mitgxXj<TzzE&dAref|j5 z-EU?7J`Y^MWhg7D(?Q(-k^ie!r|vU1yHM@;;$Pz5D6IB6JNRPShA3VO(AzCNYKkR( z!|tN|*nrS^8>*iFM8oKBuna(aQn7nntII*b;a4hkYa}U+r&3~^!3P)mjqHOMx)H>9 zzVJ`Gtl!(7t|cWACz96@eqQ%+-2Vh!p9geaOO+1pqyC<Mw-%;zr;+A=0HESmBR6j0 zeXk{}B8MG3f5Z8yjueUH8G`^<oW$=f<%pv8D?{<0*BMLX(ov9a_bb*r^4ch?7W*sL z`~MH2OnZx0-nDCYHtRP6{?Dn`{|RN>kV6rU2j^e8xpZ9l%S9IMsv|^ysn-7zP-jVn zlH(PHe@5T>)4J-dBa3O21;i@~q@hb#rCA5<I+wf~&lb(bD-n5U*FPWpdLZ2rS`If$ z(WfV2Dlx~0tHgXWu1I|4_=8Q%MEz}+D;rN0X*<mIxrbtE^)b|K+5VuekNNtC4h)`2 zy+;3eAhZmW49&Exp2;QI&9qkN)g3FmeEjO<jXdhA;L_FDow_B;G`C;7@b;rOiI;K> z!Jb#&W>rUWIDy+J<^nU@II);EAn`onq+pr<>2})Y_I60~%dXEUIwmiQ0o8-w`CAJO z%#^OiS^_mCC+K0tX_luVU{STr_?)TWdAwUQ&g+T+ZntC*o^97_j(?Kkcm4aJC2d9% z4lO)05q*?f-c{mZS7Zkpe$t|a5aF8Hr$2#+Jh%rx+5A-9^6*__tBdv-@VPLtS*FUy zG*Kn=;g=8c%(_)Jr=m{WFuVu0Pi?nRyokE-TkTZIq!XL%mY9{ZZ386T<JcJ(K7YLa zD@Orb+U6W3^weaRuO(W32DiG(sZ#UpSep@Y?IWk<+=;__eT&jmZjXb}lB!<jgTZ_C zJc+~Cow~X3)#721U9+vjKYmVTKKP=Qa9$c#UuLAb*%^@P{FfdsJb8^6#msvb%Zpv; z<tW~1#x2IuJ}Y>_mGUCVp1JTN>PA!cfL;km5rM?@T)DeM+o!h6h(mhsxu9L{;KNhV ny=L5Zwa&@?@86BsVA|wD6x$9yTL!=Bvda`=X;fl(_5ObYrW(;l literal 30542 zcmcG$by$>L6fZhMHzLvkN+^;73etn5Al)reN(v0kkfNk0-3%Zg-5{+XH6V;q5(CmB zDKOMf_eDSbzH{$6&vVat?)^j5*S*)Swbov1?cdJ`O?4%55=Igb2t@usSza3i!bJnW zt`Xq^UtV^<rv!m;PBkCuDgb}Q!NI}B#l^$J!^g)bARr(lBqSmtA|@s#At50pB_$&x zBPS=PprD|nq`Y+L5)~EI<;$0;sj07Axk5ujLrY6bM@L6bPtU->z{tqR#Kd&<>Q!cD zW)>C}R#w()*REZ^ew~etjh&sHgM)*Ulk>)n8#iy><l^GGb?X*4H#ZLt4=*n-A0Hn- zKfi#0!0p?&1qB6#goK2Jh40+CbNB9D5fKqlQPF$%?um(siHnO%NJvOZN=ivdNlQ!1 z$jCq-5LsE-`}gn5$;rvf%PS}-C@Lx{DJdx{D?fPfKt)AGRaN!j!-r~WYU=9h8X6j! znwnZ#TH4y$IyySKy1IILdiwhM1_lO(hK5E)M#jd*CMG7Prlw|QX6EMR78Vwd9zC+O zw0!*dv6YpTwY9a4jg76Xt(~3SlP6D}K7DF$Z|~sX;OOY+<mBY+?Cj#=^6c3&S65dz zH#c{8cPJF<;o;%w>FMR=<?Ze5<KyG&>+9#|=kM<igTVp<0s;dAgMxyBgM&jtLPA4B zpFe*d78ds6#fz6OUxtT=M?^$KMn*<OMMXzP$Hc_Edi5$cHum-F*KgjuiHnPikB?7C zNJvafOiD^hPEJlqNl8sjO-oBlPfyRt$jHph%*x8j&d!Fz;W;@uxw*M{d3kT&zRl0i zFDNJ|EG&HY?p;w)QE_o`Nl8g*X=zzmS$TPRMMXtrWo1=WRdscBO-)U0ZEam$U44Ch zLqo&+_wO4U8=IP%nwy(jT3T9LTie>&+S}VdeE9J3<HwGUj?T``PoF+@b#--jclY%4 z^!E1l_4OeTi2nZmfq{X~pFa-{4h{_sA(6;0U%q_(`gM4Ccw}T`baZrVY;1gd{M)y0 z6B84Ylao_ZQ`6JaC=_aDW@dJF_WSqmb8~a^^YaS}3yX`3OG``3%gZY(D`+(O$B!SY ztE+2kYwPRl8yg#&o10r(TR(sP+}_^a+1c6M-Nj%qdwYBP`}+q62Zx7;M@L6kEcW>L z_~hi|^z`)X?Ckvf-21chTM+0P=z+Yfu6M?IE8%+@y|jbo&4Emp&2iy|74<{$m8>7m zPR)LkKd$P|_YqQvaWUW}mWGprt6au&{Dh|*$8k9U_=Lx?lCj$-w(jmJ<yEGq*X-6R zkqJ6J56p7+)Z1t{E7@-zbJ>b?5X29GADfr9iGo7llXQ@maJ)<ZFMlAVaD6u!cc+W= z%d}a?3jECW$jwE|L_|xjbC2or9qRN@G~WDb;w2zVxVo0b^qqI%OuO_*N7!1`5Vn>D zKAj$}$ETJ--oTNa{r(A83h0tj+;U*6o+=bjE~aT=;MDL+Xp_AwO(w?Cp<~sAoCzZ+ zWwOVZAH6wq^tFC7kiN3)J*?uU^Ai^{)e`-iiYNQUj9x8YNKpdYrOpfjp-c9`(ydv_ zWXOdSZVGI|trL%qII^?AZx6Q1O+}H-UJXxIj44#wO3OWEZoN8+XQ<|Uc$a=`VY}q( z#IvH3q45v0v$Hwy;r*2aVcbPgJGiNlpn_k?xVcpL7F*x@1tUaGj?l4?VR`+g-ck%6 zSSp~%_-CPQ87a+@r?F3V^Dd8xOJc!e)SgvB(eR)RcYncbU_07%AUo~4Nyl@XSo{Py zRK#?3b@%QwWUN;S`i`&_<UCj?CtrXJeHe)y96Zap_6A1f)zTh^8tZK{-8k(+6cXe3 zo@m!SviC^tS20@CK^!~8i+bs+bu|fnGTdm>zPel3W~rJo(xEQpo?kPkZIE3t7dqoQ z<#5uh+4$8&E2%Aymzi&t$a;nx|LMe*x1_>^@S?tk%|~wrUIUZ!7dz_nK`e#3+hzx5 za-prFJcy4XCF1SnX$aPLRfu%sW||U|L8R3g8B&dyy0}QIfbY%V3so~%n;!^j1C8Bu zPh5R<T5rC2VjOt)I6YvZ-f4<a(8KkEi{EG6-K6P~QHA^=RotHbM~%jhKJ;f&6JZYD zMuwK(Nj-o})$iiU@U6X9XyX|t$NZRrHBCA%*&oEur^uvU_v?5bIpwwSOi4W5bX6Hp zgm7rfU>oz6k3xDeT*$%3zhQAk_uJ8~nqj8MH$O^D;{Ap)9c>V*yjRC2s-|R02gb9| z6h>}RwIYykyWtZKlqh7LCx}9p_}aw}_%nzAXKX*BC}<u;b*U828)+XNw8~u^jZ++5 zY+oGBhRE&kN#;1<%hzWYHk4I-4YsKow|W2@wW@r*Q^o;{x<$IPV37F9QA<lv0lBBV z{H=@aDl!U<EyUR=6493nI5Q>b0O={-PUp@GvF`~sMdXf<V@Rq~;~023*?WW7gnHg8 zfKwOljpUg~lwV&f&%k`n<L#|j@E&15@{lvWqNEcGLh%Nv6Tzu*(DEQip+C5(9n@5T z4-qB3`1BkkP6a~VvIyD~2-<W2BR~K8bUjEN4+KTf!H-Gd$HpLN4M-gL1jqn`GI#LN zXE^9K+W*hQz5@rf$?!O}x!?LTr&U&rO1H;I6+Jv1A4v>`-ru9dih;3&ue$yWS{#Th z!9({jM}eD+K$({yH+9e9JTIpT_}uvM(Y8PjP^Km>2!+D=BLS321VSeWKo002?b9J1 zyqIdxpLVd2nOB1Q8OeXQWR)7~+9cKA|F`oVq3@#^Rxz}2D=+BJc!6-4E<p{mOaBpc z8qDuGee7rQXTRQ^L0kr9S+|^IL(TQ$3=HpcWs~&g^G1?#MP)o@gNIKKIC8@)!@1dC zR#(2g?<AJU8RkP(wp@TC6hhX^?taxzD(O=*3uoIXXO}gK=5?Dj>dh3(OW%Rm&@d+n zm(NChQyVwBg|^-s>AoNq5BT!g<~tqA;!vRw5q>P-%W{6}eC(35!y<X}CUD5W@x5b# zWZ!%>IHB|W-<5(fv7Bhy*WiDSx0n!FQjNC)BO(oQriI+ZLlTen7#)1CYI}ORL;UyK z+>g6#Ah=E&2xS36jsM-H|7`I=-JjhBjs^!LzC?;V#Ye*U{_GweG6EE&euM)J1Vb(U z$^}FJXJS3bvH2yIf0gi^r>p=fK>rJ5cD8KjDd)GSD!Y!FKci`G1eupfSG%PPifG}V z7I;4{gJvHGMS?Q9aepbhJwAp!IzJ-rCQ-Qq9@>409erFMaqibSN%D8SBxszKsuPbE z3M_G5MMOsAa&&%Et<gUQ^xU2tBYUz(8_3+_g(6fP2s|d9g~uA`hBeUq)l3>BVFU^) zsMRu;puM{(yt?HhD{u10uyTySQtut-iD6CHWA5cL(HQ!*RO-LNPsOpuiMuz(kJ>Go zHtPKh|0>d$>@3&qbF802|EIP2pMSSecMU#%LP~g|;L{gVHQv$i?>5Zhq$gvLAPQ^E zptigIe|LBk)w+2soOaHK^Y_#L*<pIfgqXKLqg360IhEK695{6(585ZX%6SktSUeCW zh0mT9L#bWU1~tlf4IiQs?Ybv(??M-R2Waaku>5r{_%APC#&oeU46=bQT#!}$I}c=G zL_KV{QS!QM``3gjgG|8N09Gq#GI9eATt;gFePV$loM90}*k|3wrtZX`3x^gvkxMJP z)n~|_RgnK9^HpO@dJx=V+t(5Q!dXe#Mj7Ltj@|xJx)Svu!^42nFUWN1GU!74SZWlU zg3V{b!^E{gXvBtksfCdmfmICZ<d;JPX2XX(P#M&*kyDjNMrqJTEYW8A2r1su>!)p+ z{H3Tpj#{K})<C<4+;A4;H>FjZIy^iW&pF2KAuSt|Q5Nk#bN|E&o#50SFFv66Gc$<` zf){RI#oRWM{Jg^2vS<VTXC^1ITi3iMXYcYofFO0`DKI?C)43Q+>De!U5Bj}~!_XO} z5b5b!4zO%CXMAO6W$*WzTju``f>GV^KV}(@!^1N(VvqY7Y`L8!ei+yMXE9;%uEE1P zFE{IXl_5xnv}e@<sk}+#rr+Iq+PVKZub(O0IZrC`W&^f=$hol{K&fn~43ISv<g8hF zYPirmvWpu+J#o&x^Wj9))70p@`O2{Z=mK`UAm;?|g>@GR`9oN9UzP0BT2fgco25lf zq?F&4a&eVk(^bY_Ou7VKQb<rxXQ0wlNqVD}sAZ~>_(sm?N>aM5E-wxZ=h33QE`*9f zmS8bn>;~V?yk4wO?8KTz6JCOGTDs%&yK^JVqtXbsi|dlE$2>L+e`brb2;Qz}ctu`6 zcC-ynhzQlGe>cmD-{^z?+^s0+DuOMWm;I)UZImPTwOHQJTEN4eg1H2GZ&J$>EOt3l z0e&?9vMywt{k$2Gj+qQuXwlMeBQRXk+O)#&P?j|F5=&6X)D6*)r=-CwBq&(#8E%q_ zk8h?<;(4NCgeuoID=IxvQB6C`)22QTYVGYH^P5d{7|Ws3RIk+P@U?P@^8xqqC&M5V z7cTrbia^S)BX5<G4ifv7tsK6^e8)e;ya@QWPHK#4>f0FhWDr!Qb!WrnT}q9e^$>}R zgxuvT7(}7bl&YoCcX4qlO+R>G${F>9t&DYd4Fj!3Q28p%7R{djWY;G-<r7|OF}X(H zmUwFlC6b9B%`}>1l3R;;ufPCzu(!rJot1PxRgpw)Vbn7vHTzX0f@Yq{vKEbfQ9F9B z=pbqlG+6gz&n#3vi^ccxkL<2Bv4hG_I8yiWHu{YUDr5+dOdRN&<I8vnG}SMRxe+lO z``;dj&dL+C{>*IPc<d6<mfQCE8?M5t924&>YQ7Ecw6mwDB(*A3I^8Id49UcXb_s7{ zsCcNX68N-@w~fpW<y7q0@y)8^USeT!n45KeG`i~>FJ;s*71VWS`9U2S6FSyJY)BG( z(baw<<v<qPq~<68LgZ(J{=-)s1*2B^E1gCM>}ZYW<s6f;tV6TPcOgR`jSdLNuqy?e zAXLFFo7{X~=|c3~#sv65Q>=~e&1arS9&A;8GLzuG#%zLG(YzZ7FsnT<wmH|0F|~mi zKD9?;HB_Rr<J$WiueIf-gO4m<=Thmh*q2D{+F6XE(mNr*kS<vudCtcia&k?NpzCk1 zy5m8-dYgU_6E_u+(IH((J`e^bo*%QEugM#MplsIy;K}Ahy>=va66_bZ@8$@XEVlkf zY9O`DT51A5;0fpPuuFGU>5uzpmIFRty||v#X^QAvM$FlgwJ+>CoD$7i@=}|!!Dehw z^)48LTnw4ndcQ~wZ&rpf+{`AK6+Y!s@1_EK;57w&UDPHW>&qWz4Ju!Q9k|>N+RPo^ z(Aj(0;p$Ts98_+G@gHHwf%AMMSN$Xk?}<4FN$HaUX8(y0lOMNRk+~lZLhsmCNY)Y7 z?hg_Lv6TFp9~i6O&2miY@-e=J4TNRhR=8(C?AZT^&(v2HjIIJMWVUOz$f3ZfQz@5E zt>RTVdgg+lR!C4=ZRIU1m-h-Fj3M8}t_h>h^I5XIenx@O2#4spQhyU`1YhQ#8IU!3 zw!!Kw9^!-q<jP^)V8&S<7RN+Pk_t188<)wERzM&j4|xbrZkJhZp&v+UJqfE&m=88* zuNKnk&%iU|O&L7T@&NAt-}AYF?CQ#;raWB;EM@1tOo=IR6D=&OmSKMyi6?b!S?Nmo z^g^^ANE`?kHQ;pR;@svt@(ZU-O@5wzu64c;c8z+Mgm(0N-)<)0S$}VlJTBH_##U)} zoZE++0p7fHh}@EQlP~zjbQyD^g&OlS{<BPS^n&~@4X1!@so!0Q*JoF_l*|o-v`D$V z(M5PW$ZT|`Q1N}>4PetpNYE9e$kE_4d0J4@@y=Ask@=7K?b+yjKToQBVz_Fi?PJ~I zcf1!{M~<F|q4yisC4e00Zr^Ss4NRy6p|$6W@tSE8FODczoWY>#II<D&%_A{I?s)02 z(ReqO*1DH?$b9MDrW_*pxFEVvwCSod4s@kD-ME5_dqqQs7MKlDw{+Ct#oa;`ES)Mk z=`bh+YLa6mYV$gOeJXG!^qtK&<vN1P_>Q-N<1bwVVGvb2H}2UnfKb=V`kBjbgd)Bq zCdiE%L{maxUp5VVa%ZBD9T!)LJP1oP@k~jgAB3i02qhtx&vqRon{gcYcEbsKN#m== zk%g}{tq5>~UXw)P=A9+fnr;ZJzNm@*gm=2`>y&Xv6j>MsIiLAtsRY=)oX_L9+8>2= z@_VlaIxaTvIlgZZT0g{NA&}2hTJLcP{|k9Zp2BHJI~9)E)qK`Mgz_DvU&O$ZVH~VX zK_%IEzx~Yb!FZ8$y*G}HM8{B0l|!HK%&e^|Gu(pOPD7)Hglh^kZ8#Xw&EDGGPlv}3 zF}b~bjzC-kF@z5#e;!?XvG@re7M@WZdRCmeEfX2Ys&mmA@#BYy+;F0jb7X<uvq?KE zAcz4}J$S(jM@5Uen?-VoA?`IiYS?MD*U}UDIio^TUq~eAIp-XtA@QTE#V13ZtrtKP z9P|y+F7~v;g5fyaVEM(Vv<(Lb8j<KU$^VKOV|(={D=8xI%^eH=4@O?a(U14mT%)Jo z0;<uZG~v=#o^|q!5N4)-;2s%QcFL<|!*K(BeL~javtj(2Djhv5C0xhm37pa(L4$s6 zuKEk$QGSA2Rmi;9HA)*p0;1cGw*->;!00){d+ObZ84dcRka3Ugd_gZ1;yb*(mud1= z9VvgZj{C2EX47v!NK@&Q!`lmIR&%l(?W~Fbg#q%&8&y%?zHyrDm>+D^xVV`sgt{}Z zUielPcz=J(@uvhHv>(U!$By=6zRN3A(XpFA8p0(DvAAwLay?lnnpC{Z=XrmS%*7bu zirK}ndn*CaX9F_yH($+yF0#z$_ZcLwO%5|B3y~rF3u47{u58xGa1aAg(@@POmZ=;r zP&>Uj=c#K`g>%uo$=t><sY+65YxKjI-84VUYxr}A93cD|zEmHM_^1QUMF84VEL=kw zp1j%?eJXwA61*H;r`0e-(1ccKiD?ci1kCO`30C8CMxif7Q9<&r$R=N8yZcUPd+~u$ z_2&F+p*=$v;`4anOV4rPLnb~mz7~97K7{CxWlSXJu@px#D1Qg_>#x3?=7*X~bLGeK zh=R<DrF42;1S&Y0P=*JRjtE7%A3}v;Rx&oPUfEvQxgh$ZXSCEi_p|m2?zxy)buK{& zYnSlYFdqo6y7&F)!9vCHXg6o`uuhZ!UlOxN2GZbtB_^&ptO4l!i175SPlXv91(KNN zYOji)h1S76vvdkbvUtyMjoica_4E4`6ej&yma(Q^5nN{(_tr=OwF{S(DHNPJR!HJ- zz=umDWjwT#g8m$Acs=x^W+PH4%Lj@}jMUT1#N1x!L5I?()4m8Rv?s<OB17g$?uY3@ zoQbJ4&E<_26R5y8)(jK%_^tE~9=5AJFI-)3S!$)p;ptukb|V#85QPWQjf^y!@<8R| z$of?#ZBYi=brM>z-*qqeroHZt=D6;s92Ql+&uYw)okm_1H1Qkk$<@~yTEA(>TAVo* zvWJxChmhcKea*dA?ORGgjah|QJ`K{OqIqqgh^L#RN!yqx55d<{)T@W}Rh{Z39epJ! zjy7exrs?RKz)(t@K&x5uhJp76P${`D)T2%c_M{xfj2~IVdp1ik7td?PCN1CadL0t1 zqJ)~zQ^9`qP)f2a?^%+h5#FGJO3`XlZI}0?-Ii?+^g?BJlsB8@vU9Y0&CkrweC&z$ zrxc`qGA^ZR8e6B-OJ>e+vM4(wbi0==$2+9Wa^GW|C}`7A2)?m{yb=9FGEaNr24wd8 z(rXA@=Ho81gXOxA;Y4_24DmIOHTRv+i@Q`3soG$n`Qm5fy#-n@(usSe2Q84cOv;^T zGSJL$^9tI|o8>AjA)Jhbj=^}{{&sRzyZAA!r|}^PU~V{2!2skHEh=^t8>mSwN}>2+ zt99fW+957pW*e$}4Vn0Xb0sE$Efk;a>E`j~d?uCrb-*l$NHOD>wPosGopzlF7vJM2 zy&m|k(ofW48v&&t?di|}T3v$oJFC(OCFy?uHLb;oRhu+I^z(3)lusCzoxFMw+A=sG zGs0YoL}28-A_!WP@zH3vj?<<@MwAfnZeJzbQ%Aw@<F(#-9(g=SECVKOou;VF%kkNd z5ji?Y@@2E^QlTXCfRn22Ne)8TO0U9>$<E}rBK%zHZzx*AAQp6g9OYp`aKFRz$M~nX z9&pa(14~d-g7WYqtEZ9Mc&f5XzS6s;x%(0ODJ>d<^RM5-$9s2F6ArpMpGw<+p!Wqs zvhuRSB<84rOA&A4hr2g_vCB_ZmkE<wI^%rr)^nr=t96G8_pbLp$3G>nv-@&+N$ky) z21R%x7)jhkgmrDToTtmIy*I4FQQYqKyh_pT$`ajU{+9x?Q#WA4L>mE>nLFD_NuQmp z$~d{vbPD3?kz~$G)EhTF{HLtkLE<)4?C4<D$l@nwj}Z68O!o7q%rZb`!|^*_CO>B$ zXi*dksI(b+i^o?9C=i`5`Hug1{SF@0Hv3I6^`l2?ces|$>`x$rppM+5KXu9Fmd}=E zy|qeL+5apEWETG(b<Kz>w%lf?ah-uje#gejybOI+zfn&)H^lsgvmf7}S)l1)gjIGe zO}zNqSDvOQFgiiu3VfyR@vWRq7Ti|)sq3uW^TN1Y6@Fh@Vx1a*Ygw(U%rIV{>qzye z?%ji|<!j{>z^V`^;sxiZ&Q71liWz&Ssgjz|bqr7*WFtX#No_dopPn>3jdD}i@}y2m z>&X;k-5IQXOpKI^o0z$~r@A|pJ{+%qnSz+43ERo%8?ADg`eobawvsHW)wy2_!-2X= z8Ai8sH}x7(U>j*mD~KdI=nxB<YYHAfVZ^Ju@d4^wD^QvPow+C6zFdx~17h(UE(o2M zeS_8H-nF$|bVE|+H{O{_X%#$VBZZ92M2#R0w8SmTsa;nlFNdozMZ`Y$o?>8l73&QR zXe~jIx{*8?_6J3u%%gQmoE2xc;If5?wDRayM>BkQN!jRjN*#ZYbK-mcD%obmk9w7Q z)4f>@hSQF>%aC<LxEH-~u=@wz9ZW?M8liz0)>L~pW~l+wI*Z4+=tv<W4!&7-5Zr`g zi<7ay@rpq5;^i#a?GI;puw}n#VO->`C?ZJq2OX)E=bj0*d{=_X-wwFHl=hJRaFchD z&*iDRj>)EL6$q_=ccn3{Vn(>s^VtXCdujtd`UH?dUN(>!=ZVbH`6huk2z5s|tr0`b z&TEfF-;6U2!bR&d4_SAH+K2k-@qcH?xV0>~iw6w%OtxYoP?70g6%3Ifv0P7vjCtBP ze3UjkucA7pNO&q}rz8YIy<Ml=YwbGlsdnzaWlo&m7p=C>4=5c8H3mV}eTF`Z!b8UV zfGfY={ywpz-s$#MfzUR|>X}iHGn!ii339?>=^p2WrqhPACZ$m~_#;5*=N-6G_jpHC zUVg&`H9cMC>t+nwxlV9C5J!8<c;|)`_gSDdP_24l-msCIJU}CsUNtSF0LrW;qQo>n z2hQypfGtd4Dc=xWUBrJ=WpBAHneCk~hw#sr=WTj)(a~B*hcIr!J-mHf*TJW+&B_HA zP)A85>!_?3^b;wvRiFaaJ2H4#C878!s*Y67qJGeDg$NIQ9be;)<+_w~l|9;}E0Zz% zr6*^mjNtNxMpZ&sOO|{qOh6DXGt<4Z4<4O{-CCd74tB4)kOr|ts#&e;y(yZokCN~h zxmPzgnfShq1#(mV&g&0H9f<Et#s%44<Y2Vt&=DzO<mS<04LANxk}FH3(hXh1IPf)J z4T(uD<9BMC&;6fAo+f7)E$9$HZZ3X-gjnfJWeeyzpM+H@u-hbA2<&_yft##!Y1!<m z_)r}l_<Q18#ZD(s;@Eft371*!g6)Kfmc;~GeA|Agk2ongyeHqkyc~>q!W2%gx;SgN zTi^XE{s0sNLgg~Ag7f{`wpx70d!G3mcs3Q{0|8#x$`7TLD>lQgXDi=Qdac!&=W^_d zzQO_(9Ack4uUB-N+Aq&fA4fb5a?Wb5bB${17gK6+)R*0=aNI`i0{a)i5~Xc3d(}#f ziW7mFld3<y@RC4Jja0qX7K!h8=x$7volp6pGgfoPS$;Ay;rmc3dbASLU0gB)q*QNZ z*_;rRhwa&EY-nFR^z#<R-qL6p#+1O#7%_-w@LAP>{ur}Qiu+y|T{##2WTxV*^6s$| z4*EK*BYsqI=D|r_f4BBUl0k~CS}ygNvv)1uO;6pG77B<WG21i-p>z7kh&_13m|F&U zR;r&C_7=9{Lgn0r2=<UdV-mHOfU8m~DC?lF%F8A#oqO%@oV{S>OG&ek$nYEpDHmU6 zGw*rzffeQ5Kuxla+X%d-p-Xf?UH2$s5_jGLy$?L2;K@^974mF04K_D<?|4a8-RJeV zhl4@BMHgqrP@?K49{MP>NOEX=@RI>NpxT~ywHabuzJ76|Qv#C@+@;}HF;oQT>TG^t zwD6~|Q=WuD&JUyTPF>hTt4$KT-ghO9nW-Nca)_gqq$yBm6@{DD+*W-#wHM*RePqd) zzwP9NuEdH^a6o90hLsUMV3f84rBoU5nd2Rm%FmEp#P^(1CYpUW-FOUASPX^jevVNj z+Y_P_!Wx5h3vS==OFm3y#n4Z0RW>1~#TNH&xo2<!vtLvYUZ#BB{QwW0WjtjgpC$S7 zD5`H>Hu8}lpHF-Ji&P=5Hku5tLBQj8Ohf7;t~kW1$9|!UTgbx&`o2tJ#US6y+D7r* zG?jkP5pvcztu&y$JF1|iU!oH@{v-PCV1n!y>1Ht?A$<Oj{OzHj4i7o1f*46Wq0#LM zYpqzi)?+OgQuWa8J~cI_KlS0X+NWYq`==6K+UbXOe(%UJ{1?gTOYYiPro3R&li;Hj zu93SEkek-hih&18(WmUk6T~)JV~^PRUMOTKm$}m;fO=SF0HS*Sh40eIHx8!<0)Fp` zE`QRYA2g@B$haTLbT9Yb$_tjv5kB)`$Q;fhzY7#Y-aQ$YP@52)G+m<KYPq#@P}|Vu zw0-uHA__<rC6w1ct_hl$2m8BKWelbKc<1QOX==K(<4*(evS`z{)}3S`=a3^jy_uvL zv{~+O+Wf9HkwU*FDgRaAz#biI6#og#u4gNP7)zRCfw`Zr$Wj<3U1i34`f=$l+XqtU zHuO@wS3M)SbsGhC<w8X9Bde^Uc>i2&jdcGjdhar06`t+(@V|&4_WhnYQ|5uRxv5x& z+16nOlGGqN^n9duafZ;{*LNl!noMDranR*WH{A9p?uN_)D4E|_fNY{C-ZDA0BLyOg zEi@6t;s)wEuJ~pnRz5CC9$$nh0JF%txxgQjTHdsxKg+xgveNM>w&9U_6pM#0FTUJ| zC3i3>kG5Pd48(q}Ia)QejIN9J=ebM^&0ltV>c2__r$&_NvMNdW+QI!#Kd^w1s5+@I zSn^6qyf)NUBhk$h)|Q*$lOphB+5ZD2{8$xuC?6ekr@IbR0X{HcRGYh+?uyEsBw;0L z@1tD3vPi+`O~eN8?w*s6j~eEwcT(Z*&p2}aY#Mm;+We=+9TTmQ1^I&5jMW0*y*4bb zdyTll0C`hzjLTqSQ#o;e5fe`V2Q@WkX;bg}5jUkUUOm+Za$`F0uai;4TCr@Q0zSsm z?(kWk2S)Y!>E#;85(RS;`UT!9TcKf9GZA8VGmoCy^eqak(S%xtRcZ0h1~x(2Y_@9s zx8c6LiR1i=zq)T3zRO6>RU5yF2HTENHnnQMyZTUc^F_^b9H>FWmw6>qVy?S{Ec!L{ zoni~0@L<8cX_Xad_CW_B&*9<uEL<2Wa)q*lU7BurWpRlbk&lzx;5|r&Rg`%tSsm@j z9aP)P&c~ZAC*`1%l$}id#l9R2eeSW3C=s~&hFUCsQoVOjW`rZ5E9#4j!qv3<&l9h; z`{K@VWKGw;<GJ3<XH~PaE8~1*^`?z~%|*f45^uyuhwA_z*3-eZLDyPb41XK>D~D3S z)1~!mtqE-Gbv|Z|IQI#lU9!<-n2?9)Jr$$U%mV8325BqS1SKUZr0+(1$!O!L#y}H( z19z?361P1XMX+VPH{I8bVw>D}!h_sYP2NZ-e=H`c-c(VZHAQszP3<18oIO;XwK(NP z8*_88zq?4nXrt^ckS^`tZom3MW1}~m9f9FLyi?M7-GX$zDcwW;#^)5#4LY;P0G`Y4 z<7T6m_+VwHA)Ae6UN#?9tvuS+s4;Px?HhKTjB2AA#Q46!y{Sbjwt_4|$12$j(Vt^2 zc`h?xe%MPq^o-_o8fztNda#G@`{k$1S?IPzEje7tTJUAn%EKD4;NwoeDrW_@Q1iWa z=cx~G0iOMR@JzoyOZhfm5KF28a>sJNO#JLDE2V+03NkBmQ0?OxACUN232@?-p=ajy zCG0l_&ey`2EMsCfOYm^b#_JyLj|Y7Axjdva3@KINLKX@N668O?-D{)zRUcqNPvN7O z$&lHUg~1@GsF_RIO;<_CK{7d1kl)@*%)@GI(!uM{K}-7An-CLv$Na5~Oz)YVQ!axK zNJ<5H%rKi250qR?A%ThyO<(_D@{|ZK9Df(2EIXSg@>z?t>pZL4KZqr!OEz@3U>tbk zcZ650JcuX!u~17UB94P3R`8E*yQ3pnn&fG|Xu!t;SF);-bI4QE-zwn-YA4m%(LCbi zGuxP?GqVX=Al<;jl=vW`N+${meI6@uGnsgZ7w{9HwG2x)Pig8RA_W<NqV{;d!+30+ zw=SU+4G=4`Vjw1M75s5&vocFQ?jJ!?Zc%kPr~4W-K`i9xk2Agm_#3Q1IyM*D_E2M; z*#t-hzs=8E+=mP3fPC-319J+t)*K;^R+uDQg<g>mu|iv;BCP#@RwMHXu59X_xTHwp z#-KJ?wb=&a-t-udIL#$#^n_&DM!UD$^SsAkUk=)?;cheLvNPJ#sq>Sx&zJjwDbgYP z?`s@HYEU4Hf}TOzcbCcD{FtWHf3ME~iT%LF)tyEo^;%eAn))E$ln9WY8O>u>`IFs8 zpu-0r-h+=Z%^GO;g+>YBA=xWmH+cKM{`jru=&^n7MPnbS5x4t`ns>;9)bXE#u{~J5 zYg98adzA-8q!+Czi=6eH^|)mP>&*FqNJlc;++?bF_vtnXFq7L53wGnTlYWUUF1rlC zP8<NGdUlbQ1qOrt>3!co9-IAlL}W98@>d3&b#vNMAh2&UT|{Ao2U3VbgdE!2XjApX zZB)*r>oh3%_88QZ{RvK+Ui)L;U%YTKf)t1=L1vruD1BjZ1E5YH!ry~#-KMfKt3@z= zOvx`C7uJx0F~W}@SXai1->MNie$Xuf2ob+dh<5unao66{j|V7yBbV;=vZgzSmR(I2 z{8;<q<Z&wz62_5L%UGsP8{ybzS>l)u=#&_wou}ISf!3$^wk@C+0W!9X8!65&x!Jle zvVG<jw3#OrZG#SbXkEHN|3=896CYSn&qY-tc_*1k<82!+KM@$peQ^VVzSHEB?{1mw z7C)(b3BXghql2dV5n9O!IMDmFNYa;g0?gYd)AL$!fac^F*ltCigG^qc1M{~_vsVV~ z9Qj{;VXr=PLt7s+f-+y@V^~_Ao+-u!<Wf15dY%G14x&R6?#d&asm!U}Ey>c!kXE%o zG{Ot-U+Cj|tg`-0nm_qsBE(Hs4{}Qmy$lC%ponqKdF(BBOHk&|7&6SAjp$}^6`L2P z1c0!cXg+zCsUp#5w?I;iNDCX@VLYH3@83x%)og-*tIk@LMx5Mv{O~BuO$<1(47ez5 zw9A!q4haA|;~{?~a*<EaWIb^4`cTzSDDYl&e%j%*g6ZR*AJN)dwK49T6}K!D2Lr?4 zVnU<Jc^oDsJv<_w6fC_I_mbwK1M-}4s_OzfX;=n<H-#cfO%-jyRizRjl#*SFpnj&Z zwQ4f`g+f9Cb=C>C+Y1oSy7lQVP&*!q;j@`7u+K&U9wn-aFa^q{3jP^JLL>}BHl*q9 zWGv|*IKfE(*8yKgDuf=#A9{^3Uf_c4<ezoMJx)b>-Hoc){$aaW5>+V+45pC^5eSez z^$?Z#H5eW0T5nF1MI`|cBaTbm-z=9YF)BT1`}p(0*%$?&2Go<GJOMapT?i3I<CAOt zlZ^Jbnj8dBji^wPx4WUFBrRyeO!@BKMq=nHzSXNa>x>Gw@k}_7N9qQ-3h|#SFG+&2 zTQokMQ<_eK0E9`3O~pLWJ|3GgdK&>eV(hSr7=h82q0&5SM>q<N?2zmtDcgLH<cs+y zTVaZxhSq7FV<UBeXEjs~^YE**WsvWcZv7>Lnbv%+s9Hbbih3Vwij_GwVTvdKtFjh- zsS_)$C3tV+fnY7Ul^Pd-l}M&{cv>ddFDoOy@-`M9IMtn(x5o^AIGv)1BZb{XN1g&` z7IC(nMbtNe-S#XMbz#*g4?0MLD#m(bEQL>x_9uRYH-HP`e$@UX-mgnTlXUp&Nd)_l zLp*RYb&7Z{1oC;-KT;=?SIR)n!JPh?aJwf`Jx{px0+s@a@Boidw)0~A!=tgyT)ztl z7ax^q`fU-}Ma(DfvkD7uU|Hp@(RbT2#*=JJ<7)B|xNL>aOC=Zk?dg6>45kDFr*ipA zOhZqCXLO8SvH1fA<EvTt5U+Qu#xX+jl4hJMBW4F20-;lPL8z_ILs|5mqV3P-BxkPy zt4d=Jnf*?_Glb}J<lK6M))fRaU5dZLMzFQtG7B-}_i+E-K|r0LRWqn5xwv9|(dS`| zwLcUXQce_8$Um0me1ib}yn+zT`txcVN1?w!MT8&eFc8m|pKE6^+KvNP8XIn_@gg)i zH%nxbx*V9*cwB(6ck67xg43}winfl%Zx4i$gXs=V2nKCd5XF9JA8<U&@}_B(^dbhp zN5{(1TI1g;um;;f4}f-08PO;}W^N0<{<`PpuLe)@pA;yavh8ckgv7h9_LuQ_$weV% zeCIj~)+-p%0GOsdw{`sVFwVgeI3cgZvZdL?W(DEa2LSx$C=gd;R2}xkihC~hT<pl` ztePHv++XGyx+}~IGZLu8U~9@wvU36dUO|fXD*VjNr@19h-+}qa)ryZ?3UtlZa)3qb z>&!%HxTLnQ+M$RWM>olAn~F`iU*H?F6on~(Q2efRuGRkbY3&B)VK0q4ghe2f=3+&` zZR+HJhg#Ru5n0JJosfM=KcIt+{86Z2EircT>}a}ynUw%e-KGu{;3n2Rwl(xGK`4g| zp`&MwDO@YPfqL~{v+V(^e_dv8|H(sGOp~qTtAE#JLF8-+@IvnzH{5m$@M|7H{LL$t z$#L|4T7$Tl+c7TX1{LyykDF@2(qQqrv@(13h(*}9h^KNyM6-Ei6YgyIQhU?<`wh=X z;8Rf5)~$EG6=OrfkXWR5vGvKy`TndsMxPSM>x8({zw1X;uv8gm99?pmdMN{xiz{$3 zjBmOHqK`fiw3c#IR6f^i^U#OBL4~zDt2>G4n8m3YU7^#oS-e6GR#>#P)d7bo=rV{U zgmKc?_ZRWr8i;ywx7>h!ttjBTT&R-yoB)1KxCZu=FOb&?wSH7-axXz?Ii>d6S{294 z9@lmm|E(Kq<v|#aG1Sb!;=ur|F2kGgJmG5;v~NbOY*OG#?wtc0sU`M@UZsZ}3G5|v zbeYC#{ohZb%_lv!AD3~sn)t=r_jq6)ue|t9m-(~qG<TM-T8U#{!1p=4eYdRQhm=?p zt7lmWSCG2$XE=Q>_h@CBRGy~vy}p@^eS}l>8CZPI0o+=^V-@d61<H*4CbAYS;(>x< z3^@xl`X2zj%F?vhtgt=r?7I!~zGhAqHHZLUI;xPiS7HLDK(_y=$#ogc{-!;6a)&z8 z7*?PTAR6sQ!_Bwk<jVYSX&<Eu#Z^HCY5Z`(=o>MY#Oge}E)aFwWO_IApLvgwaZYXc z--b{&1V4QE)ZN$uYW`NpgT!TjelYdM7-DK~ddiapbTLDzh!u=PZuzT+hwl|l2mlF! z+9&aD1ML=d0UHvtPzEG5oZ77f^~No`HQP;xrND&kQxG|-qq33+p3X7p$aIDcdrEl) z-po5aT=1D6=F23L((>*JVSKq780c7Lb~ndO_jY~nWx<EM)rkO@(1H{_&#@wMCazpy z2&B&Z{ZvQ9d%dm)(8S3DT6ePeV&E-fhd@fFi1wvo122_JYFc4Mm=Oq7N=?cikz7?U z%mzX=I;Q&;^g8-CW**HrmRv-klV*tF3Cxv)jJr}w&99^ze5&8ys-JhQ<GJ;<qZGho zX{^k8iNTMXq+`ufXc=kA;1aSGO?)YBmZIiwCYLKOzz&PrcN;g1mdB6;Meu&V?^5EL z!0UCUUbgWUf2p0MitH#{v?p?0c+AJccr^`)<2D71n@lF%HHD^j**7XaLWD2n&TUdo zbrgTXJ%tftuYRV>E4u6!#9u9{lR04IaU;VC7yaB+h>JzxRQlX6MT#0Jn|#3TQ>zr@ ztTjh!wABv|!B|}lQtzYiP22_*zsJF_Ktyi^z%^{Ix6(9?!(3Lf3OgDqxvL;s-Mb=N zcWt4UXnIrA_CprVJZnHrPQl*l*zlB;67zE~`o%n<@8<!2mj>&}I%vh<^;5o56<(ym zOp~$yO~E{Be*-XdDkRcNNo6o`^I_{fFj6j(8N(FoUy&>g<Ig;OUFgjTiKU0SO3YbD zxu_;P#L<GGYbRAupON0Z0%g?wAZP7$`1$jq%+9CWS0}VEoVb9^1BksW@><hX^61xI z_G|7n1QiW{H(hIaAGp9QMnjKCTc@<U*fkx3uH==#o$-D?LnSMzAocJz%jGf9*in+S zXo8(jDjyDfYHW@X!e9abe^Z4H)v{*poN;f;_zwOAV8=TIr`oZ!JXdy&4ytV&me3xV zPQf6QxWHgnSzdG*5HCMum_Fm%*9^`L$hKv;z^Bf@Mdu=>?QmKX1fC@JVcIOnfyJNj zXq0lqT(mb8I?Qlw2*E+-5t>OS2~OZ)3m<=6T`sL_)&CVD0A|O}P?|Po6>LcWqz9qq zWqUKyVp3M|6iw*LC^J`LASI%|nCcgDJm6v=$DpF?`1C019HrGO%e=hwhd=u7baDha zXSSct4-=NEf_#+$V6s`rOvZwg)<Z|zRTNR$pt%PCa!Hzi-)nRY7ODqK0H}9%QDKmK zITBe?oGJJ~VKQA2R%~jUR3mTc!~{7o_*}WXF?f6oz-NEJvt^uk*LYm6(T<PR+t@mP zl;e9&hEzD0h+lk9HpMn20^FID-YUbVFd!$qgUo(h>W&8Qekf(hty(S-25ty&(|eP9 z3PxSuc&MD2AYQB91WcuY1>g~<Z$N0GCo98u`KYl98j&JXJpCh8vivIv_EnORgcATF z$zJ_DQZC$i8<?_0qO~~AWoDJyy#wHwv}I-3KFmDfrJ?*OTC6)6)m+a9@6X$B1~5nE z8Ytz*YU<7MAk6CyLbKZ1kJtf;AE9n17t_;4IvlFoY4^e%Pjfn9FHJ|YjOv2g9`H^& z0vps_>&E~Tkblb3wF~7CzI!$%_4eDVFD2u!%ah%vS#vpVh@&WeiT-vcpzW_SDMLK% zH_i^Iir+rG0-?MWkUTy^bGu(kGay=e%rm;}I=fZgMH+1BNCRmoW_1ZR8z2RaXB!vx zqNS<c@IX)YJi022B^`v!yI7$~;tc_d6g<FISkYl>_~^c-3Ouu!Bk5THV*Xo(z1urK zH8pzzTstXYpt{+JVlYUY1DVa>3w8f(Q*5h0?!{}<M>g~1(>Ll+#~~{8g_f3UR4}jT z29_QdH*Y#Eg?Si&F<ig&*_z;GF{~}il@y2trdf49&-zO(mV}sNmv;%Jc>idqP$*0n zS=fjL)Zqxj&~Ti%Bs-0xN(yJP8R16z=^@A8;{$nj7JL@_d<Dt@DFn*KLiPuOqog%A z{4fz0K{6$BH!ej4K?ErT-g$(6N;b7KUPb5W4Gx3F`2ct)6)Q4y%Y)wG^aF4Jhau7s ze3;B{G4r@h??r3cZHR#lZi8<v6tVEj(EXM&kL$}cys+{TOV)yiREym0dq1o*X5V%J z))z&D(_p?|>elDUySKtTik~~oq7O!rcj00(80~=g_F@u`@nMPI#Wg3R>U3Nd-(QQp zb`eolg$A4+PrkQjeY@9gbqJ(}0D_WTA31y3lqp8zNI4{n5RgLQYnAo2VtYLCq%Ls| zC>nsFdC`WlGW}EL#XmBILv|JJDT?0mL<Z(&;xO=<%k6=Yw=SBbXmm3(3h~dLjP2CY zHedPW<YE{7&*Rz$tybUHrU_iY;bVlJbw;-{;-L}fmy;yLp?Pr+TbFR5D!);lC6+;( zo1{py1k}^RXRRjM{~}`<8PN;Tz@y3qQdTbzK4d@XU(_xD@zU}fU$HF_W!BZda>d{h z_gfq{)$pJlHUwK?G!YkPe{TVTi$5RHo5M@euO`e*8yp$uZ!)sUaMQ6wt4+ZgkKwfj z5cTwyMETa<5uG4%KsQNn4pv(1M{)J%i^QB9cJ2E{Pv=%%vyBSrU9Fm=!$5qVA!zrG z^R9_j1Zncqyp}UsjJH`_>~VafhQIj6<~41S?mHRK68=#=p4wcDTit<udAgz^5)5!w z_!HqcO)Nzq2Pp)s6qI=wi;rG|!TRFX-d<j9Vd<BWJ@+5s-L5ua6(DNfRAktm@X(lx zsgAcH_C(1=I*&`41eq_cMVWauY-lAk`fA*~I67%~)J)c~VG`gBawQdK$7_ae1p-Wz z+}6++8b}@*K!X8<;>J}BX%Z($iX+40=JPcGmemP1yCtP(V+Op~04^MWz#xMlUPN#J z^ah+Yz#!;x1m2{6{ptUhqM+G{BXy|$(lw`^_z!CV_yiNTatxBj{rd$M=d@-JSLflv zBM>z3H=usx8ARp-!}%LQ3gCV3H(5lp6AY_9yt_3pi9*u7gVTLBvG<Go09<F!K}}i! zvd<thxNYXEU2|@*IfmY9R>NoTGZ-F^gH{Bvb`~C9(b6yS11z(W;T5`s4+17e=`QYv zU;nn}&kU;}@og&kM>T;&W|aZ+$bcAIQvQ#f0U)a41Br!lGnNF&4i6r%IOjh)xCKA{ z$9@2GGAV!*13Uy!_K7+s!K^SGf2#)lq5b3#;5JFNI)kbI^8rhUgHc_@h;$PW0ABky znZ*7!++r?6)IbNOefyjWa_|oqft&!k^vT_-uE|?BNSXqE{F@t-^8(jbxod|uRStxj z{mmb_3qdB@X6*jB$Ds6^_K-sW_Km2!%rBWGVgW|~Ls_vWf@|26D#Zu`dEakFNDcwW zx9uK>#Y!~j4<`xW69pIrwS8a!|9es(sg#ehL#_Zk<sZ7q2nb@pz0$s+|K=}+<t`pl z3hz%BfVc9$Ls9|C!oTQ64Nr`XD@>B`{}oD)lzAL&W0dso$d4ebPgn$vNd8>{F7{={ z?!*4ee;HVeWay8sMX?rDyx_m9`WUo%1@gZ~tYQ3LBDQoHy7J{APbqA)44u=os)%$G z7UN5_!2j_np}xFzfKws2R=ek+R3B@K2$|lW|F%~3$D4>hM3&%&<%91l=t9f^RKL-E zqecD8U>aQpy$M~0eBk%S1O?YbA9dTx8=A(DGDX|Ul;~@?Bz;%oX#phI!oS$7rcGsw zrp;a<0FD(VXd-VId)G4Lm6?}2!#9ruI!{!Q5aW8e1Vr>V7OzMM(D8-90L_WBv}60< ze)mFuM}*OVGkf9M+tdf6zwHXRhY;`Vkbz~Axm=X*D+mDR4Oj?BTu%{po@)^(Ke12v ze{1S>HhzB4b6zvmnSz%&peaj3r74&6hE4*_zu;%6C!6#}!ziIK#9aPjEz!>Rbs<+* zt8RND^GbRccbG&dJd{EMw#KZ`rG6h@#Ndx$FX5rmNmh0k@njF)Ly-~>M%M};mvl|O z%39Pu`O<RJZ{+vitJj;n$|}l&uD#O-5QBprt8{x!w=Ut0CP9(yJBbx1yF*F;wCphQ zB2oc?p}-F#NVchk30PLM41t4{+pD7iIurv@fHFDOOv)u%q0wQ$z0kN$1!kDO=`kwm zjk+_Hu9W3yQ@ZSZAg#GkVE-{&+_P8e%}i7Cl87hjXd1LaJMB=6Iw=CEtfhnbA!lV^ zQ{0|uwvaX5m8^Qcl*3>Vs%O%ok10HJ0UV*>MgU%`r%v6Y;+>gvF&o5T<1z$g1dwj| z-WTdN9}V}`-}og&xiw)U_kYd29cIyB+@312pWap!fVO|yt+#0*ODe&~Vtns^-kNSi z)Ugx?540J6_CA;~$jpxnm9`-X%&g=HJ!kdQE-iV1*+e-_HYSC}36vBZB(3rW3dz+! zFWT?)l?>1;w@uJ;U^T!$^({M6dM9`+aRoYgmf>a_bh1pAz{6AU#AR=?gm63a8iSXl zs@NS}TML5z(XrOzmYikLhSr^Zl~K6-Y_`x$>w%u#|Ad>X5c&E)+}x;kv`3Q6U0`%# zd=^qJ*~T<2#WL75fv1ntz_v3@rmZa*T1N9mEhY@iY@szfF(Gq>qgr78)qz3v6yrp4 zNo%{|k*-{$yJv~^=(PmcaPvXSDz6@CfxcH`X+oWJeoMQEkGox6uV%QVZ*0-vPGPz2 zKSbRndqsx2eOC|U7KJ^)nE%0KOm{i^{X+ftVMC@{enT!VOvy;xxQ(i}g^A>A_3XBH zWVKyySr#H>`+|IYM!v+%cO2D&`9g+MG-+SpwdFHvRD4)|T%PlBCY@&KQMz8Hg+Wxp zb*n_tI!~U5TY{rT|0VF2FxZSt_KeDEWdHnrQrcFh|L}cWwjl#5@7OfI$*R39S)zT+ z;`3R71wiK2_f!b{&E^%3|9`M~>(=Q`pVZ|tM59||%WZuRvVI-+|48mF2M+u<w>SSs zYrc4Y_P}c?%llm$&7SYs4TKZ50K*D2%y{i!;m|b_`D%RTVD_|N@`0XaYF)L6KU3Cx zg3RpmQD^rt+4HcSC-rfr>6{8>zsbH?sa0K^ujVtk+1{YCR#o?|d$|<V@6xHHOxb(j z`VN)2w;#*E>}xoPLRQrZYBuc-MrUmSOHThPTTwLrkKvwW=Dd`&QbGBSG=#4|X$oH* zmr_pBg23qX;%#0t`x$7%LJi`j4wIYr*32JxF00vAW~|@N+}x1Z_vL-n%0k+AwK<VW z5v8@k4jE@QA`rvg4=qz-#_&)F?A44*s&87orOhLY9+?-29=V^Wb>4X+dv)D4{VEVk z{{l+uGVGCK%R~du7pOK%fN=I7wwBfqT=Wd9!!iWR>ivTSV6Ol}-7jJbfEAgSC2O%T z{0aOk5c%aWx6Hg)|LuPfNd2O5Tn}OitJ1p&-YJPM0+3(CJD~Lb>remtoMH!Y6z7`G z3J@Xx!!ngYgS$N@|3mmC!B&MO(=26xVD$@9f1BMzyWg9eN+bbeWqZ;AWX><<WU&k~ zwt%QZf*VneS>6(#CHO<^1^BpPcBq#Ba&}b*>uR3;A`D+}TUpU&V#_Yc0Ey!thBLs} z6=0s_fgb(h;{dT`4i?mEKV>=jm*lz!@_iuO-8S2^GY3>NeveOxq(js^YyA&tmg#@w z)T*Wv4`yHgcOC*a7eu5I#RsZ?uO6r}?44uMUgDGc8a*(#B+l{YzLs3g#SiCoi;M_j za}%y93hH1hj=qZMh$E?6(oMRutWcYlPSQ&xbS-wkb$!q9@|JJHp?>m}cn<>zHqJxI z+3)JApv<3g=K>wdlpx^CX_y7HIGW%%iURo=i1z^MZ86!!0($I!9rN#u^H+u9H8JKZ zC0+AJAdLm)22?CiA?7m)Rd0(c|6M}^+=ZqAu}=IiOWK17J5hoN0Au_9(@Y62@O0*6 zD4JY5>#?Kq#Z$#UhsK1gC?2Y+v_iy~{D<7@44Hi!d}e(3-*)bIk#;X!{xWFA#T8|~ z8a<nN&-EW(>!{JeC=Lcj^)F?%95TxTX=e=nm+QPPj=j{~t@QY>q)B|(9<{(WrI!Ev z5naVCPHEkXloHVKo4O@%mzCx<qewTk_4;FHo?ldGU@PiKutc)n+dM||HU`iCog*i5 z)P@vq&O%2Z$07>(S6AwA00xzD<fwg?%~PFIm&EPq|I&;@98a$tiwN%S_BQ<6g8&PG zEt~k;Ixqg0GJMJ~xKEeUxa~A1U-e(2u>)<JZ||L`!Rx>HtUw)8*Uqb30XWY8Us?P) zhaqhuB`f$qs{WhL3n+-BtW44%y~Cp81!DYyvJ2!Bu6V7wu$}3NQDFOi@qU4JS7l}1 zg={*;{SU1mkIuf)4gTQ`14ie~lIQa9fiar)AEqv_vMZuUm<Rw&xZv~xXe};vxP|=y zBl7{tKO|tl7e!K|&5q`L94h`J>0a<mNNCn7&3^>lMAEq;XZQcxkrLgLUYX6eOZFdO zlxV&qO5ebt{|GabZ9jsJ;V3+qgEwqPga10gNE(dpt3%Hvy|cg~JF>rveQiLHe2`vH zFvyA_ah+Ow!g;tI&Nn0b@6uTqPQrJOqm5rMp@Y4SVQ<xmkp`D-KG9&-)jMCV{Naw0 z!kdeuK!$vs<e0e=)P&pWlm5lzBB6Ow>^&<+AVvOV$p$gY91WWN&v4B${EsJ$0!ihz z@mqm1tiz(Nf=_^!SPl;NwT2WTF82Q3Tbwv5eL~emLXmw!Zd_3vLc19QkMw+8rRx|| zf&cm@g>T2Ip|@n$HzB6-MsyPo8_fL{sgR7;xRttJlAAnVCtq%&|ISHA9mV?g6qNId ze7DtyUxkH8UupPTF46XIO%|HhIuN`oalmViHt_VKFrf<p>nVbtTAp4N3SCmbZyd3= zCfMx*ycfk@it6$`;K?Si5LY1#I{#qy>empGa5<^vMun?>x<O}y)F&40w|Uu0dPppg zjh;gzH!9bh6N1gW0xicgxjXkt{USi8+ZKHi(`%mNnSslVSnsuelY;ja<|e&&)ZSSk z4=5O=Wh-=n>~{{|_vgCPcHs8#7Vuv=K;^RbypKx{UrFA6hEYov-SR>X#T=J>IGGbI zS)&JTg(j^R5@?qj|5t72{YdruK7Lywl3hlT>_qlnheBp$Q)VIa*qgFP%AN<wic>~7 za%7izPC~NEJelEGhr{<i^?tqH@6R9b{n7K>`*}U@`?{~|dR#(c-igI;bHb|(-H^Z( zu^%{<v1Qy-TbXuKA!^WarfnGTH28wB)q~?5@P5mcQOfsRDhxg^oElFB69?X^+080* z^r#TG`Y`27F*Od`bSd6!I^{S$N&{&qG*qf1l#HrNC;M(5hfx9^m-B^-VrY@0M@{Wm zYs~x+D|Qk8TGGz30Y-oV4KBep^ZW#D%=m=L>xviZ9V;^bbg>txpzODKuEE8w$KN+b zY%i|gd|w2klV!CoKtMgel3;i=1{L4~x~m6`4J&wv%&Oj}`zeze6vher^?mGZ27^Yn z`95JC#JJ2qc&alow82>~2i#lL-UqJn-#mzZeA_nSMEn$co0~kPweMTl+iMB^AXbMk zTfgl8loPlq(@}r*or>tMt}nE)j2}8qLs>Jx=Su>b3n)}K25~{aDako3yBBwX8ruBg z7{YS8@zf4*@ZkrM)^e4L7LwX0nKvCRzEgq-&&u9_>Vm<0g2^Hwt1o0PZ>cDAlC&zY z|3O#ju#-VL>1F<{x_;nRujZ_nvEEw!W`WIhndX)Ps-Jvb3pr@I2~HxwSMi}z2Er1? z`>70NtMh<E?mWELg6sPH=Uge3S17jU{xG3lx%%(rWcmwFo&F3ZT**R_c;A5UTHwAH z$a>OOtyJRp+h%t8^h+Kw)foUVf|m^VT!HV4L%I2;%9%;P$L&>R^Htp6i`6A$jE(4z z0qU|6tgNE?>-uME0`}Gzm`H)3H2HJBkj#+;9@<XaYV+v^TrOxI#0=l4OX37BTSCH& zIbxi159i~RP^1-;3Lg|nWtSg7F&RiO;}SY%!>XSO7V@lZ-R=8!5zZREQPxLkPcd<+ zqIJ{1``<EL<-HmCctUs(&Lw6!RP1llwiNg(XQAm87E%-*Tb|R-g15auh8;}VZ)46U zyeC@{#jvU8q+!fSP@Dj;m28yeG~+XeGE;Wlwi`e`q?9cU2~QUxqSjDVmQG8~u(d*D zT_&JXtpk&seotp2aqchABy<s3?PV{Hox7$2=R&y;r_Y-%#0F-^ZxoaewQkR+*ppuH zyXXLGaXN>sE>>-J6DP<AIhOi^9L3JKv!Q4q5%!ZRoyW&qDqjc-1D*&0O!ldzp`3pP zSWJwU#3%TPP%524@lfJ*QricE=e#cVBdAg2n0p1|%$)$A&`^<p>C+E>eE_Uw#Cj<G zIGUUxSCD6j$*ux|24@rZ<%=ll+?!Ip1hh7@Uks<RF6(lK$jVn`FA)3Z<*wLU4rT}c z(UaYrrr|b{!^BW3+*=SX0OW}cV<YS<`<xJd%ZAmV)@#xrf<mxOo-%<1(bOPAu!bz= z)e&@X2Fx7E7X>FW6_ZGzvL=NXuJ?pR<|b-|#<nupRuhVlVzXx-)R3Zyp(K2ub{7M- zXggy~j1ctPrUrA*V1xb=5Jamge|^epwVObgTb_#O`rtRp+P3l@*+0duVn$AL6%TX9 zca-C30tcHG+~%xTHoR&ceL<1ZTksPpD#5Bc@rs{M5s0Y+THxm;Mz9Ajnw{)ZHYRk( z(I7CiJF6@CA1y^RoH|m?|CAqo33F%56Yr}24#3urs%=6yWo$LI`IbNH^1f5uEK*EX zHW^@ytBCtganbG7#`CZk;y|UkX=<l?+)>v<(H3EHM(ZMQ%~!ulZUNGT^CpIOQj{#s zE0kuJZ}=sMeA-|CfmHvHn8|u{;y9cIG731ldNn2#HSZ=Wu7BZWBS<4De$B)4bqMq9 z3(`yiS2DBrGFr;<4Rs;H|6`p_YZrx;s)<p<8QfrM=hrBRAByGGd!>P>4$Oe=?Yk(U z*ZZ?DJpj8l-CH3!c~RtuY9lQmqil8D60+g1duj9mfuQXqfdz$IqhlMDCUh-TXGJeA zkB}!2AsU}GQM}XBYmT_iHR|pwY|{^$d;)jmmB0-DCKEYtwjJsZf=Q`xs;T2lwtY<4 zef!n5o9a(W3#J=Z2?+0YWJ#VEDZVlXF|wi)G#-l5)qaG6L@y{_)FJ?LoUqx4I3BVZ z7?=H(21|d1$St%YR5|i`&JvmI;$cF~N_%$U%eSDNmYXjYcqI>~at_MEdX6{VG!8NX zfwu0&c<;u9S*g<Pod%T@B@<!X<fig`?vMK@YT95$H7RcHSdj>8YoLT6&Rtd>7OXlV zGZRI<Z1{ZSaZ<!jQ{cvmBQilbiPyOSYWON9dP+sv>GKcQK5j#mu!MJ|xO!Lnx~{vw zlXu^=d^6MDt^G~dlq^-<Z6Cj80jUrW>LVg-uZFy?p2cWE#JwU*Oqb{(k9i%x_S3)b zcM<(9tr=wM53Y&CHL%lxDecJz<I^^7FgoheO|HYj%_uW9DiWMmbK{BF>9~d=Vmnj$ zk%Or!++Q*zPo#{nzp-QcN~zuPF*Dkyr5@-Naqy8)V&s*Yu;m5qr>k=~@&M2B-XjvM z1vk&4Xu@es)Eq{$g^yKwz<_y+vQ;kW-RT#LUnbpVuWs~c@sV7h!Hy3csm~$T!Jqmx zsPqy%h0T56PhASnxthxPXko^n7Cx&vT^mHbT#?5t(e&wb5=~b!i9nola)8V|ouk9? zRzIW7&igtCvjky}6RPe#Bl`?JIo0}iazt5(pt)aTeD63&n4*OGb`tuim<))UyERnU z%rx}`=nc%((o>&U+(>&@Grs&pikZ=vxl~p|CE#H;!8BU=L8j0=X#iwp9AG|tX93_! zPd2+0D^p}wy*2v^+9*0$;dKvzi%mr~r;%qDK3drt&)Z=%X0uQJ75>Se<|$R+9iNAf zC5vnYFT?ZBlnWLKdkMvbYgZ<!03vo*XY5eC7Vxkzgv^72muximGnt+PGcPZ29U0sr zxQ-f^n|`E8IuIi`t760Z#CI$W1|OuCY^(|Lu3k?89Ij$cps?P<$7OfwAgLDL?lPJw zQgo8xPo{nR&Vsa3wro;n#&c_kLT#sM4vfT<)&9|W>`FJoG}mho5~T*)yx7(*Nd&v3 z?R1?-6a=Ah>%o?VPmO+cy8MU*VXu5wRp{Q>u$-5U)FpU0Efd%-m6c~+4V3${=%fgl z5hd-f9~0aNCpM+x^LEpmrbi3nGNdXjY-##So5pD~#}3WF)S<_(1WOr%<B-1qp-Sb6 z+K*q0o3XE*A|XN)vjbZ84Zj0?-*icsfKk$cG!49O9=<7cDFv53c5EkiR8@{Y`!U$S zE1U;zLXDe4zE=wn;6^2wleDdCw2{>hPS9$c9kw9R=bRvC92lr{4E`tDlYv#okC+mO zk64toVH#f5Xr+f#X*iD*YBRdPw5F$@6TI@u?~Dq|#Yzq?1kE*S0OOge77nBcLbmbF zz@PG%&(5|i?a6lBZC8{#T%qe!bnw{k9M;c@FnwDx{QFSmA$X<{INt3=qF#pBUbu)= zSCGBwBF8#>IOQWim<>JlW#Bfx7@YA~nvMOq+mBG)c(&VhyApWcB8$@iK-9(3V{^Fr z`_<@_r0%@PY8k2?=LZdLa7PKW9SlC02q(3eiZ;>FMfr)1Ee?>EDh})KEJINo(%lu{ zVxeXj;Z`YkcxN*1bnxd*zAP8%Xb?epmS??ogwZo(o9bQO<|jX&XufM0!3nQ7UsG*w zEV#P`&Oejn{8KqRR47jJ>~0!YaI)mf{HFaMH}w-XXC`0H);*C?X%skkCDjA}s&!{# zpyvDhXjS)Wb~D2hc}*&+z~>VS-pr(6Wevlxd$_&(dGw|t5F(OQUGM0uTx$D>UF@#0 zca?xj<L$ZN6Th{h;nBm$KIy~A_rcCdp;3A1*r#Mxz9)~ygbV+oqVE0m)N1Ps!j$;H zvNV>b!)8d2=l;Ecw^?D65;_Z*ur6iEBl7f!Wn<2^_A<edhLBNdVIe?_EV`e=y^Pyg zQo09x#0mT=HZ>gSb+ez%c6Kn3$(pe1GtG8p#2^$dzYHT<yHzO<0Wg`%9(!t(*kI$E z=UFiCd#&_y!7;NyH}2!MMbX(yoz`UF$`(n2dKpMrbDgb7P$GS+eVjigj-DMDXMU3_ zFdb;SRVAzg$eP`Ug+C_?FZq6MBgi)iDyqL{jXM7d+t!Zzo>6XcNldoU(S2g5)lsi- zb?UI0myav9TG_^gWilv>>Rav;wMb-3NT!|kIb{P^3-xgxaDgVq*yiO@wMe4;Vn|Tw zcfT=DcBu|at>14-)g(f^@Jf|IKIGN**xz`-x(lM)Aw-DpN@4t`bmHLZ{*5`rVI%bF zX;~H#F^2YDtDye<4dZ2kCZvoKl3vp@yuJBy_wloUexkrE=A-e}a>utiJ*!L^M>50! z{_I{*?GDOtt}rhKf$$nCQ!1iQdcMI&rzZQ!Fh9*`&gLsyyVC)tv#ZIujjFUV8$rPI zwn$oWpAA-Q8YqddZ-vL6if0)*5DI<Cf~=@{^D=UHA71?nj8X;8oO<<u3}T({WLTXM zxi0L)MAN8SB6q^H>8%`)udhqo>$rx*67(@Weg)KNxr6fOQWl_(BzJq{@9<s+rSJZz zs~G;!Z92~mY5slU`=qGif>8i%&>v*2gTyLK^LIVqLg_kfTQk}#YCwh8@>*~1^xUPI zOY?97IZ^V@qgAxqnyk&Ly*EzyZ>a2t$ov0>#?>C5-n;p2rgaFME5^J73y+Vd%>H7Z z&e3nNq$X9Ozu)Ms@=aoLB{}cWqNX!NoyT0;8;rzmsPEI+Q&{jP*Lma+A*}9=i#Ul7 zq<i?!IAj8ezLG*At|c|7Sy5zsOL+&=#|AW*(0hZEC<BU#>-xYCJz^dD+^vT#b38V@ zT$Stf`@vJbuiDaqazxPaZi}ot`h*OHxm^s0%D;SNB<R4^t6X*V(c_ZSPtig0gffRx zLoW%&meudMTM1*q!^`xNPO>;#f@&Kk?)H$+t}&kj>rI5&=-=OUbFufx0of$S)g5Ym z`P+=^m@led1n<u<(z-6@l*D24M0R*LCCK%g>jVEQy=*k7>5QdUV(;-@uN6j>f~wGB zclRKf55jA1)4uu?YCfY$j*mwI<14^4L!azP$=Y)tMRCg|$hxff7mFvni=zd^q=f-; zfoHe&aQy5}j`;geuI<K?J=y2NYbtcfkEM%`@?%~k27x$k`t!1wB6_wiv$-!X?=5jU zX2g~zCz3!{gCTVW6k39#Wz+-A<hB)d9`50sISRiiMSS4qv@C=)F0vR=R+wl8wJxk* z=^V%WPY}m;Ov@8xwazTvZM4i90is}f%vavb14BE3%hXvC35E>rD2uhmXfa?W2G%#! z!`*VIYI5HV5rs5#swK+JT_uIB_Mwd%M=lVqm!DKbT{=wXd;K#`5Bc;zY`^d&`8*|S zKXW`fqy-1K^lN3~pilT-E(x~NDaSI%ZW$M1Du4`?{-ecAlqs9y$=eMX640|aYEhB% zJMTay_Sx414AFF&-Gp(SZ(;#El3c~4%c?{15c~6uRv34!6A0V^YSkI6e4C`NoLMvy zD3Rh3oNG_|E>o|AlP@eTFxb*wH(%tS<d>f&xN(wj6XTyYr(V>r1U63!%wm|B(_09z z$y)i=k80+^ryhrK65erYIOoO3C!0sKx?EhgAXHRbrH~==!}O=bGPf9DuT6ig;4*(e z$uVSv^Z-|M?#mZY7MN=66&+wVv=3LhGyLO}qljE3vy+_WewxGFZam>C3>rnMxxXHh zp9CKbDRwoG;z<f{dFtO=`*3x|FyyOyb#dZ(ct!FJAdNuzF&KM|W7Ih{9s_2&TjW-u zT>|1!W~<x{w@guWTnr@zCXAIKx$Kp$Wub<0a|D^7#bh3ZS|roGnsEy50;D1pQ+j&R zk*wXXZgZ3Fm)wn@BH!mK1}1n#@bch28;ifSC@>Cm{R<my!lE8;v-dz><>k7HNr&c# z3A2pwY92x@2kANW`lCq`{+Hc4JDk!So?yv;Q{aEBaKAf8IdkAA*<X0I0-ju8Pl$aF z<(8ALlxsM+BI(Uh-%A#)iDaTW8~8I`|4t^)VPwL<I;YE`HCf2d@~yy*RVHkvzqY%e zagxNWflYy*1_{(JE`sd7kPRza9={y-iC@N|>jBa9JwJaWu4cWa6rNxT4OmedfA2LM z5xUIlFk$`#p`=O|`&H`K!}>46rP+bsN7&AytkS!8b`uKP#End^;v(43j}jos5y+S} zU**J6jZWsZ=nYwvpYYuwUG!^!-8nmsMGfnTnhe$Jg!GGmpx>$n-e>BXRwmlGGR67t zGen3g$JnlW%vE-N1-~1bmQ`#T=MFnp_@5ohPROQxiglIiusf_;x{rvPdQsjNTlS1s z1#iwD2W`MCIBaFhMn$c-2o;N{hwk3*?v{2=5K3abl>gE@V;PQ~HtjFpd2O&ljZ1z8 z`7Zh&g_rw33h!3naxsQ4UXfS7ndQjru8{*-ZCv<Dtu{AI9QaUZvjDyrHR?c4FgfZ3 zE4RQQL$scw(mW8Fj#~Nm|MH9f)qa`G{7UVo7Y&&U1dt=V4pL#wyGx1U@bh=<DH1-P zkqaMq*UgQg#;vnl>JdW}e>nZ^M?OE&xaZ#|D^+R+A{J0ZZ6{?4aZuk+49Ud9JX7&q zQ$8kg@yDfHDLRT$Z~p}ufnHWsx0^pS@-XT3N*Tz0HAJo$wfP4gyoIU~8a*zZ%7o?K zwSg^T4-4CFCSUg>!p4e3GUObK@<f$@xSJ`J*{m(e9%Rk68(aX}F^o8QdMIo(B2T1f zn;^SHQ;P>KzUwd-(vio*o9-a`n42PmU1$1<a|zv@Av%{}s0}YO=F*GYw&~%t4m4lq zMXqTZDoUK7WPrPX5VF?>_{7vwBd?26Wrw1*o{gvwA!MyztuygNRV*A}e!4<BB0?d9 zGsa=~-SMtN8E~{3p++sNxm>Y|lEeA`|>FhgFA4i9PJxqis6lCr`YV%3wvc0{>8 z=>@hJMgdIBXqQsstz#1|BoT(*>}~~{>$P5Oe5>m`FcBV^?W6S2t0o?R1&n3H<|O5Z zeKgL!wx9RahZvaOmXGjiS4E!>U~Bxrd1!)C=V+?4wX{y)!$?wIr#KROZ>bO<W<4>E z&(>M-3g>*^-5n*lbGY%GRZdC0S+TzVnh(HO7ETgs@-Y~wuo6&s?xmNBU*qcJznHIt z9^g=+8=r&p<n$+~CWs%+N|A%L;FKO<4>q^o`9Ol2r^4BWXnnS>BU02&2v&^Rc>?Tq z&&NgHdfkAcCH|6$I}P`99;SfYIObF7%ahp{cqd4X_6Zodwh!a)iOA(M!in+$<%m_A zQ47is_N?ULB2u(Wuuhaj$ayXV0Lq8Yo)j1=K4lOt9ZoDZcR_|o!OQF@(L`3!8-f=2 zW*2aO648E*ynNJy4{mTV26qNmn2XhF$h*8#W3D4o0ROPU5|W!kSB*i2<ofgK-}11u ze+yBjT71dIwCkP$V)HQkB=y!x*$lX2)$(6pgZfHxULf?XXysLKR{*Iz`I(n6P*UtU z?kaeo7SLpiHSK>t+yr~zCkkA<rf>+lPusx+lfeT(u=6l3=gbB@sS$@6F`{XG+2L3E zG*WYHYri-^h`8(!!o|YKU!{ZrbbA7`VNSHM9ySe=)^0+&&W3WS>>Cb>22d<p-mw#L zYo?r9FH{1?yR#A2N;RE9?Ram|2~x1KKj;wbL5p^q``+D0PQ-}D8K)EFrCRTOU`?{Q zf})Hs@j7Yj2X3PZB7IcjU~6gqjKXT`-u-41I+T@<<IokzF9C6E^8WYj&xuZd>}v_G zam&Lk-$9OIrXS?Eofp6Lj>G5*8{Ff&xF|7h7nW}*8#r#CD(~n#XMG$<h}Js;MPCl< z0<<j<3IdC`Ypnm7D<@UkcXL4u$y^08F13px)cIb!8GIK<D?auJ_JZq#Ruw+Um8d)q z(D(4`QEv5iirsJHGhKnyob@2Kf->@6w|q-wdp&Op9T)zTG=xfj=8@++X%PkLyhNV6 zm3gjR7grl;S-9Y}l!QL1aAyw_;Rk%1-^l<#$1|r{uJv@8kOTRt_Tq`5JmJ-l%c$y8 zM<5}8I7a#_M_wmM2m7wT`<ic1!S*}o91#%%hVVZEw)qMM-TnWzieNp?`uuRzew%>6 z#hUUDze+$yzK;e5ajJAdwttASaAMr+Mm?b&W4Quz5hL>ew~s_VDw2Q`i{{q5Gx>_$ zHGeN|?Zh)^E6aMdg8*4-Gfl8n_<^Wi5;a(Y;q%$WuoQB>+E;fp2_W0<8BY@Lr;k$j zq@|tU1kWt!s{M*^K*~QAm~mQL!%_tbgtD=G)nz}pfvx$ef|l_pNH`mQ9&;Nf(9}*i z{0~>?=Qw!+BIk>40ivtWP46JL*PA4o@G+8vP}|*Ej(5aSi<ub0z-0SOQkq86-bv~z zXGjnO6~M*@wd1dL$#9eo$j!XH<zmB+nBw%4XZrefZ!&uT?H8<0-|mtwdjV?wh4ZAf zVxJf|#z7^MXt!U#>A}qIqM{g#0;lh{DNLwU+t{wB5g^6USf5|jU7xwdj<xw9t(P9E zGy#XPmQ8HWt@RdC@GZA(V?80+-QE)!!%PYz=D$U_)W_w3(51P5uEsC+K~ihE3u->I z8F%&GL5I@;H+#!!^%C20zc;Jhv?4<NdEd~3+gikOMaKmnRL8LItxIX*kLnxF0-4yP z>EUd-B6ron=AusT&DJQRa841-c1^Xj=U!~mZ?h1hur^*5zIv^z-kJVcwMma3AFe&d zJy!+s{!DDaztk}I)$SDDkz)m>YQND*?y~VwEnLWbmiec=KpWh)@OBx@THxf2A(X!- znh)1=f4I4K)fKo3DoYgJB}Hx5I@AuY%q7of0#EJ~b7HAn0q!OhNEdDwpf1nwk^uN2 z?zf6=+5KKhT*q33t=Cca6CqME?Dsyrut9p}b1f4fL>>-vcL5RE4^v@B0H#_c!L=(z zT_4wVm)3j#oBf4qp3i=8*&3{=xg>yWz9dM%P?e5TSzyzZZ<Xs076*@Ov40)CG14XS z){wPI4nCb3at;6*H>2()jnWKJVou^Iu(aHD-0cEh%piXMO=T_xE)(n_1#Fh34j#hb zJ7<Jtej1?t8J5x?K!|-YxbyJ!n65l2G%R4$t=$O!P!fO@{3@A^j<3_O57tK~V~GQ^ zdYqIItmb!bAodx1Lg5Wxgn{K!Br2dkv8>2N33H<9-E(54^KcEYB9pru1>?GQ%&-_I zVYf+~uFb0}GbIoc&;7A7FbZdi?$1TbyO7RN8mV6e0&Wr)qivvJ5K;syUjtl>8r6AI zswZ`6^--I8$N~0e^Rnky{_`Bwt^7~eCkN3PYDxhvGeaK0_DP0Xmm4r0S?v)9bl^8u zEjbU5!ZGc}D#CG<(=Nd6AnCO|mAYBZ#$H3n#SZc7)I|nRfH8e|Dcw)%^~^Cbgg_Fy z{4IZcE%8pKF)7<o%~_n)Ao$+^LPjVZQ5nHYu$ae|2c<ZHnUP_}IstQVcD-_VyYf8n zWM7*<Y7hz{y5VP3X9rA6s@{X-vPD}GXtmjsb-!HrC^*#A-qRE?u;qQl9?`W`gaM+C zUkwEr=fFRqkNY<qW_69Xpv&oCwePg8_2#)3#UDBkj>Nc(Z<oF(me`6irSE6I*AriA z7FE1fBcz&4l+hC@Ozd9y9A<l7RD+!Ed^Iai7xk?oRk9?OOTRx(n%a^vwTm!K*d;4E z!n13zlhY&kJ>4B8#gUu3?t?p3k?R|{;hnvc*~4{9#c8Jm`E?O1{$jqk9?Tjs`QxTs zj`GXtns(6QGd4nNK|>GRD_dUWX0K;m>oQx~*pRicwB7U(r}q%D(h!?|+k1m#$j-fV z8U4lYRU_q))rF-%@622KeLF7JcUdEP_yage&kCM5{H{P7&Dfep>iFd}u+sx+c9_u+ zk%9^Dpk1;KxPKzIt%H)DXb#pGH}2rU*<Rs%z8N$mTwcDZg}nf=xz@#%-`#42-5ukG zQ1WxQ%4;4ggk>^pOs5u;b|jo2)>f1k!U_+B_~k17SCSj^cvlrXEY=o(ZY|Y@ykUn> zeh3*~dQ&H(DQmeZ?j4YyuY=xw$`VrdA;p~KF!-yDDq1x{VaP7G^$|x->&J=EAAJ); z^!bcd!;>BIxm9aFwuZb;x2eZvn;EV2{n3k`Yx#~f-#0&y5|HD6R9o<XvtY8GJ2wq4 zq-0158LvahkzH^3@Yp)W^Fy-?Tqd)4Rp6&&;8Ateg+a&6Fbf0FTl?Xd7Sncr>wM0p z!q29LZYD;n6rZ)KD{v***PauPd#&$&;9}%Be;>OycKgx**Xp2ko9hXCp2$UI4?+IG zk3B5r<^}ztZFW|oU9ufKzu_bh(WL~l<>?GjJ?XPs#h!7~?e9NY#@Fz@!6fO$oVJ|Y zY!)0AeW}cM6(Tz0duxWL%_}Bj|CDWS8ao0z6K|auERWeb)^^{I^`3ipWrKR46ceZx zKO;J&FDz06MM_^S!gNo5w-1}zd}ZDqvm{Y&gWOg{Cm7DYub0nit+<U`JW*8<gJ?BB z<U%`47Vf@3nY7Gb7dFOrh(yIy_aBE?uzs)$Ty=_|XKeMbL4MYgl0Uc`syN#owNs{K zV(PMN1XX;fkr~wH)Yp~1vb4Tqvb3LF%;Ng(YS*T(ZvI?rGwps-@zEjl!+Ci9dxalv zAN|h9*H9*$m?g+q1hw9cxAx~b%x4PP`hl2!xCcr$e_H!-H|U*TXU4<$8UsnRC<~H7 zi%A`gUsX54<}MofBAcO<5R)@}2!+Rvbq)1fM@u6m!jS#)9E@@L{xnGuhWwKGGjAUl z<+Ew|_>BvC*d)ivlWQz{W*o;^#X)!PV0#DjP9Y|$Xlf~Y>Di-a4L2-ibcRbYXRpjA zPPRH{?s~uE#i_<`7`95!e|P)#s)a-jB6>b#pUm@}B3oGc{+DV-t1D@kX_i>cqg>-L zwrkAWl86q&wok>>kbQXn=Gv^Ch(o4+9#VM#X$a*SKeHd)JUNuhKisK~wyNX~>^)+O zPv&Qw_IuTv8i_Z6a*aSvyleva5jx$F(?{lY$B4~kMh*jhTwd|~I{}0aKe|vt-oMqI z8sb+`kC9#0VWsAmgBf5WB&knGqd)`7!`ArN)A1tA^>pUe5&G7+V$xzuuIY$kQbL=B zh{BDm6J<tKbpM$>bu=;KbZKPqrWxXU+W6s~`w5McR&?>v`2KM2;aM8+dk9uEH!y6h z@mX=pT|H=k#FwWKX+1xUQ+9{}bk^B77(H?sEUGrPzg3i4#sHex7Z;i$av(YomN}T! zoC`-27hsx`OccbqMTakQ3O<?T4=nIwnfmf1L`Ak6thtpSA5J+xW_z&8TEhP$O%H<= zCsm`X!QMvF7x41s*BXW_p+$$MNyxMD8%DC|I`1dw_1h(wGy2%uV-cM4xfeG@J!kn5 zoflW1a@Y)C*g0;y0hZ-Qu4l+{pD*0cOF-@;1OBXGFh?#A*-;UOb~4%Hy$EPZ6}@ah zhtZdBX5{!0HzNzFXFtQT7E`?3tQosb%mmiNsMv`?baHq~S`HNXIE-!;JCH>b<Rl=S zhPvG?sNrX~82WeW*c<reB8!9U44{e)%&l6L_#{hb=!}<5UfaGAlnk=p1H%M5i{RXT z?cMaceVttt$>UhR>@6S{{%;+1+!HHx`?%(}^UOPKv;XI}>2Kd<(8%LEZvFB5)vLry zn{$jDs$d}egmG*vUm4la*#U!7`THm<t{}+cd+i70zlE~mHu_63MDZI*s_0z<s00Os z(g0gdPpI+ltG7hAe2JKB9BS#G;fqg(VP;st#<;H0CU_92Gx_Yi_nKuq?SGyQYJ};B z+9KipM~fSLzy9qdC+_40W1!#GUd#0A*ncN|EpRqlP1j@3?LeO|?<JfoPb5dFzR?Iu z^q)o9zlcYDkJ+ZRdaM**Xqu-(5X4JS#jA1My)Jg?JX7ia(=kLp>%EjQmYW}O7DCx& zg`Fx@n|cUM@zG~2Tl#N{!hyLylHEa?MiX-+|4t<d64s{_^3^J$nNcx3nHt{ue#Gqc z-&rO>^uHAsa%JhDNin7yBF0d6n)sTsG4~#ynsk=i*6e?NXtGE1?5*ZS=yK;1Uo})U zq2tHI*g^LqB9&w`DKC~#CZprHjqocuqwPhy`XtjdzBkh2*;8XV(IH=#(Fg`d^#QBz zpPv-+I#M$4`xKq|R;WF|^X4OY>7mRUrVyEzEt5zoFL#7cDh0G8a_B1mt#&avV}6pP z)D0E#+Y*`=G!Me)s-m)?wbX>BiG?1hk}tiZ-fQ}@^GH#TuFSY;>qaAkCM@#ofSn&U zVD0)7PRe9{TumumKNtLT&}`Y*&z(leXWV$TYy)Ab<U{wJtJAb?@INa+2hp<Y(#1-% zShjF0SNvUm*~(%*r`X#^J3@?Bm)eKX_xR;5TNjI7P}7y(JP=|*s-soOASOkaaQ`?Z z2q%RT6#w;j>k4FC%K*w0PaB(xkzc;vJS=fMtoc$wW63wZCIyBO1m_<YS~)FcNxyzZ z+1^U|p2BVPE<K~wvtrUY9-LrBV_s68WmZ-x26raDMh#88MITayos@HRTb8voOm^>y z_~-F-lXc$Y>l$O|lfQ#Fhz%_{KYSo%@Xv#9v2OVfJEWN*<C=d5l!T$vLAMz+{!a<y zEH*e2J__ci`e%Al#g{g1JUbha=l*@9RLr7^ZumPwe(<qI0fF)oKyb#_B>kC=3XoO~ zMgwRX&_4cYnTYVAz`~15z+oLz?85U$SPg)@Q<5`sEP_+<e_EfNR9F_8rTC}0^Ow-l zXRcg3djC&b;oo9I8;<vKz#jkKiKOHf=x;wM!hHX`8t>>88;cL;L_;k25mJ%GRR&NM zVGZ}ws~E?BJEAR~z1)>nG<Wvjj!4n+3zHre@k0OBq0w;9EOzIwf&J62*!k%!-mZi2 zP@{hqRVKC|B-9x7zYjvm{%HWr&j4zmf;PMk8CU=Rwr1hr?2_Y$VfcXx_U{T1HZoOW zpsFxh(fzYdFtFdzp5>S818NT#B1wYM^V_d(hF=xLtPqpqJrakM3{c|<Js1@`EG!nO z$C$ceGX92^)*<i$#_?@b#odKi4Q#BNL`wGAffDz`H){{1x~57Wy0OHE<14xpI?h@e z<fL?7fea_V!f6sImc!&p&qfCtYK&jzXVtKy@q>xsezF4f+sV|bPA==xFa3_(qDX5~ zVrnCxTXPc!B`$lOVWiosa*=I`Quo-v{h0D8In%Vsl+{}RU5JG~-pQgCkSjz~+5N6) z;Z*Vg-{)S1JjA3I93%X47o)%`1b#h35HttEX$H1-@U;*m5&P-Sw*Eh>@sz$*x_n3J z$ES+$ejD7UcI9hSGN~RkG9yus5>v0ne<F)hI^)FR`<82ldJEgh*%^y9%<4^rDPQy7 zp!k6&9pX6N6}pj<eR0r^6-{-eggsvG7Fceo3h*AT_z*$(UeTTrz-ftK6z_aIaj9ux zLI%L+<4PN~5FzBgsz$y%-&A7^r`ogX-@eg&${6C{(shlGGEu8hs|nwzlQ?MTW)orc zmcNC54c3WSe_&~DdY`8w?quVA0DZ3O_BS1Rw+}m1cAL`;)8w}&k{(cUXnIoGed*Y~ zNUw*rh%T0HcDW)~som?XFM9H0>H9|Ht2NEHQM4sZQDToD1|)?SZ<^!YZY_;yrd{hl zWA_6=y2$rQ`Sq_)VY)mUk+;o2gG`#<l=J8g@k{h8^da*;I%U;@iWa<nNAjs|VHn&e z(UjtiNctke@E&bEP%t-j<mm|gz`Xl-<WgORaDbBb$9WyMsF;KM0W89OZe;S98-eny znauY$CFQQz-&OxV>jc9*{YuQ(S`gsnj(&|4v+FrFSgUFudhV=-Ud_I3K`K2v=9RV| zb(bUFTK|ob&y+*ikVK-1{^=qc*Cn`g@aC)NhG?<<1><X6_bc~KFsnK1S1K-m_Dn-Y zgL>JNIdzlqU%vT7r-o@LKVTJN_J=gz92<YMe`^<p+3i}%r9J7BX`b)il#-B0yt!Sh z?CaRs^w}DnvvhHtqk{WSJKjeXJ;TU=pAt^f7w{u?beB@}yj+k9OI0~*mnsB75f2Uh z6WAc@)Z?1dpPNfz?gh52EdIsP)m0N4sVl1JRVpp$e1lD)%>B@JxQ3T{7Iq8fd_!5A z@ksqdu?fSLkgNA1iW!W+MxnMs6tCxc-$Fx?mpx8^G5<N(T>tySEsoz|8h2WiJ}+Oe S2LBkB=r%-8qeksP<o^SnM)3jw diff --git a/reports/tables/Meta_analysis_imputation_performance.csv b/reports/tables/Meta_analysis_imputation_performance.csv index 7f0ced4..4f0f361 100644 --- a/reports/tables/Meta_analysis_imputation_performance.csv +++ b/reports/tables/Meta_analysis_imputation_performance.csv @@ -1,7 +1,7 @@ "","correlation","Max_L1_error","Mean_L1_error" -"meta_Z_imputed",0.995683217500374,1.22529552037174,0.0562121582752216 -"meta_Z_imputed_R2weight",0.990440000978169,1.59426643449452,0.454110616671281 -"meta_Z_imputed_high_samp",0.989021574001197,3.62404793904812,1.07203953928316 -"meta_Z_imputed_valid",0.996024141348819,1.22529552037174,0.05602862121342 -"meta_Z_imputed_valid_R2weight",0.996024141348819,1.22529552037174,0.05602862121342 -"meta_Z_imputed_valid_high_samp",0.995289648693123,3.36583792058164,1.0635503576084 +"meta_Z_imputed",0.997715664577698,2.10996809407546,0.0675865516980447 +"meta_Z_imputed_R2weight",0.995827954783326,5.18573176225036,0.430368355995472 +"meta_Z_imputed_high_samp",0.995382502246784,14.0988026398761,1.05536810725203 +"meta_Z_imputed_valid",0.998172604572373,2.10996809407546,0.0660996762793101 +"meta_Z_imputed_valid_R2weight",0.998172604572373,2.10996809407546,0.0660996762793099 +"meta_Z_imputed_valid_high_samp",0.997893748220566,14.0588918701583,1.09929561392649 diff --git a/reports/tables/random_masking_real_data_error.csv b/reports/tables/random_masking_real_data_error.csv index 35281fd..bd07474 100644 --- a/reports/tables/random_masking_real_data_error.csv +++ b/reports/tables/random_masking_real_data_error.csv @@ -356,3 +356,37 @@ "355","BoSCO_EUR_Z",0.985476354338373,0.0956705476111556,0.807414380574432,26 "356","Amsterdam_UMC_COVID_study_group_EUR_Z",0.968965449672844,0.128762104108945,0.798285013504903,26 "357","meta_Z",0.976898019792069,0.113298175037217,1.0469342369506,26 +"358","HOSTAGE_EUR_Z",0.991328644870294,0.141343612759392,0.160127095950146,27 +"359","GHS_Freeze_145_EUR_Z",0.99766767575791,0.0749037098892356,0.228020291817981,27 +"360","UKBB_EUR_Z",0.997094759585677,0.0515461004095718,0.160799166711033,27 +"361","GENCOVID_EUR_Z",0.909612515827131,0.136975421410257,0.435332338970176,27 +"362","BelCovid_EUR_Z",0.996923577061899,0.0584577391022102,0.0744064871095851,27 +"363","23ANDME_EUR_Z",0.998858805300257,0.0639462153774435,0.100111024149858,27 +"364","SweCovid_EUR_Z",0.984107078760595,0.122312852774997,0.288453118204419,27 +"365","BQC19_EUR_Z",0.991219808057431,0.0947491992326761,0.351113880427233,27 +"366","idipaz24genetics_EUR_Z",0.978086957636132,0.149899608436772,0.441550907259962,27 +"367","ANCESTRY_Freeze_Four_EUR_Z",0.986168329804986,0.0925800491471221,0.316496190672388,27 +"368","EstBB_EUR_Z",0.987555385124586,0.0286432287004613,0.0394833486854393,27 +"369","Generation_Scotland_EUR_Z",0.984925142802746,0.0716150287562398,0.142388724433706,27 +"370","DECODE_EUR_Z",0.999718201250732,0.0512618437741581,0.0795751464786553,27 +"371","MVP_EUR_Z",0.999129826194626,0.0443966500674167,0.130271603592105,27 +"372","BoSCO_EUR_Z",0.998222783557483,0.0359005277769725,0.0521612337125028,27 +"373","Amsterdam_UMC_COVID_study_group_EUR_Z",0.911884325479305,0.0970165938865812,0.21003961291722,27 +"374","meta_Z",0.996714601848827,0.054356164774863,0.116558759040066,27 +"375","HOSTAGE_EUR_Z",0.996746629790268,0.103698167140976,0.193617866467096,28 +"376","GHS_Freeze_145_EUR_Z",0.986706077331538,0.192115033844255,0.419428737445626,28 +"377","UKBB_EUR_Z",0.999728018960211,0.0208145007407018,0.0297093426862298,28 +"378","GENCOVID_EUR_Z",0.999505981199744,0.0673015356139549,0.108717238573582,28 +"379","BelCovid_EUR_Z",0.999706986330761,0.0776637753057837,0.117534253598823,28 +"380","23ANDME_EUR_Z",0.996021892955152,0.103335926547659,0.275330429767461,28 +"381","SweCovid_EUR_Z",0.997928662876803,0.0362859206230599,0.104399152792451,28 +"382","BQC19_EUR_Z",0.996417131989554,0.0504240096598093,0.106532510727776,28 +"383","idipaz24genetics_EUR_Z",0.991619880579089,0.064726198305787,0.159112901032828,28 +"384","ANCESTRY_Freeze_Four_EUR_Z",0.984404222546479,0.0682552263548015,0.101730447077309,28 +"385","EstBB_EUR_Z",0.999916106361675,0.0159565059447764,0.029138947868637,28 +"386","Generation_Scotland_EUR_Z",0.993150020263495,0.0692212245736399,0.166621314818111,28 +"387","DECODE_EUR_Z",0.999905560185049,0.0384715385041048,0.0738413955579982,28 +"388","MVP_EUR_Z",0.999939041634757,0.0395354830763295,0.0699736515063329,28 +"389","BoSCO_EUR_Z",0.953628275945103,0.127214268108066,0.380480757814329,28 +"390","Amsterdam_UMC_COVID_study_group_EUR_Z",0.98314126972152,0.186044808227531,0.385032944988245,28 +"391","meta_Z",0.993455976919658,0.10198685402926,0.23546573006364,28 diff --git a/src/visualization/Draw_Imputation_quality.R b/src/visualization/Draw_Imputation_quality.R index a671969..793c823 100644 --- a/src/visualization/Draw_Imputation_quality.R +++ b/src/visualization/Draw_Imputation_quality.R @@ -4,7 +4,8 @@ library(cowplot) setwd("/mnt/zeus/GGS/PROJECT_imputation_covidhg/hgcovid_imputation/src/visualization") - +cohort="" +tag="one_causal" for(cohort in c("", "_small_cohort")){ N_eff = fread(paste0("../../data/external/meta_data/N_effective",cohort,".csv")) N_eff = as.data.frame(N_eff) @@ -13,27 +14,36 @@ for(cohort in c("", "_small_cohort")){ for( tag in c("null","one_causal", 'two_opposite', 'two_causal')){ imp_file = paste0("../../data/processed/Simulated/Imputed/Imputed",cohort,"_",tag,".csv") + Zscores_file = paste0("../../data/processed/Simulated/Zscores/Zscore",cohort,"_",tag,".csv") Imputation = fread(imp_file) + Zscores = fread(Zscores_file) + Imputation$V1 = as.character(Imputation$V1) + Zscores$V1 = as.character(Zscores$V1) + setkey(Zscores, V1) + setkey(Imputation, V1) + ID = intersect(Zscores$V1, Imputation$V1) + Zscores = Zscores[ID,] + - cor_ref = data.frame(correlation = cor(Imputation[,which(grepl("Z_",names(Imputation))), with=FALSE])[,1]) + cor_ref = data.frame(correlation = diag(cor(Imputation[ID,grep("Z_",names(Zscores), value=TRUE), with=FALSE], + Zscores[Imputation$V1, grep("Z_",names(Zscores), value=TRUE),with=FALSE]))) cor_ref["N"] = N_eff[row.names(cor_ref), "N_effective"] - head(cor_ref) - p_cor = ggplot(cor_ref, aes(x=N, y=correlation)) + geom_line() + geom_point() + theme_minimal() + p_cor = ggplot(cor_ref, aes(x=N, y=correlation)) + geom_line() + geom_point() + theme_minimal() +ylim(c(0,1)) + p_cor beta_hat = melt(Imputation[,c(1,which(grepl("Beta_",names(Imputation)))), with=FALSE], id.vars=1) psig = ggplot(beta_hat, aes(x=V1, y=abs(value), group=variable)) + geom_line(alpha=0.2) + geom_line(data=beta_hat[variable=="Beta_ref",], color="midnightblue",lwd=1.1) psig =psig + theme_minimal() + xlab("snp") + ylab("beta") - beta_hat_prec = melt(Imputation[,c("V1","Beta_ref",paste0("Beta_",N_eff[N_eff$N_effective > 50,"study"])), with=FALSE], id.vars=1) + beta_hat_prec = melt(Imputation[,c("V1","Beta_ref",paste0("Beta_",N_eff[N_eff$study!="","study"])), with=FALSE], id.vars=1) psig_prec = ggplot(beta_hat_prec, aes(x=V1, y=abs(value), group=variable)) + geom_line(alpha=0.2) + geom_line(data=beta_hat_prec[variable=="Beta_ref",], color="midnightblue",lwd=1.1) psig_prec = psig_prec + theme_minimal() + xlab("snp") + ylab("beta") Beta_scatter = melt(Imputation[,c(1,which(grepl("Beta_",names(Imputation)))), with=FALSE], id.vars=c(1,2)) ggplot(Beta_scatter, aes(x=Beta_ref, y=value))+geom_point() - - Beta_scatter_prec = melt(Imputation[,c("V1","Beta_ref",paste0("Beta_",N_eff[N_eff$N_effective > 50,"study"])), with=FALSE], id.vars=c(1,2)) + Beta_scatter_prec = melt(Imputation[,c("V1","Beta_ref",paste0("Beta_",N_eff[N_eff$study!="","study"])), with=FALSE], id.vars=c(1,2)) Beta_scatter_prec = ggplot(Beta_scatter_prec, aes(x=Beta_ref, y=value, color=variable))+geom_point() + scale_colour_hue() + theme_minimal()+theme(legend.pos="none") + xlab('Beta') + ylab("imputed Beta") panel = plot_grid(p_cor, psig_prec, Beta_scatter_prec, labels=c("A", "B", "C"), nrow=1) diff --git a/src/visualization/Draw_loci_real_data.R b/src/visualization/Draw_loci_real_data.R index c803a0e..987ff0b 100644 --- a/src/visualization/Draw_loci_real_data.R +++ b/src/visualization/Draw_loci_real_data.R @@ -4,7 +4,7 @@ library(ggplot2) setwd("/pasteur/zeus/projets/p02/GGS_WKD/PROJECT_imputation_covidhg/hgcovid_imputation/src/visualization/") loci_id = 2 -masking_type = "global_masking" +masking_type = "random_masking" sample_size = fread('/pasteur/zeus/projets/p02/GGS_WKD/PROJECT_imputation_covidhg/hgcovid_imputation/data/external/meta_data/N_effective.csv') correlation_list = list() diff --git a/src/visualization/Imputation_strategy_real_data.R b/src/visualization/Imputation_strategy_real_data.R index 1fca100..4714332 100644 --- a/src/visualization/Imputation_strategy_real_data.R +++ b/src/visualization/Imputation_strategy_real_data.R @@ -31,8 +31,6 @@ z_col_var =grep("_Z_Var_imputed$", names(loci_all), value=TRUE) sapply(Z_cols, get_study) Z_score = loci_all[,Z_cols, with=FALSE] - - Z_score_var = loci_all[,z_col_var, with=FALSE] Z_score_imputed = loci_all[,Z_col_imputed, with=FALSE] @@ -44,7 +42,6 @@ compute_meta_analysis <- function(x, min_samp_size=0){ study = get_study(nm) if((!is.na(x[nm]) & !is.na(sample_size[study, N_effective]))) { - if((sample_size[study, N_effective>min_samp_size])){ meta_denominator = meta_denominator + x[nm] * sample_size[study, N_effective^0.5] meta_numerator = meta_numerator + sample_size[study, N_effective] -- GitLab