From 6494ae82008fa4143417c8b4f5df45d0950e36e4 Mon Sep 17 00:00:00 2001 From: Elias Projahn Date: Tue, 2 Nov 2021 11:00:37 +0100 Subject: [PATCH] data: Include more chromosomes --- R/data.R | 5 +++-- data/distances.rda | Bin 8459806 -> 14684474 bytes data/genes.rda | Bin 320375 -> 308641 bytes data/species.rda | Bin 1586 -> 1694 bytes scripts/ensembl.R | 30 ++++++++++++++++++++++++------ 5 files changed, 27 insertions(+), 8 deletions(-) diff --git a/R/data.R b/R/data.R index 7b8aeb3..09f8af6 100644 --- a/R/data.R +++ b/R/data.R @@ -1,6 +1,6 @@ #' Information on included species from the Ensembl database. #' -#' @format A [data.table] with 91 rows and 2 variables: +#' @format A [data.table] with 99 rows and 2 variables: #' \describe{ #' \item{id}{Unique species ID} #' \item{name}{Human readable species name} @@ -25,10 +25,11 @@ #' This dataset contains each known value for a gene's distance to the telomeres #' per species. The data is sourced from Ensembl. #' -#' @format A [data.table] with 1390730 rows and 3 variables: +#' @format A [data.table] with 1506182 rows and 4 variables: #' \describe{ #' \item{species}{Species ID} #' \item{gene}{Gene ID} +#' \item{position}{Gene start position} #' \item{distance}{Distance to nearest telomere} #' } "distances" diff --git a/data/distances.rda b/data/distances.rda index 8cf319fbc5a0ce84a4a2cc4dc5b5d0512b315ec7..0a47bfa4116ad4bd35dd3909761766266534cfce 100644 GIT binary patch literal 14684474 zcmZ6x2T)U67d9NKBmzQ65K$CD0!b($(h-FaNRS@`Z}cQ!Rb}{p|hy<3ImS1!R%c&GVSppNI`aUtf_mQj6o>Fk-9ZyXC z8xvy~(UdD|TZPXh%G!QGA7f+caY*n>W>qf{G)y3m)4JLYE1kwV7$o z?uHHHbvo0Trc9z4Z)Q4~K^)>%jy22Liev6f$lBu5?7*2ft|F^TGsGw1pnu$xTP$l@ zMuc{M3-ka}02y247d;OX`3|5kqM@L?m4-KgCu7R$Fk;q2LlFRpnN*V~Ej4MLmIlvo z`G#dnA)%IVm^}$$Kf*Ezg3~#vvh)fzOo^OoacWXr)mOj!FpLRCpy)3Sg5EnGwp{t!wpsy zZjSBafE&2j69X(Ys<{skf1YaY2PHTTu+*vMp0ZS5j=1CDZ0ex46R2Z{&jwNnKH_IE z8uQ-JB`QNrD7R>-dq4&_>oB=r_ zH?hilZ1ci_eA0b{d7&KWC=Vbr{3VyLKokh1AU%O(H4Mo+mGy;!jFm%aB9lAi(CKV3 zJ6$Za3tSh2rvKk=UP>y0&=17*FhLXKDjLVQfeCt-$~}Rh99u&WUULcC6l5Hlt`EhZ zkVENUC?a_bV;?M($uU;XP-AY z6+l>dDs=9MKg@RKv}GlwlhmE~j- z+mncAw3_gbX95`{qtZkhO98yFd1Y~?oD>ah5rUPE=GOIqxLHP6c~JqeTwZ*WCBYU) zpCxXcj88oO>8_i2%rVYEq8Q+BQHJ5ac?AKuYfPCY`*JnBH9QQA1yHy)!{-tTC2Ss| zdtzO8hTW}BI7N#|w`69dhgN|_kf1C$*F<4E&52NF#?O{`E&v(<5`Q2BCfTM`q%sF8 z!4?3$q%%V>RVUHV5V1Q8P@x%6s#_ho0aG=AK7pxvB8K*)enE?;vL$?AoE-sKOFDW4D)gI%F6Ut6m?k!Z0-zUMsL;3nZk&7k z_ae%V#^6!_S(t<=jl*To;>9G@q?=kAX_9g@)Q=G>Q522g3~9B2e6TZB9PyK86a=c1 zdMKxq&o7V1lPuv>aSt)ZJ|^|I5Pp`5uAm}pTGN0l9C6!+sZHp7Xs{4|APrb7rv#^Q z#xSXm^3V=Ue9c$$%-(*WV=}k}MY@y;t`&ew_K?h#c}eN#FywVvioFr=%3*&X#KfLZ zb?&a5(vM8guzKmU&MJzV?8QgIymGp6|GeP^h&~L>_;(vXq&2+N`3`5DIPNhnEr7Q1 z1pqXb#OU@y__IBeGzAFT=%F04V;n>DpywNyP(;7LOqm(=WEQniDE}@K{fLAeCqjLo zdvHJGN%Kbvf}6O$=y2gfF+h|q8I4MvGW5x?0oEi=iCK|r>juW9y&-I~IzdAXFi{zZ zEE-^O0`oX&NbyI=G-)70K>z`BhQTbNGicxPh;bCmbsSl-I1T~y$5sf`im82XQwz-g;xy^&iCJ;`Dqowu; zL`O`W36%`zaCr+d1x^zn#rf!yLky8TNLeLHI|=PlhZL2?7l9GRnaP=NcyRmJy1xRG z%m@gJZ`iu9!{HAZDs_->CaX46;KVe%#1iW zfp}R&Ad%2FWDz%Ef>~u4IcK7exA5A|2+=-!5efH^bq>r7F9tcAh&`;HmSU9C(O|A* zYLyL^3mU?<$|8J1G4``uuGw=cR_mIN05SW23ufC1F6Q&mFiOOJP~Erj)wsaYH~lxV&{1Q+sB5!$P1 z=ZWH$mCO?v^gVIZCm#GP%w~y~(aOsxrdpbil%fTcrYW#{hyRFQ9tYXBi9=pzMkkV* ziQ>v>C=(_2rrIP8PVr6e1&-tRdyFY@B)tQQwRclR>IlEY5%o{$Q&zg#b;S6fF@R6H82dr zV3g(Ljyq7uuIC7Bzce7~nh>qJ1`oO6KqA)~vgoP=9v^(h6K!(doq1X2KpJv2nUM*y z6NglTKq?ft5`rj=Hn$C?RWAT)4lNy&g(gC|tTcZklf$s`Xw|#$x}OXaZ9~Y8I9@xI zA&tT3c%vK<#LGjxx0xu@A>J|qse1@#;f*7f!|O(1bU!(S@_`{_gkhqa!9HbZxBPA;N{wFQP3JB zqrQg}{A*SUjSXb7!Q0%VGXy!uO3Xbh!;wOe|HQAf3aJZ}!?vN>rfhQEP#X9hDzf=n`%2_oG^^sB&$`(j} z`%n;|AWgIcNF3Ib(%7UmSyCNWhg?T{E>u^{&jr7uGut3| zEj-ZyA=CwhaF;So65v4nRFq^Bpil%zPFX;#=!_h`YJ|YLOT=H6Vqk^0lfeUQaXmK}v$~j67XYG&&@zfpSW{2A zGb~cUG#hgb#=L@I7>6?vEEHKECA7qn=RjENM2HR+?aD?{=}bfhsV?sDgOpMDF0l+P zTVDcU&13Kvm{1jPM&V&Ak>?^LHJC%Fh$M-pr^)>NZ`bnhTr3=p|S zG_M2zhVuvjfX9#o1qtBE^M-K50N#K_0vDH=M!Fo&FQ|rO-Q^&!F@Z->DB`(7Z!wV& z0?v(Lu5ur4&WI8_y^Sn5LFVf<0Viansi2Vs$P7|HnW4>}c`t)(MF$^mLlMtc=7rY@ z-p?R?0;4|@v87y|IxhYeQ_f<>P!EeK3=w1Jnb+$uX+OwxO@g^20cWKKMnA%${sKb% z@bKDI4%!=x-Vq0tra=l1J${BjX`Y(-uh$0D003br{%i8pF|YXnC&w^JmI|eOmY8BI z4Z-3=jPc2MXqw2zyL)S>6+a=!9Tk#0n&U$u3?LzEdE(g|3SkmRs?9{_A5G@5d&FZG z&LB5EG*^S;bm9ELu)HA?xnWuPbC4U@#GB!C+5&?jJj^3A%yS}xr4Hx!!E@75C{!9A z;gknzJ{*IJw$m13YI8rTIxOiZ*`WSm}>w|B0;1|o@zrKbQ*kuj8dGJ6J;Yvk1CbU z;UQXZTGf3h*a2rm_(-FhIDl~mRTQ%-cmWtenuZ>^OFluuS;|q&Mo7RUC^Q7a%_fo% zbu`>h7HkMge!yEf#6UkJ)x&TKSTcn{<{!zs5st1~V9EU|q@$skgsRpmlp2+Me~ozMAn-Ml`Tvh%++#okL{a2a zY6zF=4#91J{Fy-lkiG3BWN^Nmlbm^>eHGvyp2Zx&u;>#Is;?2!BNKa1mVVv@=EO3$ z3k4cuNr+SuA{l7x1BUKWV2zWsYB!=a54KH~u4|rmr&U~HwT?N9O$*PLy8_Es7Q@SI zGqYmA&?+#=1B9ZH$kBAXrMP(!6%5^iP{WtWkX0h;(!a!^%(KZqdCNrhA?(AMlkX3S z!`?I*S#4@KQDd6+|L7ibpKWQtLK42JPB2GAW0oo)O_Z&kyMJ_-PK_1H3eATPisM_R zaCo*ald6EsBUUEE%2G`>GUQ`Tz+)tM85m_7Zgi60WPpWJGiN*J^ZO#i21)c4-UV{r zBp3EGL%vbW6amB0=_GgFKV-H>GOU3_6h+c$Z+X0E8DT8u2?V5~qcXF`2YF-YD4bb< zKr@;V;Fmz6uQ}TZkNqD$jRNbvD`lU!3-ihweboEE`1J6IiMqST2`L)|1*80AUxucE zg0IJC0)(3lP)Hv(BDoXD<>47}v_wmaSkV^AHb6+Oi$E@cWKuN9j12;?&PW8EVL@Xx zu>s{IS`RnFBF0eh9}=w{$^Q|LU#pauu9O&%O=n;cA2>?hM#r}6;O=#)&18(z2#;Q> z*dP(;OK0@R-T)_BZQ&e*_4}3{no(sN+rU)Lo-}*0Q7@OX6SJNos73M`< zjZpo_lPLTUQ0Y^ZUm*6cMUj2g`An_O!6`>kusX*ebMy(JItN3xsi%N3E1#^Kil>!h z$tp~=PD%J|p*nkB#wnpX$4>;Z?%`km)_(s-1tvz?qlXBpWU=+1 zp+m)a1fPin0MtDo0Kf@#{u%Q4xVm6cC5df0{ak`4pS6^XDK?U`ims%)F~It-879ft zs+-RbxmupN4L=%k>mB0sJH!bB{sN9Dk3(nqpyAiy_|qO#<76XiAPO2|PO?Su#L}3i zlh`UMItvIRyh32J*Qmx0Xil=xWhw~&0&irEgI7IDCbII8J_u|fPCy&Zq6{Kn4KxTN z6GLlfLB3F6t>}D57UXxH4Vj+VH4IJ<#a1!+t9?*7dpJsooP7J=_`Ek)vJy+$4>vSg zbiU(I_Wbo4MfiWOf)%+cPpZ@>na+7vTX&xDMUM`#5fFfkpujyrM6zArhr~vraf`S~ zqd${$5y?fi@FOe-p(e-i2>J3MSQZ_Y6-pQ|w4;FY2eGX9Du_K^3V({0PBkWCu*+i5 zOmU@DSk?s6KLfNtW)?$9iqnI|!@OIC@5&is^0jKT_d&>??8s7@sG z8zj0}hVNGxXI|AQu4g0PLW&rU9HB_GDV($1|KVS$ca#7Ubhr0H;(y0rjgblnzl(#GR1)ueN!c3DG1hK8-^_} zh9SE~Wh6I{P)gD)kd8Sd1kdUuBzs^Xzeo&YUZA#Y0k2ROAvp-o>M)0RVt~w437Bj> ziZ@#f_mhoIhl9fq;0u`y3lc5yyOF{G$r)*M2@cX@LaaX#hhx8NGvZnLoyN$w4VXtuVM`Ueux5@VbkQNQKWPHq-pb;(?EsPj134!CJ$aFo9T)K zTj#-x;byHs1q?%;1ysx+v7!uV@v=B=uzZ6d%`N{lQ%-4#OLGtt&tyWp&9e<@p;YXl zgOeDdb`s>?q1P>b(jI@ZlFw%6bLgBX0rokVlxC}4bAfjnOHQJgJGj9 z@3)b;wroc_>PkQfWOBB2fyALuhynE~S@sg{qXK|MKs)c~NeUK!wPE(W{UNatK!O3t zO7>lV2!KzlpM$fJ_Gm0kf9aNX1+plOXTYT4cqW6ymi9Ud!~@7LrAY)rQhVD;%sqLY z%0{9ok;JR@h>`8TAHnmSEejwiI$%{Z7xEkGO<-LMzMQE@=_*4ohKeH~&E(yO%S^o< zwHMCwLTHr;Gb_<%kdl`a-`f)!Lnk!Ug>%&8g>D&x>2OTvO^TY^-4Dm50I8ft70fFS7&MX;xZu|M^eHuS*~Pjym`V zFTiOmIxbi#R3v;Q-T9dR_stI`XFpJKAGm(^X;NMo)qiPnE$XPqwX5&ryLT+?e6=;t zJqqpWivN{Up5S)X=a0v96+^#<{O&y^v>E4&zaG_Gm&|edR-WWM^CjyleWo((+~@vD zcEwh@LjmM%Q3uos=lzgg>rHOpCono^M^$J23-nzfvXpbnD_1Urs@!tg-M#FMyIB>B zstRQ*y`yc%tA2Tr4$DBtNPQY;DCK>e-tuwHbk!g^73<%d0y9c z$i|ob9$4$*TScYyu|})ONa|Gl&&f>T9V(7lm0LdW;l@LmI!#;VRbjod@svButiNgVdu+l;I94Bvlnk*#5?s8kAh)NdC4OiHhDUv>KJXye0} zF3IV>r6xpxL8pq9u|V}*oZ6|8C-1uBw&uI3&NEjF?Fcfb?H~@Iy_M}1wx()DMeuIC z4EU+Na4tspaoF!WzU~t9#AvfQiSY!_FDFD0;>rzqPfc{c<8p#zI8M^8kDCS9kUU?HdbN()aYp|G4yiq`r=GT z7bD>7m#yYm{%3?&%KIa}&XEyctxx1TN2n#uC;4=#8wmI>n*|g3nFjCP7A-tm(9m6zBtGR_FjuR*e?7bjoHh&&)fJ0*L_DRT|I2wf zcr{KzU+VqOZw|MBciXRA4H&Jf2soI{;c6<>t+yUAQybLOO(L`Jb#P}UTn^&7k=cT?!iJHl9=yp)kHBPa9#U?)VzRQ~`A&o$KoT;;N;8OET3XXXs zD5uf7D^R}QGlKt#Zuewf!<8$7x13FP?tWO?a~#X@aC2k9l!EuXx1Zd&KzGj_-i>Pc zdZM!cTQ15PR^DjX>q`*%dhH0sWjaw|@u^1OHRgw~+W{+ucf00UnwE7GFJD63%b?)bqmYoWAb(M@HLo z^If&!1Ouhat;p>I6~5r2n5jF|@UFm;{$u!_{Uz#1t`}%I>aq4c7vuVzS6*tK4226L zd$U_wtur%ULAo=Dcg+Kyt@46R{@CCr=q!=T(JZfH6nXSdmv&~COp0H#vw6yN+oEn2 zf)2%=WeA=PZJfTWqkhpT?C6Yrr~XO=^`%+hXrHm(tf$U9jJd>Za};S_~$7iI(DM((S7c- zITe*S_p-HNhvb!`Y43tLwJ2euvgC|WW5)>Q5`nIz60MWH^YWQpEaS7<1W01?_OgZT z)LX5|h(&efhS^u5b6aleU_F7Ib>4ipO&tDVtSI(fIe=MoAt)HE2$_l#L6IIcB1$lX)~Uh@I0*MaqD=`)?Ow54eOwE^F!ZJ74$yJKfJ{L&&g_>r%q0{<8s~HP z`nRR=qhT`UX#B&lMgrt<{t4Yce}dxk_2IT@E&rnO32oAj#>LuOT@fld#jw@;+9ma# z=bVoimiNxpIY&q<`Fz=&O<5V~bbBv2IJ^3`jBwlXO4>Z1iB$EXXHFC7zC4NdvtbLL)oBLiozV?(#}6YPklxnYk_t@A7$<*px1`?_OLax@Ss)ZaW-E&Ve{ zoin{pc?h>}9F2KQV4+S~-*t6Z*p$rpbd!)As>IvTd(U*ifV&Hs^9;9oQ>OmKr4|D} zC|vXGUSY)vO;wHbCyXq&#O^5Yk&E=#&Z%16YwHqfL?S~H^iu{GKv_GX>7`?-Dkbx% zNXM9>=$KVeqaKDFXEo;UFpcP*?Ty&MqMVc6_8Wq)j$fHuCaCC?t9k8i?`+SzRnPea zl?!b=28aizd9Cj)FT|)tk?E#-Lv_P}@6k_N?$o#V2MX&fPqaA>;x`uEJ_`m^d4ADq zx^vp4P}I+ZJ+`oNGnFxJX&sQ3U3TKKKe<@Y&ir*_)wyw(MU#y#o1@C<9Fv#afIWfS z<@_oxy}|PNCO7SPP}0=p`?bq6%$qL#&usiREBZASN=szMe2;P^M}?O8OZw_jdHSjr z6IH?DUmhH#NKN)hkCa4x5kO}uY0NT zJDR$E5kRpu&Wn!Bci-P0)m{uyS$Qc(Z8Jk@#;DFSO?52C_)w89DCbNA9^`HnBlu|3 z6=!d!$=LiB?PK2S@qEUqurRiu%4W8GgVeAeBEO$sH;1CJe{b0sJdz;tvMB*CGji37 z_I4q?8^Ij&4)rHiac4-wKZ*cx`eLaglD*41F z>i&5N@r`ps?EAp0y7otqc# z#Rec|Ldo;X1cR{PhFG2Z4kKa^;rUwAz3aq|nQ@0@87HWe@Ful=Y5v2GemMu8?7l>) z-e0_|-*Vgdl#r*Qp@mY%i6KAlBhC{~JNYb-(V-UaPY2`WWlU9arb+f!FMr`&gk&T4AkSPV+HEwr&)5ghec>kN08 zXf7}{DTvZ*zz)M#nyN@y-4V(`YT6UAd-_UMy`OS6Z-;m+H@s_7_BJnc&a<3S$*tHSi>Fn;A^1hdT9=PJEGII=%#Cid>QpG&@U?fH(ZW;$K- zW{#VI>iJH`G!zGc`cyr?aB#JjcfR;u zM@ine)aaz%h99LRGV^UkuTj~vyY%?d!HCCLsonYch;dPeCwvw5Z$joI)ZU!TkvZ$9 zB{>@uup+1=7$c2DSk42E6HFGit~Bm-uUz+{Y0+i>xC&R25HAnn=n>DGpC;o(U)%OXivpL1r^ z^EnQfrh7$Ll`qZ?3qRPy@7BZHf2UPoPG{v7r&(CJKf8X%c|JK{O8q8z?oP+C zXba8aTYc{^r@mT0;TUH3JZp4AM+@#Hz1R&9c0!%3WrmL~^^VcLo^=rXbrh^UGrk~( zH~%7fvd#C(hvVUawS;45b*@M8-OS&YX3uvTm)EtLD%K_&1^Y{M}V{mQI){ciJQCeec=ZN6ytL5PDh1j>q@DylpjF z8KS4rT4198k{vp;=F_#7|905ML)OXOoQXo6XlB{7T|?Oa+BVMWMHh`r!D!HA=0!?M^%=choWFszdqR;W?|2h7EXg z6a7u)#cfHQAoXER*^C}ZV|HJOW-#pct1N(3VAb~Hze3IR=Yv2 zQaRV#w=Fq!ncai;B(owK>Kf)n#-l?%Q50pKpZ5;C8QZPU9#-+=;?TLiXm?;g;(gfM zbXa6}Z%l#pz{vU13)`R8g@4hl`szcjydLd~b%7S{Th&7R(?vY`s+S*&TXATF8{yX#L3>8FNz|#k*rtnbID5r0e%@0Nbv#l_0!+Ena{+=YBa=`l_+K zv1ffww7hq zMEk#ZB)$bsT4!-f{(2vog&heFm)+EH2^?oFch4xy^l|7znoffCdUU3IdsYhX&RNOI z%*RElg))7fCp-P03)3c;gg(Y2mHd+(jUX%PqRKlduHUPgMVUo_D=7mV9)Ta--S5f(I-W)7L4m|0qRKQVQL)Tt!Z z>tGMQQ>&HTH*>$fOfS)iTa8lq98*ZSe>zAwW9fwI`qze+R3LE_@oz_$DV7A(Yg5765uS-cPD_tF=tIML!O$4dk^45mV=MwV#y>`b}cESN* z_u5pd?a|7&gWb^|hI9F(WbiTAV>8% z&jIYC+)Z1y$>u{WM3|3@IlN~?>}Mh~@V!q5=uNTj{n9$6v{Ti6-`qG$frS>wp`kdZ z3f8w_MQuO9F5%!>Q+|ypz>_7fba4?`81K}v)$iy!oTz`ZWL?|B*qSehLW)Hzyg+yzU(H~n*ZDI~ov3k4BN;O?k$6R%vOFK9`I`CZ4^DEG0FTPTZl0BFddws`m zHtnxP?eKQziSbDcS12XC`pp2kBxdYMWK(_*U8`_|X}Bb>Qs6nk#u=-#hE0yNnoiDN zBp;sw?gidaF)95%rxh*kG2J{|lv)%w+pz!n#P_myZ9S$$oi-#6=k zina`XJJm*tFbM;0VVyV1JXH2Cclc`9$_t6HT9;puU*O;}v<~xed5k z$IhoUfqo#By#_fk@`AqS+SA!xl(^|~_IBy~@vBpH$7JrAf4ukdXYo_s-OJQOu8C6d zMHvkO@nXZ6Y6oRQ(OK@1sGeqEp`Jn_+0~=nSxzZbr9*bU+F@1?sHm$adLdjV#wAO_ zC(|9CNpu2BLSXD$YXQo4gKZdhBY|NYXo zlKtyje{YeVce_5cC)dtvO_Fy&z4&hGwW|;{?B7b4MDmcJuYa2V30+F8`E%*o35O&0 zS0d(#at#+ltoZuWkUv?MJ38WDxc|@^d(zJI>6M68_~*dDN8`R$g@{(|i>`>c^qG*? zMUkrIFLN!J=uzML?XLmf(#IqdcI{jWy|2@r4{-S7tW7K@lgSe`p)rx+^iM&)yJsl2 zzZ2Fm&%SDYnfVs>zPrxc+~?{sAS^PhThiApX%G9(+zuRU2jllWMMAEx)dud&3mj-V zz$d*1T@a73U461GKbD`*@EB)6POkSg=z7N{xt83f{&jBdCHw68zq^mQ%%zVO}*nV-_N-sj(f~%p9Whyz=~9KXohE;!UgC{pxFwl|MbB zO|~EZ*x_6`QzT??De&zIcypOTW6 zCpn}sZ|Oj>nS7M@k=LVAF2z5cAsYviOp8yx^HL6|VdcPgea2}n(oqva^$#uX&Gg<; zwPfbX%=~#{v7PZ#ersWEXxr2CdPmgGpI^hq@6%~1)&um@ychYltJ>zEtGf%?{~;NL zFD+OeZ1oCc<^L1N?dcLvgqQG+kx@J#_KhVRZ2y^ggJ8yL?hhu)3^9PR5AN-~53qc-9AEca1)kHlGBbD% zz-I@UX}bbm#5-}an*-?Nk7&pnxjbI;3&*U1C?n0~yc`#f}PUX^)TBMm4{yw#j8 zUL|wG3+E_P|H$jW(S8Mfxjz{#FIXj(*OK7HTeju+JF|iO!m#9+w`gz-F zscUI%mykaYP2f*^#-Mk$y|A=M7^=%OJo3@*C`FZO)^lNTB+0nEr zyd^Wu{>G8J7|51NZwBHxVu@#_K*dgZMBA_AQE2zsrH!DGCk;V_6_K&`E8i}bizT7d zs}j~Vub+11+p#a!lU%oXH)kf|%*oP!oRu7KT|Su8nYGgwm#{(O@3*^lz}r2aJD#K> zQB`}*EFdy$#5a@~ed)2E^_rirM#9Gjyie}>MB8~dOWP~l*BrZjY9J<#8jMQ`eDcBS zuLIWr#E1Ugt(Jpi5@s;dx2jSg@tA7M2FwMkoOX}vC3*j!nO(p}*_(z{+S zbtIJz78Ek3kkSFE@MEf%9u3ytz8K>EJKNyFU#kZ}(9iLq@+&*r74b7d8hVfFB>MSMhJuJTity1=Hp`8W%$>+g?QPe-ymGmfIcgpKCt08SrJe z`CcwxN73C7NxPiKH_z&6zc}K{h_fwAz&$-F`t2j>X!rTXj%Ie~gxuTFb=*Q@$J|3* zHK&3%QG&pKuDn^%$CM6Sc1acfG9~xGzP=`#l8_vA+~*7AoBX^V5nr_)-)f=b$KD{5 zD@PuppM36pAi2=RTXSG%z_Dm}W#BZH*le7=|J6OxEb+Mx5bu{8vI&H=I+`kj$iIZF2>bORC~UcWEDp#+!lXi-F*Xp zxsF<3LB3bIb?&-wPxA*Wqtg?^2aom}BT|g5LfiTG9G>B)r<2ty7t{#rH`{+~nRVp< zO5ZKYa!^10%Ff`&yYCvWC%YGXM`1idR%cTW?k8MaKge;L`bUuV-G>W~XJ+DNO=<5+ z9ByZRO<=xSt~TP8=);x}TW|mE@o2$WA{h6Dg81zdHH`weo(-(OlO`>Rc-Be zAvoxb;v2-XOR+VtZZ27cTRBnd;N{A8(`h!QeYY|vxb?sPtR8rgY%;%TW}mt6c>LU_ zms*9S8s%TzH8VNL(H-(&IVtT{purm_p=NhilsBOeqW;dsZI3e_Cl+aQBI+(daot)o z?^U1tYKo?!d(ZrnqbUL5xQCl72UlZh(Tc_iXS&JDR7HMmf)RDs`LKZ8Z9 zXf*S>h@R36eJ|aQIJ(Br5fI^iIG%C2?n~?7S5e_VY>`V|b&7Cp=d0yv6wu~vgfOdD z)addhP`^pS)On5T8#7T0FcvfprIB&%{ap8kcC=#A6^^Fgu&wTSZ%RYYhpquZ)v@Zx z`i^p16t~4I@Zk1QRl?vaHgx%$Ya5z}P|9EAovr3uwDdE?-+A^s?Yqkom*+!J^5sW? z>9u|;=raQRllB_B&6Tbupmfq^lgSr9_IQG>#7!j^qA{?A31!$MHpLeZtb3;lI;w0Vb0P zmJ&GL_~m+@wqa@hqz9oB^FCcC?G-9%1i^lV{$#fT{5I zK1@ps7HiP>|H=PHBE?(UusFRq@Xz1M1h+E5LSVl_l9q^Uw5DJ(<>w0UJPZ%ieYxym z_)_^~75v_ETklWI{V9X+vZ@-9aI0irAqm4pfcL3W(t?Spg8enX2aC79&hEGs9}+NW z!Cx?b8h@Adz|pVI?Bt$x?He~E)%YX-`3;Mkx*c)+XJe)AscWD=w_QX4dYy*pM4p59 zcOK-D(E|lW^`6&qjw-4Dxnu}i zh(!_FkqiM(eyrAoF}|%vOV~}%S7GxBB~%`hy`L3|S&NPW04fGaFgN6_=Zf|7@kgpu zXcvvy-+!V+N5|z8_+k9Ky>wA0%bvji(!FzWonvsA`F9=JXPK#d%7l62k6*)^f;bN* zFLZAA_CWN`fBWZmcUgV@$uo%uqGV*G_Ka&0)#t&D=M~Q-7vnn*!UAGILp#2ya#shp z$M#$PVH?_43Nfv0w=rjOE*-D-j4C1>fU*(x@2p1RT3j2+zS^DjB{Ls3^Fywx1V-+~ zEyNqZy17+`SBQl3y7#zI4989)>|p2d)4?MqJiP08JwsfIgDhO&Iwf#@+d@#ZQ@Tox zbJB1D5F{Wtv*eq8RGrW->ks}2a7)^6+KjT->XjvD8x}@3+1)#Ozqfv748KOU+qO;x zrB3$4m^eV{fRLcbjqcmszdCe^49(S7!{Lc%COJm3)~S`tnId<;q#>CrTl7;+f4w() zBm+@3&hM9htNt_YwsYIoaqA^=lVZHnL3`hy_k0hSjd_1{ksp6GxmK}S?WVkD`t+KB z6G=+{N(iB!60m6J-tBs9M*H*C<-;4V7q55_v~kY%)*ctRF*A~0%JYx23jCK25|)4F z<#b|r&Q1)7rmm`cx+)bm12PW~tUEDZ{_G!c;*E~%{G?T3_wF0lJr7Uc_11jTlf3E1udBPchAMXPPMZdEBC=^P zB7l{X!o$b24e-&8XK3z5>MK!pLOxj^TEI&U6?$4S9Dg1525dtC8dzT{vh zOsMHJ)QPMh;*rCJh@H>xuRI=D=M5HoFB(_0Y>^@;0|6VpvuLTq2nvb8MeYXh8nSZI zrUP#nkc$ylcAEZh>LoLU#I$}iYht?g5bk)E8(li_1{}}TQ|et_m}LUM!m@7!-XJw; zx@|>b3q?uSaa{z zg|Om|>M|ro(!=uWL5MS({zNx4^AKU3j8NcGSNILK_$Piu)8=2EziBw@`4z8JuEIcUq zZsDFS#-Gn8@Fhe@3Q)z*O9K@ql-->Z`sWXcDVA!m;X(T9S6A(1!$j#60)O5)AV4-K1ut^fC(&vDN~IL9VgG@9Rfd(21>6m& z*rG752N6TXL>Be-ZW?D%37ftZfW+gveWQj=rim0XLqs{h1~^TvNt6)0t-QV&=A45& zY4U_A@si13bRN}qtyu?=94rWQD@X?L_g7ksjFyiYCQ3Le>{sJt&}GS1MT_;`QhaWU z+CE)kZUxY%#)h=P=S!?>zco9dPEb#DI8y!)(A6&EoM@k0cg1Z#zf|I>YBl&AJHGYB zFdwwHuoiho(|9;P#j|3I5!n>s&+#`(kvQq}3qN>_9A|G}Njlwf58&S77ol`LiejmQ zeA5>(;((1g8PD>cPDKAvGw9;rtaUJ@|4q0*?`Y%b0TeiLxwFlY@0oF{ zd(A!Tt|f!W5HtmOP1jy8DA?n3)XVR-fwg+ylWgA_$g9t}ii-&~L=;^?gcHSf8m(5c zSHl@DjXULj3g1GrP%a?~Q7%CWdPnXzI5mY;hMk&HE8TjbIbMT)GxVnC>ReQ}jKdhj zW$1k9E6m!74*QHvm0h)TNn!I^)z6t5K{clu%=bS;{&3E#x$yQ`WO9tY-Kk9*efjFZ zw_-NC`yj8Z#9O|CH(+m@_4?VG}2)N9R}V zh+MneWuP&_ii*`y606@{G7bwhD2a)$G1-Zt3G#jo)zf6Oer`{Si?crYIAb0LGIj|w zz6aQ9X@4$*{sg;or?>HoH~6ELJ$KcHarU0FNA8n%hV?~(-xXs0vxnZxo!&eP^(UFZ zW9}(DQ<;2uRPb))X=hTE!{dV9j?r~xX=Jq_<}W?s$=>nEJ$a64pT}=BHF>1a1YXhG zINchJR(3ZWyc+xAGnCd`W%HEH~zkPb}7;=WCC|gk?T3ZK%O4^p3{`1jbPhBFi{pqMchSYQQeJlhmhm z@}(bREl9w{OS9W{HZgrA0zE^p!{m2v9~bS@xKYfl?tmuq27!%&KlPF?Q<9X>r|CJ!^8VN&Jp}2!mjnojQl)nY5E77H z6QxDrdOwf9@3VjG?Ck8Eo!Qylb9U#wubnW>TApQh_z}uLZpvmc)lgwDN6we>X8ShG zWVhF_%e*X)N0yEvcKMwNxX_h6g)N<0)}9A|eFZp#8FuDH(lso3WMo#Qy^TCc2haud z$HrWHtmiVZUt7P34bE|!*~6%ED;lxbn+gBb*~vg`&C=mgVLg!{$bw+_K1x2chuOAdd&l_LW2<`>-BU~lP%X|i;mUD zI;$Gbl+Y(-Vy8PAzs*=_v63(l7*!a`p9yK6eT z7!MVpjSOn2sDX8ZSA*zkFH(o{yy@y7T#t4>z~ip>L^q%=TSD}$qpEkZo9!M`z zu|s)E!S)~Md)s7A876rXC|I^lMf<6`kt==U4t3=46ooy=@*~CMU%rEALjtMuj;~); zECt!A(E>Hjgj%_T;O52_x|+wt>g{2AUz1dUu1FjA>thD?%H+Hu72n0eUwE-_FG=`; zb(D;20gPLn(mlXAVYMWS+yM&uE69q|{KdPJf>qe-!{%Bs5#3B{(No>u994EBeIbRP zK8ikK&-g^xZ0+cNlrM0H{|Oz^UP2kL63>~KF&