From 5f924c46e05fa6e84b60fc64607aaba9a150cee3 Mon Sep 17 00:00:00 2001 From: aaron Date: Wed, 13 May 2009 15:25:51 +0000 Subject: [PATCH] Added documentation for calling the GATK from Matlab. This is to document the extreme basic and experimental support for using Matlab to call the GATK, and is more of a placeholder for when we have time to revisit supporting this. git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@687 348d0f76-0448-11de-a6fe-93d51630548a --- doc/GATK_in_matlab.pdf | Bin 0 -> 53593 bytes doc/GATK_in_matlab.tex | 79 +++++++++++++++++++++++++++++++++++++++++ 2 files changed, 79 insertions(+) create mode 100644 doc/GATK_in_matlab.pdf create mode 100644 doc/GATK_in_matlab.tex diff --git a/doc/GATK_in_matlab.pdf b/doc/GATK_in_matlab.pdf new file mode 100644 index 0000000000000000000000000000000000000000..c46466c2ac9001c25b8485068380e6fe31a33392 GIT binary patch literal 53593 zcmbrm1z6TkvoK6bm$Wp}4L5PqEeH}K-QC>{N`o}g-QC?F2uMmxBi-Hg-he)T&pFTg zUC;Smuj^j-H#KQ^8{(6!$J$v-WTV-0FIq ze3ew7lYhzaHVbQKP$)-Ko$V~MIH}(FxVur!{qw^1$j3gd4kpIIv0erR)?+zDA~CmFJhFoPlx7adM7-;R z0nB9ltyXNb0pingws)Rq?kYwON%14N!>y&b)`tM>bMlVForcP zt;A6m4jyw87$mG7!f19Xf-FrmmbAGFen9;2iH(-~bEVA*G^Ypsttgra-}MRoW|l6W zrVo2(fy;X_Zlqd847ce9>d%5G1S?Vp*Ak;=eH1!H^O32L&V_GTkCv>Hx9BB6i%_QIULl7^H& zp7yMxwspNf3)hD9g(vsPzqJZx&BRkp~(Qz<4B~)-WwB+FzN?Go}6l@c*)7pggL< zSADB%{qoBLdC9IsEoeqd3bVySy`!=DT;%q+Zs z!$|2L#Dc&=l;l=2)_+m#Sc|TMld5u!iMgSzBPG6`7y&zqREDzl>e#{grEk^F=fg0g z5*+FnV7o}pAW2Mf1{*YjO#4Ju;6~$u3{*WP8|SKZG%O@S>G}=yTTYN2ZR2FNG1oIH978nSTKzumDqF>eDhj=gC(eo2KVL2+V)pK{qGJis& zV+zBtj?_U!_d2!PPKpLnVp4)f$J^j=c`v?ReW7XQS#3EJ?O^S*NM-l(rR5nc z7H{GPrAj<@ymcyCEcW)a^T{=n4v7X$NHDv%a;YdP4d#3W;)&FOIFY`IOX~LLK$mjC z?C=k>iIEwu1+y@xyl3*(6uxd^W=udPscS=5>ofsB8E91*;FYL^_Wi=<`)M2_>KUIt zhJ-Iet&m|_NYi4C$+a>|o707ksv2XS(#9UWYxoE>>MIUK` z=eybb4AX>T1?WL$TE(4x{gVm>-k`il|CNkEpp~2Tx@Ld{zV`an1@(Q-*=BL`05YGb z{)#7)vbc@t9eipQXQkq3k+##lG^y0hNs>AfmA?YU&@hgf0s9T~$r+N|=%&?P*=KKM zBsah>eMn-p4>#RGxy#I0^BR*7dz7!!>_zn#TrRl;L#M5IwZSWLVkm1ralq*bh z<9D@MNP8s38mQ~9sml@7)`~W>1=Sm?5V!V8@oR<*nqpb=O(X6Svco2EBg*z<0Hmwk=A&erE*gveJXIgpN;w`>J;_=`2B=lwv>jd6lR zs{F!b0$bj2IbB#ksql7Oa{b4KcurWZcehN#uM{E7GA8-R1PT7{+gy=`FaQ<{pT!_? zK%*#(xroum>~q^Y`oHyL|*Lx(SBhxCNS;7ef8UtnD&zkaut8$^4J+-+ANDK6~J9GC16 zfI{0aK|%rX5xj1b$(``MZ&P&H*$Ya#z>!JoECasl-S6_JFRSJ>B+X4|Pa~UT6>#+^ z1oTQP&F>3@Jc+Jg>gWAN;SDVHe@}v-9B2ZCXZdR)m3Of=AZC=+F*1N>6t=RoH?Xv~ zBYvD|<=`3R4D75NZ1oK6h?$=OQU>}aIzm>?#A-|+06Q}aF$X7$20W;`EvOzb@Cl2L zkN7vA$69|T>^C6VpqU#qJu|Zrvpym|W;ux2{$zo~?0>SH#2imq(Cp19Yip(V!N8ta z4b(`HH^hwc2F~`x8lZ{&Urhey@+T}{1_J)mu%KI6@-ZxEy;0g!e3g&eO4Q;ILLo2% zCN;%N&S`8>GWGE`zr@0c@oq0AAV)k@&gs1~B;<)nFA3*KBz{zDqKszqcT$LtUYli5NLv4c?rU?#S z4n&0lrT_*YjQN7q=*E#HYixdneN=77paa(UAKrcAz}GEDJtN zIng*9A`_FWH;Hmm!3g=0w}5IbVVs}SCWkoA+$vpbM|;U7BqE|~LI&k&J8m$F_Hemj zJ`NLI4$n_zzdaXS*`=vOq~yzyU6T#eNit(iL2 zh|=#rC`h)5CDTHlLL4XDyupqN(A-0sApC(X!z04`-KUWXGGcqg2+2*m$9mXWyIz+H zHw+OIa^Epe)q=Ua7)I8?V`ePd>x*^`+^-3HEMo(oxv06wIr}y6yqY5Vx07TTJXdYJSYs=rcj?$!|0X znXOQvRb?A|PTeYxlQ4>ns#w!U>{*n7kc@cbpz23)u@f$b!R% zCoHogfHL#`et015TFlayHc%#x7g?L>oiuCZk8*y#4QeQt)ErcQvz}v9=B|Fq!SZ!n zR&0sZ_)t02(1DRNQ6~)jIy3VgY$$IwNQYQsuaD=m@oxvgJZ2(#qCR=nW|wc0`%UKt z8s5r=Q>3N*$hYu%14*&&<_*Q)o+I<+>!Ri%d(HFaF;~S=* zU9qtbJT&0`9-)SPNwaCbqFZ_YCHrN3n-d7uC4r&SH7o!VdYeJ;bc6SP0i=#us) z=T$A=GT6N1=mf6$8kiwKva8yNRn6P>?IurPvTU1HQ%aZ zF)UjbeihmM%HVk2oS|l>jveKa?i;v3Sc^c-7iu1AwM4;Mc^~twA%p%Se%0onQ6II3 zgJI$#M7F-e-ir!{bJoqphePbN|LMD#A3g6Mr}}T5H}h|2{g3zNeDtR9jQ>mT{aEV{ zqkrqYS(u1HGr+SaXJH{`ewsCi|1+ok=DRCkUC*Ozr{;WE@q zC!-||+-28_T;G=!I`G~Z+M0asI6}ui{HV#4GpcoxL~M|J>v*=(tyRXbpqVf9OXdu# zT4?8&q32P*Rb{;HVz?j1E4PK5$P=Mi;^Z3VMGPy3dxs;9UkljZVJ+|x>p9*7;~J)V zDI8<0w$EH^>xB)Q>9tytMLuZWb(NSp)Fn)X1SzgHk`rD`kcLRM8_*rbJSazbA1Ke9 zdDN^oxlE$ZFUljdn#4}9q9o@Zw+_uqP1HJ=Vs>)9fIN>zOdc}rF~ z*>6snY;bp1X|j!{8*$2|#?YGk6Wkzqe#G0&^USdCYthfl9LoGRqc{$0Tfz;6i`r&=$E40pS$KMAbL>ha_H0bGz|ExG?AV^PqeL5&-z}zToeY-ACFDvYYK7S( zrw|P2C0!0tWV5}>vO`^3Qb_p(c3`+lr1|sMV83ut=_DX`sD*wkKq0t%(P)gTO{{gL z|Jt^z2rulLmGa<{U*xM5OkV|;t{s&m14?`eDg*9e?1PQLv2y({*slcF8|tA7%Sl}i zaAH5h3Suk10sjV&S!5aV-bG>!YyhW}fYg+VIbM@&_CKbXe~8igct6L7OWAf(6l1He zNn5-y6cX|4$HvjC{wDjjw^J^za(OqRjiq9fW2zzsVZ|H1!L@k8%&(Xd1G_4h6)D)F z@%Gj1%-EfEA~z0b$Gs!U4So&pmi{1F&uYq?ay_XID`7G%V*CX4_EpD@`@1IZhPC=? zY0Db75IYqa7H=tCxuSKf?s|B>AUuyz@TlOoO|La$y%N1?+cJ+TZ;4(Gu2j={E!MKW zfcjWsuP-PYqw_MOo4J#USzJ=M`h`8&z+~|5!fVWE6p?iI!4ll)SU6nhja$<$#OYAn z(Sx6YD)db7;jcO~s)V&Hf+up-<&sEiLW{RXPA$t<^^73!vntfy6E!3M4aJU*UBg=XP zfJqN3neK}*x8r!?&v2C5vCOSRGt3FWy}^B-`v%jkI>Dn2G^906VsnsG& zilS3r@>Xx%qyNQgo?7(~k$quXZeyeytap|#P}EY%I}#ba1R#>VBLHaiN6Ttt)@31#$j=05EP-v1kfFTmjq4!));=6ELL$89%KCu459Nl{0s3!XJkQ}wxF78+)!;g z#X78^&ZKU5hR%l2bY|It#vYp8ShVUH@B)MJT1U0sbD!yMPhCE^bT&lUUVxAvv+wa| z*=1BAokB>Hv{2cv!E127LJNCYh*ogxGH;Ji)9aKJ!&15uN++y$AnQ~u2OB3hV-wtO zr?E5|_Hb0gIox;QWgsdzmYGoT+SraSq9x-TCm8lFe~O1WUfJI{HBN4+=IpSB;ufWs z@4pe|>31Qn1o}IeG%0>98FSbY1#?i3yN^~-k^R=w&Y$xwv`LoqL236Pcm+*u>fI~P z_<(m}T|U!{u&EfsH+%7HGLFOgU!-vwF$*m-M3CtuxdW};RU&57uem3t*7(Kbw`_b| zx~P3BJg~UGfl`j{Sz^@+0DMI^uirCanb+`CaND{r;7h z@HoSFsuXZBqe^I|Z11*hF>c15fmE!q~>$!m}uT7J}{wi1xUk!o6 z-k*YvDDr0PqE0#fT=oRRzGHr4=%wY{OJ~~QI-&<0N%w`Q-hJewP2n_(`6Y+c(y1-g zyd6`-!{U+in1R`FD2e>K*U>I>whb#4mi`qCoRe4fncNK#_&dRR8ev#Hu#=V-8Y`Pe zfv+H8!4Hd1*NL{K@u!Z6D^dh|?yI(OuuL9|sj zRxp-ap_}pa#7!tw9~Xi`%71wNSTBM|Mak(L^7v*R!A4M&%_f9&TFw z5+y{{#{-g85f2@sh!#74tXZC%UCfZH_6cEf{hTM80pbSp9J!EaM3OfNAHQ#=?a95R zpI0^OzS**usH4Ev;LlY{@ly}kDJ;3REk-UIvxqf7O%rNIXXy2QA!y@|Y;3znM4-RGw$2t2uQdGq$uaP}9$EiJ&S6b+AENpSzy&Lhi?B8vn zoH31_*_7mbbIkmO71gtGGYyTD)Z82yDv?YlbU4?eX;>at_i~*&7kL)=m6e2>c_;pC zU};^HNuSjm%3oKM|SpN=Ps={Yx={a*>Le3SNAud9l(=^{6h)nBd1L=yMm% z4VPho2x6$$3So(p9|epx)9q}A`CyqBLU!yC0B(ZzVyAnZ5OE({tlp*u5(X)Xs2zuq zd5KlOo=GkUk95n#SBY*(auPsY96;L3kbXGIT$1HkIZ;4QzjL>?X09-=oRqP4iZ+XICQ!;BA&=csqsBK#8-2=7HS_RFlo1 zwos~~`JLWqch`-G&1iEAZnepBtK@W-eWwPq5+9{XXy_($C#P#dtm@O^>pdUfG8{~R zE`#Nf-<0Gv=i{?_*wzWLU@6&lP(KrBqq&&xG&&?(k$M+Bq3Pbe;<;2>5SDjDg+K&YXeXRhcbX1TP4h@qXWb;~@WDw*FfO`Q&5%>+a_B#^wKZ zfz13^?LT&q|8k1|(E=IdPoHSRt8uaeh}i)^22Ns7PUBHv5KDN_di$U9f1tmQNa5Ak znVA_tE}4nxahvvESxx{O0}C;LjhO*R%nD>?0PVy+rOzn;rIeYSf#WykbLsy}`pog~ z^qChM2N1;NIemrzIoTNi#6UJCki0++CI(hwPG&X`4v-lLYAxur7bF!YEAu1d36qtT zg@K)zorU!g&hki&gAD*mnVC45{)P0L^b?8Sc_0A9hMj}$Uu=Feej>`u!t%)ZZx9FT zBe`b@o_ped=KqE~lV)Rml;-a|0Lc0Xc~Z$U{CBTB!&w0!DW4(Fx_fHj-{|a5D1Vgq zM~%;?|B0 znI828(*HB1CqB>he{*2r1nHEG{Zabg9G)?sQcliCI=|~aLD-oeTlA+VSV49F`UuD|>Gi5|>myOjt=3~*qE69f&SiX1L>8S1!P~3y#C+PXF9*rrx6I`08s=m zKeg~Nef$Qpaz2h(PEL?~z$eZ?5ND8#ftjwImY@j_rAdm@E zir9s?WvrV?XgE#L3BU^_aCgMvgaOp!uv;ak0JG=yH{g-vX$o)2za#2=ahxx(Ru(( zPvZg9+Q%0E{o`N*StY2eK$Z_G;dmPR98Yx~zd4y6M;0>^NIX!#K7pSJKBtc;&usq5 z{fz|xJ{jgy`b7RO3V^3Z{zK=n{6G2q&Hp+7%#oGj(fps&XUL;y9H5~AGN@;<9)|)k zJIHuI6T-9002$Zggup}$V1FF6k81oQ=c9T!pC-8{I!`TsuKSxhXrO^=JazkT`j7lT zKV~M7w*mR1-=EK(^dCY0d+wj_XJWs#@kIDJ|CjcjpwBgaL!Wp(Lw>8`*%BVR^Evl~ z0pNV}6adaAKlQATze=A4V0*%RO8>F2zhju^1s5o6U}k6fbHxhea#(C zjgyn*KV6rVL>s8b?9~RCWD7$>`3iCckYq>hziHI*$2Zi*iP+z_%oc8hh9)Buk&xJm zl9Ld4#dFenv#T_`pf!AOb*xVNaM3t5IFO#xPmkoUCXBJfM8**aLk1tA5XTI41TMhS zd^s>s2^)wx3XU9IB_t#QPK<_zDljvT+}?u##*4Xzjf#oUZuFwJQ==Ib+3{Q)EQ#($ zOHc^R#0a*t;~V@7**c*a;pF#h=Edn?8ug)&P+Q%Q5&lBpTwVtC;P%9psJ<{z13oAa z;O)rl;@)gsL*V`>w~~=816atY*KxSJ`^Q#sFi^|>FQK)CgY;q2$v!PJLd zw!d#7rXqhJ%Dk_t|HO?^sIzvfbk&x3s^v**1d~7Y@TNUQV=fXEN`L^_pIG|=*A(n3 z4u}gpZJVKbU6PK>KFY*R1=W7T810kv$%hiEc+I=I0%lrOaGuy}mW1ihOF|%Y&k9^X z0UDE27;@;LWw-#%47eh!?-y)7Q-CZP>>)TfzP5HCCD^YZImFeE9g+DMV5VfAU}`{!hq4`a z8dxxi2ZL+uATaOahX_T9Wdo+K|= zUoLdiV-)yfaPL}cr!FE-d7UPHRE3}SC?aodFCqQF1vZHSP1kPm1e?K#cZo0rIG34v zi7e%<;uiPwsDw&r@NdPS3a^Mp1@hM|y@+WbuWrcl?MU{}BKU{=O~c*6A6|pAVPEj0 zk}s^Fh-|qX;k_>-*u3JnzbNeMZRho z?{~*B%auNf2XW*(y)-=7Y0)J_08ss^b=8k8Bd=?h-ap+UJ#M_6x=jXq73GUekEiIc zrfkXDT=>0JrrEDW^x>EPH~Pk}N^bI?Mzb;7ty}hgwWXbHO8IsORGT*&TkCV_b&9Gh z8Ow9~w3!&)Po$Oe>dWZ{0g^l8xNfvtYiLui_oUvhA9K;@*YdrxBzT3`Xp_C2!|B3_ zTh}yl#JXG1CZ;Ttqt1bCdR;B-Lw-bu$WF&azucgyov`AU-h1@hP59U|W+JT=bw%18 z`wLVBZz5eloY7@nCgPbzi(acoz?}lwjS|m%5lj*lX4hD>xbG^vF6B-ZUVoJ6rotp4 z4cO#JAxCsIMhbGE179hN3+CFT4`mzLd=qSBvZsLQRU_gWvRB3~gVxjL0@F*T4^<)S zL0PKTuVZwb(^>Y4J&p;F!b@!oKSFdPc3Pl$C0|w*+hkViS|r585Rk0q(y-RQ&{NHs z5aD^7_wt-PsMPj4TkIEJtn;z{59xWSFUQo`ukvBIQHwI2KVJOQ+MC`m8By7JTVvaN z`ysKPR<%gWsl!61dY^`FCGR`aE-=bRlY<7f&Q=4LFE7u5Y9cyrR))iSi+mT*%T-%Uvu8IK(klSdCJMFLyHz%Y|DRlrkvp9Vix;lLSZc2UfP45A-9i zM9TD(mhMywc5cP1NnG+5zyfZ0-`Wp>d-$+J@TSq46~0ulXW$z#K9Lf7`-jCg9R%Q;1#(clfw|$ zghExTguC-dj+k-`%g7F@LuaYU)bvNaVy$lWH1`3jJD2y>((9t%W$j!D$>UYGWJXzdQ zt;xnIeI)%pHW_td3IxGlk1yc^vX1iPB*ly@zh&J77-5b|IwZ%He;Cr^u=|b`5|nL< zdK{+}W%9$rC6R0{VmxmGla^c=sp+kcd*b2?E6-)}iN;2MB#qN51v&xg71yh)HY4BG zXg;b3*CWm~I+oN_&(&+C1_?;>*?^uZMp1qqv^QOP#RbbpI(~C?eM~n$6Gk?iCeW1N zPR%CRh5?#KNjEb@AqQ6RIefETwfw3v<&q*zR-Qjg_<>3`X?m@5@*H3!b^K_yQlcYH zT$ETzc+QOr4NKMq@%WjSr+x)hWMAI0pz5WqI@XzmmJQBcD=qzyyf9FlYciOkNA4Ax zDSHVa%yTG1es2gJ!epk`SYVUWr3@R;UY`5*pxWruJ6bv}=X(dK542o1Ob}N-Obrkh z{e2=)ulGtwU$|>>9-2})K?E`e2Iyf4tqGw>hmU5*sqAo)=8kIrJR|~xTo&#K4kZs` z!BOkci}b&}_?logUfkFS4}tD8_pN_P$N8g-Vj-#}a67CRZ7PPRrnfux75o>$g8~jm zW0mj^Z1PC;d?OE;VR%d7o6UQ-om=#3t?RV9%ebnwbf=L%YZ1CDJJBN7#{D`;*%|q7 zCSUIIs7&Ko&CTc`t)4rca1fm$h2!s)3`VI}$xkCRTpWqqa&akp9>|rXHGnIxx0@2d zoHKG-(zO`{A}%xmvZ2%(en^H4tn;5trl(u1;Qg58RYx3NQaw=S$ULD7PR8KDUmWIv zb^TUVQMhJlH=Uez6808-(4oeTnr~-qA9&J#9U0Le9-`A79#$AElBOBDLbpNi4PU$( zs}h@8AUG+KZ2_1S379+2R^&x(Fwaha_;Pu z1p@ja)TL>Tc*Z4ZWT82@QTB@}s1LhzKp%ssMw~jO$g?+RGR#Hv;!r#WOFvgDwb%R0 z-$~73khE%guOnyyENoOt^e*0}>>t5x8h&yHXTo3L%zQ0>NT~P%i>i*g=Mdd$yzBE? zngL_hoFRVQ=WIv3*qq~qSXb$ZlW*f?8=YT#kKgkd*3uaX1qr@)b#9}Ln_Y<&X2$o* z+sK-uO+%w~3&|~knglPLlc7n1>ujP`OoY`A@X(O&Y}Qqw zdX``|u%3K>(_OW$C~{C79^}Bxm9(|yfs`#Ga~bA;Ft7JH^Q4c`0x)$v*v8~{&a0{( z7j-&rGH1%udfz7uhRNTS`C7 zl=!YMq4@=Q03EGTn+Kf&7u)<<%<;srD)LAdI^p3Uc;(lqNi@4(`8ea*>Zj43=IlXg z%VC#jow?~T)EDW8yI^syPJJ^{;yG|Qgx#`~>n=$1&9C0=>bb}BDs8GIt#fauG~JZq z;db-3JM6y2*`l?GILA^Fy+n&^*rzSSOwxB0z&@&*&)z8f#J^?8E0m`&dvHju$}$ne zrm`=C9L`zjG$klr?Oi4B%?E7RNN&0a|3v~7nyAn%wf@&WC$_!#;dmb<{F7lj%p z2B+!}8SFQh()qh&dwo(D|pD%|#MMPW6o@FVUz?w7MhnlIV1Q8?-3tlrCpz(a1C@H?iFG zJ`PCD(RStuJ>pqnuMDJToV@o1P%y+IlKza)SqQ>Qbq3Vk)fT)X$X|c-sbu!7CZ2#W1K%M^?#kcG>D)4YLk8|dyc8^N7<^}PW(5(`-|)Qy z^RYvFalI13r*Z2j9=8WZr3+sunkXGN4Mm1%R1ns z&(B-WE7$^-#gKF}&B8w-P4rmpyuQ@(5}$ofcSr|aupDm_9tNe$qt#DT`0G3D+vqVG z$NTT;8lU_Q?91}~T|kWxz37Iq{b0LV-^lN{R~b=*7P&iC4G_btsQ#$<+&%8>N@&1Zezu!AS*JZ7) zxf}&87F8uuihW!~RhwK?9a^hDI<<8&7vU1N^ou>jwy#3U4oQCsw!biYEzcI8V?q^p z620Jt`V&!ZX-HSUnSJ6x^ZE?4@N|n?s?C;CGSqAc@jK|oVI93@6P(~lIC`X#w6A6> z`4BHqD7Wed5l7(gK-0i1o$IZNR|(6$=E|^J4Nvul&G$qxERJ;-XZ#iQ%sN?*StT!8$%-a2YV^{Srq@jSp2GLgOfzC{ zsY@>LI?H)WzvfuQL=Ck)Ned|WwJg`Da#^$~FT%#S*}auE9VfOuLH|x5GQ3FSR-lOE zi;h`J2>@)_j3(C*kW?|TFx4Q3R*(_j89O~xLgUA)xiMh7Nw27! z{BOP#mQ;U9XkpkhQK^w1XdR+Q(@fTj3Jq(#Er)b>z46eI6OM95Jf>2x>0~N@trXnE z>VKoZupim$`lH#|`aVaigCPGbCuq##flTtlcZT?|o=*xHBtR21))V3MuN(>|Z~9wl zxiZ|fw|=TtOoVcCZ(M1HI689!aa<3W%{AU$>kP5-CKR0N5PgNe{y~NEHcFl7Et9+8Te2%OimJM)2N|>5$f9Vue zdJOoK#3d*tV+k82so~W~7r$R72jT!IZ=dSppj6TI7GG0C|43uC1gx2tHya6U57C&s)WzwyX0Mes_lq0 zaLPbW)^`m~5N2{JxvW~)mDh68r07_BLmnEC0ZY3uT;H?Y=<2C5eT6$AqcPwTUbU^C zUNF?p-*2%lYuQrQ9$jCS4443CYZ><7cVYLVfvkev(ylO)%!(t4huIr{%i!xbuK30F z>y_~qxX*W$5@brPujXyO#phUU-@i$15NZs9|xLr`=t}~)yU$2?g%4bu=O)zU*;?h}+9ycV%6WF0ATxCYQZizjs z=OzzWN2)j+lf*28UY>x_r(tfVANuKWA&`~b^Kl(sB#%hF+-#a`lJG|T)sNWGqCE_n zfQ{>|ia9|Th0Qq(2A8wKrk;Gi@R$9!-e;R5M=Ww#EHC+^V}@6H49%$p#075iHb*o{ zFLhVU6iOR4(#Zoccs84p|pL*F& zoQ9Adt9+p=tA43m{=t)yXp#F$xXL{zzw7HK1&tgc6P zaU1hqgjb?2f7OaeN4^WzIYJ3DuRrH^K+qc?rVJz>vs?y(A;TAMgp zU|uJ_l+rpe&(wt6hrL`vGfXBx?;><+ktHn`&4Q~l=S#qUk0yrRp3$3*P=@d+VB(vn ztFFfQOqyuS^!qQK=5K{FeA95Y|Q} zs$*jAeHBE+fK)b_w%FG-5p|jO?qaAR{VJ7f*#Ajk^xTi>9h9X2rTvXd`i83nl4C&7 zEQFvE4T{N_EK?ph08XT=F7)iazpgBvoiO8LcRb80@+fy;maw&lR!DTWouUJFO9=Vk6}=L+mFa{lM$uML^{KqKj+*aK%5r~Z+?6-3+mW_R z5SpxF@r)QG82bBAE9iTv)h?OR%WKhC4(gSD4vG6Zv{l<3gyr^qAMJ)^R&gByU~5Qt zBsm#VwBO;jdk5j|D(6zF&$ARg6;4N^m4W+_eO{aj4uIvyo3tjEte0p_)I-Iwl4UFS z;{;b-Lhc+wBO{y6*vI}>FP!qDd&7prM1X9uMy9L)?0drNspM5(l(@IrFEuWrGa^q} z;VS*DSUK4@^nc1DroXa5vLHPh%2J;(j+PquT=@lhU-SN0u%$TR*E!ycgii?IT)dSd z1!M)8hEf$c44K}3^1=P9!Jlk=hH2v>wQh7Es9WZk9-#8l!?%Vg!9S+lZ%dmd(plR0 zXi0=I2pXv0R3FsP7o6$6!khb`o(MNvte{+i_k+-C4mM0}EEhFybrq^uwr0Eyr6`&4 zWEx+F@B)cylwL1ShccO3Z5ruT@F(7x>&*bENm7y^551Qoyx95=Ak}-a?un|hPb0H1 z=y*x|B`bEU58#P>gXYb1vf(iz;H2$#)p@AJM|J{{WeZ#z65Mmeqr1E*^kSwaQuc;! z^bZtCt=4-+=Ati*PQG8|Hp{JdE8k=~)c8=K9tY5`8ZuKVG5HoK*}mg#e}LIz^?!S? z625PDwch`8Pz#g&F{z#o7GVYozz+Zhxw3ZCk5jLs4uR&JUtyQeE23hDeJ}lr*YM_C|3OP{2JjdPcW?gTiW?7=%whbT6p zl8YX5>n{lXEYw`wPjF84(Td?13&?o6c8yz%@p=pAWUH$YjFb;BV*7fkR6`!kDtkW;9_mXS86nWVEt0V6Sy-!KD5SE(mbKve&XI$*+KLpvT8N?MgxrUjmtlF^X2wqjpKqXuVdC zNe!*?1)$^wNx?}GX_Amo4f|?wa|{&vZ|*kT3Z{7e!!7-^|7%z=ZHWLep@ODkPRih_sfc4}z4d5Q5&Y4kSTHi8va( z(tdGrl{xr+XM@Nw3;{bU-ZhDV5Mt8mS3}~Rn$#PSbk*|cZtaTc?1M-@JXLGU`E$@&_tZ*SKAi@VArF*i8LQ zRMQM!)7^G*zJlX1g#d|DJOsiR1{=knYsr!=Ui5oboP3`*=J^!&y4GSXtj~b-p%h~M zbXEQwaR4;xx~baWwDAZk-CT#!SFvi`X?>Upp%SSF!MSsTZFc=a2~O3Emo{3rYiC*E z_!O=$OX&-yS$fw~Y6$LKZTi}{a<#TDvTw7d>RIt@E}26bx!VjSDL&(Ga{f5I7G=`( z#tYdmo%_Ol)S;HOSJG4B))>3tOwf;=XrD|pi>n^Vxh*fpc%MPr#9%^7dfTE<>qDqK z!Rp~tw^li-T;#>i>Y5^(g9DVf5b=mx($|wACNH^0;zZw?R3;k@}=R#HFIS? znD`xf0vO0GIU8QvY6eb5vgTqmzX1Y-cICx$}%(CnE=xQLf{CN=F|}TlIR$ z-Q?5#0^a-5ULS1x|xE@L}Ie7)<#}@zsr#OxO@!SqJGh~oa_NjQQtlm^FlG=n!ogy-Z_t^DkOqjfh+ctVvlD6S?gls-v z-}yk=GW2kJh8QC8{o3pJLh(bR6pOvJ@dImMh^}jGydwkGvCi_Bp$Ek83y2K;j`^9P zFyTwX5Z|mWUDkCzUSyWKDwni&s43jE=a?WdTHh)!4pW=I!xNf~i`8UWZE??E!wa0t zr*1M|G4R{o)rf_;vn(%wX9;*>R?a958&sq zEbQM+8AxjY(?6(hS`EmfxbJh1x{Z3052|_ra#OCpQI7h05>qwkfADA;Bj3mibyI_D zFoio9uZm_ZatoL#J64Ln7Mv+EAi|tdl4Y-czvBM=yCt)XpIFI^{c6(`k>k0-PmV&^ zha!uwlaZA4c>oKoAoY{2-|c3J|2;f@b%B6CZ${ z9)iMy$D59ryfSwr?2s{S}K(@c)0|4{i@qrQO54(W=;tQj>ft}so zQ~cA6F89|Halr4(l7Ic>7*Ghq&JGgd`O3rN91mm#{`vg-mIE^j2axSQynS90=OC|k zH+o=0y^q{6KVMSc)YJs+M5_po;p8;+vKgc9*S_Lh&GZ`QMr)qO@ZCkJzQNw`t>HR1 zz|BZe=xt%~H!9P8eN&RHRakv}1Kl;qK-&DuN-*UFW#ukDS-IdQ)p#(AWPJWaF#1cd z;^t;^u(t?f7a+h!9ZLfoi8M zTC#C=mNpT18-1T9c*oyEfy{os@G#ZU*YT6Uco|;Snx)R#;nqrtKE!J(S6jz&W&dz7 zk3|684G-{MjMfbStiMf)-7<<@6KnlDy8ytZEBP#6(NVL-ZNZh8Ob=k!M<1Q7E=R_?ky^@j=)@*+T=%B_lnY7=Y3a zPmC4}j}p7RYht47=IhPaiF7iUmZEoPR)ZgQd*@DzwjqTW<&qohk_W{-m?Y&C%EO#2 zgEUyxCV#FMI<~jUE{Y9geD?z?iucuuS6~u!JKa@)GWdX;S0GY5-Ce*m_yCJnAY?n; zO~As$+WkmV<1LY^lG3H{!FI4$Aa;B4DbX*NsSkI;2d<`Gn;J5AC@r+c_Zd#mCH$wQ z>30DOnD+O=2Qmy*eMGJotRS>h!v~U~3jkf+$rgEA0!UjYQ(1{oQ+5T#yeUGnjU zv>mn(^z&b${JibB*r(fj=)JSsUfsV|I&Jcq(7G(OyvvKda8YyPR#OAsuN_(5(DL`+ z1e!l|UJUKscV3+Df^rK4_e(c`p<7D8hFgAC=!#ZK3nPwU9)EW;4EJ5x%7tO_1AZH& z9lK49b#+}i)O|(ez1B_}cE8uafdY418CG7~x9j12^?RS)15<{5>4TjmFX+n3R+P=M zql#9UuX~K{!#6z%CmQ(1d!|QuhC>&mA5*2hpqn}N?y*~k`&Ss<89PtP)NTm(&ahka zF7C%h-SWfzcG~P7+HOSO4fOTLMeem^U!cT#rBZ8xF=CfJJRGyyIe8Fa9%}8}cu3u$ zsLeDpvLPTaB2=MBo=6+s>eLp{r+tScgQuiT-)iqJZtbRJKOe9s^<|*&o{jt5+b~kX zVf;`!&@QkLut0g_PGKFgfA|46jt9jtNG8c`{(@;?_z)b;Sj&LQ`G7NeK7cT*x^7@0 z+4>W$R9Sc|b95^yx`JeFh$cd z$J@@3CW7@LipSN#Z$|p93tzKjaN?wyR(Twff6isFrXVB5!AXk|z@9kv!iN}+B6(XzW{I2bibGcz-2dd@IQ zuQ%&txhS;9Qn4>ujoV>ZjFz)Gw45RlA>Ti|tV3bv&Tq5mHco?#d50!@{Q@cm#m?F` zF=%b=THZ!IA*d}--7?7I$gpsJd`cnpsyI0JkOy>MX?aR2sMUWR%Al4ZQ&qKOl5_^1 zrpSk0xwrb18}qYqX^q?hI-Zo zNu_)Rpc$(9-w$OUl-0!Suo=-W=rRWVEouk594un?0r-ieS~ckC3v~b#0cf+L;?TH` zFe0fivBJM&1>3z$kQ3B&Zin2Crhct1q=H`~PeD)CagLKfX|vk&oZ71h#>XP58}Z<> zEJhX_2q z`yi?CNf~_R$$VVBd?Y4WFoxz;_yVW5G0ASfY>-+uXV$w89Qe&J)#YW6m_15JnFh&) zB;DyiiYP@iS#bdAPa6rnf$WE_DcmK zC&5(EJmvL{Oz9g!SprMF&#c(aN?i|V^fWoNR3*QtV?#jk4@j-)47LQk2Xy?sHgdh$ z<|AbTCyih0L1pDFsbyu-=?-n8gH$`4C4_QXdg0!1bkZygD|k!QMzWy-CLc6dbdl!C z$N_KOT|cOhTj)OCoabV(CIi7HFI#!;l)}Z1-Iq*JPFQ!L^@Qg-j-RG1qit5-+d$^8 zykxUged%(OSqW11`1 zAD*X`#hE8yQuS_Y0C}(8X zS^OCzSF>;tYBU2wA)wH~qi)B{8|znQB@11Y3*|k?n-<<`C6&P`T#jSdNT5C8>Q19Y z?zT>P){e|5p#M(kofi#$7V#eUR03tB{~;^EjuA^Lu+bzg`X1Nf;ccexr{-W!uAd~0 z|CH)Pj_z8GrnZFDn3Wm`yhDq^LCyaCZ!B$sm2j7Y(KiDduLFNZ$Hb;191OQka3C6O zX3-aeH1!mW&L!AR(&vC(Hp;5)+3Rp>(G(wARa8|;XO2273?h@UvVrMLAd_-`t8T+K zq?MxhGn!aYq5e>P%snpZ=!7&2~ zLqBi@n$|*;lPT|1jv^_wTL-?5GZ3FuUca|ZVgb8rA@Uc1q(etF6(r6KN%t%<=mtm? zx-U)QKg{EAV>O&YV?13oRr?H-bHWX`<{GlPrY90M*R?x9wq@0ts*Svm?Q1 zndc_JzTlsw!2iZZ?2Ts`uk+EpwzV?N;4WRxuV)`dU7G@3xLn`~Yvk7b%@&#y$zV22 z+)PEnR+?fot0lu=&C3t^aH14j zmF1YuuswlK^AhCV%@*NNg1=~w19)Wl0o!h8``hYDIeLQrltVG>TD`&*DIMUzJ_AOa zq>kOX8gaC2n1tswy1vQ4VW<9#v%FzIv;{0v)1ZkqFyp-_A5(ub{<)_AoPPn9gTnOp zz%HjUsa=ha6t*Up9~qLjVnd+BlhOJtem`Y>LPNE4cZgj4!@1&o-s0~$Dsq4)mfe12 zh{+Hb-)R^N35S0VN&vhS*0`l3+MhaaV+$)f^;h{@(#E3ipsrG2+wk_#eS$!QWJ7Hi z-^vP1Z(Jf#r^#rz?CFy|i8R1lIjV#`xu9ZNx~b_OXdk^-Db}q5ewb;e#X30PGcnIl z5<*=@+^2#m?EO}~0Ooy|j*0e-jYSogxOYTgO3tpGz}{;QvV;6|h4fGB)3;58tVv2y zZ}x>Y)OT8Mqaktpv|qm;a&z07h!gdT$X-=>wp6^b?}Fb^LYW*)*e`4rHU~FrP|4kG z9+kIbjQqx%$}gIVV2Cg&UfiF|8AQ@M3EvjXRCizy>FqM`Jb2o`Q z+Z@>iriM5vB&>>p-c3|T^X}V%t+G70PI3OkD(1lj%%ITq1OWK){4-uP3@_K#- zv%9qyp`&v$Qo?zU!f*N;=R*7fB?ixgYSCL7f#p&vioEprM09H#6jJQO(0{ zEdytjs5-zY4DAB@`+H}_sX%>2u58~|}yN@*3_?=&A zaPF3AXs3?&F2Hkc%dj^S6_%~V)v)R|CGs_q(}yWVk2c5K#gcr|5t@Oj$>sV+T*z+!}orq2v%l1M$eq)KJjGAiV;$?~Ac#p&e#QK&v_wqDa!7>dZy zAUIGRX3U`W#;hfC!5iUMz^7X)&B|uSUWmbg-9MO@d(^;rrG0b=?i-3h++CGm1VlvWR<|JpzxUe zHK=`WVM|MtB_xaNds37c{fDMm{zdAWqv@JN!sa2s46BZs;~Qr{GpAIaqoxNoExsxA z5tY4GUK8g)173nZe>FcFE$2k%NEN#b&34VsGTZ6Fls<@z6kv8 zetdBTx?(-z{W7XvjF$S86uqt$);V(riG;lkSb#(l&~dMoRK1CWi)|0y_2ivBWEbS- zQ8QA*zXJFDXtr69^abC#WAd6|__YW`+Y`4!>~H8~>susE?b#R=Ap`k|DEfm)_`Px6#b`BkJ0sbXRV5wA9|w zKYiUcr{Y5`8KAZpHSuo`YV6KnwC=WW4BEokK)Ts6$K@4J?;b9((KhAu*|GFX?=K@l zujG;%_#430gOAYg8e7=1e$fX*0xC?yhEv>=Q0I-i%t>eKbhXidXpCdGp-F zE$lA;;lZ--3319%w9D=qr{`8ixQMca8<`{UOIGj^pjiBkB-)(eL<<^O)opX3s5yyDt)Cjlmq02^}~+-<_jIL3!NS6yd4}>d}V-1^WRqm^&=TU-KPdaI)|! zb$FSNE<^HqANT@<#2~6PaxisAG78)iX)?fBG&XY|LyG(}zoxknA_D4#Pu1l{! zXv{nrd=IRP7kO&W>}*7t*eVFo-zUFQpbec&T_($WVzTvxShmqzcnXRrT`w>PzTE8~ z2D7ls5BK<+H}n0OvOICzH76$jWSsz^b=8eezkdvziuF7-I?bZYDVrk=o*H!#{o7!l zuh!3r4P#G>3v@ErRVnMlbuP2J+B>b?$O7)Ofnjn68=)c}R0Og0VE&lGty~*!CkZZq zxn~zY>boQgIiLD)3UQSv+qEP7){$YVE^iI9Va&v`jWj2CT%I&;d{c*1_ka^-gfDM+ zzCt7xd>7!gzca&aw)`y)+5wdK&YT?w0HqY;y3l{|^4OgrECfVOJ%MM3YPk|aoKM&m z)~e(FN>fD9s5Zt^8#XD$bYoC09CZIhKZkkF@%PIQk`cBknF}|jv;VSL zmebg-1hd}RLJB^j%cd-fzLTNhxS9#Y4RMH>@!4LI)X#%8B-S|PqE#z<&bSvAjcYW> zCAHci$kE0Mv&2*&2NIjZtbHj76nTaHYnLwPq=waTU=^mTA#a*JCiBssi0+nrrRfnq|!(Eq!%Cr~LJjAUIhXbz{xNJRiLf7-$ z@Q-6EWO^~q>4s9p=yS@Ix(Tz zHrvz=C&g$wxu+sKE!FielY__@`w&LBd}CNXVghl65Owp!!1+4geLHK^m$@VpRA-ke z8vA{1&bXH~lC-F7lY?VQhysf%fR^Gdxc#v(Vt35t(gSS<+iXF=sKR97ziC#*Md$wV zj748`Jo+*mdI@?XhuIywWocun)pl~?;O+I%fRybG$CTR^yN;r;V)7nm#MDN>9+%P4 zX#m3N%Rj|D_L5AWw zWWjC-FFac9lBU8_bR0){mJE8zXZpJe5ko9xkj^Vm7Bj+KeG)!oB~ z;2a(qksUB{Y0DfEec~jod&!*1AnHbnxHXZI$%`B;4SOezCbtA2E4fvb&ft1lvj2Xv ztI?5_kC>7B974z$p=k~*eKfK-P;vxKp37@l&|8|>*_S{Zv=v)4;-}jp*Y|>`e0cH+ z584A|h>b?)1z3K0ew zyrU|S*xcWFGyy%YW;?l=mD-WP{K*Ga&0W8WWsq!c$4D{LN7g~vNAKV5hh)&jExVH0B0FluL&6;Ibe>y4qm?UA}v&AdW-3N9tDXp{ZWzDOCHK z?|F}&AKxcEdGu(Z(kgStu`TvvyqOoMZg}UY?O)G7Wr;NB%Tk2blayGqN=#2}R*Iq# zZE#Lf;v<(a^Z7=gJi`Q@WHl?PIzkNKm}XLZ@MUMa+isO^XRF$Hi3wvH z#vknGwKBf(a4(|}T_{IqY`c{H zTc9)Ou!oN3p4N;-H*&|7MGxL^)(iGI?36%YS%{--yQ zcYZW0l0cH_1J1{{`7bc|A`*{!-tRk06(iC6Pq%UqPzT%2`vf$)!lRfgpEfH)GfjNc z*)cjs15|9g2|S~OlZbGp6NQ;f!y%+yEtY;cOLVmX(D09No*H45_$%@Wf(F>VGvln8 z-Wy?S5nJ)z(9TsR85cf?sGx&9rvf)ig+FOZLaqv9N|ij)f~x`Ch8T$yZ&#TGHHJO~ zH<@`X3!FUCJBX?2(z2si3?c?w8jvnhfES9ENknDp^M#N{WhDVnC2aGVd;1VJDv>)SCbCsT+yWjd{9#hVI ztYMTyg?=}4oT_Tq&b3lJ8{g7(VbFfzz_YbK2$$}7IiqsiQ@S%(C#g>7%}@4Pwf|}7 zG09zU=ul?w54<#&t^2$U7w%+^jpI>h?~t9!OO1(6|6RV=CH6ZN2O=j& zh}k>07>C3U?jOYj6`R@)pu4%uzFH&RXs26vLbnycwo#-jCiQL}8FOTG*5+!6iaIL~ z%IA~=)5un(Mg0K<>V_4Vosg{%W2Fr9N7TS~{5L#*QFI}UeC4vj^z0(MFT2x-;oW9N zn&?TM*CSBWlae{U7{!ngT=<}X@q-Y_gf$Ks4gpcDG|`Geo$eX#SI~zb{Zo8h>6`{~ zyJp9d)uoB7p&46>*+$%afoY#XBma%k-!6L+Ph+UJMmak!QWQk6g*)F!rFt+ z;$exA9c=J!N)AL@*2CMwZg)QoI$XA)WJ{jEe<$Gd^Q@5GuH|y}=ufxgV00}F$)Gh} zsX9n>XHrwwk#%<7$yd7DaKjU_0gzR?hc9GFKyC)hb~Sg)nhxPT5X(=JRQm!5<)&bUNM^FAguo8*@VO1=7K;sy+SwAVe531Xl!qYF<;}b_feH(b} zcVXMv#m0Lc}gPV;{!%Kt8*5TBsX5h;=?fzd)K2UYT^Rp9PZEF*=EbmNJ4Ko#KXVXt>nCCtIZ zgv;7hHMD6Y^^UxaWO5_o`|aD|_Al9WDkF~G*dlK@lZj4#dCp7U%xpm3cxmxD$@VC$ zBkH!Xjuw!=a3+XDesT7gT$+}VWoELm>^XEwL{SQlCPaYDik!dL|K!Dmj)q!GSMTCJ(qBirzV+l+Q-E$>r%{>Z<1rG~(#bYqEQvaxRI#>taNoT_jzZ`z}q;!}O+ zne}{BS+ZEZ0f1(Iu4i4L_A6>wljw?EGc#XN+RC`IaqUJdxra3R}EHsUZU2bE*>slEtKEL7(#a+ zc2XT3YkvQaCE8}xtt=->!bxmPfPM2e4FhEEj4T_*qB!EYA)Pd zo^!pWta+!`_5pCmk?In;hJv+3QNu~!;(<$U(ajmm{D+%7Ql;R*0jz-^o7%-;uD>@6 ze^8LNHA}YAAMIQMZ(rXO{n*X>QDh@1QjQkq0kKo2KR7a5Z?06k=jvLaW^h8`law{L zNE|7As3?eNUzuJYrsz{XN?XL-x}v{SzcEJtx}o+}w!ZB{l2x%{A!#p9r>D4bJI>Tj z7v-qXEs;5$dRd@Bpbx{ek^}!ZJ&d1^Kuq;tPHaMCzSijldQyhdoza{zd^1(cJk$Ku zLDm2XwBcuX5wY5{z!cw;QRRB}gaPRC6}otGb;s!ZKzEgSSh{;Pjj7U*%@egm;AWL^ z^>7^Lg1!LRxE`!Q;|G!~4d_E?OAaM;1}i)|ARUiZ-Z5wnj&yh;V}FLszzxc$(YcZC zLDihm9^ze~I~#v#*X0d$1G?Lhg4gyqf(w~_-_F0V{EAk~7}WU|Zu1x2hBv|UV6rq= z;@}ypC(K2;N)!+O$+Sf0NH4Fa!nbdS=EwgRE=?vfi^2jLVyF2b%m z_F1Zf_a0>ihF2R!O~EkX9DkI!muVe6#d4!o;N*#X5cn+OMhu8c5fJ{y9U5tm!w5CR z87uVV1j;mELZ9Z*bP|2mmm>iTZHb@|!1yJV3LmE;!jxl#E@#Fg%(!Q$O~TKJLj)sp zZ%wARDo67X74IvbXx@x;tD^^t#&<{xyG=5d77@(8K^GF!Ew-G;(_)vCAGOSjdcXCAcK`3fYWX)nG^r0H~gzhOJSc5d`nrLB}5B}C=S!I%t1EjUS|X=G(GzgY`~ z^Lm$(T=o0Cs1xr(%+n`Z9N~0HS`-?D(Fq z@*ZRnh%b<5C#naT2eKvMUkq_$9-584;-;2dIKZqbWYt6mI46ja)+Xyx%jj=UX|Kw$ z+Bs*iEL_VUi*K%ZY%fqR+!75YVrqh>_iMK73f&6SqzjgNE$eQ^#-2t7`R#QoY$0*x zKc6Yui^B)9p#f{$%vG4BUHD^6ys`L2G z`sM0&3125IF*8-rU@7qpf^*2Whvm-xSS%BsfGxVYm>hT>U-ubR^_u!Oa&yNYskgGS z)=~TE8vTjXuZ||7Wf3(opFfAs;zYm*(M2TX zj#w6_njfbIR# zVuch`5eo`p1o4i<7}XxghB&V+{Bev4+SE8hzQERL)_iuvr7nk6Y5H2ZpA{JjNpf}y z;WBok24yiscQ#=_n~OBeTKqJv)f%42UVSu^#`l&khOFCM1)!;F;oIl^HXd;xJu5wA z%$6uq&EVfC$}ti}kSJliNQox^kx-CoyB<7 zI;y$#0woHoSFF@@LOS1z>a)7hrGtf^`{&SebX2l=N`1xCT<7*+=%GxiB1RPWGYLkK zW&{yN6tEWUBn;ZU%MTlR&XpnA#y=MLzipoPkyvMKV?kxFla0RiDgqAM#7s;S!LkRa z4EU5%QL~hpA1;#E+?Eka34a?}h4IxjBaj&8y>?DYs=t7WQ=cg&c2E#!PqmoXZwDwo$iUh8 zy;vJG-7)@H^#H5Hx$GM0K}Vm?05wF9Rc3*VCpZPV#alvAD&!Y2B^1YZ?5oP7cqDs$m{TK!c|K{`R}_gn3{dA+4MNKkxF zpC$bl#>%R}Q}^`aZ0D-8(MUM|7bWkh)=wWGuV2bblT3tm8iyhyqAm;$f2skd*(R9Z z-;=D*$dep>af&)9%b}aaMhdexvhCjeI&LG4nD^vEfN64Do7jheXBUw%U%|*>8Y|fH z)Z=dpNOo@UXvW5)V}Ed<^ZjJ>z=5Js>cFX&8M)UO!XFOo82L2m;@*2%G&Spw8e zSieE~w{Pi+JJx1Xx(J(?^zw5?suO)(?^#>&Sp9+A#2IT4on+8(?hUUiKME>tDVBYX zhx@*LUF{LC)z})_1gt_f-7Kj`d^2Ov?c#rXGYvrSD(34`x99W zq(7~p`}mh??ZDccJ_={jRF@L5!f|M^`eCye#@gDGhtLKciLdIbca8>i0qtD`i9@*1 zBF|TM&Qwj2UFgakIF8I9aq_Mipkd$|d1?% zPW#*gt`iC-?Y4XoY%b3)HPP=C)HL9r(mP+ef$5_QI4>Ru-Ymm*lG@sjhx(F(*TYrO(2s+>27QE$+7iYCay~p4=A}ga**a^sUJ_Q59-ZZ@64T-~Hi`ypPjnljk z$Xh7>h}BdWwe@8E9Fp1n4wi^Sel~t>SiqLj=WpR6(xCnel;2VQ{otN_z5}{vpKx0z z3w4lx)z>WRVhS~TW_-$pzT5#Q2m`C77{DfPsSgriaR@QQ+kvqq5_MXe4lQ-9{c2}={344KOzUYJ+(BYN+R58g0gW^54?q1LuVGHh;8Y7s;17N?Q~s`PrcN0 zW|-_&*v}97#CPC-o6^}2S&4=(a*X-cFiL% z=w0?s3RN!KeIFH7&Xyrr@JgcPsp~GKVg%IPU3xQWI#St-R&Tr>KQM1~@qHAb_@PKk z2?Rzm?*XYN2|#>_3c|&^=AF)xM-dxy8qiF0m}>sybxChD4K$Op`0vFdT&0hsQiVP~ zI4ptRkr&kykokk>P~K@_>J%_4AhX7)M&G%iTZzO2Hc4CT7g9Qu^M@!m3R9fSrH!YA zvVUknDj_;nY&6^0={fva1Mj#pWn<0nP=0>jYeg$9#v*Ucr*)9Uurt5%%*(iY5N~lW zwI!_-bXb;mcR0ms2z4(tDOoIL&w(7ilDIv6cZ1}Nc=3>w8wDoOFo_rKcm}Y~ zj#7MdZ6=nxS5-#<%7`Wjd2spM#jqsRq0Bob#ktcfV52qA?)NcnSc=oH)!-;b|9Iej zN%2Eo;J>ZDk8tP&Bftc396n0+)hRw^x!UKwt1CEqw z`cz%d_K64>5`|I#@Q3fxEJ=X|i*yZL_^%!>vM71sP^2}(5g*yp*rtL@Y&Tr$j_mnU zY06adli_Q|1FQY4E{p*t^Y_`U4Neo8aA6~9c@)gBGM$5=WmBF77;Gu?LT=_}3uV)A z%te$=2GAxh^U2FjL!=?OIFMT4baog)^?H*6n1ZPGrY|S4kM2mb)6JG4)P4>g0n_vVJyIH7ZSo?Yc5lWAz_k&1&>9s`CEX*B#rWj<<1Q!UVL8EG2?x9LQg@BN;T$edT3 z8Cu6diGpf=2iaZs;2;i`N_--vUIU*cF>sh&IzWy2CmUMsgP4;rJg(tW(zNaT}DQ^9@9v80&}Id@F;v(?jXz}q(*`d#v3)sSQUC=z7a`=9(ErMr_;wl((xAbtyqEoAk};en{$7E zB!m98uW$;#Z?u^+oMekW7<)>OD2rNiV`E^3fndb}W0aBRZ@S+dLNVuj9l-O^L#QFBN+mQS7`W2!s>X;O^7@+2 znuHPrxG}#k0?Qb@cLf{HQY~xnenW6O0qN1PykcW0fxTs8{zz*37MM_hu{k0nCaByO z@Ktg3K!;f}+upot7?V&Y>Q}b7?ox_r79Uxxf{x(b+xHH(o}^LJrk@kGBoE4KMwtYy zaD)cVvPLZzy`58};DsFz;LlJDKeN=LKXjp+S^WX#ZW!$3yPf~g{z{3$Nd+mih_0Nb zJ)9|J;xi~p-B?J7zFsZ-@f0!JF!`=7OsC4#j9B);7kv!2I6}5fd(QqNglcYN^^Ajv z9EEC4w|NpY*^r*gUMe|pWDubIV?*QL1q0b@_xU^+jhpEndvB4ydm|m2hxh0% zx~?(0HSqhMyzQ}`4ywQaSE}sT|BfW1W>k$OHi$*l1j_Y+-!8@uulDFSEa;~hrxJxG zKxVe7wXN}x%B?nQ*`EwCKbPHn3N6RS{(N+4=(1Bd;ZVz1TJ|oLS4rubU&L?lW^zZ& zUG&#w%xR&@agZs13%=tGaM9}z?DQctxl#fQM~%j?T24k>4st;2&q+^z5R0p=F*#5(3@KSKSI~s!`|G)_8-)>q__H?Sl8C#pVqbeUuO3| zrIFSqrq1;C|2aAS#|{54$m>Aw;9~dB_Wako|BdRU_n`MQakTqy^5;L$=fBCH|96MQ z{~d42&hifx{Xg>OzjX5dc>M43hmrk%6q5b>I{&9`=~ff<%URh7Rl~(?UcEe4SKS-UyS=;d{G2CgM%t^G4?Siv@`;s}4B;99 z)BX$kVrt}UrffcbK?z*!e$nY@!5Nw$AS_N`{&-~5+WwzFYi=0crX9)#y$E zF?xCVz4K$y0lDPCARfRs8H*Yi{ZbF0(FO+mv5F&sp5*sA!maKD%JHSw>frpc8^eQi z@Xr956O>s1&O@i|k9pPwlFCDaj%32D01M)J`x#vNp{2VrwX3xT+{$^JMt$d&rT-e+ zx$r&EhoE3rBLvL@Z2=X;{kc)a+v}~S0;UzPh5upyePU^f(17UF)B*vb{?iTW7!U3N zLaQQBgU$v-&FqD-*Sq2$lk&3)>TLbf&__&CLS0YE4Ke=Lby9OihjBO!D~9yBqZ#o0L=twuXEJ>5bdg25gELmW~1a*zSjn(H9B_ z-aT}a;L>OY4!1c6*F^hSQr?4UP-0ItSY1P4Jhnm%}V zU)RQZ!D-F!VENX*AHv5u_3tmtZ?lG0u8m*S-&28!SH7x5Kdj`kWmXiYDj?rhP1bWQ**}4}b2Zu9 zeKmiioIZYLj`u!t$9(A<>jjvuW3yAF* z+5zrY+t%}Bf5PkqHXWM%fqpVTKKOC|Jg4|(Q#vEp;P3TX{+aCf3Z4^k_<_k^;rXtB zn123*$&b+Wg%wol_K)*j{r0t^b#SOFML=l+LE&BkZDN9WBr|#XsEtOqMLczA`@Vnx z_I&xG*YeL}Cqq1iYWnlm1Jr^D*(_rRMgNleJEP)HUMNzlt&FnmW%iXIoQ7mTJw29+ zeTkmkoOjvOfatEs8rg$8v19o4qe^BRADpA?RGRO)3;Xu;y)DT?_n6}ShC@Lc+}l%2 zp~>}m7>uZisdp8R$jzeu_4-^bh>E7fH1Ckgy_z2ER%VZJsN=(_q}jTaNDyn+38C8P zm_1$WQIz^apZFV4?=$VtE*r=h^WX$$+alaWJ|q?Y+wkqVzUUGVy;;uYn1iFwylH>` z&Ik>;YKng%O$#MoXZ6=vK3AxSaYW&Zkp1v6S=JWCH%iBvPkm#KHR|Tp2baH_RkvQ^F}P_ zjLj?cNj|G76g+b3*QM-9cg`5kgE4))(pNld4uaXG*3MRhtV9EHE1Gq}ni4AOWfxu9 z0p)kQ`KnUp?+Omf7ZfoZ%2u4;Cn3MpQ{{3$d*mQJ2}brS`d<|cQAh}kD2ahsmJMw9 zcp*JCN}s&+FNab3NN`XzEFao9lcsT6_ilUx>uHU5kc8TaX=&L)6%@%{n9bvmBjY|w zEIfYd;dXc5cp*s3MY|>x(rfvoCw_&;%@!1Ki3iD>lM?VN6MUftSt1Aa;HP%ZPaT0% zhKK+q+ny_?##FN37X0K)(P1+=4rEG{qo^a~?hVekXB%NjBKXIX8yAJ9`m{1iy)%M` zG$(_2RJb+_bz>rHElT%?4EG^Iryre+hL23MXbK^Rd$m5%l0rZU^@okmpDl*X?Qm9t z{u0`27;fr@tqI@D9zeGRb-Mpg}0t=q0Y9z-D5)Vp6;+LwZ?G zZF?YR4Q)c>+ze7P)>(5n3ydx(eM8X9)O9(~pi-o3BV7;+8q5#ytGE%;#)bxy(r{^K zoRem!=g@3U>G;C;5={4yd>4kV?eP<7$Qk5>Ir8ns``FKF^pWJG-|wc; z(hDdRG^Ujc66JcW>v+^2d=&IS591X}Yz?*Mpgq!6X|J!N6anJvqd(tCO^o4mtT4rW9sgR@r|CLXOsZjhveJ;T+k{Xyx-54VLM-(dd&-egj~iq0SFgOoHAT0p*7FBr{YnHQg;c zxR?`@woA60+Qw&fP-lrXs+H`(>b1wCa1EKKMARFed*f|TT@xXU=trPw{uh{=uRA_< zmMlG1@H0+4O9VDe6!gqRf#qwguM&C(-c`CrfUh(b9#+n)LSQtDeL9^n-)>fALt_gf zR#MHGwCP0$=e&>hQo&Jo1E`Vw=+pNfsunHAl62=-(y#@u_sOt-r0^JS4z;X+&0O%P zRoE2F$2(g-Ed%95*MC*4+Xxk%?roKu;kb{FfVS2S@{i}yeb`)zYcfvu(s)pe*Jq(c zuMy{7J=BsZ_JW2f44g>Q6*MSIJCuRQ>2RH4gb28H3vk3e{{v_H%D=~==W7iCDoukv_2u}vzp3gzayp3Xk&N4tZ@9VMy(o> zd^JdSjRZkn>4fmJDaL6Xk1IPa2aD)+D>YyY)CJBySfzTydYewd$pTC4z0ls) z0Ga%3E!`Boqsutvi?)DaVdl`P62I3$4%XD^vt(Tp`h?UHHiqW&f}pS4q;}i5n`!$f zK6ThrSOYN|{X--}SQ*_%dF$dj=P$zBeq?91LgjtOtNioRZV`jbizfzzP=2O}N8ORN z=E*X!4T9}qESrIz9f>PxqPp6YZzrl`N{=1US{@`PM>d+&A_c=*DoY!KKnHzvfJk6%$ltSD=2 zsyH<4KK~3}96pkn{E=B{4Dn{ylPi*j7#^vL*K;* ziS26vx#?`iPG1+*+w;j*jfqoUa)(rpXAda}C$9g1Pm9n>WmxnVmcT>o-jO2);toCS zWHjMI*WNi=#U5)eA%}&SMdOSPiv|9eV?ZIyOHHMj4O5$be?8Qxj0j!qeEiWR-7;sr zgs&RKwplB-;SzuY;;-#luh9S3qqQt~z$2OjlNCCub2NuJv2L=+8XJ$tue640i8kam zDSERX(Y*Me*nlTAc zsqrMu+~~8ZH-&ZB>w!vUUf?J9BoqLR$m250dP8Ayx491}8bk05wc>uX@>B@&+-J&& za66x48IR~5oZR`w&e1x%yawa--M;-vP@^X1-ffDG-|dLGLGxGR#YO2e$6@N8h++_yFakJ*gLA!W$ys6TuF@VLoo#)N})aXq2H{}i5)XbnGv=cQECVz>PSK7J|o zr`D@6iPDGZ7huv)iuI0~C|1OfZ7!v1;m^s32P<>QG@XVnc{<*T26yE1(jy-LIxhe5 z=wkx;wjmu#6b^eQfButV7^dmr+^)EVXM}Ecx-sX-*(<&LuzZ5L8jdM1VnDD{k!Q7P z&XDsYOT{6qLKEM~(wFtTUJq)$BHj8d>Bcw`gpO^pY?3{U|9rj=LG}s6TDh)uEf7GW z(bRy>w}19v^Fn1!FVWm@bzG&=2IJ^$UpUY7fqrRi-hVf!TrzQ#2Vq8KU8q95PUVA9 z6D3l`(v^^CJ92Sk&>HrpSM!d#YAcfEGOg>>s~3!fB`p1mtI`S9th_LwCBv8y5Dp0H>x~T7r^r785_sG7{SD!$M*xBs0w0&y0NsA=E6+K1=0+P0P z$KHCOunIOSb7&q(4a{iwFqLE8_n)grl2o|3Sve`LAQz-Pw>`iwIuu3?hB8(k+(YPrTxC z!D=o5Ox|e`*r;~XT`Gu^u@chA;YP|5_2_CMDDcD${9m|hn*E&N0Ia2F6$L5aLsJxi zZAckAV&w+2{A}D7hsUNBI(NF%=oqkI z>WqnEclIF}aY@Ik-hn^nctc}*u3t$zN|SGrbQ&q|Cy~Qb4NKDpKqlM`#QUZWJiC!R zYh%2E^Qp2h}4L2VJ;_D>g2<|lHOo_%u-F$50=v_ucYg_8zH7=_T2mq%m8oke=$ zXqVNnH&$9!*gNnF@iD49oXzR__ocF}DS6 z?aELW2hR3gOkn}Uo}A-lQPq@(_pf7=?XgR4&OS4mS#R1Wz6(H#5X!t~1zZXx<~S1x z31ORbgjtIUxbKM?8#D;`@wdjTJ`T+*8l+v^8Cc&;y1_0X^7~KnPcy;JF3)O_Pg$^r z-ub$oipE}!k4Z4f6k{^v^}K!W%Axi2yS%zwXu1e6F3DT8D&>r)ZC$YaKw~p$1+mHltOrH zOZ7>dv2n%;&O9->Yql^tH!rkB@3EpwbYAHu$aT!Q-FGt zTzeY+;V1DRATF#c3K;xy1S8HV=vjI0YB04uSe{P&U4fAHpjys}>BQh>vqgS5<}C^A0l zT&UZxry9G$y8H4kvb+k3JUh{5rue!Ih5Z-%vWV=zS+K-CPc4U2ik}3YxpyReGq&Mx zd_mdQDE28)W6}O`b1;ZkB(;x{-gxl+N#`$Q+c7ewb#_NM9y{j|f$Bnpk1q~$T#`+*5`Cf8w8or;_XOy+vqEk2R5H&Ta_G0D^2kDq_# z@Pvj%sm`&hk~6k+dWRdeZ<~Q&etp+lys^S{GAFM(z8BE9>bS(or!zM5qNBt>+kG4Z z4ryA^X5GePo6OOT@%Zg4Y5MCUHl}xLBKZ<*#MI?V#?k!jIif-uW}YMAhnj`20IDs> zZ|~~m^3dA?54o*)M%8T@P>KInZC@D{*S2h%1PBBT!8MILG<4Hg0>M3aaCZsr8X$z= z?(P=c-8}?%_YmA2UT5#K&x!1F@44^3$H$*FzqxvjuBtWW`fB#5(l$HlCQ5nv@cDMK ze#2*}EU-2vIC^23S_$7%V%%G8LB$jgptyFgRTDzD;U3_duMt(J66n`RCbu5e$39yy z+q#_@=a_iLQtOUUGu}*9@9>z|+94ikL9^glHdZ%Ubru8c3PP@qsp~ zs|n|B3N$U){nV^(PS18so7y=DNoAVptau1%g++AK3#>?^NTrjN;v<-egJ2W9qA%X6 zt9|p6J%D9o4VXfT!1s!ilkJdrlTsg24P27b-BW%g@LFuEOS$k#=ljs5p}>ixP9scp zGBf0lk9PPo?eruzHmEK;bIqF%%J?`P4Dl)Mm5=B%M!kwcGjh>hRCyn3Gg%=fGS=&I-qn3TnOR9+&|JaOO%;rcT8#da)8HFrh}r!-=Ou(MF>x%naL}4f!zsbmQw9WJn}<#?CGQG71OzY2P$y zdD`*0gyY6nY;^N5R>)a$UJw=#M`tbq89N}{xE8A`uMb;leA!s(qBzUnuz(M8q8K=1 z;777!bh(E$A3bZ*A93g1NQr|G?x<41gHkw0-7d^@+_%}*g72RyJ^vcXE= zV|EedVanb8)@i&^kwBEuOioJhY!VAW&uOU`Us6`don+#I*P$xYCzQuo?=x;nTQ#8O z>}_P{NX+Y(Lm{C^Bi@XTSp;l}VMGoN$Sh9OA9GMFx)3KfkY3k0g*E1arJsye>d|hM z>{w?6GQMpvd(HdAZr|T+kzfpSPwz3rRzid`BwBYa-j1emT7}7pi);FsjZQG?g!U(u zCGqUTtjYv%I}ODqJJbBh1WsRYsQKhXrYs#ssIAC()PX24+yK>>yuWJhU?BO__P|VM z-)19qojR3I=0c5$w%`Z`r%O6Msae>|uw&RLKs^btkh`a`e%|^%`s;k@0xraZp2l~O zGs)c0ezpw#c2FpOm?=(&AQ(;_50hIT~@Nbnr?(GedR(bDnD7>6aXNxwQ(o?LwEj`J4C#$;|2so37t z5Nb9SprgaS6ftf@1CA>$)l*CQu&@ zx0oC_xNSq7x6|`+hU_r`dbPA{Po|^v(hQwRJ-c|OC9WQF#ZIR`iFiDfE3RB4$DZ;a z5_d%DdeDL1I>;YB{agjR7`5`L2cwoOLX~}Tb6mo`^($H-XpbAi)vGOx zeXWNvUW-Wz97zso<`G&OjIl&6&R8L02=5J?>c@Edz^UH8cD=ew#(8QGO(i?eTe`T1 z*sBa0TEe0c77bn_E9b{GSbC`Euc8giUN4I~;(TuHOqUb(IVwc=8L^*CmT}^pCS7#i zAbH#^c<$=f7#Jn%;#ptvph6n8#)6*+UteD})L8g(6Fjg%ccEao1udx9=C-GZp5o+= z?knaYFIu*3roaQ^=0WB)cZjl^ow&YR^q88BGJ4O=K|@?t=6%!fOxtk@lWpXR6} zmpz5F=kQGQ1Cu^YfQcyv&o&B4))FRYRN%Sf$Y~OiEC!`Ogi3`?0-#p(#iD(W+juwnO?Hd#n%3??qt4zXC#1R zh%hlOawFvoZ_RA=`ElT2V(YBrGhxFL56h@A9#i`-T`#l0p?+qxt+Y8~Xg?wYV#7(I z781R!0~W4S`z?7kJC^ibdNG^}J?yX?Lx0V0o0zbbI4bGlhA0+t+gd8+o=B<963KlR zkSADV{t0tolt#+UG2EaD@xvo@p2frY0h~&r*J0J4N)l2y`W78K>Pkw6vQ?V882C!% z%wrCcc?#Om51z0IQp;%0cwM%>dg3Ca?r;0Ci7IWOfX9r>q<_cejV(1xqcH((wv3q{ zRQIf1p&Abol)8NUWeSZdFP}Gy?eaXe=48{~9Q&oD9BTWO_!*IzdjH$7OdT{QDP$?t zuE*7%d2dFXJ+I7bc@f-EO(}|+QyzYujL8Fs6jr%M5Of8W*eIhqFk;$mNq)hC-G?kGci#Ji`%nTy~;muLvwd z@3s31lpE96l>=;E-}(kS>`0ym!eG|gV0LHTs(pA|Ca^ET6Lutz)`~x0D{teE2C(r| zxpR~jouBHbK6GH28wR9pPh1d zx<>QxEay-e9q{^tGk)&`EG=1t>ngZ}i?70P+l@>lyWF~tI};`Y#$q64&sMzyq?7BeM`>N_L%Q=A6I z3z*Q|Mv0GfvzMkB;`Q#W^`y2L3@!%6wVpnV(q7@~M~(|GuewCZ4x0`V#Ur^DR)u>e zv-I{-hi!A|kd)zNR%iS8aJob_&wJ+dfkM9H3ij!*c(R0}Y(-eN51C!Z(pUOEvxin$ zvm{Ol($CQ;b=Q?aN+Jn~yKMPzph+Y2lCE=&K|J%zfED|EmxD2}w4OJ%C=6&4iuSjh<86sEo&AW}F5K)SjYtk2av0p3CaMDEIA8zVI zlyHDcjuZA)1UuRMkX92hy0sy!J5g2Ta;h()a_4Bv?_6?%x5T2fRr_EEW)$D0xW(w$ z9<#cZ(#))FU*L&-d+yGxyUD_dn1KS|;-xR;Ae+=wZ3p5Fol+EGR+!-_h-Zee-Ib{k z8z*Tu3A^V`e@ocODL7`~G*q_iQ15~~_=8esih-BjWm?b-<5r@8BwZ!)E|nODI?$xC zyN!DKS$NLlcjgQ4w%Xs9uXhR!&dZo4)$#>ZS#|-ka~VQ9`E7;C@CvrM!93FvwiV(* zw_usU%RvmMbvIdyUvV^c#O(XWRWqafFhf9av1g3#vx6D~w6RieM2~QJjjl-Eury z6EZlYD$eShA}eFgJSW=%oQol6v3A$w>(*WljNs!(X5hZP(8Wpd4}SK(Ft<%Fm%Uo*;T}BGLQ1sZ2CylcwFUOtHuQ^F{VWkD`qMiA@TKCEe-!uHx|k%Kxj38Nnc%(2@kd4oolB-F|nt{cI$CH*&|RBY9LzU8*1`TKK@&% zfb?L^&dX;9`}L=lgqTX`PMXOs;>xPS?*z?u#^_r237>$=TT!V{9#tS91@RQytupM0 zUr<37?%kQVF1cEW8QqtZnS{iS@*r!{jS$Ub z0udHIl`lT?_*^4G_@AeQunqt~h4q%(uxkVLyW~x|3vr7$F3aRSnhN2PTT?Y*wt8LC z&yo`^-JC@kH4ZS%!8m89b77XF(kV!2YBJz(@r(p)EB-#z=ObK0=(&67^A1lmT!yrh zJ9R}ZxY7w5#fa;2mGuObn);a*WQgDdxhzCmJWm zmVobf6c4d#M#*aj$BR8`9(BDmu0No5hN@X{^=)clONw-!#p35O6~^_aUzB4J*x1G*1Y5^2)wadFxhBot-?10r`@69Rl^ zm@395a4hD9Nnu+&)fAwObjZHQdM(d+%AO1<}c;Atvj ze2m#$CTtt4;@}=QzSx+ydS-A^08jIC)>n8%pEj6)Y0`DsmuK}n^YlDVj;qan&FLXZdnjYy7xPtADlxrB~2 z{o%D60x=34(;FJoITiwEN6(&dtghD8qZ@j~3Ll}o@x89D6_s0Q#GURWA(4d{Y#8P0Ji_f0NOI?&KmnDhuWJ#>)!&_c4^B6*9}M#Scsy=%$VoK@_b z{q>g+&BOj}@z7_!bmVdjCCSUMM((LhSiL-M(eK2!=aHU9!E%R`Tc0(y<*`9DR_#>z zuPo6UOqxw9Tl1X50&9uy{UHZtt7o+~^f9-|rA2Dmj>@HDok!szTIb=y>91-EzMZt1 z%BlB~`Hr`xUr%z;`fg73E$^mPV{nh#eDi^UyDWcjIZ1MKYcURo1UE~DhSHO1z@Usx z+6QyF_~42Arg}h*3goHACct}<4$qk~mcfCZ`R?G-qte)-e)~9ON z(i0-j3O3JeL67#$bO%lK+2gGYc^W(wQLsCxpplMUREv~^$bZFLBqS@&C@6~9YuYa5kl7z?P$PL2}FLoD)uSkAcxz}o{L$7yDx5v1S^lZ zN6F3#^ZUsghlFic5GWcc))C4}PcPH9L@``s&G&~mt{kYAdyHlVa4ynrt`^Hgr53w7 z<{s+qGnQGGR8f8;Wwd$J_!ynjtQ)VQ`pahQ*M9fNSo?Q#!+R<9Z7>f&F^fdm?S_+v z>E?w^dT1JFb?1#u&ywHEt1n2col{3&d9!$G$Yp1)XWYKIV>xu)TN!27<`GMx;>u9v zSr|Ikn9FR~+Hs8LKc=c&s!4lEP9VYmqB)Es{J8xpDrMqynhk)3W(vOC4HxRZ&0~!U zB5%gW!sz#rP5YQG{Q*Db^lVdcE=BXvi$)8$%LSMkkE+1_3_Zfzn3{l z&yPm5(EtgE;*PK!)!QX}OPEh8XccyFv3y#M$7z9+5p+~PiH@`biNp=@Y=SqyM}B;) zBU&5}?eXfjO<(Sg0-2F3m`y}u#;^DiUXh^I7lQCm%n#6C4lqc45dM;U$FzfsF!0Kf z~??I-1B{SefB(lFmGXjLxv#CG|{YYAmEvkk3YroPm}TnRj>FcVP)@4fD+_ z@B(@2;QSsj5(m;{X)M&QrOYU3k8UMBd>rf-Vi7@xA_(xzd>A3uUOoK;Hyzht0T%0- zNd|YKGkQknfBT<}GEjPx40H}V?f**J+ zM{5Ja%GProG^?=G``JSn81%mdGvVi4#f=d>KVeRO6-N_B>CQ9mbyBsDGBpc#Ktw8f zI56mmw26XNP~6msa45i?^xq#e8$lG%RRf{nbv1ua{L8l7hE6nu9P#0P9B zN4KF^UoT?_7l6LhtE>~fIjeCu5u0a1-RN|E!- zZ;`C@2*1Eg6srW1x2}*FL^caq7T2(%-2V-Nk%aOXa^md417w zm)ly`xNacl`$`&*Vp`)FY@Jk19TZ8? zYZI_5K?!?WSh(j~7KFv?R{K5AL`to)Z;$9y#x5asl2)AaXt2XkkQq~GOBJ19?J@lf z67p8dDzKE!ez5`*$7`>-K&=e7b9;H-{a!0#s}P=f@gb`OwO6bgpIl2QnTVwO%8FI) zB@g8u)mk?%usH=7zMtjJ$<3e05aHMCcGi21t_<1Vy(4_~0n0uERg$R?-2MbDj(;qbeXEm7x*rRv68&fGVrrIzcV-wC z;yHKL^?qnn%WDDMY__nFwLPln6*B=eJ4nf*$Y{}PJo9sz?4o821Ps(rH*ekMbPVJ7 z@8g`{QsmUp4L47Z5pEIHk!gqQ7xG8_^^gm)u>zBmO&eXfj|vhwbhTSU?v<6zPFG3W zId*nI0X)EM)^zqVmIue(U-S?^po`MxQ^X4HVAn zaJ-CrDRotMMQ3&ZR3b@#n%cSS1qxvOP{V;RM61T5%2lJv1Xa2nNtll z$Y)uyP*(#J;3t!eChuYHj%>$7W97FLXURbW!96D}0=lA{+b&Hf$i!z&wfh>zBXt|( zEeQQjN!MacRP6=xi}M&B=QY+`39p#h(IWmD@66iV23TI#h#CJ1==%p798}(a3FjT|QRTE2Zy>R@x&7nPt4hQ7^yy1svyfxOiHH<}Dqs3yQd$K*fKWC-&r2`bBd z*JS*CgA)c8;6FN?Fmr(ZHyuv?smNf53OD}Q|LpG-8G`@RWc;bfkYbW%lKF=wLxJhL zD&sFj#$Q^DKYFk+>08-pL-iO?U55Qnjnn?`^%#0OW`9aDeo<%G=)X1oQx9Te?8sze zrfp;NuRd(XW_mw$WBWVZ#&?;WUS;s&XvX*Bk!CdC*3bUbIpcSWy9% zx1^U;1cDWSz-ktj5*F4?Sp~~4Ou6)Z5`)w-LzWY;G%SjipCuXkBwl}(JAmFSG5{AJ zmKdAyG%+PHPDM5^XDK!!<0)aot&H)?tPk zDotJ-0_<584=5F)RpgOvBfTCEg~cPjqRL*|T$iWK%*?jhx?YFb6@#mLj~~p(HdhFy z>Skib(R;xj30%Xl0ol-|(H_|#iHXYwV7ea5$p#z^MSaNOBiu}g4_A&4VBA&?kBx}i ziH%8)$WJvjN;lLpLk#Pox9|TXBKf3;Ui6E!i1=rk*Kh$j8(FDbaJa5ozF+I>IX=MG zn)0t=(c^{2#wk9NRSwP9dMAytK|p%ZB2C-(N%|vX2m}Pe(!%oL(qRv(5n#eqOELq5>;*O3(t>W1$ z%JT%bB98)K#V8TN&atj2zm2a<`&A|WAqoTxmIJb}6$;A4s;S6F#w$L0P1myY?r@aF zlQ%YAB`p54-LZ;$nVtB=$=kPd+iEzPN}c(Gm&YA5F>KRl0z0`0NpZ3SW}HmY%#C`z z1wJ8FXGGY8P1MU#W-k?Li~{tEHKpca1lvrj>=~4>nbTTrn-_w#4*NEWj2rGN#xeIcJoexw}jePb-LN-A8V0#cgq4eHD8IF+MYkT;JUW?pj@l`LM5k zkg4BeMuz9M=Iv2^d9~{cn2vz?x-q8gX)7u>(VMgPbdbWG6~WyW!(jO021j|uup{A? zWpqSGOp8 z*LrSxuU|mETrk-2&055Z)V03U&6XUr7fJr&{H)HZQTzHHXGd0i+nrwgb8@+x|EKP_x6Qu(=M&$YEE3Q!-Fx@i=O*m~@DXP+piObrUq3Crn<)i!@T}xhKd>jxrIn8fADE8vK zBs7V&Fpi3|oS$u$_bB?JYJ;UY)OIhZ;dLkF@D{Y?GtYIi*vx>&cGcL;pqgB*U!jI^ zc+y*tK7)I2jQ@Tp`?{vR>m!Wf%2LOX7#f8AhZWwhOi(0W=moFft~(4nV+x_ucE9lbW4#_{nHd6LV`Kf7%+x{e*`_8sM(eL*rnSbq> z|D@mV-8uifU;(hQLY0JnOTX-YG|T$m(68Ko#J&|Qkmdg`>>K*OVPEP{VCJ9fi}*L} zJGb$F$-d4|_H|)*s&j@l;n%^U$N7rCJzfS z59<~D_V)UE%!5km2mIM@NEn(wB4Oq|M!c$xbIR%j$S_K1{>i@a)KKzOqQ=E1P^OJl zjU7+V6=!6i6&s*?NQI##-Ev_Z*RmuoStug8QRl?KGhIhvdR0{wWmgoy^hsxY^TtVM zw?APAp45J9;EO~a3M2_9H&&x;ZM;ob^%^z6P&4vCtSXsStwOMDvY2YGpLCg3wH}*3 zAa+QrPIF(V#8Ka761?Cs^R4&=8qSA-*wHJI8557;Rpq2UPPS(1$&{tIn4>ltlf$O{ z^CBxGV>rnX*>INTTCZEij9AC@`Ax-p5jtzgU;t`e^GlK{XX6VzXs_-?jT|wy7^K*p zmC}a!?Z~WWI_DHd$$M*l`Kp5A@hn%AHN>*R=U3~i%nr`Tqa8=(G^3hHQK9#^HRRkP zeci2>p$+QO*G3tBbVOu3WBbq_ypdBsO(X7o;dHuRM1d1w)`?M$Iv2x=Je{g>%U3FW zeQB^du_+F_autN+4%xe48{9aC`S30IvM%CQq3HDsvn3f9P1i38Y`o~g-x`9vu5)bP zYSG;}jUVeRf1p&IJT5+Rbe-9JIa{&a)nUZk&67ZszE7w&gngf;;SRaR4#VsxHigTZ zh*iy5sD!I{UN)ZWvOqRT)yf^}F&n1d8uAhYXR-QY?wm3~>l<<~nJK@%;Ei`=Tu=V`Ey8P^iA==Cr#TvK>916q^M0$V)f<|Y$9P*o%CGF*Kw zCN``Yym%Gg9`PQES`K!k3rlLl+-^$rs~av=aNic~i9qh-cs@=@~j>Ft2))O}@g zylDazaA!PN$YQrhVOEfrI7Hz{LL2t2v~7}fM0zzhF@ zyt?{BAnZI1QeQ#Zr}T#k+UW2nR!yyr6`T5`fIOyC@l3?8If`?WmUKv>a#%0CU` zckob)ogAP1asR_QZti!d#fXozP@cVVeef|{*OTudCcY@_9Iq5O8t$k+;>iRId@i(E zz~%K_9O8pRSbv&VN7a%v?{uE|v|qu(P!D;->3C)&K9+5`@tC?g-}rFH&5e!lxp>#Z z;hV=10`d&iMKdB{*JKL9%0nAM#17qihD!K*pdo7d&|R#Dq2xz|=U`c!?Jbx`WzrNJ zPFXxFMqVVG2m;D!em(u~;M;s3+y#3pjWs7pU!BjwuZ4$2h50m|cnDu}8T-7_^@FrF z0RZ>aS8V+E24WEhUjZ$m$s_vIf+9|LYPC1pCC5tSs6JiJj}ZH4A3O>XhAo$7ME8mG zl%{^ymEHn)fSk_$p-063e&be56LB5!awt;U zTkGjZJ-p|wB&->Y=968~A5xg`sWW0o5%i-4noWZ9zKlhx#t*OIw@seTzN4Dp&}DBN zS|`M#%Z_>)>3c-zQMeZslLi#|jJqlRhB(iQ9(Pr7EgY1XO;=nVL==J1% z`J}XP`klcQZTX<|INx2=wry4h84p`qBiYqPSI8W%pIeWXRj(;f9DbL3fbx8gfneF+ zaGOH%t9|DbLWdp&Hf~q>#vVq} zNg|+56o)};zKZP1xflP`b|KEoYLoDuGJkH}y(2Y;T1mUqL~*IL7U|A@NfA30Jaw6V zxK{KJ>z7Ko;x=c48V^-wq`XtsD6YU~;lapGFD}N~H+;WUyXBX#Wnzj-;L?5_=a_VP zw%5x<}z<77+%>bP_wp>@B!s7(nmp!1_jX{Q04%-hK4tz6<`Blujoy5zL z)b6YdPtJN&%SD=*IqyzI08W|t7@h9Z5S z_wm!)&j(MpN;VprM{KYcR!Pd2qtoLhHJ^Kpz2we65+qrd6HHr}12I~)a@0i*A*CRe zr?YWTK4V&6&l{EamCa#X<6&mHD5ZPYh2l`2V(Q)xc8jK-$&hcgc-8d zEUdhnU@3!>lTeE656G>dW#!>Y!WHcfRooVQ3vUcstnSr}WFbS&oj_BgOKz%i+M_r8 z6u07#JYXg%es6x&nCZGqDra~)#~7Yvx978u=+Mp)ecZ@Yb?F0b(|y(9AK~l_@qK79 z-prPKupwAmRY{lqUQp)oa*>=ts*L1w$w&)@Ys+!ff*!H8-SZ8%uchgCS!$gfmBYd= zj?2I=5<+y$c?^dbGyZyP!*5JI7Lymo#wHv;FyuB?%bD|t^|l}%v9=~ zO%*iXWd`(|)bAOHz)FRUZZ|2`Eo!Dem!@QzYrj^p#(G9P(Zs7&>68#HdT{)8?5MT8 z8jr+nzzCy&LMV~0X*YaWe9`E#b2@QLkUDtrVjk`B1V=AnO|lHCc!nyKl6ay-xk2G@ zhGyK2?5jdoO+yAQW7)e}acvYi36(uxtbQ1>`HZ2?;r!s$oUt9W;z4lhm6T0qumJ^? z)ZmBds9X{28W|rLkenRtitWo;jly0Om#ykR%wDOrgmNs>&Z*~g_ByfrE2p5D?WLYR zHo3;`1;wESH?*DmPz9Rdr!$%X1U=}O=_B5$7+Ai<57LhwcCVO<`X#|pIaxpj6jjc| z1B*w&_)O0oIeHu%RI$cu`A0k&zRCrTw-JS4X;`6Y*1$I(t#Q0and8e+IywwJTpOvA zI%A8kRGkTqi~d-{xk}YU$`2Wtrl7;f$BRG?h#xOS+(ft0a8JLU!BPA$&o%zUo@p;= zAg*i7R-{B7Ikr}d5R~lSiregz`~^gR{nauQK<+ZWmOGl4F_vcrytj#JC$H>~sp^j> zhI`Jg_pIpHLEQ;105SCNW!51U!RP=m_PBkz+>~A-__z(-l)m~+yrwgQ)TYp0VLSg2 z0s|Ix1uXjWiHvTwN?d9U<%n;Cv#VjK-N>0#yPuYkqfO9@5%RFyO`<4uQsvk+Lrl49 zLo++;$tm<+>qW_;5veGZ>ZZ#DVre3f9=t5pTI7j15uhuKb%KfVat54erN08wFbYm| z1$VoyFb}kFu2!jYMXjPcdaD=oGw5!mkI!=Q54T;4&R32V9~&Ie+PZq%lEoA}1Z zq;QtPs#{OFD*;pcF}1jJedf`_9bcw1m-fcnfB>Tzx{;U_Uz2(+_W2t+(!%AHXkyRH&zCqxxi@1o zW{9ZnWI3y+lRXl$|02|t1jLppTZ6!Eh`nk)jv$YGBmJJ$d3W7?Zcn#Az2Ck~ zJ;2@~I<0U~oi|xzts#$BD>=G?5{&RLNsv38Ltl8l;FPPvy3Lh`i#XLn855jmrlIX$ z`@t*Likz5zt;d~OtcbOcLNmWK&o-B7c#=f@=qLJ^ z>E6Y#LQTbYhWY|c?EUrfQHxmWBPIQePk9+~!08{%F&jSaKZq4wsmpSXeO6wo<%%Qp%E?D6@N}A`87mbHU5LyqVW^e zWvdfa@*rCM9SzStWz#?*4teHCZbBrzsd(y!x`Unfl&d}2dec#Ci`{%TF8PDvN2ai? z3&TMQiSBP?uODzFQ`Vbn>9FhK-3I2yuTN-`9}^3?I1la;l3aCP-A#p898e8c94a_7 z#$a#?Fv#b}Aykx7wSDQUYFsTtL=TDqc7m(mVsYs0;>p1RfDalt?}=*T#R)V+%>&A1$0ppw(OZX{3~ z2uj++80J|_^5ck>Rl8x>aQtepN>L#Cu;=#BRBGA*=e3gTqBd3;C*>MRA!%y&Fj`J7j+K>)WoVL>y%i$u%n6d0VWYgvb-XA;HwSmGuAyB6&l5EDq8{6sBl2)_c-yCx zvm%AZuZ2N^^dzt1=3nU2lG$i=CqAs6ZYbO5xX}aHI1%=)C)j?nX9y`>8w*|@;voZg3h-kvh!MoXuOy_b#z1^(A@(;A6R%8^fvgTY=(_LiNjPHUaiMSs(^neB5!$ zPo5p5BoPXNZk)t+kdGLP2z8W2lLSxI2hB;ED=i&XO%U3jHrR~F9ZiOFX-Exq?;A{Q z#Ha?(MNpq_APxcyL9snafx1bQOmOf5{&WmI6DZ9R`IE{z3OhB!*fQUg1{dtrb?P_V z)^@u=Rsuvd3BD8Z;vBt0&MKIi5wATWdPS|f*|?LitJ+r+ayrWb!roZmeGqPN(`kOG zPHI@dsWr!~nq`SkvSf0MDeh9u@Q5%&a$+*6NplK3UEOqZ6C2WkJ{;|K)<5A>Cr(`? zdAjSR&2nTF*E`ym{#DSXd*3c5(yyK^nn`Y-b$X-dtbhxp63K7m+q}T|^Ib}fEMg!- z0p`l*tFxsCr7a`#?Sp~^1nyEYw`AuR9uFX+MVh~7Il#<6q+$QNJhO_hwD4<1XqH1> z^xx%~nIWLRPA|cb|4lu)9}>+F_J4bRsST|XX9TVPX2xXsS6w$7rhgq3+A=vXIsTZ2 z{v+-2V;cJJPi_&Jihj>K{fq_l-M2qsfqXyR{0j>NivBMsf6nCQ2I!mXL6h=` zKi0JSS5*90NPa0b_zRNnMGJpI^4)a*M8eDhweVk1z8m334>5l)xcApV!GDJ3=OgxC zpnUiEPf$SLlU{!rg_#9v$G@QbiTYnw`SHZ~KSc7&$;>ZU{`{^#1Pk-`l-obCKrjA+ z@+azlL1F%m?SE#KUviYcK>00he+-o$g7(+)hyU}K{gp+22J#0X{p6p2{{a1Zr2pfQ z(x3j3);83)L1f~$u!BCjP(O)8NQ|M(67=I<#`Spf4NrUfwnVE}c+f0-Ekd*JbF z+iIIx7$P$L@wtHClBz#CMcUdz&rVn08bEESXCSYy#K6o51~PI$*@)j-U;8^(P(w%m zx-_;h7tprV2T%)evOp_w0$G_s%*-q-U{(eomwV|w~uMM#O=RM!Ye=wM%wY~u&5CA>H1OD?50AXbXu>uSL z{}^NDV266``yaskpJVLoP*?ms1_XZh%THsh%wRUClYSY4TI;tlHdbg<{ydKb2!cB2 zmoa9xU&nwTFzavlKp<%1{^xsvEI{_(@j;k>oyW?|&H-hOpVkFJfXu(G!O8*r6*CaT z3f)3~z84x6-?y-zKaUN>^4q;^EYP13{W1^24Ek-14g5XufB8HxD|A2jb?kRru>--s zt;-JP_ttgRN|>wT;d6tr1_m05HiJJL`Xc2AO0nETEkC2iGAoiJBW&0Knh*3CgXa z0sz$)Aa)@^Ar5{P5Cp=?!7m8pfIg`lg1iEJg1jI;W(Ws2;J=@u4&dfSw6WE;w*BMT RXJuh!g& + + /humgen/gsa-scr1/GATK_Data/Validation_Data/index_test.bam + + strict + /seq/references/Homo_sapiens_assembly18/v0/Homo_sapiens_assembly18.fasta + CountReads + true + true + false + count.file + 1 + +\end{lstlisting} +\pagebreak +\section{Example script} +Now upon startup Matlab should see and load all the GATK classes. You can now setup walkers, and run them using standard Matlab syntax. Below is a sample m script to generate a read count for a given input: +\lstset{caption=example.m, frame=lines, label=MthodNames,basicstyle=\tiny} +\lstset{tabsize=4} +\begin{lstlisting} +import org.broadinstitute.sting.gatk.GATKArgumentCollection; +import org.broadinstitute.sting.gatk.walkers.CountReadsWalker; +import org.broadinstitute.sting.gatk.GenomeAnalysisEngine; + +myWalker = CountReadsWalker(); +mycoll = GATKArgumentCollection.unmarshal('default.xml'); +engine = GenomeAnalysisEngine(mycoll,myWalker); +\end{lstlisting} +\end{document} \ No newline at end of file