From aff202664083e1476174225327bc0c230b953ef4 Mon Sep 17 00:00:00 2001
From: Venkat Malladi <Venkat.Malladi@utsouthwestern.edu>
Date: Tue, 17 Jan 2017 15:36:52 -0600
Subject: [PATCH] Show cutoffs for H3K4me3 and H3K27ac marks.

---
 H3K27ac_distribution.png | Bin 0 -> 27493 bytes
 H3K4me3_distribution.png | Bin 0 -> 25101 bytes
 cutoff_analysis.py       |  60 +++++++++++++++++++++++++++++++++++++++
 h3k27ac_processing.sh    |   6 ++++
 h3k4me3_processing.sh    |   4 +++
 5 files changed, 70 insertions(+)
 create mode 100644 H3K27ac_distribution.png
 create mode 100644 H3K4me3_distribution.png
 create mode 100755 cutoff_analysis.py

diff --git a/H3K27ac_distribution.png b/H3K27ac_distribution.png
new file mode 100644
index 0000000000000000000000000000000000000000..b98db091c11c8c082fc22ff5933da5be57a29e19
GIT binary patch
literal 27493
zcmeFZby$_#);9joDk6dk3KHT%q*Y33x6~q)ZbU>nq?=_CBDExxmXuaPX;2gdq)Qq^
zTDn2t8xQ-u=RM!|p7)FE`u+Ly+Iu5=J!{68bB=M3d))I0ysa!tag6R5f*=%FIVn{H
zA!0@l!b?X;;s0bCE_{L?#EySq)sMg*uOs&Z;P<2UaypI(LTQZtOYl}A!vaB=5v<fr
zb+@?RL*7mi6Y&Sjt(ceh35ls!6eUldTzPu*Def{QA(fL2$w%&-vgbp|jN}33xf0|}
z`-H07H;c=yXsw*JI*)yj&eYJNHI*hwlqB>EQ<NkiAbJ@lNzHi8Jz<2&nDpM~Tg<OL
zClfa3;^O!mM>Ypn;~lRnikIdsI_v23L@U79kWrGGGPLMlnC(wWqCedz2#DbauJ`}{
z_5Uktyb{Op^Yfc8|NPq6($e$v=&6Y4XxE0`h^Q#bU*E#}tGqp9$L*^3KG>E{(3n1W
zV7dQiJ1Z+|Hm`sHMxb`h8__AYjQs6V-jkz~GB~K8YgnmEFJxOWZ5731kX2SzCRl%r
zoV@sc12wnyJIVa2oh(aB%Odmk<OMDUXy6Ki8rjK=&dl|>UaJN?tGu$Zrm}L#WJ`Q@
zu~o0(l7L<1dX%M=)c}L<M5LC}fX&C55RrwpcgD4AU1};8tCP*b4!`7Feb7x&qprvx
zx$(O%WpeTX4AdjU)YLTd-Me?SbJ<y0YjZ{ix=qo%(M=PUj(JH*l6`%B2~$Up9+g*A
z?8-AJFUxR4n|KOMaBy&(K6j2iyL<Kq4Ff~&&d$z>T_Mdp{Y;U?3bk8)$7hxdA3l6&
z6Y}-Tm$x}NIdkM>@J+?WE`l#JP8|IFMa@??SN;6_1nVVsfAe{7Pekv`sQC7eRnE@N
z3Rm>HE)v1xqn~r4tG@U?d-klD{K&%8*w|Qt$@=<wV{>zNN1DRGvlj|OO%Dp4v(ioK
zLo$<-lWQFnqPXtJ$psjCuXigCoq{HoXwy$HFE1}=<>hHAC<HxztV`y}!O3Y_*w<C=
zYS-!vZF6vQcRwZ}?_cu~U8*69x#5_Ou1QJBbq!6;52e0fpK=SG=mOgzVOX`D$5XBO
z<jAPO3IAKI5(gEG9@7_CFI}QDIDx#vXA;TD%RkT2%szJP7{6uLWvr@?PgO@}=YW26
zczC<5#ju`jskDmtc4}QPJsY<1`}eL81`&S8Nkwd~kt+1l`0d-Yx>BK8Yfl8bny~2R
zQZ1$DV``*k_xuU}C~m!Uk%f{c2q7}`1+4g>(4y|eIn)%-&v2?68s=40i00~-s;bT%
zk&%*m`uX!`{<}3em;BJ6AhnSyZ&s}9kaP4x`J%Y#{ri{2#0)q!bJ~72yo`*9=*o=w
z`ST~wbxUe$>W&7LBVe^H+RlSg@AhPAc@>r3xEqc|UfT}b`X#q-+_-^HT^_Hus`aB_
zsA{?QzQ{e{_3KQw1h>Ch#dmWSs&;LqWn?VBK0h;LC;{!<y+qUNDpY|pPFIfS!1fne
zy!;f2NZ#!)b<8LzC@3<iKhfG{fJpXU+dG1Gi;Ubv(DO1XkZ$`I8Ir`Os;a84uAcho
z2^qb?G2|5uxHVki=i1uHsHg`{=a9RWik}D%pC>+ve%<iaKfeCo|2oR$U;8FaA-ZwQ
z{_KA?`L}&#5CWX7)c>&wbbSBP!6;X3b9z^odY`pmiUI-BSwLq1L#bjB5YYJg^{Y%Y
z0g_(#_w4d7;LuSn(Gea__W$Bn`=}vXEnVFha2)Mv2+o2^5<#wK=jAE8lp8oZJIh29
zBI(lq&jZf+@<nQ5a+2ff)%Iysq^5|AC%VbhiOO8Sc~%|Uk**YbJP-90(~K{0@Z<k>
zyjQ&u3v=_fpI?I_<KuIWp1O3Imv`}TZ`FE_UW~YR*|lreoO%shQCHKv-JT-L$;v9-
z(9nRoGB`7F%4V*<#E2wY<wkdG*}>5E&PG?ba%Ibm6CA`>2Ll92^g6*jN#(h^#quzB
ztDc-Y78aR|zJDIFva(vh#<Iz4YHjVsvnmwb`+5xg+H9njVLK6}Xi}*7jvP1+4qo2g
z*L-Gc`@4(aGTaSm85tMmLy3`+fCDc56ta;6Hf+am*4~cWhHn#{le>{CwA<L&*nu8U
z%dWQzT&EDis7zVOIOU9<e4~6;g{aFnZs@>yh!Jtk>&Y`1Sm<wTly8m|QW@DD-t4OO
z^+o?KARy5H(Jm1)Q}c-Ed2q0XswzG?B}EL5#>ld2J7>bu1*{&!wr=sB5UD0Hvc;oD
z7AEGo@wT*7=k@E?2_`x^I`iq8*pA$f!&&}}9&&hfBcuG0gZ<6tfr0I9Z*H+-cQ@uc
zlWsn#P!&RwR#%?yzb&OdxPy;q<5N!gnHP=W-uuHoe|g08z6F9?eI<;8!*u!DHC4^I
z(l?va)8^g1z2Nk-YJ0u*x2NcCIHdaQ4vbhrj7wnG+}&J^?wskW+#Jq^13&kPjHwHb
zVde5w1c_eq|0yxYM$Z(f!1Hx8pqotMPj=7He&X!8j$2CHN@cFKDihYJ`9z4<qfDZJ
z&KY%cbMsh}g((UOa2gM`mc~5x)_d2#T@0OGu4meO5q18$y(F7TLRW?R>e`@VYr})W
z8CWJQ8ylNVpCd?(bng4=ET-zs0DTAcJlP!at#+eGYB9xIDWr2oY~0b9$ha-}APMh|
zQF#Ic-w@kY`}uP}3B~3rxj^XV{D7AC?uLegLjkztwG4?pY_E}TrJ=_RtALOYCl{9@
zer>H&#i}C}y93VFd3XM!k&9A&2!oomw6y;e9X<V-)Awf;k01uKcbkp(YU`dkG^G7l
zI_teu+Z#zIHt4wjh-z~-cW-ev-fmy(pla3K+ir<psq#VcqccPZM9H7y;~EAAxlvIx
zt-KqNk&y*`*0V$9uCZ3Gi^I92qj%HN(>+2c&Rvs2&5hB??GaQf{p;4PZrDGO05G!B
zo#%Qi)@QqEMGuC`5QKJvY*$#In7@G~Xn)f`pPl<sX|sFK*W`nk4LkXRS}BHnkvF>|
z21!GkHg{5ULvNyk*qhpjh&Zj2t9j=1X}6GxsVPn|vC5(1x!(NVeBZq}!TRSxL1IeN
zu$4W=gM~z=07p1aw_TuM6fIE`XigB*%dg(|@RxuXI@?!hma#yFJZhY9m@(bW#;0cX
zDyqv^?;Z@WvcEB_NN9D}#jAVly!BabaLDgIf|$2<U1S!FRpNW_{kE-bo@S2ro0S!Z
ziHQk|iN;90uIw!q7M9aAH0jVToM;{0!Uw^Zl;hqCjd-b*yDs*qsfe*+V?^Bxd{$#e
z-WpZ=Jb(VYgRfOI1^ky%41cDY061879-h2f#|sxOXl*_sV8(el+|NGIbv^<!l(^wa
zTe(42Ta~X<`(U~|f4Hlx($0K6B9dW&0*}W7Tu4?vcz6`C-&PbOWWff}aI1QFludqr
zoifoBZ8i4!Sz|*(I3}1$!jN9jIt&9VHQQ&F*lm?xnSPG%!Sr~jL{AWJ)!QFGw6Pr~
zp<?MXof&ge5(lLZnqyd_o6=|v-EPtZeSOzTiah$g`PDZHBmH!vBmUOGB))g#<;pFr
zs<U_ShiLD<F#&K-@(KzjFLjV<85)3I*UF9*#RxlRVJte+l}NP+k*eGe?=_FMX3VJ9
zTI}8GAM$v(!x!z<oaxs&yE`z~)zy(YNYB{WL4OcR%Rz#8MbF4eBFN-wulGTU8bt=H
zfyYjyo|HN7P?y~c%uMAiCKo*$1Q&D>y!)RJ&z<g3@hAM^Gk#4zLu)r*vSN)Q=G8R}
zoIiGvA|=-$RMq5rdCmq1G7Pg|%XpnNGdz~lDXb}&x$|gR#Sml|jXztF#jN=aX+0ZX
zcU8@tXNxb6%Ul)|mBFtK-%K|qz)>AxMv&?;htf6wpU(JH?4Dg$u{Wv9?E1|p=Amov
zoVV9DPk3MMF(yD34+m8zSbt(Zx1S9g^HO;C`NG%&(e_3c)%=ZR=U7L;T?)}Wk+owb
zWq%PMv}m~A3YPexCoAEh?K$mVyHJwssYvfuIJNBQFK+JGGG=eyS6~Y9dRhklgEM+-
zmDQQJxvzm=doov*5zuC+<&^j|VYx9t>7^9Y5#;P609Z@p2dmAyS$X;4m}YEw!Hkw|
zCz*0M#)fB4-&x67%YYskT|q;!&)Ab%CIx)5jkkb%%~p0<Z_!O51$>;6vbV$if$k_$
zWSSVYAFr3k8acK@32SoterZ{L{_<X3mVdE71~tt`v>ow9SlBa~{j6A_b$Z9{P4Ys&
znMnd6c}s&3T4YoV9e3-%>^H5TY-!bA;*;1I9~Z5)p-2}7F@qk~O=3i!4uH|<s<>E(
z@E!;Mb#<NPa|RBrGf&dp80|zwj@>Kli-?Kgh;8m+T_ixvL(v(!FYJf5FP|&WDZh07
zco1Fkqe?P4;a)O)Oib3#pZ8$1yEJ$rNI$w9Ub<=uFFy#c9FA_lbhB8SE5Snhz|*|Y
z>^Uux<as!Q+_wc%8T74kvt|2~LYsNh^>fN^-zAqQEK(rTXVKAh4tC~aX6k(U3_KRG
zW?G-|sg%`wMgRIa#dK}5-{lXoiqXR&mqc|h-qYO`6X#&XAvhs4j}>ngBkbr-VX}Qm
zHi*IY>+igd(TgdcfGs$6_;L1H*Bu6N4UWvYU5UX<lxKOdS(y@<kDn*ox>F)GHx7-s
zvo<Z*I-aLbFz0pZWgU1?%*@3dNhE3Da1KpF2kpHyX=C30E=`_m6rrjz)-jlw+}7bc
zhY-y3V+AdOLyc!VwdqT&Q?|(2Bu2q`@GZJvU_L85d+x&~Ax=^j-H@rB6dA>HzlYBc
zbsvnmtoF&Q^~;!jtfKt;>vskkU9$SR%Q!S#$-)A$<7LiINf7%J2LptN7ry-9BIdGG
z$~a&7dL@gEH~%$TsgG5Fh0mZyd^yN@*YvB|Wl>DUz{uj{zK}bkkbsc$?k<635(~Oi
z%8m+!;;z3GX7^#!E>Ik0&HVL6%g|i~IU;pv0Cr)+WYi_N^Wqfc=-7hNoF&5MCyYEM
z#0Z#A9vWe)Fl9h2)b4lmeB|bD2AVf>`3x77$hGd&8s=B9A;c``hMp=R`quW%>dXjV
zi`AmY-cDx6T0*~JxgJF`g`uYras+eOSKs813)r#?4~w)$dWXnV2h_9YO%-biaSM{D
zi6Tz5aCjSdOXB$jjU%*1cGFYeANxS8fEVZ3>&k{*^XM^ZY4#~0T<Ed`5_n`j?f27g
zbko3wkaT56C)-GCmLdL()Su&`VDsGKVuIiLOgOLli2CCCCK*zAxOE$8O+{F-uiwuc
zEA(^U+Z&>%?h)v<(Ca&qUoC;WVnA)m2Q-Qq@%`_8nYV{s)H=@UJjubQhK5=Tq>>_6
z_|WN^ofqED4e8dDTk$#Y$J~flX81FkBqNiqfd3HTv#?Hw)ci&rT;f=Q59_&8lJfHM
z5eW%-zv@2WbvG|sSX!)jFd%my9NNTQrTY2txs<tNQhLq$w<T2g)Xrctb^Dxe?}U+>
z7k{_HNYJ2sHj_E%q0yb!4=>&im&O-O9O-<AByk=t(3{LOtVkm{>$zUsLXQa7Uo%(P
zVY~_!Qe?Qi`#VTSTU*lCAX-il5xrBF6mJ_E4wK}2&*gsP_ug@W)NeYkd_mRN*f@6e
z->?-6D+J7UPcx0A^r!1$>(+d|dcUO7V~X>hcK@`zT0)0hq3k3J)5z9ru>EM&93y}x
zpR?WB#Sh2pTzUampwxu!yStBu-KQ8(2>ytSF|2bY0&|9j#-pWrXEr-4GZS*kpo+bR
zkabkq$5j)Ry}!te=YM~29ow;7c8cT1e&nZa1m_M=^?JN=u2Hoin$QD!b8v9rcU{oW
zH7NfuR17#OM#wIqYO``N)AwK>{X5E**o}x|J2>+zk_T+dm0SlO`ukDPr)MO1$j1uV
z^@LxPt=M~y)I<#KG4)EVY5mpe_9Q#}QXB(8_zy~AqzrQRZwX?aU4w%;-@hxiz6Z*}
zMJ*dM)1znCeTz~IgU!m%?_Kmg@B}`?e7fxoKK1cqLQ^v{&pZ`5x!d>dCA-Y`cQ+>z
zB7{%c=}-D--ODGY2(Hki-hG5fYC%dU>6-QS?d|*bQ*(M8nfQc-b%DnK=#8eg<%jTl
zr$NPs1?#ahlW_)8lTs~xef?M~lY95nH8p#J6M>(so$>4B5H@6K_^}@ffbPthGZssu
zJ*+skOTnF0vvx%j-=-W&^$WDQ{JKez2gb8hC}Gsr)+Sg#*h)`FC+yX{G*-L5GnYRX
zF3Z>s@MYl~cP4nPB}3XWN(gB4qMFkdYmBkx`p8jAN;K=73RO9^^HI50*rET;6E;eE
zO<LqtAml*6E%cKRZ#r&9H@U9=TZFwwvIq=|utw*5lX~?#ciP8F>*{1kNL~w{B9!!p
zK=7FBCfym=%1_;Kd{-^vrk0~RTAghYACYPjA(MA!5o#u;{MlK{;-VrqpZEDWIT6hd
z)E0*;2lUH5H*HDC$T$TAx*e^*g<a^c_N~SM0#t@AVPmMKqtj=r+mn`~_fa`pE5B-V
zd$ETEc~momOQ4C*F?>^bJa012YWKK$b;40$k90*QroOBj7QSXqSY&h)z&fxE_ltEk
zZr@J2lcMXtwz|3+nUFBN&=0Izc2SW!Z`D@vl=j1i@2C1<i_HIce~ZtoB?5Ef{0Zot
z{o=)xzu~f4y5h{zU@fJ@_ajGexT5BEI!UiIf}!Sn@6GQWJlt6wuI>mpff)Q>ayXtF
z)^O4k03Nr02|Tqv8Ea>`xiGZTo$tenEt|_LcN+2DguZXh%*^P%zxV8tN<vO*YN~Sg
zX<Ay>Pj_C1g@tAB)*z$)jV%##OPQ+Js`mtAuBAlZdb`QgUD5~HqK=Rvz7KXJKJw}o
z$w4DAd}bk0yheG?$Z0v)*ktkiK7ZV7e!RPh2LzF2mgudSU%3(U;QOmLp%VMrs72Vg
z?!!?JWJ_o(Y!%NgU#!v#4-ZF&Vw^~he81<Z^JwaQIK?#WXCG#yt@z(k4HZ*H;(NY3
z6CZpvN(A)TrkKJ6zI`N}($nXD?#-jaHk-vSvC6g+%FY?ga}@a#bBJ!MYzF#B*K)Sz
zC@7i!2?>|ZlNS)X%iKImK!D3e1M%6q2?vXEGcu>~3?o=4ZoLmVN=SH%QY}Z$yX>;O
zx_}O(h9C91wZF5)lV*2%jHbP9NOD${evgzf1t{=%dpO0S4#^h{CdVCph&)^>7D>ax
zA71R@OZCTC@HH(H)@Y!#&hf&^3g;b-qTGCy#cQ?{1+#TNCYW|ILPRSF;@FiJ<p(BC
zUVdCOL7%$y2fq}0GjC6j2zR>Gjani#zNlaKGq<$pXfY{g#rBC09*j*c5JBYV(5yf4
zB%1&^KtcA|<gG2n8?kZtKNoQ?uxHcSMsH&~XaK&9K0j;#VhcX?_m(0N=Z*RYh3vpT
z3C8B`-N;rlY2`Vi2;3Iup-Dbh5^5e8?nV6I{4{E*3rG(b7VXT;OojwkmD>K<A8&Tj
zJ?sFxrlVZkS$H~1Lb2K*_Ei_1<ZR4bgpr|rAdH9sopeP@-*nePBQ;jl<*}V&b=mpx
z9aGL9)LelRZUm$9=%5U>Zr@2wkGh1ds)T)qW<YOC=2#Ey@*IutWvldUvb%CbI};9@
zkGve5Vgb`p<dGtODACaDJkYwKf}&G7J*FqwGv{5A8{PK0<s8qC|K}qMR4097d?o~|
zMzmYUxn*DZks$G?wcIq6_ONIY(NIm-Hnig>YayjLcdb%S9u1T&C|Ihnl-=fOH8IkA
z+-v)hlbc(~rF<{BS4pw#rmEG&vZjd@84Kht9U8W^N;8Z|^d3i^i?y{Eu<7pZh5!gG
z+(L41>8@l*Z~X0s!2F6kbjT}GbkrJ<Kt*7J1F&~O!x%}m98C|Sj5a`dVRSDnacmyb
zt`gkK;-mItH8=lNjE4=_ri1p=8<BV;-A;ue{8u99-FUlus&o7~+_KX|1h^Iy@Kn@^
z?OfQNI;K74xcTV(jrleq0sliwkfGfQERZ_sy;HJiZ-73X2MMI^NkGzb90f2|g9P{Q
z4<30&?0;6jx%*t#F5UWKw+gj4^n31bekJ-mt*3LYiudCZUcGt?^$DP~4Oeu7zWzK%
zvi)|gaDj&cQADZdrRTTUsx{J9OE4n_1e{l|-ho`KG6!M79)j|n^vSRjJ?9r%{^MBX
z|1i+`=`MQB_UWTYr+9XvOKyz5OY`-lnuN&y-_0eD=MrE5@b*fd5G;wwyL#OrjC8#U
zGRIG2Du;6;Lg#iyI?@w^txfUU@g~~YPCpLdzJG0=W&<qKgrZ3|>E-1f1kX4=(zNK)
z+hui09^cq<zBA#J5xqVIvLy~iNE!g%B=gsA-YJXrT8g}Mt5W)QR*TLt?!&>LS>LqF
zZZ>7sn^N#~j4$um3m9;|*LtkFbK)h&dx;xtWzP<_lKjXiY)8Un#7OW^IG%%tC#rMi
zch_9HVP_d#kL-v~)@0P=!sUM$Yue>&{iZzL`ja5V0^heE_jr<1Q>DPia<H?Pu6><)
z$Xl<cVU!CQ%D;SLfz?edIbY5G_PbBrPL;8c>I88C>OD~ixvBnfutnSa;6>({_Q}};
z`=A18OZ7kNcyrj%3g~$()41jMq$;bQx8dWiL|&99T~EF22Y0#!0%RUNj6;CH_AD(e
z@b|04BwMCzrSfgEx*^p@J<JGtW}G-TPfX(@Xl#a&YvL9;e?8W9bu+tRh?HDG4X9V@
zD!F*3L#=mh!OmV_D8o!nWhm?wYoLr1QgZh12V^BxQ~Wu<-CM|KaAo>772*|dA!=eT
zfQ-g|-&2AzR-EV2VJC@KK9Ps~3rkw6exId*THZNctfjR@SB7lWbR)WT6qzLH=;(^v
zmhS<=vikh=Xk$}T`Dihd_geUambgs7a+Xm0*cp97*&oLw=1~NJJ0dPF=hv_MX{o7B
z9v}4~1!<AV*wA8%z@%;8B=T@`E9vOyIL*lDBb{Fj8Oi6|-29<0{9MR%Sobtx*+7;?
z_Wh3o#Sfb4pt3LvtfUq&wLpaTOqKy3YyRmm$<9|^pErT0*esi(xIKP_%dSON|B<fs
zr_^$HFKvmtk+M1BYa|%kyjsS&hleciw!=@((qH&rOc$B8I3XiHg|D+5e-pa)=_Bl+
z>45(|k?cC`E<0Wq+$pv>m_k5AI-TIT(Ej7SNwIM)DYoMg5!q>a`s~zHjI^w*rw?#~
zlf$^%_zsirrANnHQR)$60e%grfd{DJJ^y{RZ>ShL$?(}IOe_}h*;%u|c2s%q@|(9M
z*=C<XzJD9SF`ny@4HO%v$g2}sIKLx=$3*LUyje0m+$tusi0Yj?Z@zt#)5)EB5Ch}L
zuBtNB(9p2SV-R(F=SL~=@zx^nthu>*A1wz<9^zB?b~l@5oMwkT7rVEX$Jb%uW=B6g
zG2A|njOGt2ETz|Dn(n<ns~_yDPb(k=TuT&2ta_L7_wy6Ti194)(@FRY9bNzAt6le7
zwa=<gD#u-aLqbk-nUhmNEkil7OPFL1Dp_b60G7G|b(kUFgT1WW+<8d25udziYvk%G
z=~IT+5)OHo%-`Ry?m3T$kZKVjRbRW5rXd%af9$ur{Vb-Qlaq6OoasQDCP=zfboIs5
zmI7~ODwIi~Mu^2H`0Q8!>4;_W*?MDJy|>YnHv?;u^_ClX+>0Q8YXRO)FX@|@q&$AC
z7=M(2I1CbP=jmS`@q{S^dtM(YVzK~xc7e|>DjSt<X7tHuXC=V3UU5bOLqoNkor`sK
zbw#UxXlZITJ}Artet15Q7m60%n?o5u6zA#`%G|kgC&5uOySQ0BudVI&`r?QYP^s}*
z7DHvux|JU5B(1LMw8+-QU3ZG<{cD6j?XR!>S#6zr<Ikvi@7~WdBW%|%Yf~EeR`0jg
zH?1{<+i+a}pa%(N**jHa$tZo_(et0ROVRq3e_&r-$RSD{r#LVeDXZQ*4Mj!8*wz1F
zid(Jzw@fjc1S_(2$@Rd$tDop6<=+*yOa8SFyEk|I*UC-FZ4Ex`^*xTx+wno}(xfvI
zyzjG)xoYC);OJNkB^+SCIeB?iIJNR3x`gSdsi@w(dv|4T!L?cq>OO`(+oq5lye%l`
za11r9^srd@)v(iS=ab^Ix5LiCA&;&e^{G<=79D>h;f#`!k`t}Zo<7wdA|pkZei{;R
z0z3I$>>kg<;%zDa&F-NF-|Xn;2XD1D{zzBHyvn|;m)dShfE)WiFvG)zt@TjK)zH$)
z&?$Ui0gfsmkuA|n6~FewB-CPOZDwb?)i)15Y5>Cg8Li=wR~iOLQele`0q5=Jraxp@
zd|^Xo#ZQTIXpkZeYMO93XF#l5sx|nJN@-2Ef``x3kOcP5Ki?ClCDL>F%i@{}9H*3_
zR%`)P3OWHx_2|LXI*GktaaZeJ;BnX1C9KuLzOMQ165)b=`rbq&qqDQJd)oPO%W-6l
zC#Gqy0hGq45($v4F9`j4bQA&Xf#Y!nV5_^f6r^wABXNGi8)`a{7_su82TlkIJIdn9
z2@F;eS>4JN%<bPYo_&~)v8WY_Rd@Kq`Gk`YF+oF%K@0GQkXHE=ue^k0&C)W<T2L=~
zc=OHN)Mbe7Jf%(tp5$FW27jahHRJmi8f!}Nr;HTvs(aQu8t^z2AtGa9`qMV&nc7uU
z2DeHuGeHIMXihNS0dtW!3~zCs#25jOSDzxcLh5B&d!GV1fvq0g$lv)h?Q%F2ln4Yl
zMP6MYp{BeMVSpaQQNyZro%IvJw}CwMtsct|?zZorE<kGL4ri|6@KgV#z(r;oido0t
zg!lq^)wzO-P@k<<QY~U+7{$v@?3-(x4F?{Q6=MkzOqLr<ukoqqR_U$w85n(z)RZ5#
zV=5s+NKMa}@@mQaN%pcB>Oqve<@Y?-Ut+|Tst;0AGMbL(9YZSoQJkU5pJb#P!L90B
z;#|yCFfC;ex!|_CI73Eo^2$FND{p7)%;&Lf6gXAn_B#ci$}myB_R(%AAE5b@Y+%q1
z2lg8EFI(!F81`;Tts3AUp@{Hu$ab@Fm%I95_+)BJ6sFSGt#GnJ$OyRl{0brD>@&3c
zYo0B(l9#uqj>s7p#6x8n6c-gMJ}~41_onKZ4J*5FOzyz0AB|6mHWY?A7by{$!)>5A
zaOT|~JK`@#E_}M+vF>@0M!##SfY+q0c+1wu&VUO%+y`|-2UKThGPf8Xl4^}hcsyNe
z^KEsL22hGm_1(5k)AU(>aR~8`BV~*znlb&&){=`eeXo5yudW(6*N$@Ea%CbGJO&#i
zzMWRLZ&BtFv3RZI8ZxbmqOWLUlAP>LiR47Vy(M0hLP+iDeqZ<JVGQHmQld^eEzaU_
z3nW9mWTfjnUh1H}yM)Iw;KR^f!v5Z*QI?jRN^yzq+V^PQ%M+VIvm>ti%8N|>IkjvG
zI^5*Ak;B@3&G`+@`Iid0*+OX>ro_K2wV&kZUZ1_NH(Ith7fwk#m%lsQ9^QFkZ(R1^
zXK5Xcbth1vkyNqio@&0SaaNbuuV4OAlSMGPa)l6#O}#Lqg+EaK1Ogg$D1m`j3g(ek
zor}onUEHc1j(M*Mp!UNJtzP|%<Z)IH(^ik6ogHfmpNB<J-un?J61fE4>LMkgG_*m%
zyDu`;GUV}sHvRYQS&=sPgFg)~b(g<ZbFvggTTU`09j&B~%aKU)?j~E)=q5m^RZo@Q
zo#kW_ZNRK0SaB7Y_t(U(+X)@0vMW9#8h!HWF$oMWfpo2^_9)X(s;l->nXhK%>Rx_R
z+TLoF!-Lf~d1)$>55$UecuLoX-%vVtCffCixbEK>Z&~1ln$%-<ZEzhA_o*wE60TUi
zdO(qwx#_bMbK(!iD5#E0J&aXOW4$}3Y@e)usd-l>`Uql!E{#@hRi+`}Z26UueM-Kw
zT~l<qXRD?1eGfJZNCtjcJIsa%Eq3efTFFF{B8;%&H5$8RJGdBw+;W!;SdA~{e15~}
z!4hs7-k&|PfeGEL+Ev>d2pKtNK!eDkSzQ@!Z_BXNgf7v$_UZn3k!g>f{Pa0C>`q4W
zLGqb_{pQgu1m*Cj2^UD$?5#>jHyH2Q3)q*73o0e<4ILcjW2pzbeH>72D*mKFaK!0w
z5RO7g8EciwH)#&`EdL4o`O6}B9D%m#-+Iz}n&y1^1Fb^Ffr74SJVJc^U@r<}G{;-1
zD}07=hN`N`YimxRwK;1bA6@pLRZ>iN_9nJtfJ|_9?^Y~D-VsS=lwrt|H*ris@p~Yv
zn|I&Os|RnP){O0F6`aZ{P3r@1L(Bm<do=5J{*AYFQ8L)v*pr`y)=vbrt;G6Vz;fhp
zX&Npwp{*Sj1M(WTIqzO2{5jx*_X+JthRU7-&@8C`Lt8?L{L!ouA~|jQS9$Y%ZKIa<
zS9_SV_nkXV;Zs-okLCXPc@BK)$pF+PWvhA#;*N*Xaae_Q!-MRdq)I(DPtxV_K9%D)
z7=RuBbcWyv|DmP&H12R=<wLQ8CDi(P8Eo8Y?^w0ho7Z(=x5`8lG1nY!f*d%Xp4l5%
zyi8s2l@dDrysqgnkLjIV4AysWmDXF$@f$isYIKOt=kvJkfbYiNAVkGg@FxO~zvsk0
zv_NY@<H1;6>igtTB=UFXXLeV9*tU*;|KUW$IYak%Syc_Ne|wr!IRO>6@p)?)rt&@1
zziK?uS#eIDU<Oxnnqjzi)Ur*7&wua|c4zCGCWC7N4C?O+M9IAc<euM*)+6(}KBVKP
zWjm}HojwIl5@kp9!-3j>@j;pR1uJeVW$o2z$n@m!Ym@t~!Q7;hHK^ar`!;-PSG`xC
zEfAOraCH18-*WM(u9QAW{dKyig^&Edq&vPKQE+bJ{&<;>Pi=LoRSQtJYcwSyo~=aW
z5qBVzM3(jXmR^Q2&_o^wOEdu251ZG_)nzRo$`i2axhy241>9&)UGV%==|W$jEIAEV
zTd2?4#Y$ACBVZN9=>CgperTYq*3Q`YO%#_-G9c~gZY?7XC8b~*Zrxw=p+tyyN>Z{M
z{xtY`E0!I|u=Bct_v$Y9_mtskVsrdRYRDn!_}>U+{@eWK)|_u$e#Odj+wytAp@;Ph
zbb{8e5<F(vZH6m!WJAs`Z`2baiS|B3Wnb@I1+rkw{(6zqj4BX&DGJMUR$&?=OxKj~
zYa?WW|2QtM;{OY)&jrPZ59J#zB8y$Ekp8&Tt0anZV>vEgzU|`jVan~VzaE2&22o|6
z?c4K+S1c6cCyr?Y!3b#|`CEh~UXLRCvrc9!!vL^UxW3o_DXGcscI5h-)}P-IBBf-^
z8LRd8Oz;JzOMjVjHlW0*zd8NhFMkUn=_Q~DmR_47uc2YY-Ef}%Bx#uik)38qvGc05
zAS-d((BC!$G04DqM+rM~^ODiI|0VIk02AfxmLD`khJ4Q&)efI}|5`?pPy<+}GSg8E
zwu63mUlt-b9Hm#U*E^!O@<Ds$Lhxx%R8+L0wigNuv-J}D?}3v{aG!jwF#NwrDT9N9
zj~%<_e2Mrj8pA%^jrH4Pu(TJLJ_hF`{o_V;_MQNlsYI2Q#nlFxXyTew81H|hb}650
zYHD&uiR$auwNdGox_ZPE)HmWlh&iBd^j}!}kXb^+Xw_8O3OPshlEo@4vs#|8Zj)3+
z4Ppu)nacMEzHofYf-|cH{*|(q!Qv0omSLZ&0#GgumO7rHrcRDF^3eeakd31_oEQ<8
z-)-Mt^RX5Sazy=ioIO%Kn(+AR8s}=jJXK2UPH;PLD6`Gpi?fOK`!6A7yTT4yT#)>M
z2n0P$prgDE9Kg!TN|D2aJT~{@St=?j(DHUJ`fQ4zY{Epe5vn{jE2@RGPAlK=|ADn9
zsgGxKLJ=NeK;QSX4zc3m1sVfQ)m`CgCs?I1bOzMOscU<Q1GL<_sSsY8U-QX=cDMKQ
zGxErnFE3*0Bz!8&5<`ROg+TnOT=<~517yph5~%((wTCL~--tH-ym^o-&pnC$YVbk@
zra&PVckX#VxC$OCp`QLkmEB!d<ooxZ_cF;7#|$4@{X?BXyWzB7=waQYN}}`F$#P$4
zma3z68sGMebbm??3i4UIo%_n}od>nk4#=ypL-O7LD(zQAEiaGNKDC_{0<y2~jv$7(
zO$>F)$`A8#rXOas%$-fDt+5?i8oL^I^xT1RPFh7jZqgjk&K6^<ijqREX2;^|+!*<e
z5hjVCQ6<{O)ih97pBtb0{l|tvOzi<sup0Sc2`Il>hEk;VecvkxX>_5$OBef!6^n`i
zz(L+cbs;sTyDQpO2lcoBjI{q|<(KXytp%C2KPv`Vc3DCY3m%}#Qhh;J#CWZCa4>ip
zB3i$D4g^i)6{@wj*P)B0*U&bOzCXNA368d?h`(>%B}L(C7e!!-X}B~>wI1SKdmR_8
zehml)Sh|;2Iz}zaWdDi5U+{dr4WNco?2<HoZKj<U*iT&F;k-=s`jlJ8uO3Y~o<lrr
zZ+q)A?=Oc^^yAluA5R1o68!2HCpt555y}Z|LxbMT6UAqQZ#%aAxQDztb~wBbH*|@<
zY0F3uKK3M*N~R>1>_xE>PiFdlsu|hZ`#|C_Nqayt&ljmU1I14w1}ITok~(i{)R@*I
zj(dMp2lx`ZciLQ2>yH4|Ce0moT$F3=pW0(k&_{%^A9NGMoz{qY3%M6+6cV&oENDz#
z{NV4)0W;5t2L=X0!B_)IrR`~3#G6N;yrOc11ec8xkNOcnx98wYH%Qk&+r@2A7J)eQ
zd@F%-gtD(5PthTvd~wNJ)L|OxUBJ@C${oluA1(A_7rCH@Uwi%JKOeCPaa6N2m!uRG
zy7WOADjMIkK)(W2oGU~qo=6@PNs{$jZ+E`)1Y9^&{8g=x8=czOAg4h$11Y*0q)TNb
z>o70W!=B{_OY-qNESkECE+F9u2EmbATHAIV)@e)>hjAexUpE5U-@UznPi<w`1|!uR
zM|UYXTBtR9sYXpP{pdS`_-Ye5NNRECP^;1Y8a}vn?SQ%h|I0~OLtmdOqtn8PU7S^<
z%9%AMKWn4#!UbG5D3KBUFwfOev|4E$dKT>G1A<ggV{2itfndLgcC+Z8k*V?Z$0&0C
z*)Su)(WmW{`q~ey%N}!o!lTOvPzCWl%og9mQg6g<>4M%O9HIUOnb~b!rveO=kMzaB
zE=jkbAmz}&8hGX%b8raYq={%qyQq=<5ELwjc^*(ge7$HQdsd&~MO6G1RlL*EFJp}^
z&BG4R;W~Y*B=P{k%;#I5!LBm&48D0(-ef7vNFE{Vg|-k_;1{?y+zBil)yw3tq?6dS
zM@We<SX)xWUnMfD;=^Cm^e2-wUI{3kdv5__hflOnR^GfGiXJ4IL>~VN{=^$VcvJlX
zAslT)v(9;IF?z}5&nyC5mI0r--BkXA?p1D8kS%Dnm2+EDd>xX&!v_a>A1ynX*E2}g
zT?$c!^*vT3Tul`Zn@8)p-5ErZ9waPquy5m?ZvBKz{t+)Wm@&`ZYHs=@2r%^7Q9}hA
zocWq%2jbcxNh0AE6N&=M`}$=V8L;Dx!;f=AyjsU4EhXzG9zs3M(2yC#)xgC5p+e2z
z2X#QVuG`XmEMTOnI=Ds=UytVSusxS#QDe|AAfMMboiHNCsZ+}v1-$6;R-<QC6a*;K
zlFc`<kNR*6{|bC`P}LOVfQL(Y4u();5^It%^^iFPcm1&B`2umE0->GhS5rCsN|P%}
zHu<8mIkyt}HcxHO%VK^*HSJehL#8K4EFTub0wcaA|4v83YV<NOF+VOY&ZS;9lra}H
z2&`;uoS=ecXJ>a1jvu-s2<5_|gYft;_QNz`G+TzPC<NElmDsgEYlr%)>Z^MnoBms_
zeQe!xbD=vWR6GlCieNp8879A<t324>1=VcAIKcTKT?xcZ0VGaUm#yL-p9$#9x8E)p
zZxu+MieBt+zH{q%4%)~15zriW?X1tC*J;opLLJ+2s-?TX*h-{gztwlYaG`9D2Xq`o
zBCew~1Yhgwpdf>7Z*NavLavbSE>{W<xyisWd4xPBFNa)hy;vH1kCh9FH=9KY%vu>y
z<uD&#rz02e&p-_`7#=jRB~EKq??i%lv~ywbA58iJ#~I}K6i9-wCg%t~No#Y_1c$f^
zB^OfD;99OaBbon{gYY$jCR_uQ-Jl0?>CHE+EUO)+#G=v$z|aeon}(<cl~Q~w=>((b
z4}&x!TyH6TTYlhlzE34JShcEUD4=Rp`_12#Q=9~}4H@b&3JM#??|T(*-n>!U6t*A5
zDJO_#y01>oluji8L)W>tIg;Ps-_Mgxi@Yc*enD}>Mle=YO>JOp=Z%SunXF_Be-BzM
zLOojItFohFTwFN@ZT1q2)~o-)imN#8bgD??VqadG0kK2Z{KpZWzRephzhA-4AYswo
zA3xqEF{_#Vg@|W?=aCQbWo5r2#O$P?q@)27igwd25?3S~CAK1{#lF{AB3q|G{Jyp2
z24re>TH4puAa7ICG;6_F50E=#dj1z>^sM+?Pi`luk9ggx3Y`fgwSE|Gr0F|gkNu`P
zsb2+Y31@$g7KjipS64{-ChlBJQ6L;;1Bm)R14hfC@r$5_?%V{eZ~LF^RRQ!aj;f7K
z?g<8A^wx>%toZ+Q8;6gOu;wI0HSk}-FX)Z!1(wKwH`7_}wSwYK(fC=BWTFx;xN!7;
z&^-6gZ?}gxykwgbs@@KUjhgZH?Pamm?|kSjt5obuVGT6^$8aY{Sd<p7<Yaon^Z*qJ
zKpE9qbBD>^z2GNb(Q=#Y4rIFmANe_fa01R@i{Rhh*=Z~6Kx-~gJo0NQYyEG9%4s4^
zWH?(;5Q!2zoza>;zA3NGpG-7;{69!(rFdZte33=xr6-&OC*?LY7}}iOO<1w{gBv14
z@he16I3HJEcUsY-*vao<Z?Ur7Lkp%nL5FCM)^(rs=?h#EAK)fD6tnI7AVi)VMZqlo
zbnivY-i>qSU`e(G8K*#ob^?90T1#BQEh`?4fnxlH(cOL^DWQA^6}n)auXZJyL=T0M
ziE7dynlj@ep<<50L5$L@Z_20cv-$sWIp(KO!OtI543elmlo<y#d@S)v{nzbBUnkUy
zG|s?4EWqO)xf`n>@e%`v9MtrEqn4wgP)<E8#~-<?C$A?fc?mO=#7u5<bh=`dzWp2l
zf&GoG4+J>G`Lj;iP-Jd%&#B!|2rOHCt`~P&OrC+!N-*dI-l~My0*UC^-f{3K`{qUA
z!xFwE&SBG6>gk3mgy>J*HHm;TP!(_qKl+W!uE@C&yPU^Q9z)el2y^KY2c#E(nq;DX
zY5`3p32ii?SCkN^g8hgp(GPmmo%t=%P_KvUOw^LXXdN~lSP2{g=+Oo|xm+><gTWKj
zqC9BxN!)3P3V0Ann>%?;vgVQUG!F&%lr0|*o<?}kph#f8+J`s*=suJQP<+Qhd(kKk
zBuEFg2GWj*FA$jGK`8t-pIe?9{Thl52qqn6t73FxFxZzzn5S9bIHa_OJO+?BdJkRZ
zyP&r-{$vbjw{g_+?j_Vb=TWO)S=f`PyYvFzYzf92fePP0p^y(ck!Td{5v1m2D(f+b
zVA>Db!|{0Pw2Ta>dw!W#{nIi1)R8l*boT`^it_I#B|%7^xijg3wU!LXb5)Y^{pwdQ
zXSw&0^iZcFi9A%PE+Ewr$QwapdLOuA(?kZ6K0{5Q{tz!zSxfj$){N-e^p@iK@fizy
zfh#G$w?J(A_&JLJzYrP`jUZs+kUNljfHE`^GojDd4%5+GpbCRt&uDZg5GA1~ia4M-
zEzy0SMa*WjxsHaYvr+l%jC)08oQ3wGlv)z!wm9mKPYq}@EV~fCPw5m|!uY=IQgSth
zuARUM;0!`=4;%>wt46ALTh2b<{nwf<*p2pAQHvoX3Bt)xGr{81DPyed1**^+RPaY0
zi@jQQx_Z)9idVc!wm7&L_1WmnUuMKwRAf&+`C$Cs-&f|+mzL~-DY04*4+E+acF+Jo
z^)2)^8ZmD3L@{aD1#eOdxnF1lg-y>}mqZ`PI{hsScxh#8(X;FSK@OjOb?=yA8(Ypt
zKUs=M5X|E?!2Wj~06b=8<-;jtl0zQ3mSh;ehcLGYxCPfxs&#04_Zl~opqn+{LhKBV
z&dtpdW#wg~lG3WDFzZX3!=ETI*l}vvY@ne(bu-_NtG-|?PESpE9n}FT=0F7bAnYP{
zZ`;643yx0CO-(H9By)6kUfTFH7TW#W68cU#7{{Y^MjitW|5{@wb{}3*jQ$Li9?%nw
z-W^o@UirrsBw8|G4M4$A)%uUS0bx9U!PyD{lhvex7<gEuoWQuHCh*vaV}d0t#B`e>
z(`gx@S-%xaQIvKD*5XPgu<R%R#&T7Yo7`u+@4?7NL5`zGyJ)QZlZP3)w)CKOR&rE6
zNI8o8$E846|CS423@kXIG-<YMUnz+v**BSw9X-(wSq@(NDEXm4{mL7^NmB2vd8p7?
zyzr@!T(Hm{Xgc<fWRDlfMBGsu=U;P9<LAJWoo#~iTe0vogYBr~UA}*J`3|4o(Auh{
ztgKxA@2)0?;0Z?I2uXZ?T#E{DO@fv$Lf8T!7@#eg=zDZiivh|YI$?*i@t=?lt_@nk
zb)^21hvOvqBjv97jNWVQF123(2k92v52WGJPSP^+)_{9&x<%%9P-ar={`jX*%*??@
z*O9S1k%0sC)h~rFMBHV*fKz{lAwF_cLsCl8NsHt@l{AT~QAA?)<Lm-;iU7jr1vKi~
zlBe%H|E7B~>eBOwGcVh2MSW0LOM7<Zb;Rq<(6I0$$39_H9ys_`I9zpXF)ne;d-ZxZ
zZt8}A-ny{r*N@NpuGnI(UBiQX<K($>7Q)Z#OWr>?_4tn!vXmG>Tx%hf%Ca;~2sb0N
z%beK9>qFT@MdLvKd(CP3HcZ04ZyrdMF94EM;Km`3&StwT6Y!C)_M1blojE#uAPuPl
zywe1Nkj9;XlT_^CiPmg^Dv~1zAMy_^0B6OoOh2PQ3qu-S;jTaebxEnzFj?9#(5S#o
zw1ANPZQ)<<-@pGoSh`v;U^iG24l4a;7jF6K7udp`z9OsMmh}qdc;Tndo`Hlmq=S-&
zWkn=Mq>O;=nb|`-E2p5oD=d=r?>-=y^u;d;c0zG7Qi`k7g!u6Ao40SD6c-oIGyU}G
z(~tR$rk0UcS54|nr+iLVYRbvUy*Dkk?1sCk#ZI?^0{$pVuS)!t-IZ?2O~FdND%{P@
z#k@f%GQj6TuY8du-Q5W-%5QCNpV8A?6BQErD3uSyl4D;l$IM(nKm=<>=L`pPCEf7j
z=&vch%|Wr|zQ#vrP8+^alHw3b<h1))^=EI7m;H!zD-j`~{6w+5M_L5XRB#cLr<{kE
zmy(y)Tl-|VRkW9bM18@rsdS5dfA!M>xi_WsdtX$tZN6t^O`tl?$AFbOzA{ZnO-I)V
zH>@6BKc7dr0hc4Sb2YPoOiS)PcJ$~4KEA>q1>e7aSMW16Xn+f-uVKdddAt@pRN+Cn
z;swL&C1y`zDK1g5=b318d8#aQN;y$2|0V=g!K2Vn%7Kq|Y?soP*6h5zL_uP~>XQps
zBjGOe^uR}61ah0y^z@E0=Q--i7P!S(3N$%y`NI2p#?EsmskRNC`xzBybw7TW<^0``
zs#Fe^<oiHj*hRC@a0ea@iA0vjSF*nS#k#p`GoAi$r5se^Um*{94O*2s4lS6+k9jRu
zfsozB>O{Cr9tV*e9~bxi?VXhA;YtxWXq<P_!@^F&S$f#VYh+|pQCX?*$NPBz@?&-0
z%h8%IX-umZl9xar=EEW)?!zvhaT=4?PHI|yuI|EUAc=rX>Xy3tX}FJX40jw~oW1(l
zIW3#VpzNK)L}Ru`JWSNkZH#2tXV-?5f=&vy&Pz5Go~CIQ7S@98dZz1%zKm@XuK60u
z$|tZ`?5|(HUWNR8$nqpstvR0q(c`BwNL4Mo?(5!UcQqYSLIi}@E%ww@4`ofNV!w1T
z4S$)nnfPV2Nj0ih5M*-%5Dhgwy#i3|R?~TKLqcJ<aABIXtY_|<-y@@#kc%vSZhye0
z;f=}W`?4Q*;?>1d|C&8ZLoVsne`4j>%JB>;Dm6GpQlgS|+lthku7vh_DGPp`0s4mR
zq9^sLg5zJum|axW$;1T&xZ#Rrg-RH$2+P#f8F=Cd2!5=LQ<A<62*7eMv#>BGP46=q
zKYCMi<OCd?Un`a*#EvGnPf2Qhrq2A1_AnUNIdgw~%ntXtkO}VaSu=-^v@L(W!eww&
zA$3KSfQss&V>#dEiiryc@=>7cPw2=+F0L1_x*&{KlwH=2pQWJ*yJdFbrH^>6Jss86
zn<P?FUz)i^md+sH4Lk&XKe!(3Nls2~0)i0^X4ALh(yNkHDp>wkQc{X{4qI2)JTt8F
zW_R7-twCDU7(YvOu4%4Z)|NWKU4`eDl7rM`&}rCmZU6FdIVhsVNesC~!&6$5;}f^a
zh1+rHBbtYg{O;U)W+uU<J?&b{1O_Y>mfj@li;0n|x;O$OW4AkeD93!OyLnxNitg;a
z>%rf;*~_RkU54E^{zO<`4-SYBI&+eWO2~P%>%s@K-=zG1<Odb%?+IbLCfuYs=np?z
zpV^<ttKS@S%lL$0VLe56_OEd|3W6UwA}|c$!yYVWo^3P5&wa?S?H2aybR&&zt^;TL
zuFpC>#jS$JSp->XURHt|AE;jzMHjH-%|ZN)zr6Zgko_*j+I;z@;Y!cI=g;rKEdcZ@
z3%#zSquM=3PsHD_7JLy9!Ff0ir{$&W{g{Ibw$nYN1_q(MOatau0{>$I>Q3v~A;vbP
zRk~u(i+rB(`kcY4WDnE8eRPs+hYdKc{eD=b!d>wwHqVE3A;)v?v$3`j)2{NWq!b6$
z;S+_wYrfXzJFRm*G#biZB$Tvt5k3%eu0d_ddpOf{^DT{mHgmde^SfZnVona`BzI<(
z$iQO(0Ri+k!XrPYGkqkagVys6i?76VTOWIi3sApw_|;#X<eU2lMKyKnZ6}&vu_mz(
zuQR~*mUF+2?um;`3wy-4S7Dliw+jQu_!HPzSRTVo;!{*q#=`Pxcbp``j_Y0Wy!MhJ
zgDQL_iI9cm=Ow()PX65FlP4D!1>)-XNQjBslBLNB2?-TtGx=>x?OD0nu5g7n{)O&N
zC9hXS#0V2F(*k96y**uc2W*DFf=2BnH;Je+xem{S<f9ZEiDxvSy=)5OQ}s-n>)~=U
z>~wzX^UKYruOEZW6<&~zc1sF2{kH5VfzC2vs%Ng8(7w?phT&W}I<kK}GHG<y=gcpo
zCnB+fa2DidW;G8kacLtIwbYr%(atzZ980IRoP5WITFC?#gSzW2(ca2<k}l<1eMC^d
za}700=4ced*7McvDT1@qH`-v>1*jb;*gk|UenUfIBOc5yN2105XVx#N&gr<7)%9RA
z1uWJW4j$Z+com{LGuWSG&nTz4aD%$jg${zADHhC;N+=TL{5U_<hpcAw+apL9Y6C^}
z<mHc@VpCQpxSC{D8RDjMOd`t5`?m<XW4X{BtA3Y!!{_lSHiS~<PXjs>wc2Xo_ZM;r
znu(~EozMk76VW;2L{VK<YRE>{Yi<rR>MJY$^vs>P{@@l20BtD;c0+*r?PZGrtavFC
z?gBUew!Wy*NV+e=oWblr;jf9tWq4N*4Ds>K-BuPPXIOLJX?zH_Nkrud<M;e;HJNa%
z2YZ9)hJ?Kbo!*9~Cl|TT`*bh*N!%W)jv13;F>Sr~{z{ULai5{P8-v-rj^C+3Ih{j;
z;wIU+U4kF^itg~DZ6-cy-X*-t`)+mOjRv|E|9<FhoH85VhgOm#B{>D#tbt+s<i=Xd
zV(N3}0=u)yn`1i!`D0|d>5+H$6a9<>-B{q@gdck_-)?8{vNiDqW0t3HL<(p};T<%$
z;6BvM<;bY7_r``Fcq3bzWBE7Y9qa5_sMu`|2mf7fU=lfLT^X_)v&k1Hllo_9pyD0w
z7tODguWx)&Tcjx+F4qmGop*C9DWOGPVaNzl8z|Liq$tN3(`8OPZK=(YAb0qXyIH+x
z-1X%Gw`J`)iP(^f_oT+I<b2LH&J|yw$Z$}mSH`fUHH?$X+1l~|<VrIrcX68Uj|9$a
zZ}tO3zd_H%N-0@czPmLD3|-CId*FjjU_byf07uBH;3hk`Gt)ip06Et1_~_=rAtr*h
zk24JXZw*IaK90}*j3;iaZ2cLPN7tF(-HX5--g$$Sm;VXm8HR@^1j>o;Jl4}(wTk<!
zMR%rC00M@XHbqe}Ff_w^e7-d_5HNFd2mk(U4KJH$0{r=-$l?^dK?uVvS$k0`)yb=9
zQ%!k*oNZ9NZ;}eS`2IX^2v7f*j)=~W^{0hs8SaO*AS6bOKn20=t_-u*gyRgt5m4LL
zFPI0!E-fjET$Bv~^6AnGdZCb-8p-{=#m^0P{wEmFn{8d$TA?69?TrqQv#=MA&uf^v
z^NGqS+xI1GrUwRI0yzS*)f1`n5%8<3E_8-rV^h0*Thhs?$V$iBnoGoWq2m5Pyof8r
ze?4m*MaAP_JEBp5Q|P-j6jz_iO`vY%_r3)$qVl0V%~g7I%2<0l#7kVH<TJXJxdF=n
zZai}A*ysBC`eLr@*Uy#@wPdP&4NXnW$T&q!t!k}fZ-2$Y!U8!j>IT)>I=HhZ70qMd
zw$?6t!)wd_wJ#-$>WzwG8$smezJMEZa`hO0B#6AJxpj+}+1VSMP7~;MTiz^;Q6c~*
zw2XoNtOe<Pe0;t}oe^toZsz_V%EmTTFhFDD3hzc}Nf2v-3v?}T+bO8qa$?M%X%fgo
z(=CPoIn#F{j=YHmukP<2Br`I%3AMrD00P<{^xD0DB(Bc=G}g8XwPoJ(@<VM!=ap9x
zxVy5sS)FQrTI^ur*Rs(-k7Y685-6qb+`f&ghFhxeGAIsE_?qnPtp8aK^{tpygewI1
z@9Xc=nj0E2!aqYtrGH(kv-LXa-}$I}pi%8=tztb@=56rg{<PJJ%8v?hedXf!_5CxH
zWU0eP8dAvR%*`)ns3g{n*_Vj_FnjReUa@W_yxE9#V87rjT;@_Ocd<ddyu6-MOh7#y
z)zs!)(ABbn*MfY1b#+{;_uyJ|fZK;)w*A^@e@%ibER$5TTG>|M&c0Hg8;pa{Y`phA
zxT4UT(bsMlQ#D+bZW8!yakSlk;35ml((<w~<TrVy@Hz=8X=#H38=!FwPPGI9ID|LP
z&?6ojzZ>A1*`%)P{vY>{^Edu_c}Yn=!UFaRm)?ip1xCPU6O&p_Ka(>Za}qwcK17am
ztx?`nFCPORrmlFgU-Vc_XkVsRnQtioq7_f~wR0!4yX6zby#s=ROh6`_nwECeZderF
zM|A9T!q{8QT)O@J{k4UmSV%7T+*j_Co;d#v3h0pc4EpTOBTkFM@$fo|mw^9O3m+&A
zyNy!-$eZrT)eF3a>WQ!QRR-kehfG5ZH)7QNI?3dV+*RDA>M9~Vr;Jg(I@I?|5Pbc&
zgO>AS%{GOlycLNva*)Zn!LV?Uq?|l+=4;r6Td(*ns4vRANP<&|IL&lCi{vC{^8w9H
zB)qW(?n0^@?C<pzThRg)?b1we?uG-V(sOeXUPdti+vMwv(+boy)S??thg;znznOin
zqa}E2RLr9ybszzrkb}pL*RJ>AqMoYHsDGD@))qaBYPb$m7CN1duy>`%QX()Jv+ID9
ziZ@++zZ~f8sgI2h43D<nd2*bM`to@X`a{f4c=hVl^zm=f%2RP~_J3>MpNgM)v?TF<
zlyWg91jEw%ZM+qYT>KwI=^w@#5gcRtaU$eAgt*g3$;mH*cLVWu=bc})k8=dvuCItA
zQNMTNA{;Ye&At_*c$(b|;>pBTDHFk~4pY%#!(qQu*l_sHCB;`sgsL~iz^5TA=iR>#
zZa;!$b;w0ALgO7v7do{MRRoe#W%H%kSUZ^PpzP|L?dXQgN5U(PQ)G}N9rl1mHg>i?
z_rAgHKBBmSMAFfTyw<tDXjoX)ov_1etLM0eH~Y?fpb39|7%0yF5?Lu-%$CK<4G`4)
z*idONH6*Lelp1Wc2`=l*d7y?24+_#7D9C=<L!))X)r`)DrF9b>A9DfZ$>-hGdsFZI
z)B;P`c!OwbdB2lDE1%Mts)gd@q~zt%IP^j(DDhLd$RvTi!d^>->i|8S2kVN+@xb$P
z6bG}2;*V?1sK;YRqvo}u%iBgTb{0sImhk0G|D(6-42o*o(hZ6z0>TwZ(tvOURFEu4
z7Ex5vl^{r#C>cotlJp`;NuuPWq>JPXf+T?!39T2&p=pAEq^8Mq1GA3r&6|2PHC6Lt
zYP$4~LZ5whSmFED`u0iyS0Gm)G7=-bruTAZN&2V(^my~2YLhtf<wzWt{nF3J_008q
z@)SyG%?$T8jBaT-#Nc&_MqtdIY%R6eR3t_5=!@{!uvSq%qhPJEEC2B<tQHAvW3_*H
z$oTjyrflqz&$*p+JF!vJs_(Cu3t_T`YTTuv`p3eM8!wI54|*`wJ+lwcW9tJsxx?Gx
za|MS*qVW9UcqG=kP(C~C?K*4uazgq9za}JB73DY{3uU8}D88vjoDQ2Jo9~$?B;3Y2
zI;B(TVZmnXf6d{m)AuioHOx3GXdW;H-K#Pre9G+iD8hBRiZNoDCeJpzqcbSXZGUho
z3~$$g-E@;gHo%>hf!8V{{$nq8Dwk0k_4sJV4=z_D?IlD~)!QoebK1EP|KGfxV#bGy
zb-5A~W*iG6{Jgrc`}EBi@wS<Ec^@ilQTkQR<BzH^pCNJK-U?<eW_%v(Cz=~|sFN}(
zhrn^y{f0~q4qxO{v9R6F;JAAQLaz9+tDz^P_R+#G0zO8(PPm5O52tmeZ*CJRvB1_K
zaYWON)fgV0!DhS%VX5X3t(oG0<S==WqKoxwVf4a7o;Oi72~u&+>~ziF3a4RHHx4WD
z5L14=lK{VB=+mxHqG78nYK-A~qR@6xrp~w*j84~d-+71E>3S9gYm*&~YA+|U>H(3+
zte;qt^^9>=+bi`rOh9?iIAG#Q`fk^I_(I}t+d<x84iuV-ib`S>&v~md1MfV{C2ymK
zs|g7`X7?h$TKn9%D1f|U!KAB<xA<CyO5jjhQGPgT%>5JB{p~z0lgF0K+bj$WN}!Mf
zf|w%c8^zb^Wt>>Y1Ef)q97WxiRTP1dTvmnsH$y|8arI8L*2;r6sp-<$QNkwC+AZOc
zV4eGAipIUni`ubrgLpQ>frz8#Ga)B9Wc-I5xf$z%k##A7xT1S_^PJEeu3HKtMZsVq
z8?i$@Ad<VU<Nd_vtfl}@)6D8Lmh=(VhUp4jYFH;*SA|%F>q;PjV!~RJxN7ZkAB9}`
z54+x3>uIO8A9eIt!hWV+<@@ng4>`tNI?=a9Y=Z|bD8DCu;^Ngk$ez87BrI->p?xhz
zqc(kR2R()IdqODu;P~esq<aMAxo<l@BsO@mzNKSmURW`js(C(Axcu<@HQ2e?VaTYd
z4&{A<oYhurcJ|7^CSpdXdcb0oEw_}|aWPxUq=&r@%BbO5#T_s=nl=s*o#bkXB|}6<
zaEZr}aYM>$vrG;e;F@UdIE_2<++9Av+gwuElPwfUdESN@zZj^>3R~O|L>cLvW~-b!
z6mVr?)jZDwzpI@;J2NNZv8hAGjo6m446UtiO}F=24#_-}bW%1<isGLYJb(T<$gssg
zxvUHn&5AB(SSVH0Q|v5F(xmVBN`MD&y%_ri37{_)W!Xmed9Il>)XZ6;oN4Wq$(fQb
z7YP#Jy+iuKPXGv}o+^Eps_jqj`wfYKq)`a{PNT#hH)7RrT2iTZvfXqj7-R3)jenzS
z3-}&~oX4^84KH{XIw1#Dc1=E>oSvVSvlaDo^(&`{V?Xj6YDm0b$Sw-Q#>>~`(4+=C
z-s_Y>&)v6mnpsf&gPp|T2htJy9epXZmdck~F;&Mdk%Xr9>^D)~9zEiIo)4?=x;$^<
zU<Gv971Pyn2LDzR^_o@#n*s<gL9Z#A4Y$VDAh_o#FtzYg@8mRYTQD~_kb$a~@|Ymk
zmpSTt$x-LE7v|CWDi@y%pDKRWPj*I<In>DW>T~#~Rb$_3^YB(Y+UWpwty|-$3_cdY
z+HH_5j+gNFesA-Lb<2j0!zl!36>k=eU%_(C<fX=;$RC&E92;D&YCL04qT}be=}|oR
zw1iqP%S%ZpReE<+3M=Zj&q_q^=4OTK(!df<#X=(0b*d=0=b4d@1BnCm^JqwE?xD76
z<F3#6H{<MAi@G)F!@9upJq{C~4zlcySfJk7ba60n>qyRtL7Yob#>R9pmXjO6Ak@K7
zEsHGIUVJ)TJ)0)~nWTRE^I~GZ73fK26w)Rim9$v*zhs;+Zvw7LO1LiNdfq8OMW^tB
znsmN}BEB2*6wMIl(DE}r(gO2hxzv%cBfYCLSpMOur+0prL5RvkC=}i%YN}kk_@puI
zt;0dbvyT}lhIyms3*&rix4J?%o?OqS{VPb98(wbk$NccL`taVd)#cs9v1i7p_q%%5
zy7m0KM%{;Z0iVYQ^JsqqEtOnH|Fk)&CiSNksG!la<6lxuOLR&M6Y(BA={KmdDFF|=
zko3tmO3ZrkEP|}csdAzwm|a4u*wu~mVt?TfZ^LR+f-BLq2iw!X?{Ao)U<-&$59HI*
ztd>-q<H=F^$G6nUbn<Vn8<J%kGAciN_OpX~hD<$t4W0g@BUyuvpNa8td|bvS_mo`3
z_DWwMCvIRmK6+P>tW4y!vYGz*R!zNlkJ6un!G_o!t$XPEtA;CJY2JROc<FE`<^0!h
zo4GRMngkAwYyDIZ;y&N~0QST_;wd2%K1$I@Ve=`xR`g5^jS7{@jQ0Jyxx|dh#G2|J
zd{Z8MZ9K_$m-t5;s>91WBKKM(9Ul{;me=KNj<^Hyw1x21552Wl^1p@Ee9-!FO2Zdc
zcPd?1=_11!hFE3#yg%tMH6kupjNo`pD!?Hoi^}i6|0uh6YCW(h^bdZK@8)+H^VPQw
zDx)z?ob7v-XHB}CT|6#5LaY8==4|w5r^;9~%`K{Tlss?7V=FK_=2(jfEyFelS}#;z
z=F2|AYr$!4UeUNi+BthifagsTfTQ0658K6Y3R-u(k{x_XLTU(*UMv_k#kF|wtX%&!
z;jv%lJ$l)|VN(Gv(i(EwD59ufML!Y}`Dv`sO5SyY!eS*}>&pOMLyIc`DrkPuG*P;L
z03Ols+aa@yQ_82|F~LEi(a*;=9whw6LsxMHQ-x!vBXk`(1Nk?gNxdTqtyb@~JU>!(
z2_$Yw{Z4VOUe#_-JPD7AvS|*YSQ@U#8`=Xgx4r^>VJNZjnF{BZ28ql27?!nLwoA76
zs{^kdX(XTYw5gi>R;5+{C6wx<Jo{MpGuesr)<e!Vo4%1zo}=&mZq-5deGUein#$34
zH=D9x#VnFVHhPxuRfW%EPZsh?)Xt9`K4-HlF&`=n#T<LY3+0JL4ijNv<ZcLy)srQ{
zgM*Jk#_<J&6xk&0nLr@~Up?x4sXDF4CJRUmgc~^6CsnxOl|ewa6ZnBL!y;gY;-CU$
zz4fFp11D!MWah_dX|15bC?{UeR;f*Aczv17e?A>>?$WLKW7WBQQ?4=*PNyH4h;v48
zG-nbHBnksQzRlJihF1Q6!0!JD;lom-9k5f&X!oRakvR2Ht@eR_n&Y`zeWMG2!O09d
zX#sZYsdVF_eDq}~V~GU)OED1M8L1XZ<=z7M*mCPG<{#J$GsjWjD#RQH1VOiN+He50
zbUB|aiSBJJcL8}pAQCOUyuN@!n3IoB9s^>b2K}jb&X>@H)mDHl&`edse;%G*QuGYl
z&pTCi`28Gx>H!x?i>uYt9K3w+z5Jei{5Am&!$@%D-znv+d6Ya&q4K!ipS-wIE8Lvl
zBLc2H#P-aIBJGr3&rx@s@hy+zFWv$iMd^QpF8hoqCW~5t@82;9sDJ+O-!O<Q7hIzw
zt}zXYW9kO_^L0WY=kQt@x@%9RlKt*WCl=pU_xhw+elZ495inyg5Y5peqiOH;VKKy{
zsRKJ24aoJ^3VhFo<JU+4W+6;T=zDyQN4uKrPp&Jlb6K$;8P&Uq9D#!9LjGmHV<@{h
zaRaZ(OMqI^2sDJ-1Hd{mXa|tJ)5H6ho}K+O`HxH72GrJVmcYDqe*VlYAaI;3BcrJY
zpqC{t**wj(cOWiUJPY~11eDPa;ZD)h{{RrB7$LzfR5wD#2s0j!f6VQ+44%a&O3sbA
zJi|0&-X%=_(s{VMw(EBmjLzpT@)$T-<!;zV9Rcx`>ulLJd2zl9X|I<gdWW#_JGhcF
z4Be5}G-}WNB`P)GsEjE!K=#EM+(m1`ctn`4AWH%x<UC*=9;`<zmM&n591qs0y{!50
zZQ-95%=K-4s~3w)tS9jf;!5ylY4BREz^<_0E5UjWBpch(y~~vhTP&2!w=23v*Yxu{
zp+wWO$zCRr*SS6GWLS^QCoHyPsQlerU0q#!{M^ipvDfIk<>JalG1)in$%pBJg7z$2
z6vDn_qbzbn)QZ8*W|myLY2QLyc;99+i(;c85H_2uCVn3I^HqTX#@>n%YK@o4DL1e$
zCX{HJQrpv6q)%S<UCL~}J89^(woYkFu^UA>xbCM);S<$x(rkC;3=_amrExmum~<D6
zY0MLRf$M$E#<Wy&Y6cg}jIkK@T_uFiCr;AlnD&dX>y}_P`+t7RH2jvaSqBq#ZGhdH
z=<*`0_Uwyx<+#k75-052{){~m)at_<z974M&hsopa=Yi8WWmF_@281xRkY<f@hdoy
z^!Cr{`5(#0n;G1qd2w6vJxMPsB+!%6zkp%V@_aM@WNXGf@Fj#bV#{C`3zx8qpP<5y
z_wJ<iTI=Z4b4Vv9;#*5vw|g>hioyY0LVE|s(H9OHIkj_`AsX55l^00wfZ1(oIy_JH
zcHT!xsaem-p~|qte!m}ziurkYr^^?i#L5-u)jIAK$Xi?`7FCk`zrK03R!`7Fn>E*{
z1P(r<5`F=HSWKB+{2Cfg;(`SyV9!rmu62xd&9xHW@jCgC+a3f=iRf=LU(hatT7kd+
zx(m1l>!Z$V{ITWVUnfhr4G72)Jn)%1jal<p%0Y?}oYz5XV1D*%wVws24RfwM<@DN;
z1u_kQ$;2ppH9(2MAt7!cP{2V?ZEE^lpQE?CyUo5z8_bDX^}e;}=kD%oz$a0`!ND;}
zNj{x>X7pVjUV;9ZxMjynGZU9G0|5bn=Fw4k%^ctR4G(8->I*`z54Wsr3P7AS!<ZEn
zQq7VbSAvj>i;I2#`y79xkPUJ2VuO&W_Y;!kX5&^$k^N!(N6`9$-kWQ_zP{m#)OXH-
zX@n+^ot>Rhf0EZ=zdH22067Kb1)vmh$g?fFemiZ4QiRH6nrY;YfYZC6t4;%Os02u<
zZ!C|p1VQmER9hmQ`{81A9MUmHpx|D<2nC*>KgaU=K#!Ll2&ILa>K6ml_@s6&{fYPC
z{&rP#SJ0NraG%0uqmd?SEmEzc9gOXrzhB*Mi<>9R6_Dc4bbj5~m;-2BK7fk>ycPwI
zmr0<#3A8_9dpi$w1=KSX(xCWI+;vik`LZ#bxaqA?7wSon+_fc4wp2QgCqp?R2uZ0!
z+NkD7RR{$1erW{ck+*VY*B&c#iQCWu6hO-+77S>ys72dx!0b*a0oZC#=XM?i@^IGR
zrP?r<J`k}*_-^pSChl^G1I9mR@rHL}lX%E*?b<c<bh+W7?el<~^Grp_?a>}*N~cN^
zvyNCi_KZbMQ}gmx1K{5X>Zfm30Cu|goB{HqJlbxYfyW@7h0jHL$RCi;1t}z(LjKda
zfnPmQq{ZX-?tOD}?&Z-Mei4zc9!$|W{>IgnfnDiakB0v~5PC9#u~z%F@`nn{d3`=P
zqxQ%P0-8fZ*8*h$jEu2vS(L~Z^6<>wgbWG!d_i*UUqAGV_Wbw9BFQ^^vd0iS13tIv
z?o#9jiHQTI*vW4<Sd2y>Tv(i3UZ2#O3EvmM!GnM^H6mC6iU%61L<1m~W<aY!4E<@^
z%Gz4gc}!a7@rpVaeO8!oP=e02tgHHOL611PGYf$<0pgbNm{+QC9!~_}*nsW$x2jzJ
zuK{6%Q}hTPfJ{B8)*S&1e!#ZUA)8X0B9c=+xv)?WN(b7YuG65}k&Vlb@N~g|({g33
zPSChE20=t^dpQ~p51t?e0)1t2pd-NpC_t-48Bzz39(drhwLFReeR%KHkK5AT1a7SS
zcItwG^D^{2@j<|JN|77`v^dCxwLqn=28PHAiYXZ9nsrcE5eGd&P}0Vpb!m7Y6R!&D
zfF(nv5B^+W<+S_xJ>PMnAqjLRh$9XqYY-cjLM3-!*;o`rsjp*Vnllc`w#iqV?f{Zp
z2O4NY<<>7?+%15F0w5|1!uumNE?UsB$EM&X4MM_K&^^;Lv$KpWEUGD2oPeR1hQJ1T
zAqgOl0sQGT6_ukv_*gHGRLQ*aiv*|)M2)P;J8QwAp?wv$S|Rl4FK|DAy6XCIphGO?
zdehXC#GV5M17TGcXzJZ>^s8HV$T<g%W$GzcBEbtM!R#O$0mzYiueGy)ZFmO{AkPX1
zJ^<B*+^jk@S0BZUHpr39fP`@7^-y}4!;sZRQWVhKDN(bGTdO!MbbHT#%L;WEPJ2Aj
z34iTG-q=EeCNt=Y=iSY^={WwO--9E1YIe3A9t@n$V<C5$G$tk{R{+k*WHN28yB6rH
z!#^;f(aVO+FxJBR6|}iT4(M(K6`iZba3MgkwPv2qr(<2sUBXRZ7W{8%X-T2mr+@$b
z_dkjZB}b~A77e#Rr<+xMrqSbAHYLnCJ1{rAe0(9!;~yBIBnMa!XIMRHuO(JV`|sSM
zqOm|1Jc9mI5x~;Xl9#QiPjOm-DSL^Q9e`6fv%0DUhJvuU=_!<=7^nd6!Y?53*<5V<
z&3{%dDRL#8q^Iw2u5thQ{lcG}a@WCd&5zYpEY?BV0YW@u9(@|ebbQib@<7}!6H9A7
z13D}SK?OGuyhLtD0E<b4{?17c84Xo;S(zgiLF3N+LPE>si}r)XJ)lyD1f$;2n0Z@E
z%XX?QygOY!1p>e~?%LZMMkJ|AmreqI1Ti_Atd}4a6_wl}c^w)lI-<GEdT);pz6l7K
zo|<A6wNNR`O4An8AHh&iP!!KW8?AZ&d%f;7*_wtaAC8nBug;{bcXKjIclS(vX3xBe
zi`xaEaShcqR)}lBLOq3ZkwzuN3m0r~DEK9Up1cScYT&rMSL){e($1lQvs1d@?R<l(
zEeI{^CE`t(?9%M^|Cyh6OI_Uxu3tDxO?~9K%S*&+fSzMZml}Uj@l4IuD(F++gmN7y
zX?vWejD~CjxtM|G)(P@}LI8+MLID*mt4Onw7<SEw<PGIP#Sel(SBX*h<Hm2GIg|y0
zI<x3L%N-KQOj{x%-OPAvW*ZI-JILOg10ir3jg%{rXpCh8ltBT1vrwqbM%>!s1AhwE
z0)7{b0zcUT1d5*0om((;un%#mseZ$r&?cFQ1g}I2$n$P{d0hblYJhRJfM7)04&@4H
zwkZ3Wef#uaIW)>z){I{DA!cu_>g^Tz%Xt&by+8{m*)KgIA@mcT8}{CrOZRF6?7Cb?
zMgGQCV62c+2R90=gRD@zju@8dze>f9$p5QU%vm~Z4&hl|un60LyvoMXa336>a*&!q
zrp2P^2?ZS78_=EB`H++M_ur2O1|C|hT0-+AABatwU)%<h*#ZJdL4A2zEF=*>!sNC0
zy7});VRFFnP(o=Itgwk{IC$_|=H{6ly$P>?RXtNbzT$fT()pzhgO}k>5<O*QWm>*<
z!Na}5ChyIDeKslQ)RhWt8t^jzB*vgQ7{vgL4|K@lt<oWeVG}e3#Em@7wehpR3kF<a
z-a-kioC!8hBjs2#1PdBBh|-|5gdCWe&y3J9eVi#E<S+e_biCDyDpmF~#Q1<ue#tHw
z0aoJ?appEyuH3AYoNS19U@HV6P8_x#1Xv*L$bRG)?N>szHnHm!BNJ0F@Q;XI{#%vw
zU)BvWS!ix?gF;9MwiXP*fDA6!?y~<+VMpILxmHMFeV}NP?TRC-(@5ooEpelQdqQxm
zRSgUb%I2OveTv!}H;Q1Fa#jc5#VYSB0hh*DSy_pI?xn_F(=q0}+!cLkaDr3qEM`{F
zaH^ZH{anF+kI;5nYjSznNhVqKGZt$L3+lHuY>OHLwIuTXYM`*!jdK^|L3Z98=*8{~
zMI6|kGFWP)(b=L4awFy!&p}!LZjH0a_0BOPI4^LHtl<3W7Tk+`dp-2$lJ<pN5OeQT
zJXnji=!yuifmp~v;EnC24=4BFZWdbpH9_zL)5xGep7ZF@qi}P~DFUeeBt^}7{VNuj
vm;3lD7<m1^3>c8W<NqY`_<#JOl*3a5hZAEXMUna__}spssZy+L{_H;h*93>V

literal 0
HcmV?d00001

diff --git a/H3K4me3_distribution.png b/H3K4me3_distribution.png
new file mode 100644
index 0000000000000000000000000000000000000000..6ac03853fd8fca931b5a86c8f1386c7655b8726e
GIT binary patch
literal 25101
zcmeFZbySpX*FHRe2!fOW(uf#H-3UmBh*APdHz?8#Lr4xFf&$VYEmG1T>Chse(nAa#
zl0zdM1Hbbc?&p4==Y79-{noeE_vhzY%3|h<GxoXnaqMGn-#%5AB_p9Dfx%#8@^Uh&
zFc<*~42FM;_&oR|-*Bx7{5a<<Ew4@t{&*3;3<kenaFEk+hQTPFLx1r;OXgU@U@S0s
znFs3b$t&ZYUNJi<r&zZrd9!c4l*AW5qTE?--*mDV@Q(@5&9Ca&rS@68i&S(eR9&r{
zp#C*+VMpuS^4iJD?}yejrJwS(M?9@vDhd{yIqt9~ND;=6kYAs^m>@5GZ;yFWKjGZB
z>wls3gKAU5iNX)e<`b<YET3(4EqMNrW4m+#4~$KU?xlzU_z6R-1;FsZPZnEZ7U*}>
zoB#jU|BqIqbq_DIwDj55d}rTKsfmEku{)Pu<;n-N#NL=Q^7lIqey3Tv&SYV=<-x+`
zjl->O$>MhInkRE}bI6(ZqUGk@Of~I9U`M=K_V8YA{%(0E!>U{F_4V%EyX7x`UMkeB
zSo*xIS!`H%u(e>_l`QNr9Y~{YXjruFz3%z0<waDc%kp4F-N{N>x7B35ccQ499kPKD
zJg%|t-lexP-!QBzmqJ?Ly^m-tF{ms$KC@6Z%yba*pyzsB&h1q9Pt6rPb=y&ucH?<^
z_H`Oj@QtYlVYGaY8?dYPE++>I%c-7g&&<rS*J{=!<WcV#Mfn|n%hVnJX{Qr`TX;0^
zLNChj8pj(hc^;6Ml%!#3II%mn*qafasC!d>@_1fw^qcI%hXG8!$Btyz{*qcqc7D%h
zxjIroC+cc5x)VS|p((^L*v9N@v*iYv>%MFw27g6ePcOf=R*Vq;r#+vfq!Aez*+`jr
z_~Fj7TfgUo=ejtyS@Jl54dX||^wiGIZnTbsghWJN)VLA0I#k-+x(=QL+n_NbiqPiF
zPEJlH@76IAyno;IFqP`Y<Hfy((>g<6tjoA_w>aH>n-AR8Q}``2o674h1xL`0LW#u$
za(`VDyUU_&q?ZzQV`qKR2z#`r_3T;XV39$7PMVJ{jlM$@@A_4eK%{zKSBj{Hp57?N
zdTXIu%hk0q)q8iCd;Bsq0unHS2Lb0T*C*<%EG?gzW$Ub|YH8&>y^|^7eF(4KjaHF(
zrl}dHprA04Dtfe1R`F`fOkVZD1HYxEC6U^kqN26WmO8h;yl5fe()lv`p_o~R2Ij8A
zPd8BPv%@j&v0O+J)T|s&0CT_tFZ%ua1tC6Op+Svi#^bQ4mzASj92}xcBe5J>R-=`+
zx=6U?sI!KKhAlcSKHj=BQ6TZYZ6ezi^8ll$tHtb(NP|4f)==6+eoNJy;&IneTiXX0
zV8<+UVOCB4gh=34h*_k}D~qE)QCTh$qCQfon#J{9#jD(Mfn5lcP84y;gN)T-T6&<r
z-(vWSscyMNBwJt8;&@G!f73^^jt}XWz8vd=KXb@!l;Td1<ne<vuRm|s8cu!WQ5^|<
zE;C_*R$vf!IXL>h1(FIl{{GHw-=z2SsY6GRysPUj#=99<ndon#OjF<60u1~{L5FhY
zVUcfqHAMF4Q6IM$xF||O)_un%+}8ACV`D7`3v`fs6P_*$DQ;H3V%VcH9cNlGEi^_k
zY_Tyhi|Ij3D0|-{J*5;84g0#CxHxWu5o+pR?d>~dY1k(e;@c5LHESZ^bXrb-3w-Zw
z1YSM!N#p0nb@@#!s#PmXk}HY{se1VP#)1`r&Fpdjah=(p*`HCFT86J$rWIf;0usj?
z55t(mb*%>rER#q`KMvF=X3_X^o4%AoOIYBy4N_$Cpp(!I!YYS{Hs;pnY^IjkwSm19
zT`jyHjtx>=FtR<7L-EveuZAle)x!(;8<p$R{0U*ODTUs<Ti3W(_q`@4x2oe|_`aVA
z@L?rNj*g{+#YV;7zDegO;*DkHKFAw?**YEdzIexWocBFBvDZ7`=}VOC?1~f=6;-^h
z_l&O`=l7#*T4f0%2FZWaTX<OGqh902fT0d!Ym+@sX-zViZlJ0r<}H9{DGdV$Q%=wt
zE4%(cwHS}lXnAxkv=j?}FeQ1862ukwXfKjLvwD9KUBIq=iL`fG@$6$OOFJ4V5uF>R
zT+;E{!Q@7)2r$XazM?`^lR8?8wM>GjQLSH?=dEz5sfTBGjP%;!9Y!dNu+Y-3x_)Ct
zH2OlJYE&lUnTDq3T{_e2Ts^J22&R>rb1A(C0_Q0Aui>4887Eor=Mj=pQ%WIr3%J~g
zXcQCL_Me6J^-P<bJyo793NN#aaFIHn;ESZy1i?NrxC_Q&K)%9HbWtJ2zR9{2UrtRw
zZ_>oqK#Gn^3XBq!g~i{Vagq8nl0x`RgmVAl*&X#ZWGp(vs7B}a-UdfyQrblLY1a|p
z^30G5g26`TEUB*2;$=P|;cS9`zHq@mjEI^sq9InOJD<plcTVCVSk0kZS4eU*GiB1B
z{=lfh7FzZcYpN1Ar!JqPd|T)~t%sLG3LK%9y!_j|yu4&?0*4?*`ly)=`EwEV4>N6<
z`|$kMd+?w=39_hREIH3?PHLz)YPrZTd%^c?H~|%~x>Ft1&wFo+tFnz5$JiO|%8R*K
zI&py?)8mSNDF37(DvB{eMq1>liHZ_}0QZfn=-8bByL&Uy<SfaA_`bOP8nWHf%>DW^
zl;~>8nCGDroyhT92sRj38R$jN92X~#d3wy>rr^jVz-Pw28(hCI&wJzHlHm62M08}Y
z-EB^M;$sMTlqj)XdS+>v%NqIPEzkL5o4v^nlgaDzF<a4(V15+PvRA%sn*p)!UE7Ru
zJibu72p*V(fu|9Zx!;)~z+B5jSPaNmG>KlWyOuHXOc6jE3=<U)l=rMka6Dl@AJO_y
z^;pVk8J_`Xv?YBLtE!4+v}sFMIORvG#9ajN$j@f=GtuX!Jhppf#Pi63n)o;fH>cay
zq=uL)g9ryO%fCNaru;TRsHz-ZR3~^_YpHiK6oY8IjWf#E&CfN|)$cHnr|rLJEFN)f
zT>8Du5mG_^8k{%QDLB}kQ{bNS3JT!Ji;X%t9_NFz)+sUVJ1Cznd_PIt{*@>yhLkT{
z>4ueVXA$6^?cLB^BRR$)gEq|yRtmQ6b7Ue#PBwSr1Kdn&v`KriikPJ>N<~<>I#Eq2
zIYyVK{=z2)htNPI8P6#X?=`R=p1-d(N6gO4fAv;bW%3?nd|Vw5PK5X0U#XQSyBx-{
z?a)Oqpw}hpJYX~|z>aM$hqS0`Xk=+~&coReeBd-=l{YU^%Kz$!@~`U<@|oJU#%Fj8
z*>5%yscUL#npafhioC*tz4UgLmgS&XW!*4+VxG7DbbF1U8U`JiM)IJg%+`lyUU{TF
zT<H??kzjEUxU~)<7d1|4y~6nxd=hx}iFz{Li}Zi&s7d!UvQKSEcxoS=Di=Y1W<xMz
zA`Zm1q}s0^LUo>VC_F!3N?<pzl^$v-W#_}#{vrZtnBN}m2(K7(NlT9JFss&(|L`fR
z<C%K?_p$m9Z;9flQ~?LSHFpJ3&JdUQTUB`tGkMto@({*JV`KWUad9FB=OcV-MpqHB
z?_e)h{ck34gKwF@(7~#!Tl=1uyPIPyUqIlufP1ZKrw9RC#El&$ufx5vY-4$w7Siq;
zU~OTUBdLQh)cO;4K0-;miZuTLDWY=|&|<>mDh$}i{I#N_S4CKOMOZfOLmMiUS0VZO
zudDC%Uz58-UJk|*OShawAlfN(m9*z3yXYO$xevS{J=2b0M-jH4dOw5Mc{S{D3OAu#
z+ODW^mK8}{O!AvYMCX_mj$o*vq0g_j<}Fu-O6eqhYU8-{Tsy3d00vS@y01iI=$@Xf
zoRX(y<XzFcY|<L~46M+*ckk5I)iVL8XnDa~(?AH@DEOC`?Wrr5Bw+J(Us|&!q@+%7
zp7HXfoc`Ssa`mlDxz~Y{?P#SEfW6^w$Qkq6c^Yb)UFh%I^dt*AZ%3yc%fuUc<W6{P
zDPgCGBzx9<_C&+3@k$E_2#CNJBMf~GOarCUEBEH+ohZS<wlGHnSHj1#!<dLX(p{F)
zeoY<@0LR;w7rI}avY@CvSn)yJvoh7|&(93}AEZ3BA47QRwb4etb!?c!<zjSxaeW?*
zk=L^uk3)OFocSH5q{tW(zxhUj?KlNT2>V7}hyD0jI|&O2L5}t0EMTNhot?`76c|mV
z7xyUl*~vFzbXl!xz}kW_GVUBV1rXhoPr7f5+L&&(&A$q~^y>jU`ME=i+;KWeeB#<p
zw6&C#wRLQKe2x@LLfN0Kt%07NkK&%Y@aoEvlBZxG7DtDDDv|?7>kY-wV57<b<^xMS
z?$As|Cuo;!_jQii8=!8>0|{VWLE7esOebtm%nE|?Y=ouSPlSbq8B@#Z?QCqs=8s(0
zCklm~=PkAudlSW8ePOd~kK-C?@NKYbz`BYbZ#0txN+U7ZN$N&M_LG8O+O=8i+^hEE
z@auDUuMJA3DC$%FuU<SKPC_?`ZqpSr1n4x$zFvoxmeymViD)#{ccrwge{8J4_^&(i
z9#eiq{LjCgtJ@v1MlR*$D<iX0MBVZ#D=YiZiHV6Z37F+4#4tod#UviboLWf%IXEU0
z^UJV%7CkJh<271bTwJ!2L+F74b`bnj(j;UJjMUUJ=<REt>FDS{d<VdL^sX)o<>Ftl
z9Jg=FHoqZLQ&Lijx6&<r5g5m*<J5r`%_oJ`-o?ev%;q&eMDeOs17Wd}wG=oS^Ub*i
zQ{<VE#i#SsG&CHP#1;_A5U`42v)ox40H*_hVL%6l%gi}MMfFs~HeVm^txrzxu&ZS+
zz878VDYxunli%6d89}<_laP|usP~1^^6^`uGL?;LRx3vr5wMNesHjT-LJI?jAfZXb
z#;5p;K5oP_XapPeWg5|+F<DwnsNaFcW4l6PxzwMV_xW=V0J9oklBVr(mc5@IMLiMB
zbX%VwIUnFL8!20EJ9=My^ZR9p8W<QDjEZ~gEM>jDe0wRr;W+EfMcPk3*h5Sc5wkkr
zN3SnjMyjg;W&jZ#Z!S^MqwV_V01$Fo6BX1xSJixXS~DSm_<n(9gWHIRQXX)?B5F?U
zk2+;$x!{G!tu9dk(_dHQQPqw!<rclH(6P_gEL0U^LT_XGL~DIV$}N%W^+)i%aW}AX
zi*1r8rQ{4k9~HU3%7jZCe!h3_Uh;k=Fi85<lQFSu4`7kW{x?g>&!xsif;@vV4PNtz
z<8rK)p=0$W7|cO|>DukpBh|wZGSmR6s^>c60P0~O{p8qjjJms{)jCDy1yXWgi>m3l
zP4!&vSvRtcd-43-&fMqo?nyT57jIY9N`yiwCm!MaW8kjd#x8T};nBi&R6#`!IIPd1
zIEGL!q5f$8QUR-4OIXZiYNpU+X!Zo=7jy<CWLWQVacQt!aztKf^s&|3Adkft-t~H1
z;zCEQR)Tu~!uI!OHEU!r__r?$wp&Jiz1##F)S?@#47}!<XC}-ic+p0#9jDKM==s+W
zmC@7Zolk<;=0=C`o|90Y#WM%wu8nmS;kv_-L|Y>8?XbV+qZ}`QRrumY66ZoK7b*a@
z2lDNW7^_i#!RAyhcP=m&9Ed^mnKDWbxdmS#9@!HHwpA}FJV&~54P4!Lf7(I&>eZ__
zu7^L=p6GA7KwYc)g#d-(6c`*m(m^QUE7JhBS?nNG@avaCPCIPHJO%H?!XX9wFkWN_
ze}S=zrzAxQ9wW45eU8=7Mry`G`)H*ABoD`dq2%}CqnR=xv^mx!ao?rk?$l7YnSyJJ
z8p~r=5#IsgNy5RyoMsS`Bd_m>_Dz6H1Ku=0-2KWwm+sJVnvPZz7RwNwtCxayO^C4u
zRZbB1vg=xD7cB_uxOe?0!tjO=Mo`Ye)(G=sVs}bh%w1Zt0t_w}i5vmQ*pvSZ+0}Q7
z8tW9eu$-g(`;F?W2j7wvxQD^%&nhJNc7@aZGHxX519tN7iLPWX@QsQ<SO4?rz3>Db
zOPm+A#hH$pM%A#`uB0N#Y~})=4hSr`7@d4KpWxm`^0xOaMoM>vvnR4$Qnj@vdO0<W
zIT~XDAJd<G%<wsqjy%y(M3m)gAzmiCStZ_iF({$-4;ura6R<JWZU8j1z7n8r3!BkV
ziv{DZI3r@QNlD9acBv(uzG5~@DLR7S0ldf$ZwOV+_Jky#042yt=$XCl9I#423|>re
z_Tq9bDU}(pw_i!)z*?uu;MV#oFz{?gfavDFoc|kDEanYKV2o2}b7}}Ky(A;94TC^d
zj|QV_qk=Q+Cx+d${Dc?{VGxLrz!vZ3R8GZac06GpY7aF5Rw@2JaRrw(o=3QkYWMNB
zZ+!)23OVuL(*X$ye?n4XmJQ|{>W_NTg}6vfftv<AR!wh17lZ<E<ZuucQDpH>5rFkl
z6GqBx<E9#%VHrHDdZfVCF;q$}w|s>YARrbT0@4!tN?05SUK5HN`qjb1@D=@Bm7xT>
z5$be6g7$FO5BUj8uAM-Ry))oL!)YyYZHfvva_RI4kXN}v5||7S5hK{v#&p6()q4mw
zU4Tcs2Jn*o5q*GCALj!a27J15JYJ^q)^kg=x!FoTw7TbI?io-nzl1ju0M8FY^jk=9
z+Inet*ww4=eCit*GU=YooWI5AYDH=Yoie}#yMc|vV9F|;k1-jO^3134@RBaEszXc5
z`Qpsq>|KL#Ruy_tCTO@6SSV^}enJACRf^pwQis1YGK>`2GAAe4I}NXsu`Bx8x|RaR
zc0TjCEIEJX2fH>R=r2OR2$K^R=5II4YF4T{1;Vk#Sq;9rB(ri#yqKTs@<*S#1``xq
z3Q2p8(oE~kbHEK}F3g~6P8`!G*x|l?5p(*=DDj{~onPnh7dq2^(%w+?sfn~)fQqXX
z6h@={!GWQHj0~>CdL64`-MH9mbT}`TE8FlMt?v-VkKyK<po;YIjPFM!Pl%{-MatPV
zEzcfO{DgnN7dBgocL=h>Wtfq-{zR3pn(&>pt}^BjRP2QbqkdSL)oxz39Kd37ySa(1
zGT#K4t-xWCnaP8NZuHA3g#?SNI9C@bCzF+x1vsE`yuZ4o1ZUkbk6E(d6FSfF(&Ek^
z5!vInQv_PQuxqN8lkvM~i+GEghPHq%=OS{&lAKq7gN8#K)aG({8vvtM!rK5E&i&J?
zkXz)!G4fL~2ES>>h?!Vn5=(RW@g*+Q(pBgYt7#wFq~5DL6qQxTsVZB;EX5XW7^P^S
zAkz5LGE}};!lE_c@dp>C+ZmIzDB<|T&BUue1CDFw+XrFys=p62mT-dO@(l7S4o(xk
zV+8>bneA(F_sWz{>P(6S!--rK=na)$Syn2t0m7ZWAgUK~6Oc=L=zJh9=TCCf?F7?O
zoaVtNz8BT`ewWF$9|}GFqi$ub7(^g|YSQIn`W*#(;UYx`g8P}>luGibq*i+$h?^<N
zYP~d*H!(}xjB|@#maF}?<RiC+L6j>sy#Oa3N2`<zmKr)|w7BECsyF)Acs-CYw13S|
zobU;_WQAOn`r@30M-KSpb9%!w&&n`a)Sa0aW~T2tmphr`=BMtJK<m3Ln)~x%O@<C>
zPf(144xYJnX$dhS)7<|nUePUiTFYItat}b`We9GiihHhbGsDEV7|8QKeUhGb7{aJJ
zh~%BTWa7J*`MVuZNh4c>P^Qt_g}&3<Yq8vKUzdBQ(AJn@gTCLd-Qb)gif9R;;kfkH
z0wfh7q!o2jvD|7PFQaG_<e}#~5;Cld>omK?_mZ<!(u$amr-(+XdyiYq$B{0HlD_rO
z1?Fkqnn$~<ePdsp`>?ApL=MF0H8n3u)fOQI5Xrp1;JnIhVB!8Wyl6Gl39WalQJ*|i
zl&vZ_jmR{n*mAG#t<j^;(EIZbBbo0`EB^F}RVm@#qcZdEL6F^j)nCC=yW!7WzyHF(
z7rP@e6Uqn2?x2^gjgArxy*3}oqtsNyvp9-tRnhv+VO^r@V%Xm_4Ltx8r>~Ay<s6OT
z!&tm`z1eM#hC&sKgs?{-g=l0eF77?5i)qj6YseWdSPNn!``Y+b<o(OSv+%L$0{ic&
z_zJsmWJVEox2mC#8M`PeEAG4T(bBpIq>4@2-rIWmY)mx)baA*p6K?yN0S1c?E>9AO
zcn4#ZpPN)l{kRs`vItzSF%P?Mbw(E47`$ib9Zl-*Wti$eS^Ab22L}gM78XKQ*vHSG
zr#}w8*1Nj{GUa)DlfFiPf;mDddXdT$qKnr)f`s$XKF}woO5~Qe^RJuE%WSDHM&3^G
zt)MMF3J`ui-dVs{@+L+;f%~vG*Bb`hDurXO<fB}gkd~O}=)P&;<>l7aR+D=By|wW~
zKJ)jF!<h<yM&DM9`n~C%vnMj)4X{vofD3AnqN7^q(Kto>XX{QAle%p*8bDEX1-6k@
z#M~#47ieSG1_$>#_podiQ<L$=XN0kLN>b7&#uH@U;ai=895-$}q>$L3y71o<{eYDG
zg7?{L!eHa0Z{o8c<YaXiJJ-Rhy8dAP_2v$boav;(Dpt*OZ*Ueb^r1od$wQG(1mBdk
zH^~1(5U@OIsxK;1PJb%reYh>WbeNb(H|{<Y_UqR#)StQbfnuYEN;JUx>RMXp=|dzE
znUQUp^5flA4j>9enbc3{098-uTcVG^21D*ta{q)VUtco&pNW}l0+hq&8H#?p>DJ3V
zT*Bf}CLqv0^kS5H7ZP%n_*}W`s%{vafI?Kc#|{Lw7+q{Rsry!j%VJYf3V|LZ4YZE+
z4Ubo^zJe@zhGk#2d=$toxULMPN8XTi&CRCeGqXTt%BR!=k)Wy<4{wSCa8kWWn`h6T
zPxFd)JLC@eShm-gB`>W7Wvy6UmP-5$<;+P>cswN*AWS*8t29pcLH)u#Leh!ZW!NQ>
zKz^4cb&yK40NM7u+}z$TrtJ_%k<Sk3t#Pv_40s^^=f{naO53DV_um&P%FD}LL_u(e
z$Wj8I`BE#regBWh8{pt(2a?iQ{r>(I(xu|9?%ug0WQJz@n<+UY09^U6l>E^nO0gF+
z`9I5Uu8^=a|H`U=ZfVs2q}j;?X4!CZ<3{ta?)#W#R=9ig|3y)#NZBg28o1js<EWG-
zVKCDcZhy43m3FEhg+Uyzm`nCOq!7<sE~-&}_3D*wwL`keBtYqv?`mH%ZM~77R`ry)
z!KZe&=2xhsxA@n=Z(O>G(F-xvwgAJK5=M>E0(=v+k0fACAaC{F8_zic^4n+heVd`4
z`<4DbHlHNjFSI40hQaCww0Uej`ha&R*%Lf9)skO74ifJ4Z01n?8ge44x2N-m$g^q{
zVR~UVid+y-n4JB32&Nzx0VF_P{8^9Dic6<qLu(fJHcfOd<?$HQ!p2`HH3MCN0t%&w
z_xJ0qm^{fMx=LQHDnEYgL^`2lpV+c)_k0yGz-P21`W-A*%gBWgTleYfPza2Ok&{cB
zxIg2}t<?Q<a!@2<{#LF>EH=&UO9A<hr!?GBRoZ$qd;px_^0V?4NgdOODEMymwwGD8
z`jT%P@bGk+ThG}?kg)N>Wy#E0&DtkT{y`AZlS0i`MbDHqq^5Rv{E%~faKGWkz`A_*
z0+UK->#WGNsuaY8A^G1p**FQn??R&n!<Lp7041U;4j8wDEME6}8XrpR9fLfa(-31{
z$^H5lOnhaTd<I@`lp-E2ZZP36$k(r50dt(Y0x!0(i;_4h1H9-k*7Rv4|C9~i$#=rO
zlDyM4vD6)nTv<tTYhk8AL}LJiVjK5=XUy88y5*gL8PN}#ugeoz^&OdQf4c#Loz(&y
ztgeoy>YqOIN^o{b4rKs0dg5*@$u=^b1A3Oy=|Xt26_Z|KmW!|4>4F~c?X!Nc7q-oc
zfK5UxfOyU&CF!jGhiDW|_w-JRV>$@UM*xhskg!avvfF<KkG*J3qUoF=ti-edgD~b_
zL@+xLJNhb)G|YcBm@rx2m+O@~lx*+`FW=!YTiIO;iWz&O^#U*x9CGVGhCtrT6VJnQ
z!HA2{sbfDQ7IS(WR-C|Y93_Y!TtQaMCNL<fMbtxx+w5eY*ikJy)p>o!uijZyY8~$5
z%$;5On3uMAW9t{o_~J(iM+qa|1OT`V#{<9;r_=z7`Y7XBj$+eW8!IxF!1xbkJjwh^
zHS?Q2+?Znf%ynma6I8ZfXZ8L)=m8*6+5#x)xqsb_(|N}1y6;==Y4W)`0k^xdsaibv
zkP>~u?VVP0kTQ!|Xo=5RcUF-->crOuh>!)CanHEU=(<7|u)ddsQF?#uyaqTgsst@b
zCbThOv-tGA4qtf!D<gt1-{M-q4l}kZ<}eLo!$oZ^STKJJ6D>J%`l>TXX(%%Dll-l7
zuO5&CIs*|T7_uym8uE5t?@5u&aZQir8OY@i?~wJCoK)T#wq*bFW19mL%&xuSPiT)D
zN6Mr8C`tF4=U7=!no8FPmA<U2$NdxCUtAl0b;Wi$4xf^@MPH5h^3bdUP+Dlxh_#MH
zpV&xC*G8vk(xS+FdO6a=!-PQ5%1KM(lQ<x#)6zEggV5Lmz~VlNmQ_wNdejRCG%I<`
zD=v69hf1FC05w58ZL1{+HDPPE5<6ZzQie}#E%o1M(yJB0U*nHJl&l&s^V)R>h8_2;
zuXDyP1r^kteR2ZOmsKe(k*EA4HPusES5lC9uYxSHlAWFXnTAFzpjv=ur5qf0|BP<!
zcCm!f^2*bV9qujc+?+3IRL?tMwzRe;uGzlNUZrBSo<^<-drn3y<_~zCl-mwQZ=QEZ
zkxp~XhN~?qf!$!&gioSG?QrXIZjGC>+LxnVK~4iN=cTn2Wm{{hArLm8{KChYk$yTU
z5}wE}Us4L^+`Cd(05b@k^pJVCuF~m^E{JgTHfxRcovmwGh%;W+VJX1#lZ8gS_<q+e
ze(8Qv>1b67#3;&LAG|oqtU8_9ULSUN>ap{Hk4c|3EZOyN4^9GhnHRUqH>~$K<V6-#
z7nX!{nj5S{u19T8n4p@FeKnJM;R->n37dkxFS0#e#kyR(4s!BTOeZ;HfUkW4zOwDa
z9Oz8&%D8A#mj@tCDab6bbF?yfy>4{mk+7Tp5xMb79s>x_dXOuC3j$SV-mPPwA%4wg
zAmMtJ>ad2{aePh({};8*Uk&m~j?VkP=k_v-Byvi~cyRd!DYr6j1E=d%@%X`PN{8te
z$UDPO)(tC%M>q1La`KG!4J;j=SuH#JUzJKs%4{a<8pKaN_dlrCZozo*Jpz;#w^L<Y
z5uLV$e*l~Jz8<E~D53eBN^9Hc&BAH}h4K$68MjV*>P7l~CW7SORUR=vuzw{=VoX{<
z*b`%NLtZI2Cye0Zj*+edTcZBtB)#7=g0@(fv#$K9m8I~Br-Zxi;dpq*cG*}xN~t|d
zB4D*xfU$r)_`3+DJ3h|X+3Jl(DhEUBN`SvQyclPv!%lXG->9Lc|8+=nSc=i<nkwdN
zTB5Riw-;vpRE_TbJ@1{7!iTWu9ER0Kpl$>N2|C+X-Ecg%y|hzDu~!eXhPxG1P67ZM
zwH#;8JQPQ56Ls=u22{*<$L%S>BA+;X#mMu8wW~Od04d{BGfoDD-yYnp5kHiF_$#?^
z4!*O(GkzSi-LZ<wvH!4cI~|x}G`ju-Ge4)=z=zT^5t#){djce&fG-G<kv=SU!OBjf
zvpP7<m*oSKnNU07Os8yHU`@;<XSz48^AeJ#J|u#U^^3g)=s<cEwN*}vrxw^?L@DI5
zCZl$JuJ7E9@Sx!L^&sPC^bft5lV7fK1pCIVedf@|11pYVXR&<^;V}v9S5+L5B<3aO
z&oMQF0U&I@O_zW&j4R|pfcHjiTP`1l`z<JVADn&{yjPaFqR(eWLM?_M|2w4!pdUzj
z$XmTr)SD2P*{YwmU@v)~Op(jp3ZM@@lq^IBbf4Nt2eXBJShV4it217<T!g$N?i7)q
z93o#`+>PpDSOxd_;bc!SNZCV7-n@irMR)A4IGUmF?1xlbrUo$?Ck8g;A3Yh_667<(
z21ycvtfDAn^84dC#W&l|QTqM@9-$VBLWoJ7Wl4bl>qr^62=8Cm{VN>6bUbhl5w<T%
zKLADiN=QpGch-@W&>im+8a|0V`uB18^c7sv3k;Oi%(y`D?@2^tap$^1r;=I6Hpl4;
zKyz{8zvR)j%hpgv8oy0uMponFRhsnMm!P9@Hckb$sAb8nuj;=#oXQ!7et>+gleU!|
zypnqA9B#-1@5xXnjL;K(fglZ4Ja|BE;&`dJaCF`5rNV;Yk_}f$JcyRSva6ZEkYEj8
z27Qwq({%{&Prp=7#|~W}Es_OZtTB9PPr>)7Cv?~ced_eVpKu)t6H~HD6H%AA%sLJ!
zPDvr0Oyz(C(ckwf4UIMSaOW;3Cmb2G)~0*fTO&mhFQiAh(gybIKQ+t8dwpI)U$JYi
zE@y>alx=@NSDxk%)Ijbw$W)K#LnyL_4EdONhlh0|n1sRnFH23?;U-}=lRwW(8Pj!w
zEURD<KQo+Lh{ieG%sby`b(l=wR8<Dq5*BZekqxtHCUt{yr&G6UB%*;|b?m)jZvDVR
zXIW|y+~mb9MMqCBU_T)i#we-@ATk=Dd*%JqfqX6NUojb{S>k3L*TA&$KYx~c^ytxD
zZf;FZ>aq+QkVp#~z(!<JMJMMIV9ae>Mzz;%7#lC_`_L5p7w?Bi!P?gaCNf45^>`zn
zS2<|=vFc}T<1X?3O~yeB)>c;V#h&X>ts0QQY`YicQ6uMJ)xi}fCUOb=<`G~Dbsz(s
z{F*#|%Y}%II>d7BS1gnr0OBV<kV`sK#r0Gq4qoP<CAT{G5`~?!&tJUe&~FXMPH&?h
zfdG8*>D<CX6&4T1pvO<AVQQKMh*Dn`oTAhlBA*||i;>LMB-EoD-{t`O0@(fI{Ja^+
ze^iMJdTd)1>Qxo>El>I!*^c|3ViN_d2cH$&j+Dzc2T@eztHU<3{-w8~=K2IP<IMQD
zSrKqiD2(lkAo>Rf@AB}d_(?fs>-hN8)6vm^D#Fx^BBL)a<#Vhh*(^aRA}9?+PQRh&
zw;K~0cWUQ4Jl^VFy~_wg6qOD7ZqsrG3d6pPUbrMRMFxb=|Bw||P<53y!=C^TT;3hE
zM>kVQ7J!nqBs(OK$2hm~VbxO3(HMlB*`{z#nwumrfbfkS)@|dz+|_$Z>erTjZ?-Xc
zKI!c2v<>NFcAvsi`ta9-;=0|$A|M(@*LVL@F9*@+j%}A)=U_u~_STvEK73|}8hibf
z=U*GW|27U*>^vD=W}W9=9>@r;5*;CtB7xG^cR(GP;vUI4e8z02x!GBZp_1o-Gh$HQ
zng!bFHx=U;wswjJ31NuL+3oxleGCV5I+CB;Hjg^UDX=a2H&y&k9xObQQBY9OzL$8<
z^xazZ3^^o@MFGLI80=79NlA&~UIMSln_EgtUPim4b~RzZfqq=RbiuN@9Uvtkc2W1^
z<!P0{D+O`QsOiTiOpVYh06~Q$^>^<`xb&-^0_|A0(k6kLOPB2@{GrF*ns<K!b4r4c
z<0qekb|b4F5!aW{#c51l8xI1TpM&K}R}{E<<;v%eA3ye?zm%2Lz<FRWHhCAl%yGTW
z>fmOg8ZMyJyd3lcny_m9sml?_U_tyBsC^hv966BUHZfe;HvqK0?cX6hsr}w-^+$FF
zjxARYH^Y4kzyfka$}RTFGHaY)Re}j+gT0N4iXtbf{&t>zG{wh<2`_WhidhP{DDHo$
zdnSKgmRW27eNVV476vQdT6D|ek%)Q|7GsA!O)fgGJp3oMe8q$B<b+hUdT7`KVY1;3
zf*j0fbSv-zbCkO4pPGB$g8XHlG~w$*rD%wF??2K``ny~UONcEekDG0vNP-6EL`>7f
z?3udhZtO7rw&8Nni~x$+d5FgtYSieAQ1Q-e6}wha3l)QiPeMHw0zFDgXFNRgN$$3S
zxwYTc7f-gfFc2*aAx~C&CsOWw#J7+vLxPk-qOSIiCw49ES2(M2)f7B=!z=@~_f4{m
znl3HbGR+Vr7ZM1O|94>YTT*lH28Z6=`ftO!Q%p^RDJ6hW;70oIIp!9l5tRlJ^V<F!
z-NPPJOqx*2>TEg!qv!+z;+xq<={Eo`BfvYi1__#t))EI9uP&a;`UGSFd^a$av&b1=
zUe2Z>DQW-W50SFe3n~#DzCy$w+g89owuUD71l|BLDK7kOM7j9>6)r5UM0C!gm<y67
zJt6M0@owJ-GcyQNbqb7$UbYc;NUMCbNYrU(i6F&kq8C6t(&x`lg&085j_>}W8&rqq
zyLuk1oTv^|3uZOT53g3k;i!|V<UtxbIv6<ZNe_y(?({V+1xNgwFnA?XwOWar)0jhq
z63LVT9D1q0lD!;MzV=8^M2~cgqn{h`T_>2_4a#p%5o+#bpkcY&dSIe^OFekcuZxU9
z^<#+90ozCkz)#vbV(^8#1_}18B<J3N+9g$J=f>8$+dECnqz;#MpamkWAZ1~AOdOAU
zq8S-c=bBcOx<^)jO(R%tyMI9vhq7>=v}6V6UQRHdqq*&fd%C{083g#a+x|m$oajff
zV!yg<9z@*3gpvZ=Sc{MBUE_xXLY@LexlxdD%P%N+GB7Zp9RFJ>Y4Sd(xWRl%>Y#{k
zd<)K;6K;%4<5}~i=JoIs4KT*_5`0MFIvc&kvXy0@p4st0Lp1&a^JPCb?I0>Zm=vKs
zaLCX)4S#iUCz87!gx)Dg|2=!${8ob;1jLHF<D8C!^eAe9@TLDAw^Z8W;bUTz6OYh@
z7g!twq2plp1&BIjd!75sWzh|t=anHv`|L5Skr9SiLB6$KBjn;^@NPkXIJs&5Cw#bK
zKlJDQFjgU|)tl!a$>t+RL4#me5C}8ueC^->-@uW(pPii=pkS}h)cqmr^*R)CH^aPs
zyMjKskZ5ASXT-4i`h8x+GESNj9NNy_ZV56IJJDp%2`bxZtQm3Zq64bf9+X9zDwmP=
zm^13EMoU?)fGs7kx1NMdK_8Ms5q+#Y@}rWxGm7@nPEKGqCI-_9zDdKYaeC)%0p>VE
zfbzQN?}nvdY{WMt09V3c5h2Uf^(-WpnHYwfof6okyHOA>p(_FHPBN;WrEd2lErpmC
zop2{G6N;wY8mP+P2sg2}k264Fo;p9B3-XHF;gr0k(emkiHJ|UV$mGJLp8ooE-vrn~
zAkG#*?J1Xb2?B;ncoTr4r%2i?amwAl2Ugbi<x4iSADBHiB*8XbSm%sga!cAYbxUHV
z^HVRy*@HAR7St-pY^0q2HB;6mA4nsrM{(SyO9U80H6(Z=;GSDlxtkHRM!P1!ZY2W&
zAOSXkn*<7DT)5q&W~T3v#lMXYUMX4O0@OnXX2=9E+g+-;>~^h_A7*U0Hj05VT@7h>
z2zc(~4Qy7$(KD~F{MA3e%#8-GaQ-u!*V5NP1@xUOB*F&oL~HllomSgkk1}!*N&o;^
zt$lX#HB?nU<>lSYX?F+$?gZP`xj%LR)w9LOMI-^hWA!zN8la}Jb3d8S>m1?L1nP>i
z1kq<Uhk;|j05L+uI32nCy0}?uAdRB=d#6gK{WJkWMUZ3ChiWPj;}ZKSWiW%NTVvBU
zT#HO#zL)UM6+t3&qq34CCdDXeGa)0ZSbx?=e33~StOd@Quauhw0zoyj!7(T%avJ(-
z#0{!oh==UUI6k|2d;W%**tQoYSrM4fI1yM5oS(1p+K7R?=j{{38WAc8oGD05{;yKh
z;J&O(qvrINdo!V=SzEV+^(c_><RLxapwZxh!KS-D&?GnA$1tj>b?l@Z#Pw%=KOOmP
z1L@W8Bb26XhfGSYV#o)xD%-iwK_edsKnDf#c0TDhI?8IOau}iR>6UwC(jMd^<n{=<
zN~k{pw0*Chf+Z|0EK$%twe>w?{^m@ZH7K0F?bW?BkPnJ=Vlq;FK^>UOuxWg5ZLPlB
zm}AiH%CL#-PZ*3xv|P`lYNvn8`8LQJx>Pd3jWmfs`L+Y?|E7wW7&bW0JRWyJ8z`_{
zVoWuvvP*%)Do{AKIBcFaQ15-j4=U%OCKUJ1(Qt*ZpbqyxlEQkVNI?-62X-J@fWj%L
zk>KR`01!-647(b?=dO+T!EAJj;b{k7w*}DC=>9Tpd{!*J)31j7^NYi!QvkhKyVcdl
z<FiRxx7kWRQtw?eK6xJ2$U(jmb8`Hemd41g>jaj3sT>~^H&lld6a08XMi14R{rU5!
zasd)`KwFE82X@rHp%7Hd8I_rJiY!%BSfj00hF;Qp|CeT^W+Vh_Ofo!mbKW<8{06Ah
zIr$IV(Dw;L$ce{V(0G}Zm6{*}4<Ms$H-PK;>FMcxJM@ZioEZ3OkV~<@zdw05)R-db
zx*})W3!A#}Z{;>G8anM2eP8_wZ7+&DPXmmj5gY*ZsqgXbTAb8s%F4<?3kAjW^z=3G
zVt26#OdZmL8Mdel(Qpoq6z*o99nZ{3>`<uZc68BkZrKh@wZ_&x!Jp6#Xml6O2gnkh
zTV5U;8;cp^xI*HBkq~uTQ@6FX%_#<Wk)5v(t;7!7NCP>A;;(=)vnlCu?r;m8cp~5y
zlX&uML*#Tt)Bk&ModHto?>-VSEiEm#{_=v$I{%c>)#%c|8lUH<Mds<bX)$=FsE-|}
z8?vkxkjpto!>b?sS{va}!m9x^<H~;)(mA)vf>OQsb<pEe1o(;ts2YpP1T^JeEhNb9
z|6VeOi0`+HGpt325+!qaN9SreP;mg_KnE&}xWVB6@7hn(<q2zGTS#C)0-%Qq<Sd|`
zj6{H(L3!y&Rw9tc-VY1>?CM%W)xluoB4c+7kJE>D?sDiMcTeSn4U#5Ehkk)<RYEhd
zCmFnm`Xl4yDbjmw>?_!yz5<;x_<G&m63|_3Mln`n)s7yUEi@`i5}@_taJbtiJBZ0c
z4sGC;0xS~hZ4p%@f>jT5oZj{RIuXiKqerK@WFy`fDFTPE@>*o1(4QanLjTbU0&A<l
zaA5=s%^#;=T`;Yaa!=9&M^}be#J9Ql>JQ{|H!hhk<LPeNEB<G0<E9JwpnStoB<<ef
zZSMuVZ$x7ls0u`J!~O}lFfTtJLi-@ck)tRd7DtWQb4#LRm7n}-w8ZygI1yoTp^}sZ
zprR1SFA3?a^F%{Hiw7xB$2;8XgqR`Oeexb;Umii>Ejg42t(d<iIMdXEn%e_z%YijP
zM^%}E4|FilK(5t+l`JUJe4U>NX}X}$=gV%I=1%I?=-A|%lTN@5<(#Kv438&KMcp0y
zd8b-*?v=pGZ?})UiwEwA3W`&0W!@>kchkoEDAEax187lH;`+^s1<@!Dy`Mo50E>Q)
zMVpd%xJw}9*_EWU0G!3)RiU^g4b*}!@66`uogKPG8*wQu9lX~9EDn)>m+n+oYcHdV
zwSeFZSoiq!V~k2jM>Qz2k7{tMWSUafRx1HH2;4kJXXoeh%FFvP_mn${ROBJJYAPA`
z<srZ#$t1wXUC%bggE?d#)`8q9|4cfK+=9^z*xOu^QZk+~C~P%GPH3~X>?Ul+O2*zh
z=)NhI=S{>G=}<}i@#^VgLfFSTgEU27L%Y_ST98XFXP7G2CZf6pYK=V1t?_gr<Y!?&
zfSv!o1u}hqJOO7C5(jblZJE>*9fT1$)9QF08Qz;IK1qp@90R-5XQtE?l_|6=HAN16
z&z75|bGYu=SuHDR`I`a&Fu6b=b>S0_ouy#^E&Gr>lxUvCIAA7r5Tt5n1k6h{-5VDi
zF5u!>uZ21(HK&6w?CokigsYhNrRlH@&;w(Er~$wQ5wvN=r#0#gqoEBu-Kx|J+uv@3
z%4-@3eIu(zxKKw8yfeSsv%V=pcV(?GsXl%~TXy~&&KltMJXn1zVNRfDO(rkE1*EJ@
z&2Dw$zg=~MRtki|9&hi%H$aN-*pEueIx!c~qc9F)!ggc`bW%D%0lP=MWmmpiE77so
z4+S9;q7tQUuPt`~6y2cA7&2|X8TC#SBEo2`<IEYm2Z*%o6?-NSug~TO#PpjtZ-R7>
zay+V6ku<Aw{ln(aa^@d<fV+%yD%rI?0nxK)kq7`VKPD#Tb6(yr%sr9ADN7vC0l`~K
zLr9n0xTcc)2e<h$Ocw~;FgSA>R>z}Wz+qcd?a+uqUT3&#@Z3<RF2SWyj*4=SV1UZi
zKGqCv1#Wljb;i{+8KNe!BEkZM8+^;_uN2F!m~$8kg##nnXketzY_sQgo&q&47zzc#
z1rdX*5loOF1TH!dLK`$WD1KY~aD4W?R{W`qi2bKB0X$59h5;<soRK&gfHW9F5_G%x
zB2_%r9k*OG-@|wCEjvJd4<A+Zkmz+v$&?uM6kpKVwBTi9U|YW?`khmV`95kG3A+tn
zL;$$j1rhCIQ}hDr%W{_Yv#p(Bqxd54k{NqHzcsd;pew`*LG+-$qo|P%%D95^SkN;D
zs;N+0ejovn#TPwoaxAGJD51PBwHFr|x#I6|?YjP1Qu!aObV89o&MEGi0xuc#@^;Ar
zlwtK$zme2qH-1-kn2@vx(r?_2I*qd0PL)Np&FkyueHBfLohxrODtQ2+<Ow;hiAs}z
z?R_#cJ5LbPUDEjzQ;HDe1<k}DgNo}_sm=EZ-%ir`xxs5&%D;5#ioDzCeO70P4-O2%
z4O`0kb$({FA%gq(J@83(12-ncKq{5Hky{%Ev)4VHL>(@4f}*%@L=Ix&dO5}dEt#z!
z$|^pNLd#RT0FZZ^@BJCjZ)M-*T*=h93aXL0A<cT|*1rmtgj@ub*fJMRgoR^gMfVib
zCvs~fr({>VzddIpW@yE&-Wd{5^Xzl+jmXnV*hJe!knqT_V(}z_z)k}d3~s!~3ZS-P
z&BjhSfnQ+YBHFkJ+!5qC8|g1yz));~yja=XS{03=uy}G2$;O+pdez6v;obljLw!Gu
z^!Mv-Vs`}i>i5ko_q71@-H7VoWEbp)`jns|X~a=Qxg6x+4GA)To&rpbpm78n1}Se2
z4Js#hR)armerNt`heYVACT#jIX(0$6J5)83q1sPhv|_vO1(D+~BAE(<sF+6PY#0}J
zair+`3tmt&$xbaA79B6?t~*lG3B|)-zW{)jasygK^_gD#Pa6KP-;Pp}rLO*aXizFo
ztYHH3H)1VYw*J*i5^7vOQ8ilhG;-|T_8-rV-5Qv(QuiRcc4?jJ)^m_bu7WDDtNZ+P
zG2>ZgUR_zq?&QfuyE*A?i~fzm&<sWEPS&AJhLLMkSq03e2h=6w6DLE2qPlxx6}4pE
zvF~r4c|8mG`uzt-<it7Y%}^;S;;6Pv9dcDK6{6RikaUgAAdmwGNDjQv6nBuTytQr(
zv*~2ipCD0P@+o?-;FqAgYHCkk_g@j5*{lNic<T1>1NNL8hy?`vU4;F2-9cbr;O)IP
zfHFef65#ad7f0`Ye`|M!{ap({fXqn6X0<*9eRi^tj~TipD5wPzoxeIe#lH#`_&k2M
zE;yrz@qw|UUFi;iE0Q`}hvw<?W&8R-cHv497##M1NiOwzT)PhFKMTPAt*wtclLS>j
zQP;r!=1hgn@D~$9`nNLtpxGEIyBo1C2D_#agl-F^%>d&BMdi65?+!{NK&sf)mlQbP
z0xS<kDvZ|XFg)3Vd??-gY>Q@(1Xx07cr@8A^@75r3edNJ>zx*Xf97sD{zHz2`d2`P
z8EOLtInGR|L%M3hBmb<Q=Q?N!UhGMKXltaW?4b()!;o%bR+X=5E8&9tf)b#Ijl~5i
zP@sNLzlH#YRjB_}6>6mh?bg{(5<Y=`jpd)W)5MlWK%deBP@91Sr7#-W+Mm--wsV$-
z%W@1oFgN8<|FrjrsDny0TVEKAo$?^MeA$^V`v`9cl;mqP=qa!z7ZQ%g{Xa7e(B^}5
zS?2fHUL@~^%A=UQw<Fm=^^+2vgjaRp4!%?ssBEzpf6@A?IrB?CbkO*$sXUychh4}W
z5TIe)*0vD9#4E&X=80!4;{j;;Nj|D}UQh<PT(~<oXk!-;Fq8N*-<fn3wBt9E30eTe
ziJXhm+df!G8$A(+A=W{K@u5>4ZTL}Y%-qp%6T%ttI#U2B|3~Q?d3PIVS1+z88vSA#
zvb#Q63@DNTh2(Jw$U560F}DWL=uuE$ho}N|762*Ur3-WQd7v1w{i`9tF|CB(%H!}I
zWPYpo!~kAU8l5JxU%SY3mA{}kTwlzmz`BrdAzXa-1L)-kok2I{Da1Aqx1N0Xh&J?8
z(lT%vkVlt;ISmXBx*qDPsYT!8^G>;d2j@Ej7Ci!J{R{^ouWiD}_)JDBzSIqGXH`Df
z(Ke{fT{sXyT%_Zd^W0n0P*hYLrPy_GaRHHfsV}Gd8JdaqBB;BATD0U<|GfpsI_Z<k
zm;Xxc-&n?nO<GyT8G84Gf-<U%B0zsXiY`7NI@bkS!`Oj-#6sCKQeiD%)qh7G1*kMw
z0N4JC8?x{Fi+!=n<lWx;(?N$PdksAD|7*Tqpz{Oo!{}Qs^>wqhRpXC|&PnZU{@hYs
z%`p>me1GMK@W~EL+30gsOsp+6B%4tJehhNxxoSCippcc6hFe`pDGb<S_M;GL`Fo7-
zC+a*ou3cPu$5GVNaI$<lIgui(p~2(B?8_`gt3mmwmp@41!Noun!>=c(6)npEOI=;;
z6Rqfpg*deE(~C+mU$(Q^buER}=li(Y%SVw1AIRSEslL%tynOD_Be9N!0}~<z7JtmJ
z`F`!7$nr{EH?_n{;KIgmv)k>+h=|6irU0F<PB%64V{6usqzEvUEzlVUs1F|?z8?hb
znE)N$2Fs+GAT2F@*I>mE$QB~Lt3=dPj?_0>Hg=e0CH?gt%wn@!sT}F4LwK*)uQ`r~
zK7BItr)KoFa<ZI^%sJ51#v&n+T58gE>GEX)sqR$q10Xfwf&U%F7;rV4(aJC)X0Hg)
z16Sg{`TLKxRBIT$Y=aNh9dNBzU^7aGZEAN0^|bjl>w@K0R}c4(y(K-udw*!L?y9@;
zMATlMpQo%cEggo_(8>1Zs4&Co50^Niql1I>X+iOcjI3;#_tEYjwdd^IuV3*xED`h9
zGo|WLFs06)zv_6q{Y_I-p~FLach2Dz{@-u(Pl$tqgM~iCP=x-$b7<3!rJ(ZMBE6sL
zLK62urmgmo%&b^aXSJx49d)uex%s_m8|riXZ%XQ)qW4fnGx!n`60d`T;?uf*eP!71
z#y3G)UVA<166gz~`pS0tMXdx*h5Yc4P)0`PB2h!2O7^$6r<z!#J;CBJC7hO|snB}+
zW^;4%aFu-;ubMh+J^T$HZL3#WL%`#c%cTup(%oH{DS6+LfbrFE{0a&yd%iq6P;4$v
zHhGIayiz=Ubf#5!R)97%qx{O?`+(-r&;vcDpk42i7Y<ZkIZXv`a0-#O5Kd*$Bi?!E
zB7%df%Q9-Fn`u}?$3N_m`WtpjXAHJlYVAme)@23xzGl{JJmCqWC(b^VoasGLu6CpA
zx7Wq)4aX3BN(y8*A4;G8)Y@CP6d6*pVfK2#i5q*WHA|v~M@8i*3lv_iR&?=0%hJ-7
zHxqw3WIq!*DNquealCc8viJ@2q~O|(h)SxHnyZIOh{fv~X{tLjM<uT!dA*g`#J2Oe
z#b|h=PXpZ^YT;2Gojov#rb$b_6>i|WC`tH<xqcg?zLRuqZ*}pcp|#<3&31<Wy-Idt
z(J3h}&$S%qA%(CP(O9E>{=xmK2VBNOzjcj%aI<O+1rp9CuMx8D=HRok{<cF^HD(L9
z!hewcwZdgF^u^Vvd6aoAm4tCU(B<;fes`;-?$3vLhLcrzVwD&1>97?mtbfo3?j&)K
zzg@a)u}?)s<>s_>r-0U)A$`LZ7S3JdeNg>0M0sQTcFzLgMuFarh{TTZBCEyN{GW#K
zV0~d9zl9tZ9!g%n_t3(Je8Sh>3FY0+Zov<jr_FhU)i89K-Zh`j4hmB%5)*4`Zg;PA
zmu{w>V%$JOx6`J&GzwA{<b-;j1P>*gJ)--?^dghV;rmkp^rWk4gO8Mn#Q>r1V|tn-
z!4MXKw#PxAi!!sb4Tslnvid6Z2()hfH0tcUZYK|n?>nQ79^okZTHoqnqiy%l-01wm
zLUdxH0#Jz!>eGVzpE_{N9L!$lQwvT~g*McsVLd=)CKNNJxzVFP?x+8KT)@UuEKTU;
z#dfY2vWMp`Q8}6hYVJ?jKUF1m@tTelY|5!0&XxIYV=5~}(j{$JnfvTJ_A@l`FFfWc
z4ebZxYU%Q#&RfC3CU-bkV7JMWs2`3+LuQ?ScwoKn^Obk}L^+tr#*T@~k)a8lnkM1+
zqx{<MRCOsC8U5a28|@<s6OoTp7oSvb8+1MT?!7O-zt&4|m*s`qT^3f>qE@p2m&@1J
zU-SB&+H6qV7mN-Tq7Iomw$uh!nB%T68&%lHaF`&M9jd2Y`#!yX^X3vgy`qK1XG>a{
zS=L{Vzlr2pOMS10--c6BX+<6Id7XUW9p4vO6m04j0?zPzOG`rH0~zk$-)^wnPg3r?
zQqxOC39fYHCYqope5AQiKtNGP;^Nhze3SKkUIF6DJS{rx?GWzZa@=9}ALxG2Ror?k
z-uF?5$eHhLut|->9m8=dgB3YwbFO+t!icKq`V5!@lHPnsl9gdFdPvs7G5-|W3fEE1
zX1DcgeT6j-PTA+@si+pFXwFIOgjrstEWL3R48q)XUv_`{risXB!ix#^`lr`fo!IH#
zX+8M_S$%;rgMV!IHMDR4PSA6ET2KB*GFExD2#Y(0rOddSV_Ybsg02bDmurVo6(^SS
ziE2P$0XWa+6}p}8Sy?sMYshD3ubbR4q!PuA=YAVSBkFpD=h%nSG$G;;UIA3g%49=Y
za4<0d9rmtQsB#q?tY^M%FeG>06AY;)!@VWy;4&wphfl}Jk>UNvb(g5VGSHB><S$Eo
zSB7>txuq2Q_~4G{s?`){^#N6x+S_10(!r;-FMOfB>zb-!EV23ESmhi-xRyrNOOceK
zAd{S8M7e+euWXIKAHgtQ-yL>?c~H%G6G0UHaVx+q=Mt4xdb)c{7xV^(*v|8Mag9q0
zAFE+Y`;8HVl){ZGFIj&p&{f_`Rj~rALBUf-5F9LbRYB1h3x^DHo%v!j^}`V5`T0v!
zXf;>HA3vL0-&2INjEVjJgL|Y>`qjZuzqI;$^Yj&dZE@1gakxv>T}d>jg-zVo50##T
z-6&FlHaK!N&2R3Mw82hrUeCo(mfGLw<JYg;YTDY@0PM*Er3GfPRFNU<H8jN6(0-p0
zKQywfp+P^IxS%bv;&rg-P7&(fX4y`JC@$t<6oEIQ(dd;stWLaSJlEds&0g<&5qyN(
zn{zqgO;Oh!Y|?mrw{C3o<WJ4n+Hq4kdi(fI9ksnD4|{#n^&_x~7nPTTjDI#qq}X-r
zF%a%`!6;PJOx_05_ktytg_hh*#0M64@r{~n*ue*y^vQE!(l=Olv+!A6GedeF$#g)w
z>DgbsOE8sc*|Wn78MLR>XRy{j@N;0c7ildQ+Jjo~_qsI?D5<t$81dMD1&uvX%LrbE
zoO@|iDAxBb^IOWD2*_om#9pQBvD9`NfCAJ|U#aA;wy)99F3w9Y5OpZK3^0hghE`V#
zf&XL546-ieKzyem;2W5rO<0-x!R>zZOeo$YCY)$ND#y<(2bl-1Yz2MW;CqG?dx8yX
ziznupiAj(2Qz3sNzH)FeDoFpS#JqH-bvy-?<GsQfx^By)P}1H<MVElFM&|!C3n(@U
z66k*Dh}XaA*R*xr2{HUrM&>Tby17`I#QBDGUSQY!WX&9`P8C<KHnyE2#82<dF9$Z3
zRDC=MgNz{e7g<$+IU_ju{bU}+QTdN$!B~slL^pi<lD^O`Ye94kTPo4>=)2=Vf6)E$
zFibyBXD2JS5-^ong}Y)zN%P<6y{D@fTfdal|K%jWPb?sq6#OIW^^QPlS>JVz@iETd
zU!d6+ol7Ef^~K+MF(pa*v6@n>c4~Lw`47qTqO{=)k7T&HIOmTIdej6zfVsyNft5et
zP<S`JQ2f08g;4G4@t<Y8fSMfXnaX)nXV@8vnVFd}f!6pfqP6GVh`Z$I>CJ7A))krG
z{?O~VWxMQ4Yz0>KDV@U8j3zd8L2Mow|3BTGX*`wd*T)Uv$P^)EDvHdJP>7r|&q@+y
zC_|Y^WRp`IGL|7z;TS@hg*dcD**c~&RK{&g!zNQjrvGoB7tjCmd0sv*pZk?J_r0(C
zy02@kb*=SVU)~{dIs7af9L}JOIgxhU8~p9-RqnfpGnT)-$xa)%^YRtd5A^j7xH%Hm
z8?#RQx*t{6&+-T2lsS{olv;@eyhr4G$o#2W@{Kn|hZUN!!F;Gc_c+n_r3%}=YxL{&
zf^PXU7pnrAzBi6N+K0exzVH0w$Y-J1yl)g)6@KvN+82ek@jNmy;JnAe!_zb_Jl1q`
ztjVeEm}%fcZTIZS#vSUq1!9qrSGQA#)}ql=Cos>VpnC7w4t#HK&-??e5flW3^XFBT
z)w!q7UMM@aE$sVYW24{dx;BrVI!Tde5GJU12gDn$ZO;+<pZnyL_%CgkS9n>O;Eqz^
zr*n+&$=hdtY1;X&IrDJJ8zlO-X73T_;4k>=F-Dh%eOD?s^P9l%qR!_J+nuW<%N#2-
zd%d2O7RZNrw7bkdzPeQ`A`;i+AaC3<SjiVMz{(IfHtXpx>CK5PIdNmrcbrJxEu2nt
zQ~EJ>=D-xAOD~%>nXaTB#w6G`DpI9#USqqi_aENG+_;~q4|Q$Au-sLM6Iv>jy(;;V
z!O`-z*5-kYO;@wsue<b)92BCBoz3d+M8rTvF`%P~zPY8PGfq^b+S|yv@clNQ!{b%1
zZ`QS4KJ-vdPte&s59Q|N<=VB2gh+j1VPPD{@2NZDnT0+heSN1VCMH&JUai~w_4Q74
zbnfpF?9*j8e|;|K$WV?iE|$LWgQ7<usGT^mhezw>oPN}pm88Jt+k*lKI-7ly-!Hk_
zCB-6|b>|57O5dxYk%cENE<Jfxe(|J0Mg-HaI*+9WQ=s_R2x+{Cqe<ff<<i6y#?yd`
zB8FcHyhxnD(Pkt!Tt^#@!gZCDXZByY+*fGE%6BLf`lkw?X?>IqJ0BTk%r7k1qeILv
zO(7hTr{3qmK@*dc&|Hb&gr@jee5bc)p!oRD#@1k!kP$}WY?#Q^>F;V^HKV(<$Dhqv
zI!;M)n5@n4aC-A_UT;a0m|ta8@&3&+CQN+vFyogUb%wIGn#Gn<!XYdBUcY|*)RUa~
zV|M5;N>n(iniYi-e&9g%>z4wwcQj%qnH|bxLJQj9UO+Y&h?3JMGEoDl1!0;yckV3D
zjb!BN#P>ARX_Git+w~1%qgzY_!xM>GPJ1rp?e4B4$!&hqQeF2>%4M<%N<7Ye>SXED
zIg{oz)(@52;wZO{IF8*5uU~y^qm;7Qsb^r2o0k`cGhDJ<bX=ShpcBdOMjC^Z-@bjz
za~}#zoYL4C*#+95!rQHW!*BH}T+#tW#U$Im{b7}DVmh53?71dOQX$-&+y9_*vsaYG
zc!S+JZg+v@8PXG3#dyUi?RFOz&U_XUcFSMiwy7jC2W>KFZ18U$;bGj*mwWk6UxLA|
zyI2Wy^z;jh$F>|9nV2k0O*K(04gabfC#=!MraWvqPG(|aS^^#>>1N5@yLY#jDPQS}
zjERYPva&NWweTr-)KV;+*`*_yEG4FX?2Fx6_CnlKdwLn_?9}PsS+q_LyOw|CD=_BS
zWi$HoXAzUp;P1VMTn0x|+1cpjzDccNkF@0n1~k!E?UNhg>msaOnLH}p%2!%@Dk&~f
zH@0rsGpheezQ;(^<Vs5CE?!<stIVvdi_`U&zuJDeR;S-}i7$v$v3p5rX!Y}+I=gPe
z)eGtAl#;r!pbu`zDg9xcmVu*|H(s<b)^Iu2FPoGl-(O!ypI!HJWlrMdf0h08*<!<N
zce_D$9-n@6%uIsAK?%hRQW6T<agq{>nb9np0bR=fXq)}9FxWqFZWrCHec7BV%D0-F
zU_-sDh?O(HjH;cD5J?0CwupLaEt<CJhCjGCK9t|uBG>lpN9l=!D}Bn0v<t3w4aAGR
zdf#Tu|1j9K!_)JdKEq#bqwn9>YJdH;veNhId-}1Eb(M4H&J1s3YwPI9d-^m9)WN<l
zSL=%tv?jgEQuF@w<ltqolb8Kb6Y$_>K$20Ar8-;FfQC>+=|3iWKT1hU=bd?C86&7#
z(_<{a>20pW&mGRE^D262>%;hcmg2G-!DP+J!WtH4Hnz*JqIjH0Ue?$YAP6=i6>?BU
z=2c6}%58B~ReoV%Vdxt#x{he<x=o?D?NL^m+qr{_OUvKiKmYv}!TMc(@rA~8z%4H>
zUyFXFrKiW7nYpKqsk?qh-$cR_=vD5fq(l&j#G3PLt{JvM?zMwU%gcFYx!Qvls`JRz
zgS|?BtEE(rNb?qc^;8=+#r-<*Leta?=%%VEE5_8^+&sUy=rGgu@Z_mewB^WOiin6v
zmhkW1{TkdNDQN>ez40{XOxkRV0%JGw^M`-XGtqfbKF9-mf<;+;<`J+zz9hlxr&q0S
z+&DJb8tVusBF~(LNdHi+r!W#^G@itbn*R=uI-o(1>H>HFz^5Zu%G@%}TnjGmi@|5Y
zR6s+l>E<TQ%*=er$*Iz9K*6aqoq<p15RxzvwYAF7cZ*6$+zAWY#>gSm@lUc;g>}g=
zgc~$!=z?SwKK`JsO&6or)z!@g*`->DR%-wLlSKFv906+7WE-P=GhpLzhK=EqzqPUw
z0BwnZp&=^`@-I^Wk23AuyEhb#l9iYQs^Aowk`~NCd)dOLzM`GI1t!jTHCDl6Nli_f
zTwGj+YjwzHAa1{BAWTa;d9V$z#<ZipPWkV9dU_-~u2VkTgCgU?@>MngmGX6NY5XoH
zXOVh}1|4-kaB~0VL+}8KHbNfiMQJ@OD=XVY+2NB(Bx=~!`@JBO)k*K3+x&OJOok|n
zE`RrD)(lL)^?#!7UPWwms2Xr{iog*F*UxN=@_!%o(t{d2;QGJ+Az<q7uY}fgH5ll_
zjW@oFNJzY<=o*j5DvxGEQ0+H6n2g12Ze``-zrR4ieNa&-h$3I*sL}v580@-XypZ~*
zcP0Eg#f|qXE5D~$u@{w>-_=VL^G5%$@7I!A>i4Y@J=)OGVF<eB@9#@(KleY`1u08U
z&^eiyJx8Lz`$-pUk$}ThK$gDVe^wulT&~rIq}?x4Jd2P}Sa|T)^5+)SZw-zZLLS`&
zH5sev;UNotkBGcHb^5>TTIWhgOr+-d^z`(k&sBzfBW`7)Q!Z7o-8C;SuN+OP(!9|F
zW&)!Ia*Z$Xx*h12x2pG(`26{^;lpG7824w0d-okVl04m>vXh^`eSLLlb+KKpGj93~
zb=)UT(BhHD@pJ;(*Wv2BDOSegRo%xYL5q@C{j)dq$W>lQ*=67hJ07a0plv!70ulWF
zxH`C|+ZY%sJVuU*8Ku8OSAHv;b4T2?-)T*OH*lJP6L-X2hH7N66vm-mz6?gL6UoO)
zDI+%DAVmQv-Fm;zTBJBZ^c}O}40S{`-Qu96q==HzefY2p!B^wwN-UK;hTHIJCvnWO
zx3@P;SG)ta<|_9+%c{hZ3gbi}K?|+^@V54#j22r}G*evIzalqE8Ewo!7jaLG`}v&Q
zT<T9`h5RZScq=f9m_g6ba0>0kimFi-D!w0CLG!f7WMI8{wae7bPbHdEUVSfHm3#DQ
zqPR&yW@Z<-UPi!g;Z6>eKcvXPKOYXfeGs=jlK#tw3Q+)0nNgacZh9cUz&o?^myWnL
z)<o;$<MRm{Z?!ULt}bF~>e5b&qwz=`zb&;@t-t=w!=F21W_<hwNK&qdyH^!5I<;Zp
zV9O56r%nEEe~*Ai?2(UKEWfhT>x2V}ZXe1CX`N2J5@-_tN2&}8EgJMpT=nNJC=Xuq
zhs7ZoDUY0fPAMwLoS06zs3e~}$->H-2V50l`K=DrKfZtgo12^e#=)Vin%-+h>#d17
z1V!#PmL)o{PbG>O1t8SHIw=oXn^oSU=%T5wza35g8nxsQRoJ*eQCoVfvSRD%C~S)~
z!B-LiK|#G$&L#k4Sv;qVTqoxv^w`G8n4a3PaOv`8PCmZKl$5ipm1ARe5cjmbel3sJ
ze0;V{Nd3n0!K~T=H+(%URVG<`tP|L3+pNUudDK7NHU#><D><*TwYh;2j%RJHiz{Vx
z;wyKywA_YksY2Jgh^S~Cibr*zc?E!F))Dv7uK@<r#+8AA0iVTk)f%Qf5NT;@aEXqE
zJ*0SHXK0+HFEpgS&D9w?S~OYckum-8qeHc|p1!^Vy5J!pnC~%=Gqbh^X4Zo8ps&K_
z%XJy7O~hH2A11%d$yRp~f`i*<&K6s}fWFP|X~}jZT~KMVw6atxdq{ueFobllr7NFu
z<2&Le%0V7OSS_lgRMYb<!^Y;q1)&IbK?mS^JJOFtBdS7rJ}N6qI9_A(*}#H;Wk|?z
zRQvMJ*LtWycqgFZ6$+D&Rnw*ZbJ)>etssZa;TGmsp0F}6DYvpQMBdcK_aPIO*j~JN
z5xQ{VI80jRjV}^vYQsy}XRZyoK^iqRGjkaF!^M`@=+~J6HP>&A&V0mwZH>HT?CL6U
zn~_anA3ljkMrW1=as4zApxt<M-=u-T&|7j|If&#JFgiM#9Hl+59)7_1@%_})HdGA*
zXRF`)S(_ah4VYwuWAIBVysxIN{<^)L>g`|==)9IcnOT<F6w#T-V)??JwSY-bbiP&=
zmT<4A=xtbN3x>S=r!@`V<Qwe5lszW6!&$<%B&TqM&d+<H{HAj4OVZF>QK?7924EIe
z(!BV37SD?zB;@Ai_FA0Ms~vGN-SFo>PEdg<ks+=$6c$GI^qA-5<X90yf>$FF{x0&!
zNY*np7K3SJ-YGVZ&oU^H<_ioQbAj>09MA^{BoK*}aEk6LyGWk$i<3d^1MN(Ssj0e#
zhA{))0)yb#9R(vh{HJ^<8^6-sgaijC#Ry(5u3NRW*Ta9VVwoo1yJznpNlLkQ?*z=a
zq<}jM3yULSb?jYZ$-V54umK;SnA$Adg}t+~x?1;LxTl|)g(Vk(eO`Wk@;3$p6B81d
zd~s&_tRslvc)O%-A}DW|vIwX;5nVYrD3E?Ekc1CAcL~5#6|X%yEh;U|DkmrBKR2v*
z`ZPBUsL-M!BF$5m;@qL@q=sN1KqkY{^iovVJrU7i|EF||i;2-4!ZD^w<Y+iT59~=5
z!pb2O17++89!j30!??c$;}tkgNX6N3@#_$~kXjqcfQ~*2Gn1J!ODg3^7<x1;Xt%}-
zcgqpin^CFMba2S|QE4fPSOk1*PGRA0z(_z&XY!rty8I^lgh$Xn)<$7rVQ+FyMK9{u
z+3m%G0uw@OXja5H>$njPnDk8Sc{jH<FF8A>Xh5LS88>4}AQ1Qt+0j&>Fo9Ws1E_%W
z1&PGoUoj#X<{7T)b)i}ACv?k%ljVx6U0hu7iEiG!xxTRhI$t=fzBP;r)j{tRGu$g9
z6N_ScRltG+q=r}UjCY_>(9zMQDqP-<Jvl+%+Ek9_m2p^^A9vF6J;s~yr<#`v*o+;i
z^0z^MADfvWh0@Xopg(PXd0A04y#>u})K!JWNb5z;l%b!Yp`neY{*Qp*pXke3&uw-&
zW!?x=wS<PNdEg$>1z*_%VUS7Xb2h{5|M|*c97M|;9CO7E_!`2OfC3SV6=qs?L0vme
zpeR?nKL9aD!1w;MQqodu8&nRYUeN=C0oymuT@PG+xJcK3P|g9y{epV1|Ni|lny1t@
zmV4oD19aj9>=9>BO*_ogB9{G^;8_J~K8JdZd{SfiH}<OFrzUxE{eSx-V~d09)0#=E
SJ^um?F1=I6Crh>LZv6+iNUXX5

literal 0
HcmV?d00001

diff --git a/cutoff_analysis.py b/cutoff_analysis.py
new file mode 100755
index 0000000..fe65ead
--- /dev/null
+++ b/cutoff_analysis.py
@@ -0,0 +1,60 @@
+#!/usr/bin/env python
+
+# -*- coding: latin-1 -*-
+'''Take an TSV file make a plot graph'''
+
+EPILOG = '''
+For more details:
+        %(prog)s --help
+'''
+
+
+import pandas as pd
+import matplotlib.pyplot as plt
+import numpy as np
+import argparse
+import seaborn as sns
+
+def get_args():
+    parser = argparse.ArgumentParser(
+            description=__doc__, epilog=EPILOG,
+            formatter_class=argparse.RawDescriptionHelpFormatter,
+    )
+    parser.add_argument('-r', '--rpkm',
+    	help="The file with RPKM values.",
+        required = True)
+    parser.add_argument('-c','--color',
+    	help="The hex color to make the graph",
+        required = True)
+    parser.add_argument('-f','--factor',
+    	help="Factor that is being analyzed.",
+        required = True)
+    parser.add_argument('-l','--limit',
+    	help="The RPKM limit to plot",
+        type=int,
+        required = True)
+    args = parser.parse_args()
+    return args
+
+def main():
+    sns.set_style("white")
+    sns.set_style("ticks")
+
+    args = get_args()
+    rpkm_file = pd.read_csv(args.rpkm, sep='\t')
+
+    locations = np.array(rpkm_file['ES_D0'])
+    locations = np.append(locations, np.array(rpkm_file['ES_D2']))
+    locations = np.append(locations, np.array(rpkm_file['ES_D5']))
+    locations = np.append(locations, np.array(rpkm_file['ES_D7']))
+    locations = np.append(locations, np.array(rpkm_file['ES_D10']))
+
+    sns.kdeplot(np.log2(locations[locations !=0 ] + 0.00001),color=args.color)
+    sns.despine()
+    plt.axvline(np.log2(args.limit), color='black', linestyle='dashed', linewidth=1)
+    plt.savefig(args.factor + '_distribution.png')
+    plt.clf()
+
+
+if __name__ == '__main__':
+    main()
diff --git a/h3k27ac_processing.sh b/h3k27ac_processing.sh
index 589f794..8d8e5b8 100644
--- a/h3k27ac_processing.sh
+++ b/h3k27ac_processing.sh
@@ -42,3 +42,9 @@ bedtools intersect -a /Volumes/project/GCRB/Lee_Lab/s163035/Matrix_analysis_PMIT
 
 # Get RPKM
 ./rpkm.py --peaks universe_enhancer_H3K27ac.bed --experiments h3k27ac_list.csv -f H3K27ac
+
+./rpkm.py --peaks universe_enhancer_H3K27ac.bed --experiments h3k27ac_list.csv -f H3K27ac_all --minimum 0
+
+
+# Graph cutoffs
+./cutoff_analysis.py --rpkm H3K27ac_all_filtered_peaks.tsv --color '#008D14' --factor H3K27ac -l 1
diff --git a/h3k4me3_processing.sh b/h3k4me3_processing.sh
index 8d9ee3d..181c4e1 100644
--- a/h3k4me3_processing.sh
+++ b/h3k4me3_processing.sh
@@ -35,3 +35,7 @@ rm /project/GCRB/Lee_Lab/s163035/Matrix_analysis_PMIT_25842977/ChIP-seq/universe
 
 # Get RPKM
 ./rpkm.py --peaks /Volumes/project/GCRB/Lee_Lab/s163035/Matrix_analysis_PMIT_25842977/ChIP-seq/universe_h3k4me3/universe_peaks.merge.bed --experiments h3k4me3_list.csv -f H3K4me3
+./rpkm.py --peaks /Volumes/project/GCRB/Lee_Lab/s163035/Matrix_analysis_PMIT_25842977/ChIP-seq/universe_h3k4me3/universe_peaks.merge.bed --experiments h3k4me3_list.csv -f H3K4me3_all --minimum 0
+
+# Graph cutoffs
+./cutoff_analysis.py --rpkm H3K4me3_all_filtered_peaks.tsv --color '#FF7C21' --factor H3K4me3 -l 1
-- 
GitLab