From 15152248a448defd590b98fc1f49df0e92d1cd0b Mon Sep 17 00:00:00 2001
From: piyalkarum Piyal
Karunarathne
- Last updated: 29 April,
+
Last updated: 26 July,
2024
Source: vignettes/rCNV.Rmd
@@ -111,7 +111,12 @@
starting with raw (unfiltered) VCF as in the number 1 of the workflow
chart. However, if one believes that they have filtered VCFs for the
parameters we highlight (see Fig.1 and sub-sections 1.2 and 1.3), they
-can start from the number 2 in the workflow shown above.
sig.hets()
and dupGet()
use use a
+sample based inbreeding coefficient (Fis) to improve the
+accuracy of the deviant detection. This Fis is calculated using
+the h.zygosity()
function; see below for more
+information.
# Start by installing the package if you haven't already done so.
# To install the CRAN version
@@ -703,7 +708,7 @@ 2.2 Deviants detection
# Run this code for a demonstration of the detection
-dvs<-dupGet(alleleINF,test = c("z.05","chi.05"))
+dvs<-dupGet(alleleINF,test = c("z.05","chi.05"),Fis=0.11)
# z score and chi-square values given p=0.05 is used here because the data is RADseq generated and probe-biase is neglegible
head(dvs)
-deviants<-dupGet(alleleINF,test = c("z.all","chi.all"),plot=TRUE,verbose = TRUE)
+# in the new version of the package dupGet function requires Fis value for a more accurate detection of the deviants
+# this is obtained by using the h.zygosity() function as below
+hz<-h.zygosity(parrot,verbose = FALSE) #parrot is the filtered vcf file
+Fis<-mean(hz$Fis,na.rm = TRUE)
+deviants<-dupGet(alleleINF,Fis=Fis,test = c("z.all","chi.all"),plot=TRUE,verbose = TRUE)
head(deviants)
#> CHROM POS ID NHet propHet medRatio NHomRef NHomAlt propHomAlt
#> 7452 un 705339 67808 9 0.08653846 0.5000000 95 0 0.00000000
diff --git a/docs/exVcf.vcf.gz b/docs/exVcf.vcf.gz
new file mode 100644
index 0000000000000000000000000000000000000000..75f5f563eed84432682d45b391e8d5fa1c64d7b4
GIT binary patch
literal 807892
zcmV(}Di`EpWAxmS^X8qv@?r;Yj_&abTQnG$(EyZG-^&U6?Uw+HK{kMPr5C8bf-~aQ!
z{FlG`^LPKt55N2WhL`!D{w@CFZ~yMkzkKJ63nshN{_S7?`M>?kfBfaUfAhz;fBE6B
z{}=o6uWF3>H|>$X{OjNS&3FHgzx?fYfBmQTzy5Fk<$wLNKm8y7^_TDd@}K_v
z%m4Pvzr6qUmw)|_|MCz2_0RwC&;Rt@zv)kT`{%!Z`_te3m;d^Yzx?al|N9^Q{GZ?2
z=l$Kk`G3l{`~G+QcCY?+zx&63{<}Y~Z}$B^{^KwI_{-bx{^ghd&wu;npZ??jf$#9U
zzxm@|eg6}_$B#d4U+GVO{+ECFr~hqzlOKQir+@x)dvJg8Z~yQ$f6pI(`Mdx0_S;|n
z>!1Jg{`&8K`!%npKS6zmAAkSxPr%>)?f3uizyI=&|L~vx`Jew2yv+~);g5gwSD*gy
zH-Gx{m%sh=$G`a9r|*CN)2F}vpTGb8r{DeM?|=G>KYoJW`BVR_Pv8Ia^B=zd*`%MX
z`}D(q!2j-kHvZFZ{{sHIHSq7h`%C=qxjpphcYlQc?(EOT&i`)nr@#Iie2)Lwxle!m
zpYgxf_O|rjU4rNTslR9Ytnxqo6}JZ{V#v`%in(a!|y+te{a6z
z_xAt4IQ!+*{;zp=_T8n5kJ$JQ|M4X_dMvp2(M5@D>(c
zRqGe?f42V}e&V0M|M`b6xrBG~ZljC%Ccgb|?LWR$ev{fB^pQTH!Y9-(CD4QIwbt0M
zzOwQw>$?r5zimkG@E_}gDW(3l2Hw_uiTt)PycgR)GQb}d{CuG`9$x82u%u#b3zi=~
z{nUcRdzZ_r{W}Lsv=A+Ab%LWv32UUZhn#^Ran6(^f+WAPZ)=S~(vCJYNJ8Wwq5r}k
z1rdTFp+rb{)$)jt^oj-|#8HH#lGX_6kzot0fBE7Z#D#m0IpO-+wh-xyH}E@!+Me|8
zi|@+X+DDXj1>l=NJoqF(!nJ5|b&nQ$bcvYTABl$~hx76UD*lZn$Ri_%7DYe7
z{o-?oF3BirM`4em(AMa5X7nQ?+9ZQ0wA-Xb#x;nXH+*FH8AQP+4x*$5QM;$25=4ei
zi+qWpd24w3-TZq?7ANr#nE`~25jz=T>lI(@U_B^_6D!7^7`;BBP-bjWYPB}w7ds-`r$Ki5SEp-LBFe0QN8k>jQCh2Cq%;gAQdp9YcAINg
zi-jPu@XNt3`0j1rlMc3h`^_?3Ic+cJNDCO>EBt%=A)RTjeUG+p@c`F(Upm^}kunsq
z`RSIRQ$g(PzN>At%ZIm|M&vW?pWq&JPzKptV=W2MUu{yI9I
z`ZbuG;WZB}WP8F1ImBCvp}t}nYu*;4U+gyY1-Ip~;4N1e3m)3HUlu$X&)H^8hpVO1
zg4;3^iI$_@6N!)>tL@%eQ}m#j^%diT3a9|0T2JXuJ?PQAXfkN!?sI
zXY(^7z$oBjBe(zBwm&@L@wjm&s*zT0>HiL&(K4Pr2p->+lFf14!AG}rB9@3pe|x)j
z_Qz28?0*Y2v~)4ev)}p9*V@*8_S>BFvp=<8Ml3S?bIT*!BK!7}KX-L|yo^ol(##ee
zW(~Ndy`DbAN46+V{?OudTi(~dhcmwYCV1AH_w~nW8CUzobGrfhj|Dpq_i_8zFqRv(
zSAhxys}FDKhlTwM=Q5;bp}Z$rAW@ZwiliO3@9qmOC#C?|;XfuwXmBI8q61ey9KQ%W
zJGB!W4;hQGV&z(<7PYT-P=60=*%FQX3SuqWRZ45wv%!}6>T4N;ciUPTT+5ah#AkWf
z-qDJ+G-53yt)yWNnG5$^Ul=&8GfzyA8_2#rQd$|@Yp!?So?p`
ziah;*efTOt_e8dx5V}tg??~ta5w5m}TZ}BVW|+*cJ+QT-446XP2e6wUu~nsog8K
zK+%@kaKR6E^SGz)1%G+P=RFNal`Qz7B)O&V^XR3;0STmP$VJ25!fFOR>xzXQFDPwt)x_+WJd~z@x49mCde@
z{gMFgf3pac-~#7#+4_>o5-Ke%@4-?Nr}1VB7Rb~^53#gY$Y5z(pxn*{)3OT=7F+!Y
zmsV`t!o|mze}+p-;@5CVmvCuY6>}_46`v1qTzaaG94B&E|D*^hmEp%yJ7wGA_sYqn
zmhPAFV9Vxa`EPyO^pMvb9yM`&6mOeCkD40aOWhu%2(>0x2vUBqMSK|HV@paceEfLS
z;$awi=P=at4k5*VjE8n8MP-c%;?l?fnP*8W269eYfEW%CI!){$6|^l6sY__2;Nmqj
z+-MxGp;0U)5LK)19GrH~OT&&{Qn~k%Y^&HTW(DwzxPq4`hlRykpL(#l-6%2a3*>2S
z1D?!>g@*+a3%^BFs}yBKEeqd%a=GxWL}>pnW-9UXd(88Nm3o@^!U(5^osYI_rI&!c
zu@ImTR3FQ&5QN8OP$ZKA1_e}tEfw%Bofoa4!rkPeSAruXT66^w684p%!OoCo#nBhC=6b1cjyM`aE#InCj6s3e3II7x$fh6o7Lp2B8$Nv7
zZ9p<@_R-mAu5X&jTwGWvzV%((Er
zP+XL#9Te$S*$wBH1C1($IZXsn>#+9gpn*Rqr^DQ$0%AY*kQ
ziBe*Icr{cBRQ1Dszy@BWKxv8N^EnnLA~yC`bxWMocd~GGv|87T-++T3P9ypY0vw^w
z7e|ACLqAGd_WNd2n6nC~Pufu;VusTXlkBOEaOVT|}jWE%C#0*qQJ!slyyB;Qo@=!}UTKSf^F)?Zfxoa)Am$U~xwH?$4U7!cy
z9P(BlWQsYKcw0L8L@}d`K$T;W`%#bchnt@>CbCO?@&M;Zs?}JTb5ut$k>e>Qyvx0G
zaV?9KF;RoH)MIo-O*xY1f%Z%|OitIBsIQongd=!9NP@e?#AqXpLAJMFQL|`w!6+!;
zzJ{Ae)C_jchytu^5G%V(kjAIV0k)%#2}~~_9m#wWvR?@1PdL0CM$fBFJRXUpF~mNDXnFQ=<>TI$zD%32QN?9ycH*+&Nt4&w((Pcd
zO5*7}j~%=}UKZ!2Jpb@?knpy3;oZ?J(6=jJl!viTAosTP11vN8X!KF7TBSB8^`XTr
z#%KTY5y$Wr<%?fF2j@<)0$!y$Kc*3HZx;5MoBJKs9J)uYRsOW@X`el{f})sT{vEIW
z26Qq|wO(h9K141>Y4ztPLIi|YiytKda{H!ASIav&eF{N^vG6{^Jw!zXC!LV<#XdtC
zIH-<+s%~9|zBs;!=i0LNtmOyMBoOYFG@nRAOP+_>j5XrBPuze`kWm&2ty+|~oy*y-
z!#P}X)z@|}&*QuB9!}}hgoai!T<_jRxr9@Oo!SSy`3Oq-L75xR
zCvXCn^4$M1(>)$EU{Slr6)>WGl`2nh>+iFP?Z@@ywx9~Ql_FrYqeP1)5Grus6~xFg
z#BEr;RT=yE5thbd->8otAvxW!5aSf%VsLMTP!ZDqi22X1@Aa5ruUq&Z90Bb{u5L-<
zIhAWE0*$2~AVm&=BzuUT=RSAsPOeG}_~ORML0keYO1WF4VU}G|lew|Wt9xAHI>ta=
zi_hP--wy4bq(-c#wT>adO0Uywkp;+qDnUyO0t3$A0>{q^Ho>Po-?JQ`Be~A?bI>W`
zW?!qqY-MY+*D<{O;W9Ulvoc=B`H=({Ib~%}1Q_(3
zEvAUQDg(6BSHK0eym_=y1*{-ZZwzkB#P}{KciQiWr(}f2?>mq61Qd&J_YTKeTQmFt
z6$U=xFih;|QOYUQQRkRM?^%O{r*rAAPf%}dogOi*u8V`PU~@HIp7}ZbX4F&|F6Ob}
zy+u@V2FRCovHMrZDiEcPK~fP2^pq*-3n~t!eV&r_*~DX{TI(5mP)ALE*sr5j>vdX9
z%D;M3)~4ZxKN?qoS&w=hBvhd>jbY~yg6bDY(OwgbDa1r8!Ahp~_1|2q#c~zi>PX{O
zfI;w_r{!bUHhav`vh{VD00o;Jcj&ak+c=Z>@D8&BV{PTd)LQ*e8YnZtad-nlT*Lp(
z(a)N-3qWvyW>))!%$^C4?TN^9YFTfwm-LRB`jrmGe`uZ
zU>J;^9yNMP*!y`a^C4i3FjqPmTltg^wS|mZK;R7e4DJB>g!W9rl`bt;+jeb;JnOUC
zu65&@j&9~^!Tam=7f1X6Jmg=k%9>4F#;XGWeCL=n`Kp8Xq>xPOytDuaeL@
z{jVs&iPT24G!`=X07Ub}WIZHKeME1^8DRxO(K=2jywwkh`Zy$Deo%Ol5WljCv?AP{
zSrZ45G+0=ZcBAYyqI@bNiPu*j7&6@Uh&=&`zW068yYC|)L=m~IavdF{=wbD8J4N2!
za;wY7>UmJcki?V{xF1lDEaezGcW{S!#E=Md9PI)6!~3}i=ubAcDNw5+TQYKii#1-(%oiMTHAY>pHW{
z0zma0#&P&S##t0r839`aF>Y{0c&iZJ+XcD&JM##0&-j7%_U22NrDI!OTe9Sk0c?x5
z3;~DDlWY>O-&0ap#(Xd53ZX<}c=iZz_lC(paAJ2R?(k~6(?UN|!|`MW4!d77$hi=L
zR1{z$9*F}G-DnfS&l&rH3%S61ASRO7kE1Av#2C`vKE}dtZ!rQN@9%QJpa`5?1fm<-
z^=q`qM6CFEDv|Q#u_PiT`iUEjEty+M2r+V`FrN5S0|Zm!!ZdfW609ReuX9S{XeSB2
zp&h)WsvZTo<(H3}y|#-+HW}jHAY57
z{9NnkEwcUt&g-J
zx7X4ISa9RC(y`1AK)a26T$+qbED!mV2t?Nx*}}9?-W0`^JNw8^p@}-zI52L!9>&4x
zF$sjesKLP#cLDCx;G0Ad13LxkqwdN>>VW21kV9XGl0$xU!HWj7Y1S+RxEL;bpYUVc
zHq5-kGo%s5)6Air`{I&mEA8;$X(=xj7WLfOGS!U`2_}`PU*vqvR!)lp9|k@hD2{9c2PN?D8ufZ(HGGY|au#$fHsjnOrLK-OCr
z{OsChDf%&7J|aZgj~!O_Dl`UeM}iqgg5^@bfH=hvH9H*r4j=>!CAKwv2_K@j
zU~Yg8F}oAnDP?8}`H9cL-~sh7>)@W7C0NqtfM>!TSD$CuKbevb{iTQ39Zl
zl15&bO0IF?6IFtUcv3C1;lL+Kv3hwenDlD5%w1i^`vCBnB=OYyS2^VoY_%qVQxFIslMJx&Aw%SWl`^G;J(p@9h^g*d4TRO1rh0AeV>s
zb&tN8iTc3BgNGm}xgYDT3`_X71=EDK>NJwTh+aRr#iA#fp~*3V_qVp65Ar*S95{%$
z$SF-RyI4>rV0Ho%B@8a~+rXT7W+5$c;o2h7J#Y*QX9twP`!~@sWDJml5+JeT
z1dCE%&fL%+4~VY^_A`RN)e6XMX3{U#bPVI=Mh`Ab`!Rq1;pgvX(+IkO?hA@y$Y=c2
z?tch@q>i;BNo5&mE{ry=tq>3cda`}jn3ldGzJR)k>YrxpjDeP4#aSJ+S|Kjo;1anT
zObRvi)D2!yqD$+q1+vx*c7ucJ1WCvRCJ6;qq9(nJCTo%gSKS}agu}thx-t`~@dP4i
z9e`2Ui&-H4o=CcIj%I?8A@^@O~Ua8b0mio-KRHeTLLr(2)@=JThNv<
zX_j%w;Sdz<0FtMOg=t_8NI$JW_E4{qlsPO~emz@HwZl;J<0=6N%fDkrLOR2lAl3X8
zydrJ*xIl+xI!Ix(-_p!YI1`7oR8VRviT8?8P8Nt!*lzQlYFMlw9lcojK%Y+^XmSsM
z?6BW}1Y|rJTD{zk@um3?Smu)GC?I%oAgyomSWudgOpd}ffPyBv(R-L9Bp3!!lUl{3
zq;}I>mcC_C-_&%Z_K>ofyrbTs>w_KKORnuD5Jc|kxw$+c1mh_KLD1i&ie~nFPQ?5m
z@eB!w?~%kTYZo%(`-E?Se$@c}}J^5q$caFJwwTAV$ff&i^{2Fi%>
zZlE9Ea?9nD9Bg=2MDe7r1OPbwCNa(}-Ks>0xJyfNq#nM`K7^l;ICbi1J^(-N9f%bS
z*U6c(wy@MuI-Lt}0!N}0U{aTn_U|gEh`3XS;wdWDQK#BMEZNb*auG!j=v}AEZ+~H?
z7gB^Ip_-8lcCRTgj!GOL*ATu;jTBS`>zE5jPhvX|>Cdwjwb+7)gXpn5?KJgJy|=Oq
zv6U~6OqYj$^etO#bUd3WOCZ!3NxC3Ly&W~A3tSK?iByHSS5Z??r|^^jqGDw}5$}2E
zA)nCFw6Kj87I;`m=RI3i!Ng~>u>=up*R~Xn1fL<)Zu7P1Sz$2=fs6454tkL-R?F#F
zWZJBjwG`w>SZ|<*<)JH1xrMdKB+gyYOPDD@M9SfG?gi3lI88xBI1D(#?%H4nG0OTN80y!j<4>D`9|?6!O^#QnGK&
z(zoX6Z5|)mkx@Hz3Ky(vVaq3knLWUXgH01PiRw$=N|TmGgd%v1#kXoB`q{oTEMSrD
z7ew-P)4>U~_LRaWRURo8UU+dmm|{qspd1u(pJhfTUVT-CacEp#PXulAI2vKKUp#$D
z{XDR!Ta&gTLJ?w8-!UC^2dFQpWj&|fgFBa_9re(DFarsJiK01`+y9k7XF#eDc8r=+
zPuPR)Rtw%mU7aF8uTjG!`NDvCpXaUMQ!gkY)C?-fSyMrRFY#3u!b4f2@7w;*dyS!Y
z8h3D6Hok8}?gsTIrENrU{DMAaEAK^b*vedxR^9c|0kLOfTTGf3udK|7k{2t?K9+y)EO#M35r
zYQO>uscfp^Y1~S>Ja%q@^U31)(;PHX%hK{#Ik8qNS&H+(ss#=S9vWr|hBA+gkOMY;
zL@;n7d6{m;IH#K?JVdv{bqa1x95Epd4jUK`msGe^Q8EQcCKgt7$f~x5WW=RM9|(_c
zzlrnCp<|>6cG-K59sYNG`!qaLIBo=5ZO>1;8h)(Hzq61xLt%D-?^w8@N;;ch_mClz
zT^2c}GL>smivCrnU6b`J%4<&`hrZ!m@qX<;de8mGP!9b^$JK%rniuvCWO3~o3il;1
z#X{M$!Ocoz;ms||Qj+(#-woFek=zj{SEc?#!srX7+s_|}>RY3`T@Ze(>l@Y>&xw<}=G26LS`#0$>E5L(10#JCRO!0{G9
zeuw!2F{tZ1V9nCwX>jtP&;(g#i3s6iU5^q*%jl6JfHt0~HA9ni`BSYAB+j~p&V=Oo
zF_#^_FY8;BOa?GAmQo%UL~c*83qZ=Y6OlxC3Fxlf5IzZI6cA-zqi
zN_VTi+I#szZ`>9QFw^7VY`cj5;M%N0z79B?#`2022uEfuz&hYBl?93lRKwwOU_(O$axD9#~~6A_ACy}DMyGnAJ-6Z$y10_bu8MF
z>8aN2gj%z^3Tn3UBj3bx+WB~okTi!DwgBUq7t+4>@i+gZsAf{e@N@l1gKyNuP3hU0ahN~@y@;l%ONf;i+|S|V
zgl9<$g
z?-9!33L2lU%6Mpm@JJM45Ylb`aZ=%8;IKrH(qE<0ykl5GBo~3l&!Losk;_6InUU8!
zS-a(=!au@=bvF9)?VA2Tqlv?$2#uzE&sQ@z3tomC=E}-dE4c};g1csYfDV-yRADcE}6A@rt
zoDpMza5H^njEf5q&of)sA5@$HK9{5tP+mnrqG69w-EqWBCbs}(;gIAi3vZwGE65JX
z3McKLtWC3b0bcG@tddl-H;5+!H~=%B9cQIPU%zGw9X}C=Px%M}_B2g7s(kPl$ncU3
z8{X*OxU#VK?lR_(uVca24%D7p8cuK^%~j=cM?N43f2T0JcmNtVXPIPZ?gijFEe#zf
zLlh2h1Fo+McL8tjW++wCMSas?yNbA!EV(PCeFP_GxG?;K#Ia-j}&RlXQe!bW$JJCkAoW5qrL(&AkzE&i6o-!M_C2G&GAEMW>8
zQl6La*fy+T^p)(8(6@c*P<%~mCHr=qUlaaxl;(MWMcd^k{Bg$oAl4KK)h_1;OFwoz
zhx@U63vYQ+N)`Kg&jCTKp<-Wm8xRZ!mF0ah44|zx(HKwhZH3)@3Y9m@Q`)KC>m23V
z@Vm5pW5{;w|)&!%&%D~90$f|7!lRQsXI_(3Z}WLSWOh9W+>UTU??9!-Y{6*Oin|)N35E4E+Jw=Ua@2%D(+%OPZwmliJWrYsn1%ra<*)GSSfR)3Z&?`Z
zTA+z!9!Z@jD>>ymECbu~{My(sH{S+ihIWiN%e=O=#D-&XT@8L?2_Axkd0v}%qzIA>
zfa?(?X~R!x!RVt3*?X#=Ce4zff
zl>`EKw2Ht<#zP|Gu8b+{bozC^6J9)?6^_NE6|O}u_W|wdC;kPNg3@FJ(52(Ck(pN^
zJ;BDucS|tt7HQ{4JRFP8FQE~un;e~1$)r~liqg9|@i8)L+)7aZGIXI)qjrxjF$l?9QiWOG%KB0lcT&Ek)@Nng
zwPYVIv&++*g?dz4%0_A>ukQxJMLF^mEx;e(>7HCY
zh$QBic-L|kd008e%4vnB`6o?OXX|oISzssL5kW5U$jy0Sl~mYN!5td-&td-OOl99t
zXU|tcS^@EiwyI=Vj#CtZ`VXd;!WY&ZH6jHI)XkZs^WouwjkQWs8xVWXLY_96Zqr3+
zkgAX5>80S0Nitf9L+SE*Y}|Bl;(Q)0z9jM5;m5J{dt`iMKjyM62;9m+a*nI5aV0CGOOHo)gj
zJX=CiYi2{2M|{yaC4wXXk9TNmr4$wmX3KOXo7@LdOW8jNEBBJT~z79J}su%w)9Rv4Fr%^+-fqsTk=-K(G<4dH-p$@
zF>V$L#5x@0Npl}Q7xhiMKc)oZYxdyD~Yrq=)8u%O_so-tij*0Y`
z1m4wz^^ZMwI(RYY;-^B2eByh;{Ug9=hR=zjFw7gvPme1Jfr(|4M
zeKdznwj~3#qe7as>T26pb((O21>5AOK7RJ3TjABp*nysb!WM}OTjb`;a%`Gt{r0|V
z7loPjgGHJPcUg1Z!f`f1@JGc|&EXYqtzCEDY>E}9ey4Dn0WcguZK%|bK1RupYtICM
z_R>lR;5@@2j(yeS+@gtPnRaLylVcCjV1e8&VV&)0$mO#8Kf}#1`9o#SUO>3FZ~vkK
zJAwOroV~bqjNfsc+V2jo`{NeoJ^C&H`i*`m&*bi
z(yfzyjzah7(N=|-r`ST!2$ygbYB9^J8c7LCY
zvK=H(z@Z`hH*E=Z0r`FfM;*HAIEXPD1J%*>(egJP#TQ0zwaWXWlCt6=x&0Y#At#g-
zpYC3gtPGX{-*9LU80@f-K3t1?|63DxuVJlFaGhF8;(AjkqEJd$BTE(3I!c;wT05h^
z9`&-is3HY~bFLeD;#tx1-gfebY01kNUakU8jO6y$tb_aZ?)KhwN<+0^!F!JyQN2B!
z%2or5$bLiwf9!OAy_(|1jraP3hf0c`E7)b^e$?MwvIujsWYV=Y`T~an=q_gj1Jwmr
zUO~*)80twX`}XO(&I=srO)b-%p&8%oaD88whuG~d
zQUQ57*K=p(vu|2uD-p-bYwS_hELgGF+1NQ}w;759AyQRXQ#~Km4x4kdZmv(e(2YViK_77m8noi5{hx+j5zNHH1)L!|l-2+3*5@tQC+#(7w
zC4?1_BfpBz9o-G2_c!OW%LAzYR-J`?GXXh#riTwiQlou>u}a75fsTpR*{OWROC(NY
zWN=hVMEgpbb!WR3hlQGfWc631eiwZQ?FXi+?^48RuiDYtkIz!PFsOummUrGRQb%2B
zxh<$kKmYL4_dk9A{SQBW|9vl1Q3oxbj}i!o_CE}iiQdmk9(;5bQw`^7x(Q0AcWENC
zST40Ipw;gAelU3J$H<`0b|1wTdB;-j$e+>?VTs__dajW%|1F5Gc^DFFty)8&(IHU==3BdGQ
zEV(1kj{-aI9CLyamlMJ!*+v4ql@bp*1&)U(WXggccs+I6=vDrJoqV3NOOw6_9j$k(
z$SRz}qkKv^bbJ+>VfBg}T0}5qPGb&GyL5iBq%?`x|4HSWZL=BIaC7w@0^euwmh$(;9cP#p~7F+v4iXLtU+8*
zCT9l^^@C7hm``j~C;pv59z1A6GQeHFZS7afQ${BOo0X6o*wKl2*gS!>3^pQC2umZv
zRh~kHLnFeD3gfAO64A~@4YaECr(79h%A%5u78cRwf3W5RtSpqb+HRt`RNa|D`r-<8
zBD_|O#Fi69cy}A!ydlQBPh(@cVmByd>)%mqq!I6I;S(_2mX73Z>&QAaMK
zzjSOvC%cvvtt`A~@x#9kP5>Vssa+XHHkVZU9;YDp+Fi0yNOo7}(-3JR*KKv794!sC
zDi}yqp^a
zhfeS*7rzdw{7y)nmjl&QORrr1JvT+KEz@ZC)SZ4n1ZbnjGsI7iR{FtD08KFXwt0
zRwnMC+GAf?SR}UM#(x_I$PG<;=p$|*fFof7sN2dK`!X;<>Kbnq6`|~&QlX5{oKEsx
zE|yST)}x9OreFzSSw|1V0Vdk5O!?#=pCM#N2?6f#BY<>n9KeXTdXDbAe-E_r|fnw0d_}wT@0X
zjLZ91xZ$0BWLTwz4{zRwC$Adi{#C8Cu-ND()P7QvdP;?m+f387+jw0yyH2!q&wzlH4a}b&
z4T^`P;t@W&K|?fiV2$xWZps`ciCEa*qfA3|3we$WjljHCzX%o`&sJ0UOWfOsviNIM{*?Ck;8AR%10{pt<
zZn;2-JH?$7cX%Z5BXFP$Q1JR1p_(ByW#~050kAl!>(Y<|6&zS+aMFT^62$ELhtso1
zgs6i?1qf1`4*2@9ci@A2wVQDk{wzRF!W>~=^7`^&EO~7??%@K-Ni(Og(@)n5kdNi;
zFsl;tPBo+Pi6yDSnMr%Prv2U{p(A!$Bu|Mp?AVbZkUdO~d>n(5-qYdN;S?kfBlN|(
z3RKNo%6A5ean7&1DzBb4T|b)cG4ig_O3xBRG}P{&ir{f?UCw5aY2k9A5M{tTVWA)l
zj&=hsWIc~Ol?u_9(dILD)`QgzC0JFUJSfBuzy~ai`(!>KR}0ru_z?)~ox15Mkzj*Z
zkC+^+$N*Sjc0s3pX#tQ%Kk2~2SG9fHJzjH_tP@
zb+Jfsu_>80UQ3oM6e$&E#9djsXPvDK0=lQ15i^E_r*=<;FdsUyyd+69DD8|7G#6;!
zOT8tc=Yk^&C{Y+;q8=LU?UYmk6o;q}Y%YS2y=}nOr0Prz%f`djL})E(0N^n00wDx#Ef*7BWd-_Iy|NQu+&s=UE_RQmi_OBIl_6myq-~aUe4?p3)Nu3NBmj#cUfwkPN)UNW0go61p3)k=BM)a(6n2VUm027B2cEOz~jf}(li6!?ArC56JIJc6|OJY&soe9mjB)nOD
z!U2Zk&6X-n3nO4f&8!bR*9WBGeP#j2x=(}d?o%eEp#&^l0k|^I(cz_~riyrX2#YuR
z-5TJ50EoAd%^BIPV{|W4d+FY()K3_J|G>^byHCJTjODv*Jj-{Tl(||^i2$hZ_Y4Ba
zTX3x^n*o4%o_!@>CqnD#SyWPI5w_AdEwg`9Zc%_mgJak@kA>%_-8Q^eh^3(1=-9C|
zS8c;c|%o9^-uUrDfa;{gy!B5Bn-hZ|T(@qmK
zto5G{Eb}@v1*`ytRVEedw6LTUf1f?6!42f5SB?ZI5h7&&aNBmmq#voZA?^3sjVI@`zy(Lj&MQ`c~)zuVaHRSd`8358n5kQ)EoOA
z(A%h6VrbTQA5qL(ZJC(JD6whk#;2TfSipkJCWI9K5GFD8i^Q~7uVE4y&K=cRWcw4w
zEYb@InHdbTL!78~rEZmG%pcgBx19!uE;iF#xBnS~8$;p(z5R<46I6BU;_9?$?TYeJ
z!KWQ}c*fXaWCkVwh^0(!#nX{4iD0t9W6uA+juA4iWQOjhOY7D26bwJd!1A_o(|}3w
z3=n7=Rv|r*6X+}prV$?_hG2##ztd`VW06d@0#->g;PJ^a64!=faG5?>Qp<^ndZG@~
z-Q@T{EXGbo%3)O)Tc*OZ4
zEq-X{iWC>fLj;~`$!BkoASx-yj&TLbM!ids*pfAzLd8x|1!T)X1N?b502)#*zcoh%
zL%^x@pL)lsJQl|3Xp~ywLE)7o>
z7JG+a2LGbvo#FC&XQZ*J4IGz1h-S#ryd~>HJNFeUxpY4Tf>}}t>t607Ela#|I1yi<;LgOK5pqNoVA_5J`$2-6|Ar%c$3~8+#|sPIvECu{g|mZw5CeqYeGTQ
z;Qk*H3Ky;+g+f7vA~q_j-Y#529wdbrPXh<&2ckb*t>sKM%+G4#Ge{5_4u>c}$Iyr-
z(6LM!0N)U&>PS!;$)7XeS4B}#&R({ZAVeYF9$H!MJ2I>bLYB%49HBj+JPjrly{*(W
zY#lcf`}0uD{rew&`sw>0e*6K{cu_+zJ1${5`+QLt>~JTL{)XgbD&juh&4->nrulYd
zeF1&VY_OtX>qHPrQ2!eIcA+}XfSeW8YC%$LbX0zktWr(en}~yOEz-8Hgp>HRUQoq$
z%|l3Ev(Ggzi0CI;_9UEY%7UT<+uc;ztgxt|Xp0|g#04qKP6Uf@7D#Z0IzkcZ8Ku5m
znJ&p)Eh|K?dq9^G>6WQ>7UdwGxzx1|)U~D~a}R6qA3JqA099k`>7C(GYO-lkNs&Qy
zP$r9VDfDLO;OVO>OQ+@w=&eR4@2_m*07^O6LTbg*Wr`DDs$S-_tHVm5R5mA?w0k85
zbo&`puQHcALFnW7+RL5tN)jvqHqOXv*wA4^$+9c+rCs`We7`En
zaZ-@g
z{Rkhp29>fX78yS6w7$Eo1tCh-evOiaZHRZo&;XvKbI+%_fTuV2UBQDyds?UCWK#NT
zXORb#EWx`1%KKTmPujk37vUuEx3n{L-A`iRz{rF;tLuj(j7`h$7SHGxVjqBvofpT{Q~l?2+7hgE1N#kGjYMLfIn?-bOJBbK
z9v2QmjRij|WIIzbVSTC(T?Cf=*jVVjaHqV$#Z-F$M1TsWFyAQJh?onDP6EyoHwf@p
z6=WgRPXquBd}Q?`4hP22U$WPq9WhwjTlXgi{A1$bk8uf~Yz>`r9u@KH^9%QPZJ9~o
zqdQ<^d*w}Q??3q|hhdfX3&{Yp@;Ta0XPcK(I3hmKo6&NX|=+bv`$E=
zp5&32YdM@P+)(c4BrLftoz@Ysa$+4*UZ@D|vFg$k6b4zFDsv!}1pWv7Yp`S_ybZ
zcW@%qsNAZSRVuy4iJ^C^ZD8y@KFh%go=Byx?Z!kpbsQcw4O4@+5L!kj!y_>eAs%fe
zNy)ahO|p-iEX^Mj*B@w64Nrq4ItD|=Ruwh|=7OSA=LUnV?8MZ~@PbU4sb5E2Zzhpt
z_F!jb8_b`j{OA*(LL%e~+JwB`6&%WH@2E|)_sLXrMogy07iMN1&nZ^Rb2r4RTM`{G
z46MLT_AjR!DLkND!
zhjt+1ng_XVSz+sD_%>0R{RdpA0bo_ZLUtsj7Vs+DgihOx+
znp0RWnQicMcUPu@ZUfv)!A{Qm+(%S)2
z(EtaEatM{gExAXyCbHF0MchA}g{fVUO)pxce@%&$q0#PNAOgmTz5GBVro>$p#~{Xp
z9hJ)Zrm*<2Y6k=H+9u*DK3u66S7~1IRL6**{rov><`rfGqCd-HYQ;tHj@|&KR}UlW
z!pH@^7lY)1FjX&s+A3bHSXp{1=>yb((k}Nrv(wsh<>)&lg1|cr0&xno-xl-xTKjB)
zAdXKvnnPF!ggRv+CkyC4UVuLxJGs(6BE^~Hg+ZtfE=%S+0+Rg8r6K@|
z$_ybA>LAs%asZD=%ox1W$tWhUzyYI?=cj$U!EmuGV-ffM$wppVh-=R+t0r6s86m>K
z3LcWyBM%HSi&qNRK!*(Q0}LCHA_v@vD+EZc-T8>6uqV`Pz|59t7}*M0SMg$(;zx(x
znhZcdF+yhA5zCRMLlqE>&-yjJ>=VzWlqd^!~AZ|CF|;IMhAt4u~8orOdLzLf3N<
zH0EZ3JLHRyq|s2h*o~f2yzLs-+F%A==$MlvGDYL=aAv0%AmV#*sN4?<*xH4w
ze=i`HpBjd4gP@i06>>;)^R5YmW48HZe61!#EE=)>`-zvp^n+E_3IPn2MuHAI#a#|@
z#XWmy8@(rTVJiT6!2ARWPD8MgqvrusrX`NP3_*ho$m>@TTfY!O-T4pH?K1>{-1!JQ
zySiFNp2&JLw?2@VX@IZQ4bHa^xbRog-P1@z1%@Z3z5OaK_9?5AH#MW2WyU~R#@QkM
z0|%pv<_dq_6Fexvhemv1a4G60@sb^>s;z1rlQV?u0nuht9FmbkHsZ2VPpo$N>G^tCh_&QJC@aB$u5`PGUB^+C)cvukMRh+IHJwh
zBPHu5DS1UvoY6XuWCDHy)E--Mi7|`hJA^bUKMWPb<+vhflfD&Qot|WKx8S5B>cgl>
zbU>_MGkOT0Kvin6@9;H8gpe^s7@9mSXnWjR`CU3YKCz6F5d0i9%pie1|6J(DZih_u
zehBEa1G)rRicH_NlLvcv*)#>WEj)42<8#UIs8gp)6~3Q@1^Z;x>2+$-jTYk|O4U2<
zR-wK7k{JgPML}tOiweMdny)x(AX$LK;akrdJTEZ0_=eX?ljg^PA+Iwr1RouHu9Lnd
zLV`snXMo0A6WTU6%LoLa4yeYJ;^`g&r!=HX)x3^J0o_K@DnU(O^0I6q^mkW&q&3pA(-3Yc_FZ9u7$1vlp|xEQXi%}(|eWnJw}l_I0s-a
z$6NWFc5E759Y~XHt4lx<%Ni-xHBOD*4+Z(gIKe!D&42oPn9@5}HlTN01=f};K
zupW-y@!tsAqvNlLs8|gJtc;mz+;DSrAl|UlI3Lq=TgikGlVe!@zgu<9Fd4P+WrFbW
z5(i2kK&ei#3IM5iLMKdGEgEiG4Vn9hgZapEs@4bifmk$6QCg|W^b%TNV}q{@YX@1#
zW9#=f8!wT=ShQ!S0eON45>IKDq~xLW!cqhex_GFl&EWyN_vQi^kq=UK^(%(ha<9V@
zw)5=Hx{C0Cc-}{z9CO5Xi{x
zl4U{?E?Fm_jEkYJB?&;yuR=U^YO)TyZvL?Ds8kLRI8CDiJUv*%{cFH73a$0*roj#3
zG}dYkBP*&rk7GU<1i|jYxyoIA?S!$ew9*OF&@!Ohd-@w_L~`)-6Y&u0kmLcbOxaXA
z8Pc8YLGk;uhj}LFjv?8yhD3%_`xX*m81j@d`&royJ04Ae(d!m{Cgc$hmmvqNALv~T
z5Ci-OKmXV+VvpE@LiIRE_65=bW1DC8i
z2DPmmeN@
zFf#{N0mCdHPp%#THawQmAUkFo7|)R#n}%w`oQmvz?_7fBd314%+}_@D1$kA`!szaS
z*-n2~n1|~%ie08|?T&C=5GW&?Ox1}V^011(YtcRBtbe8R)Xjl&(MR^KMYVTztL}@a
zR8AhNQ@KH6ycT1^J~(slzDvAzrfn+f=d0g}ERy8cPhtTR>Vyty-*ybZFhRg@e=~m5
zLkB959CQkKm45V)r?%xQuT)hdXk({!)MNRQ(Lq<>`S?_Nv!Li8EWypL4*2_Z&-VRF
z42B`{$i`a?shbz>>fp7v8G`vh74j+08U}hNkx+s~!vW7u6jMly+#c8ibn_ugq#mT<
zE&DRJn(EEDwuNSeX-#&oI6eQYTbKFl`B(fd%hXTpR`*x6l8Am9RMF@U4
zzNfEr>3!^_M+^|WSv0P@R_jAt1-EF;VvqN?4V%J>kv4$4%P*iV7qDGh7lzmE{Vs2<
zT)^bwP-z7>%~iQ#Y*b>CYc{>Of~nPhYb_XASgXaoj%zZyi{BV@tBt>`t8kkEc)J7T
zx15|rxDTx$Rc$xiz_0#%zxm}>6xsX3y7`emb)t
zDl!P>%jzq+)pJ9Q)$~k;itm+>kdPzcZ`L8EDEpAMsg36DbWM#lIg)sk+Tw_1YcRDk
zyL}u3?N+Or9uqCce*?nr%*5N7ZqloSc_gq>tL=GL%b!m^uHJ1JhFVya|2bLwjFW(&
z(JHX=ma+2gly}{XTw_t~Ikw~4aXO5e_w@^IFA}=m&DTG^62dd%VoI!duU*l9Nvs)f
zr3ra=k3U>754+f!TBi91XeELi*{q$0rwW}O%fbnfe`?e&ZHm_6DiM_ndegO|fEZLmI?@NzS(UifQN
zeZIioPH+CVdjK%%`{!UFx|Rs472}Kn7Xd@SB!_m7`?L)33Qh1_tal&JkgY_Qlvi=L
zKL2NwR_1)$ixS#XGG-6Fx7aep8-6uio4N~S0RH$sT+Q7)yHMjd5i+z7Z_yu9{w;W~q
z@dX!Smc5iGIeLePz#fpTX;B(ivE_5}r^{&4pk!`%Z!mCF-?otR;@h^UpE$FqpGZZL
zoD(XN@Ovha>j6$9@VLk$6|045KVkyVN^TXn7=WUJRC(q>i5z{)&-agZr$&$6
zn!?mvc!;hqgce%amNkqPTnsKJztncDe3M)`g4l@(4_a_yX%kp?-sQVsv^i+6Gbv1-
zQD9c}_Kq*FDWldE#cq#n^`!@bm%bG=LpZzP0YaL_uJOj5lVF9S9~jmHzm#7Y0ce`W
z!P*L#&L*gQVEZIcTr}T}-2UGGdAw|KFXFP(l=KY^g>&=#)96*O3u9(=*l8A`;
zT5N08P&Y>GOR}RFNmEqBjcA-U?%BrAZSU{I7JGCrBR9S2Peij
z?_}t*WG?IAl`VQ4r}Yj_bK;{pxcJ4b?7LoNm=Y?Ej*tE*>i5!orkxZ(!yJzWbodCQ
za_FJ)`6qP8GQ
z+bz5f=SQw-7tZP({FW~$lD3E;<&j*DGb~c*jQ|LzwW9QSX#asbxhAca7dee^cJew=
zprAqvLb$ZELk4zmRW&oEtC}qyucWb!#+*ZmeW_4T?AY^=Od07V0={Q}w9IB8n#j1p
z-*Z?{5Wu|c@&ebDQ4U8yX>gPH06G@&fL?VJ4oT$TOp40olqx`{9n~_bTU8X=
zbwa?QsGJno8RCfWarTN`;VG?Cd&O1bz{XcvualeO6a3WfOXa2>!FK#oI>=uP2G=}+
zS@_cdD^pLKzWvzIz`f1Ul+pM#pVW=JA>dbVu|$&6AC!qo$Z#u60z>LCT!y3H8pQ3W
z8NR{eeK&6p3=c3VdTZan&52(twpXpCI<)!geSlx;?(5i!qdU%Uuv@CY4;M!jMNu6(
za*MRf9a4<@mc|naEv?eV$Ukl8gr50(&iwB{ngcuaJf6tw=`s6EF>GGtqJCtjV5P~C
z+~u>PSpVv$UhJD|Cwq_$sK}#)aFkAPrgILQTZVzwgA2D~J^rB95;i5K9aOmwb_V3BoZL#Z!sK)>
zRT{n7ZB}pqXf`utvyE(&Y3~{TG5lRh3nlDQ`zKa|ru&aCJ=f6(NNxs4pyj&@AS>~C
z&bq0`^Bm-M7dvm)$3QS6)W-<}teM3)VW=VK@Qn-=Q`c9;RM7H5HZ~zDZUBXCXBr#W2iitxSzsLEuUr
zw4gU%Q+QQj;>tlldN^Qop>O8oHX0c?Xx%P#KOmZdvC5>N&C_B8LAsO3l-AIWtfG^1
zgJ+OA;t6lHw{v*LuKjR7gyO0tj1l)kfQ05`sFqjC&w*5Q=py9ZyG_xg=ZK
zojQ^KNDNLoWE>Xy0#ELQ6~saU+rf0dD^tjHXq#M0Jg1=hiGc?h{@fZ{5k
ztaVqPVc2JwrGmnC3qV72OQ&^~$nAJ|w_v!v1LZ^V#7}OOJ3YS;ON&)=OgNv3oDroDtskxd*Zm96k
z!YVGpX{}dr&q3j5fEBH#BVX0_#U^0ee!dIf4r4o3IK!n|MuD`=Sk+|1V_Dv+ww)I)
zqs9=JzR2hDE$>uxW;WWXu&C&bPs@Ev=k{IRT=GF2&C!+EO7B^Yo0SKCJeu3vtbl3I
z;vz>xuk$q>Z)jnc2?GL};2cuzprJHU9qSi9I-8UFOA+D`90oci$8=WycK6Y}SH3OL
zit->A%W%C8wl7E-_2Zj}rOd&Kl)2r`_Lb*SI!>9{-14{~eOY`qTK(0@+71LG4!9d*1&~gX{~{B
z83~jgkhDCkk0JyD3&&&?h?gM)dv52-!+K~P-zsO$n3Swp(<{#Gpu7*&PksZmo&qa<
zw@@ZAblsWxM5dv+KuhAolpukX#4El6JMO(usyzr%LdQHm9Te8TY6%cz&vm;)ID(h%
znr96%m#Xfc4v827>YdHofp2Q!=891jKk~%)2o+O@Iry#lsxtJ;GP5u4PNCDVk<^hI
zlwJPS5_s46uL=}aCt4J#P0(Mu1d7J6b5u=V0|j)ZSj%Az6rk5bq!2`3Q?<`z3k}jH
zjwP1qO}aZvus9nuwJbqX$_N_kwh^CaJ*xJ=sbgv;=MVm`VkWBi!~rM*z8k$6aYQ>6c^UCL&9hG-y1dBgKYH#$4eDD
zGA+Q-#)~L-LQQl0sl^en476(KWI|i!c7#-nmDs%$ty^NHw5y7zJ0&EW*6EJt>5kW-
z;KIO^yyZ3pv-UCfHvZ7R^+Aq!*ip4=eIMsu;VDuD`ix7!oCT8)*Y;*#sc~a<3(Cm)
z^ABZx1CkyThQf_SM4BS&F(V?fnPB--43Wo88jHfv)D~PTS^S5
zgH_SS+g2Y%ycCxve>avA7j%hdJu?>%Aj1J1$~A5dEIzOTmBP=1qb282fJ}myOWf?&
zJ=v?m9r*0GNGai4BPFiBJ6kL5(}Dw96wa7pO%En+I3)L8DQ~?AL0unX($?~UY>s*r
zU{-J3%t9xHEQhgaN&2cCQv$otcqUX717hr?UC_=@7Ci`zSl!+udPI>U4V|&_VicIZ$YJ#
zMSn||6~F-K;T3wPUaEDemyaZRC-mO!`@w#cFo
zwh|8&T9pJs_H!H*C(`tgwLWQ}kOGgvu=^ksif|{jxm1rwsPBy-L$|ueT+lq9*K!JM
zL9xg+%9G7}Z$7aYL>iK>Z9Aab5i(`rx<_>;h#BhKu=&^iaTSY{;hpP~PztjSO;wYz
z3ktPxF4DraRMzkh@3DaM*W*KgV+epJcGgb#Hwh136f2`gGr}gO*iP(TmU#)TG;m45
zgB02tT#kdwog|ImGD3Aby8s*4xGKtYs>hWZaTP>dF>Ns)sHa}(f7!rHmv(`E6{r!n
zcCH~K$pFsa0dKT`3>B8;^J&E$RtXXi@RrrRNhPvJYYyz*dmP9e9z9F0m{f02ubV
z1Tnu$H*GIwVBG<&K4P<=duqE27PL|>6jKom38#XNZ=pU=YWPyiZpp};m9lr}IvNUN
ztFrKl7rz-Q>iF*8-b2c@!!k@OptoT1t8LZI3Mkn;)dtL=o$9Jx)dqt5ND&P(Wop1s
zI0XnNA(Lgbqe?f5Y^fea6o^26&{2)6x)`yU_4)ShAEXYm{3?_tSZ=!}L414}3cCCr
z$>_!Aj!8^uK+Ds+
z=J!NrIq+sDK-v55gPzw8pd45eu}1IutbMelJ4pPN9t;_jYKDbN%kmbtpTu*KNU&ye
zz@N0_y)LXL*DBciMz)1FauR$7oNs;?$~xTw+w(dSP>{SklUVId@$NF>d|!!%cAHjf
zhz0eR-zHzPh3&R!`K*e{h8oj69q6;3s5d6Gk3!V8Ocp#SY#*B@w}sdGEZ^<9ySynt
z7C>Msgg8W#-&)d;UtY+OFe7f{qldyI*`;?HfCb&bEG=QRPZMAqv5b5{s2mWf`Qq!7
zHNE7*u%w5y&{IvON~lqvY+ADK*#se+rym_itT;l;9hcVEY2uK(tWeshlIv(~a&KAI
z=hFgDsB)zrV1}x_x2s7krH{adC7+gd<76iZONj|f(Vym?a7Dset~WJQd7*r;!Sm2j
zmxj)~ltERqjIi-9{Pep+({n&|^jEu~j!OOd#}xUA066uLhB=4Qq(-iw>Z~#cq+=vE
zFf-fX#0@tNwiTQ@m%m+-b_x}j2pX26#Q?1%>Kjt_Zw`H@0!sUG`C)>79+Z7XVA>db
zr-@8+1YbDtTAfkdSfGmF0Jwup^&d(Zc5(;js5-j1DdCyvGitp3(2{86m2g6M1SA8q
zs0zx0x1FE9VvywOb=4!mo>)KFB4f^RMh`8FY33uIhG%+O>|4Z?m^zcSqIMVn^m&@@
zYB6>Z$Y&2XZ}eFz;EUqqHKkZqiDU1X4Fv7~g;baS6
z<9qj_J{sPpEAYpJqw`=KirPjOM$gfK6o$MHEPM8~vT8PYAkn;~YwT@xKF+D&PYL3oR
zsaBO&<}*4_(B~>t)4JI_&ex1}iGd^p4@}5j!Y6>jPI~5iDV*?)aaWHtIr3r4H%@
z(3pva0{>W=MW|W_1!fc%w*!C0&L}u#_E6JgjHm*?@jleNfSQ%7zZ(A3-GgB9SL@Y*
zC^m?rqB$=Xg3D-P8)Z#J*N?P;gqV+iY&*_GShQ5tindNm-`b{~eE|px3o3%Mty4~H
z1h<&_c>u_KR1J>Pf#YI+7C9TJH-Q7l&5@A+nT`Lbuai|2Gcj4)%2F<8ahB}(=R7z$
zZW(&qbaVs??@?H6hmV7Ib`w^2^p2#3yV93ozPD^Ti|J@Nf}{G
z@+-B*Z>^g}La%JJhMD*~?hQma{f3{ob4Z
zi^?bMI*Dv`1)odBr*dxko>*9$KDDv|uv^?%e0$4gvH6`n;zAs$L#hbKs147=mT9L{
z(fc|gLJW~C>0TB&{sY;SXa++nXuL;1lv$|7bVk{J`}jcXM?qZ&{|MwW@F|0zl6vtl
z3ASp{$|Ael8+vp3Qpym*MD^VZNu)?@&ir289?}`P{K-y=+}m66tw;;!mh|4^+cK&=}?E4rrJ<0udJE&@=f4>MeB&yFkAYy+E0RKa<+2~``#1e
z$Mh=}@o0UvEEaz}%2xT|j)Py^F3%_}dtCh7#9LpS!*g#Ffb_!KTP#)ZA(IEN&P0E3
z1CdG2%5cy{^ggVB3Q@PddpX^ix72a+OPcfO5FE@@pUjG>nhRZqE2-i!QgR=Sl(f{l
z5WzLQC&FuBfm#dB1J2E#N`P4nx;~84v^fH~pzwFm@M)K7u&-!7h>x5vc8y1+VlFfk=0lr`fp1hSpX*Lrxx}sPXPyR6V$|zJpS=Y9p|+#N8fSU-f4n
zg2fvVETy8f;`e>#kSa)|OVR%|CW9jqqf`_{2o0OWMzP72{xEr{8+Z=)_OX_Omw-UaBNK{
z5dbxt%2VhJ6#0c&l|5rn>yiBct!AheQeM=KyoGY$!j$<8=f!}OWYJMbh6#B_8$3DJ
zmRSh`JI{-(HI_Fl2B{08ZX%JT65>U&VZ=CiCgTIn04ERRE?Iez9TwYpisas@!s`e-
z3=XyLE=*^^T)$*Igj^{gV8*7$t^S8^fZ@`SpDA*t5GnyqTOCy1ZQvAP5McT&F@Y^F
zf9sHFGHEl=C16t87rX$fEXFw0aBRsf(&d6%W2inTIh_0py}>LYM!U}$Z@3fO^TUqH
zl$j^iFq~4fG{Ji{O_)w8{sWmzKeB{&P`nFEYSE95vG=1%;>?Q~W3(0a8dv2(07&}a
z!A{*dCr*Tr9GnT;X^p)qq)@%fwz0$jRP9@Y7Bl|v-3#~9xlEKYt0m3VkG4x;@*pWY
zkU5z0DX{eEmD0oqD-*`p@-u`TxxIau*&kciLxO~pOjBv?=Vy__f;xC`ig*f=C%I6_
z>L5umvN#xuQI@vLPX%Z&;$8;62<)?m2jDvHVVDp@$&zYYM4YGxz${aSH(l;dB24|?
zjqM8yQYg0lj-w0M7P*xba$7rVxC=-=WqxNpf+*HtJ$koR0AQD!-I-hc88<(pMVz25
z*~}8SPUYVxkY-SfaK;FsTzPD|-3@Jp0OICx$`vDZrN*xjy%&&>$DixUN`*JS6_`WdjZm}rAej?15nR9MCS~0jXK^Gfiyp41?hkNI
zjU~XRtIuCLdAZV+V0*cVXhO%S97X0sV{4aY3HfnadELK|O8p8%kC;Opc4xtcXk_B<
z*=E#tO=EDD)z=X;aI{`eJ9F&;4q!NEi$Q1_TSHA=t4)D`-mER;s(l2^?QhK9q2uDe
zjL|CEPk5#sw!(?s3v8#E9Ysux-ov$VY1v>A5U5cz`i1@Jctwy}kJe$DAt~;9zUjw@
zwqn3@zSbAIgghm6VyjkPZdxJ@L6t}vSN?z!_5QKYtZCjZHytA|w5tW&JFoh7kot
zWSR>rpbrp2r4Ig{tes>i6$JQ}z$f)iry;)ba>X}vNEAwWLfyw|*`Mz58Mf1hgR2!c
zq$^_rl
zv{87hH_;h$xUnn@=`#T;p3(#IQNY^@tAg!{Zq+hS>ZZ(nLPFrKgtJIhP-6~8#Q+RB
z+`vR4S{gZQIQ#s<9?C{5L&F&kjds0SFa(zMm(5gv)s#*ouUK>UrORdo{-+Wz13ILI
zgwb>n8q>90!^0ug_gQe)VO4VO(`tI$-m)I|O1|txQ^f*!*u{9EZv!GFuv9|=-htWx
z0{}ur0jv3_%%dDg*;^qWzp*vj>yWKUpPGb>QLBI>hMw8m0=-dOfyFl;^%mg+o(j&|
zD$c2(A6sP2>ak`xiHLDN;-lC}@QH1WB?1(W`_TSnLj5#BiqbW-G7!|4EIxUv(0x_K
z2k>|yX2{@h@~7Q+J`q5G7sb-^RR7SWqS{Fpi6^BGFsE>|$#LG|a{l>yF#R@2)TX2R
znd-wiGE)8L83eXfY^L{=v?_2^&{gr&Fo7-&&GYUD2}6CBgVDOfm+Gl}5rI;InNVS8
zvsvgvg*&$_iF85==e2DI-iYA?o^KfRK+(hj7y(g5oqok66%9Ip;#a?0jViig#Z+lg
z4J{rM^=K=JfI#RJUki#|blVcoYPMQ^nH~(WIghJL+f{sBZ}GMD{b|_26BwIy!}}vDN@x~qqk}#QHtyGJJ^4%QT`)|BisH-^xX|Z
za_N(q(9qCQDZg}gEFS{(kXnLTO&W|b8*5pZg>9|Mt;+*vL5?HClq}&%EmWf*y`m;6
zlD-03uCJY=Ew5rL>Ecku6CI07@x9M`wEF2{_b4AwLC+HZRJFfoh4XHof;Ovd$@;iB
zJ2Q<`1iOtt0HuIewR<^@js!v*czF(>mesa}H!@1H8wZ4iIqlLa_J#oh(cg)ut%@r$
zb#U#5aPa~eh4`puSS3awc~(c4SdEVg^rI)*d=tmp`6PQ4^623#u2$4QYAJa#`}|yE
z0^EoRj6;{R;YwXS#sWV`3lttvcxA@)EQqMz|$u4c-&$G>Rx-|kc87P??0bzTQ0q2F)k_T5>4Sbmh
zrVm7WcSM&vT)ncHx5I=N&sBtAw4V?m7+{Mh%~$cWmS?mjoh2c8p6*X+Ri$uV6Hb3E
zHSJD+-+*{!2k;0DQNpS4C$$p~D()p378&ml(R^y&BYcgJ$b?A>CLj`m$egICa1ba|
zR6AiIBA2wCg_vGmWg*sA(0ml@00f*Jmmiu`sU=KZZB9y3krNFMbg;{u&Lu*idE<4<
zsy^+S^)XbdiL+5ci0q;#-xv%Bn`s#X(M+5`tucfJ7*eQ*3F^3TBILt>A>_IfQ9m$D
zDcO!5)rbLcB19DvQh=v)oFlq+BwU4Oo<@n!zz8L`oBBr)8On2V8Jtq_l*NEKX5S5%
z^);hc0OSKMC4Aeq79H~&OTV2qfWY#km`z6t=QU~hwiqq13@AX=4isvJ(aZH!$G-X0*&pjU==8}I{B_B5b>v7
zYjtebYA!-Xpf-_+iGnY)=!}{+|D0LoB&pvoEn)N{had;)8;SLHz%&}D5nEcOSEhNA
z<&u(g3T-}_Cqespky0ilD%)d4t0~8@QGx_5x2ttGywePd$Bcjm{;C<&=8@vI-p=Bu
zqU!}go5^Cb4v4Vm21V_R7SU~06!gboqP^(iDqw>ect_{0#)K|&PoW#Fir@6KXC_BO*3tP#l>%xRnlOS$-
zKcZA%p8xr;gcRegbTZVPNMJ0y4XkeyiozODYJ
zY`#Qp2%dDNyxTLm`(@qof|1J`VaZbJm^NzJymxvkm&dfx9z{iVyl1rsxF^dG+%9Iw
z`h$dVW3S#DuYx^sFt8_5yP)ANommvNZrDCZLeej|#(($EakGbfX7@g6s$K
z!hP`}gN47|VJT;QhyC{ZI8Y^L*5{x?G)o9(Xwyv%fcdn9L~_R*9*q}7*%^2Zz*|UL
zm%>s64o2eI`-g=3WoZ)}5ywQx!yesXYcp+sZLI)w;@ncXmqOgMl$4=A>X1*Q>L=3I7RrN
zDR7G<d8Vb?M|Nz0(y+jwum9cQH}Pz
zG@|HiL*w=a=lZvlYPi)Y#gXu2LhrdDsPkc*#4zPcT+0U6%)`{rnH=*P
zB4F`J$z^@|a+H(uLSfA!A$_EfF5x@u@aRTFcBu5gEj+A5Fn+$e&yW~x8HxMieI4jDBB3mtu
zFtnNEaqRM$xgAe&URmwKI^1b}mMsPJEyPCY@0PdcmA0GWpuyDlDjWjlvs^{^w5?V~
zTo$Iz9=zWclXmVvgnqHv9q!(?e>a4(0UUv}v<9HZT21QGAUP%>au6|Z4e3?XjyXdY$9VBpL!>3BbNCE?rA*h(Zjnc_e
ziDQH)5pHW+e1?waphVbYcjqt=(!rAn@K#ENljG3C7o*U?I2jd41+5j$$33wce(+Z)
zCZP@OSu-!t;|!7bcESi%8g^#h;Xe)r`i`>S~N0*;WlO#Z25E-
z+t;TY?f2x`*L<^((nTyitiQ3Qpt~Zd@b~O`8R7z)$IVJvA*{8Q3zf1_rZo=$3avF@
zxwEDj535((-wiLT
zX*-UpyIhZ!XUD46YnUA2Ox#RAo;tdV2iM$JLC8^%%d>`Ou{YJXvpywS_mm|VeNHI!
zO+-$AT32w2
zYt^<6M&^`>d3TZFmHsHQ^z4^A>X)lS6gy>jv!AVy6xaqXsdYzsB@G%eu7>jvURlcUa#a5thMpvtSa~`rvW7
z=+R9N(6_dz%V?}K>dUJomQ!{6IUNH5yH2Rr7e2p|-qs3%7oIUNxxJYgVhK*iLbV%3
zTYBr+pU)=Me_1-Je1
zz1f3EMW;jSvp-~%r|))j)`^s+L)YakTR$C?f^ncmB06_Gt_E$&`>MABTV3^b?c1jU
zE+$pODoa>qq8#=asOW%hb#v6)XF|U}K>l>dpjKfExs)lDvgUP5-q7t9!x{j6@fDqo
zDsJxt{^2-NERR%;^h%*u|K@7Uqqcer8QSoYugH_UR1Z^{1!CYp!SH}Xi%v7Xbt!RO
zy|9_h2OZb?BqY_%`>L9|Wc6*Xa%*G_;+k|52m}}gvVavMUi^JgDm+Bkl92$59#(Ve
zW5E3hVkZ<3w@15Qb;5Sp{uv34uo5w6KqW~Sq!T?p@!!ZSy62xE#zX6mGvmg$Yt=Kb
zEe97?<;ud?rj(8DpYz#{>>C0n)+1HrP@3D3;uYbVmdcp#X#mESi9Szy_$n)vt$vW<
zQDqGY8~hN4^lZ&iWB2xI+MZQgQkPi4$Y%7gpBA;_X~m%@3M-?xr70{pxjZkS3%PCg
zB!9I?YMR@ZtI#O&Q$lmfR4&hR2B;IBi8l^(59kfo501t4bBJn2*fCM+ra@Pz46$6)
z=Yrnv_G`e*a3+xhVY7XC)4n)>ZDjn_S|`ssM;~748+=%x!jRMDoL6*Zu$4wdxvJ2p
z$+}ni{BS7hL|doy`5{xGM!&lsS!#1_)yUR#ut)1%Gh9&>KFT@>SYXw8^$6F_Pt^!$
zupr=tzHXoBI^G6QD0bwv-gYYww~QV(Rl1P|Sgp-wui(Nckl^MA<1fEk$o3$HAd7LK3Wuneu^hZ8X8Y6+>BS)tZUB
zc;Y5}UHV>)d~nd=5c9U0bGHqQB>p=K?$f^)cGqN81t8Q&tiX&I@&0|Jn<&YNSoEe+
z_T*#D%626x9h@^LLJ+e3Nm~%rQ`9u&m@j}SoUc1T7yQ2
z5JU<>S)v-~c;2{SHjT{CMK9VBT@(vOw8&!ga{Q)w!C4!Xc#mhxDy3tSHqX_^SU}Ud
zqH|m`O&KyswA?-Fy!BLs2vUY-|5A>h*~bP2Iq&ULB+oUqUtVafK3c@-G36xQ1ae1NU@^1Wj@-w-x^R%%!!h|
z(-w?MUnDG0cUaR<^6Tv~zDzrZF#>AK@!#IEx6?ewC0EPBku8P&x4mwT@VRaAI+8Vj
zL`u>W&if*Hs0FL4so7k)t>Ak;_jMwlhB-9S$n6u=-L80)2dC*x!
zL0zOzPmobzkdb~0*+aT9m%CM+sIc83U9q~BC5t`J4mgYzP#ti(+-j(9BgWgVZ{S*D
zruCjJT-FL&U@#;Py!dAi2P(6?IAz`yn{QS0qxXxIj_=F-pk2V$vXSJ|RZ_<2_ud*k
zV38y&6Lzgs^`#`_8GWnv+7=&3bX?Jn;f1~(l<$%3l^?ol)JEGrAr^8K{S0f%(W)0-%wAyYP9(g82
z7ef5EV8PEz#aCM&;-rn3Hs)>O7?;dN3KG-6;pi7%ACht0CH$>U#
zS+KnPjr!`gpW!Vfkm(z+?O0%PC$4`k04y_x_K`wpx9}V(l(qZ=!pPku!@`J?vOsYj
zAcOI|!Ua--=L#oAd_Bb5f4^sDCt}sQe}ym=Yz>fNr+$gpR_m??Ulfw)-gfl)5uut)
zr^`@khU~4|VQr6y^Fe+ARETb27VxVTq*fP8`Cbb}P)8_Js2?gm`^rM;4C{QbB}zhW
z`#~jSx<>+MlV!ajYr`IX@aotY&Id*b(R$BMqwbt0cLnL=1DW0KfU@|a)S^TEHU|<{
zcM7mW5UwLIpTNzpeYnZsp!p<;CJ)g~dO_xq@&ws$w(3)*a~%!vpmTr0>IC80S9hv$
zvl4et?jua6-U5j0&xuWGYa-=Q{mAReVPFWN2ci+1q9l$cu%SWA!-TX`T_|~?1jNY3
zuu!rJ!p}-?{&RknbztbPrJk06h4tjV454=|(WQZ}>t9j>KXL*Y(~JUoSNVyRoOy*C
z4<(U;_d)L&IQeYLT|7Zh9B_xnJQ`LS5`dpb&AZSad&qh)L_k=g9RMD$cMOq`T86O#
z2LTp9NN=7hKj8>QbnAVs>lI*fx~nLcQT|yYbd<|oLDz>@)h0N@6EW^B+_eByO>}yl
zRcFEMqgfzvdUh>+69ZTzeiHA?g6wS(`UgQwg<8%4se^92
zVDn#vF70}S82&BYB7TG#HxDP)vnbLo=mN_{
zTwxM)=_Hd@jDE?=lVVR}-Kl_dPJ}|+%tnm-b4_xR=h=(N9iU2GgAR^`0-XS9x)6N1
z*!sZT;b*CsA9u$oHPoe);?UJj`Fd@r6g@YHZ`R1kXFxyt;`bZ<v&TgM4Q)Z|b3V$s!J}bxVtab$*S2v6y_|c)
zXpq7d>jgH|)2$a0XsA95hCw9Hx_|D&p-x;Iownt-y}FEi9^4Uu;uXN+UBK}1(RS0f
zSPl@0TC^%5u+WBXY5tpc$%aW2sWzpTo7H@%t3m`IavqrQ*sXB4iBC%u?vvVJdEtZ%
zI$CIUqjYJO`8k53z~G4_K$AXm{}g2I->%hjS$}0>x=j;wC5I3m$5p(mTLry!+-YK@
zPIzt97lgWcz!dPX%za9%zOVYPmX;KJP;3RU=P`s2HM9~w3x~^fhVxAFbF0ck{+6v>
zn6@GK+`;|;MnJj0-akb@U`yImwMh=6Es9#t*>)7)XqqKxtGUj~S$0EdarA$mt`psH
zB5j#`*;0Iw|Mv$B#m0$E%k7@sx^@(15MyAM70v3r%YCD0hhlw~wP~krz;rakIv3fw
zi1c~{di&IJ<3*yOxWd@adaBt?_fWIZKD0Xm8QnS^cNd3|X_5YymD_Mum9@)?nCYX6
zi(0`xW*NzTL3Az5W0%ym-NCJ<&%-BRN?efSO0qi$viQ(pk2_@8Cj@*`zjXBPHs}oA
z{k>Px))R5@|5zTm9-C{QYX
z*@rinQSHj2*&5eUjfLjj*J%zj2{eVKyI#gIhwsjNm9*awiUA9WH0q3a@WL$5IT#pN
z4rMq%-+Hv;qmTX*0IyJqhpnty)h78om7^oG%={s*4j?YR^kt*}U=V8t``Q|Dix+Of
z6xD~K03H7VzzSE3yS?@%nDaMLwcLfn>yPBTa+|h4+CxQnEL&h?PPnuyT{?k@VQ9eI
zKgot}?0J4fEPRA2bI|v)ON)FCNdJ?9IPtg3VFa%TFd5r>w60wWssPh481S>{Vh0dM
zb)T>qOb;(WoO(SFv{CnQa2__8qyVyPgZKrb|47{CW1Joy2lxQtj9*6aY%^;u^w0=j
zFb%$Z;wGEm%5Uia02{^Xa|%x?IrT~2kI2+bwYoKV
zk005Wxm3P=02l7)*>vVqobeauPT1x8P;}f98%XliqPb*FCxS;8PdB-=OW#DjT25Y;
zfZ>Z|X<`3Q7~c94kvnElFve16rH4ZH079jMAUb2>!sBJdB(M@V&D&EZolm!!Pqbr&
z`mau?{a>;l1)34MSGZ8vQT5=l{LYi%v$e?q=SdfsC_EmM8j6gVN9Nd4uWT1-YOoaO
z$m#5L{uN<9DQ%D$G`E5KYd)Q&;Tr4H6jbB^Z=pAx!fnY)|N6n{w$c1wpGQ^stI
zEgj0L#6Zlz@w5c<{Wp2O@UJo!$!_S6H7FhZ0k9fX^&XXTqC_r8I_Cxk%_i-+J
zyJN&eW;_Rt)`ybB*+loOQ*a(hsdg#qJjzuHgwjqIpNh1@^C`A;Zx*TSv`u=T?V3Bp
zLD+;t5}0owJfteEChR4eQJ#t_#9LAms;ZWRQToG@%ww*!fP-C}DjjmegPyUzCuHZj
zG!eNZ;``bP#18w$e3KK{Z!3E8@!LeSS&zX1Xe5-B2SUh}nw><4CLaGr_rWz|a;nS7
zFWH`dwDKV9N}j&*09afzCf++w@#UV=0cbrf=+OF22AQP)aLp4hTFhtt%sc+%%-Fnc
za&2gt*Cv>zX51n)Qy&sbS-i=73?dfV2`apFd`nhj&XE$eonwqg=iJDxNWb0*nV@3!
zWvd1iYn2Hn9jVT?Gd81pHmbIrbz4Mmf9H<#wI!iBASsDP`(a!`yv<@Z0Lr>{f^x=_
zyz7zyE4=KPyrK<@5fanlEh-k8qw73Wrq-|1>esF})tlB?00rig91a+|do$+57+24a
zL^04e)b!L-Gf
z`2~svRO%@FGi-Bn!r_EOmCZcUNGdX|s8aQAV2;vVuugNpqOrr-Hr2vu0+7yO%T%x*
z6{*e*5VApn3nsP{W*peI9IM-Y^uoC>>W;cOfMQkXei$IZjI{HvKz3+LNVI9`nB=k3p`d&+T=d7CgG`ej`LRsCr
zOx;{{lqshn9r_>#UpCDeq+MavK0BeG5qwexL}4BZ!53a|wf_PBY}#-nQWfIQ-rT$G
z68r$>Vc-R$PT`v~pX5;AC^4P@de$IE;`wRlXA6Ef@1l(gw`$XXb36G=D(iN6!CwGnKZ4o1ha;?EGFihsUGmb&%z-NDTW
z`)GMq#P8$#1(Xc1Uu!GECi$VAj719C3~_SHUU84MYfF#1%~$k~g7<=Bf_%gc+XLBl
zzw?u3bUnU2IFY?3O@|Zwq_tyyAyKzT^#Zui^ZwihFRZSyw1H7)Q?>7Py4QYXIR9YXbhT&z-R
zt0gF*-dj%I5$}C)k;x9!qBn@!#l_nHYYUu&_$`a(#&%r$Itr1(@6$&fLBiz-yhyp`
zI1mJwYb$wz7e!8n@7VD-IHq$^u52122P{sdx<=L51{x#^}fbNz$gqQ=K4GCG$(%jXVkr`1QWP4QvY4z}7=h$^m3hKHFaFMbV+Qlfpu#4xI)
zIOcs`>dPv?>-XLTXw_Yr47-c_$7_E;(U9`n)G5Lr&{?5XG=}5a_>^Kswzq~)k|u$6
zmS8Vy3roualvfco`_N{?0LYNxwHhK!dMn6V#CA~63@@$4W4nE}bG^L6lq~eg?m}dC
zh~%nc4WCr7Ivaim?Iv$As@>$FjJetF(YZ*F5y+T5@H81I`ERr`8bVq1LpVw761;(+
z4NFcTAnZ4r1BSoT>J`+HvqQi<#c2Ix^=$BZP=OJ*`=Gvd=
zv5OGX^NFoU#5PdvE15t$=t~(ncc#3Hvg?4YsO+EbQ@Sc@C4OT2W-a3{=J`b1Ayc*x
zI3fTa7Yc`YJz+_hE!A`Fl!EE4@>VVDL-bX))GKL6P5UdSM!Un*KG)8h@{rj;`9crV
zpU%^d_NTT(qJR=D#hf!C=f(IYGD9LntJ5Jz3nI9%V{fZt_=V1vwL?igrdbt+#0@9C
zDm0ohvwa1-vLmkQ7^i6w7WkObk@^TKt8a!M0+d+;WF$k1=Vx*!SF$P+LXmR3R#psc~5%
z@-wTasVpepuBn^&`UR}1%E+w+7ygnZcW0qaGK9)&aXEtC0{m#*GZStpI|>A{R@mMK
zThz254vGkzaxG6HA7?HZ%NJC$-JU`gTeojxr1{b^T^Uar!xd%bKb{nWk*?_XcM{Q5
z*-;!ZfUQu;FF$AQ3;<=Pyr*a?tnLJcw^Mdw*11T-SF6}4F2L}`FYgWnOp5P3z=0ST
z9(RYjD)k5N;USuJl<*
z+NZ63(wh9~eQ