From 28e6573155b97cd8e69fcba7e28a075dfb3a39dc Mon Sep 17 00:00:00 2001
From: Venkat Malladi <venkat.malladi@utsouthwestern.edu>
Date: Wed, 18 Mar 2020 07:39:21 -0500
Subject: [PATCH] Update percentages.

---
 intervene_test/Enhancer_percentages.png       | Bin 17440 -> 11299 bytes
 .../Enhancer_predications_updated.ipynb       | 175 ++++++++++++++++++
 2 files changed, 175 insertions(+)
 create mode 100644 intervene_test/Enhancer_predications_updated.ipynb

diff --git a/intervene_test/Enhancer_percentages.png b/intervene_test/Enhancer_percentages.png
index a90fdab321f3941a01141c6c713e2377a307f3a0..7dc7d26f9d13cc052ec34e5a7f91a2d896122f3a 100644
GIT binary patch
literal 11299
zcmdUV30RHm-|mu4gzRld(qJe}DnqD9NM%W-dC-7JljbxJrfgENL_`{uW}4^Hh@@3C
zNP{FzvNS2xs&hZ~w)ffp|Mh+6`p$K}bFQ<m_O;Y{-{*aw=lT8a-+kY|=gJXPr6p`!
zY!nJ*3H9Is4GM*YPNB>(Uoank@!t4yH9mfE+)ve7fd4!foctZXFSI+T??|C6T1Nh#
zQ{U!2jXz2|Dd;(A+L}4J9CtWHIdj~}&f3<=+Tz68v!@&!Eo^PpZWZ4uF1cx~xs#Kf
zjD*DB7l_+BoR%0okdsWItff#7?9+0+KiKA>b;~(wcBq!)O?i$*E6aw3lE2*gRX*qx
zm%NPG^^m8^JsDw=m9eq@H@gj-MedoEetJ{)^7+{=jr$Geiu2p6tyYcHi`NMUU5Vv+
zDz<Kw%8e~wzccF`X8ocqg5G<-A7UC*XEt{*cXu$|3_C9I@fCd@+`PE#S4_h<{9!=D
z9P*FGi!4j<d1voAvJ{6S_Rht}#v>8r%jH5?B#Q6C73TOT`GrjqAIDhcap7b0x&ODf
z`Ci)d{o``xyJF_O^A|2O#A^yJ;oQ78-0Rz2j!k+IM{RN|`dul$BVV&+C=|o^Kt8#Y
z&a%q`V@)Qys;Ysi;d_#U<Yx2+-(_dZa8cf~<jh?{;dq<xlFDz<nqxmaGNO5zoyVz)
z$K}zPcUhsk-O{qG-$o^-78$Htw~iKbOH{W%_MVN#kt5ghocf}A9(bhGc`jNdu4QeV
zB;)cWNp8AxY2M7_K&DNb47L95w$s5OA#rPLa`uJoa?!!#6A}_U-3uwlR$*NqAI!FC
z``FVX82389AyGFmFVBD$eCN)iT*saSpV=81`wtf_EiL(}JiE@mDZR8vTT?TdJf2MC
z`0+@7i}5dQlY_QCiP`P$%Dh2JRe?O3`ucarGtFurB$t-%_nPRuEaN%xc)z^-*!v4h
zCugGOl+4Gj9=j0{vDRz&xuAm-ci!yG)JhJH*rcRQw?aa~$0zn(m>++-_JQT=*Z}vj
z{)hHmm3?JP<%0QFu2|8VlK-V%EW0)8>|jgZ_nz>Ovd*s$cDZ~pWDW)SXdOO$ttHzw
z6$_#|Z{ECw)YS3knJm6x2&R(#SFT(!czo>KA!X&U?hupyubpM?qjgHP9-p($ypweO
zdhF7rOTto8dX7D>8tS5z1gMvn^67+#8x>T||K*oQ_-lJrhxb&n;x*35>(_t1bn&9l
zmM!Yu-rjQhjdLh+0SHF1@_?P~j}|W9*jp`8Z0$5Yku=m*B`PYajIFGHVp=8XHhglm
z-#Y5mt5-EFEaJzTHRDv{)6(jva+E{F4GP^evCaiP{r&wbczC?NODv{X`(ad$x9O*Y
zjJ>~~vK^9j>b=e(W)S1eoD^Qae*eRV4`sY2GsY(v#k0m?GS2ndx<(HXQ3DL0+5GWb
zOv38>aQEcsYc=P|Q487eR>yV+G3rf0wMI=f*@w2pOdXn*x_V?()OtQ6VM_BNq}8Tn
zoSY=$gk@!ou;E(8-pt_p_me_Qy!Ep?yz^-^+KSby<E>J16=|a(#nTG4yGO3E%$++n
zJyI(wC|I>iS3`qFHXgGT5fM40YvPqzT2@wQdu9CDh9Ds<g-fgbW!%t^Y9fcn)X)az
z`!>EhK0P4DM)4K9FJ-?U@hu`OtlpGjlpHvhl2yus2))u>9g>)pr6Z{J;7wPMT*?Dk
zkN(Ml)NC0JN;5m&+B)jdqk~Px9?2<%-wbXA2b&LlsOVvyIzZX``tc=#CY#Ug`Uj5q
zB+fMyPRaem<$afxLdjrbUCKM^T`&Af;oP?y1F6NvHrem`IdR{rTeokUIbY-m5aY&g
z|0Ng_XmjcCVq3aY*TEspB`G1{(A~Rt4;?&s#LO%<Au&;a>UX@i^D>`~r)L3)aEYRc
z(u8I_<1uBPLH)&r%QMVsLKRKj+}vI~dv=6My|Q)dR^+fCm4wNmLx;39G}ebyKV;=~
zExG!P+_nI>m7JRk7>ZLgRZviHkO`ylBcqXljETbWcG1gakqDdk*jPp0psJp`0M}zn
z*hbZaBqVgTw1WACw{F$w?(U|=q@`_j?yIlxFhnZ7c=2M>^92~bV<cN<uhn}ysgHJ2
zhr@duaz8rc@>5r?Tv@|2z#S;KUSi5B<7}m}Qc@s2>WQCt;C83uNs<1hON=<JdOzsN
zH%9NBz4-g@3*2_|;?<sa4GHbuZAe3uC(B=EcMB2~6VrBb8g=V=R#T&GW0S%!EG9OX
zc6ufuqBE^nF>##MFvnm=BfYkSD{t+O4;BO2EW7E8fQ#;)9%UNP@>W>bAdg9oUF702
zdSmj@mE7E)+xkfXnd)NChyi|BypC0oaUXqj;=~C_uSqABdplI5rKLl6+CDg8YMOcS
z`7&z#<Ktx+mJK2mRaK&E*X|wZsm+#>#R#5PSi~&qq7-^~Opp!IrtFvH9Y4LyZM}i-
z==-1i8DIZxwP;2M5lWF7NIXA3aYKfuggy_Y^<N65#EHfv{X?p%IwbjC-4~_;-Zc<s
z%a$z*zI!*mzu(edlLh5dSBSpqD}lb6ogHt!4Zh2C>`_hDd$ND(+n}bg@dJby!Yb<8
zwZ#YoeV&zg?4JGSihO6khtI0JyXUd6uzX}NQs3F+s{pUG)zxp@SB>L;-shC?_%Z!M
z3@-n=l5KPzzo6g%l2qr^sfU(UR+(5)9vP>D7!t5Y<5y5rOeiR5{qyqn$}cR-mM<Us
zcvG!EaMzb3C?2<HwDodw#fnu@-0afafg=wXe@9!0tOTASF85`xNqPB|<zs;Vf#{$;
z8s_GY@|*`^|A;Uut*T4i)pf$Hwx*`$*1N&m?*?`B^zKNQ0tnu3jl``;OhM{o$1KxC
z%qUN*t13uWN9Q(zNPznN`}Y^GUS(}vkH>iaF2=oTi{%?3Z5^EiQiQji4&1+gzvP(~
z$?lII&D#syt_B1&&;*o|0vF<PJ~DC@pO5D&y`EQLnchtIMfLb|iKTTXKK{39&HI}R
z=1}(fZ*iIA!1lYQ8vNY;m;VKx<n`Muvtu?nCWc=?P|yne0toc_++5b+2M<yvCR_tn
zZ&Qp@3(u0-h!^#|$l^8U)TvYMV-339XY)y^Oo)u!=f6s#@uu3I_$pqH!z5m?<ezHp
zop#p3EfWjoO|I`B7;x@?Bizu`M37otUA^V;@e7JUykb;9k;tY^8Bzjd^h+4MuCC5~
zdZ@~|&ACY!MdInxr~R+vg0Qw98M9q+5mn{okCKvJ{vA<Ty1Iz~cvg09Iw+L$@JDL0
ze)^}=;fRdj+qYvgG9<ONv<xf%iom3!*{VkQXSK2GEm>z0Z){Mq9Qjm}5FfAZ<yFWZ
zpOK-3H|UX-+qqNMCi^IUjvqf>bd4PVUri?8WaJW-Wz&|>(qgdMFLG?4DN!kimut7{
zV>UK6{+}TtCH!+XeSLjgvhEognwkyir`5uw9TjtPb3?b7`Wd@*uLPJRrln~t7gWuV
z5g~V!cl1-l+D7|cZ}*xURW!A?w{I!(%2&98?~Vz4`9DZ~fg5{AEcBi{?<WABRfRhG
zzhuxRB~2zi$-%)Z7~vpiHTR`8wY9ZhynLy(fB*Td+qP*N8J+oS!&&lN7%3k<oLISP
z)xb!tobKMezJwYT6zo2A^5mnKm;<X<uYOToo&LwPFS0PVC7To*?`*ZZ0Kmn5p*RWu
z6&gHbX)fU!G2;TguF8O${NRD2h>ag~V|ds>G#X%Da@qQ5L}FrMrd3ml+1&vwqCuAR
zUGK3(ldzHdWu4c{I;9wH;vvcw78W3RnkY5tCvlr^rnPni>k_UQWO^bBaeRC{A)MP!
zy$Ie|o9{km`7kzC1aykF;<w+{1yyaJ4!iC~;g(~#VIv*#+wjA-Xpxlz0)(Pht21x1
z`H3k)p}zj1vNNM_-11%48G~#aqnkZO8Rj=!ru_JkJ;HKwCbZ!2@YK(rKW{81y2gO@
z=*r&7Nq5vhG5xdyWtW%Q-j-H-usa^JvHbY#M*9<#D2u_COppXR&<r?}Vdf@=J2O17
zm*%Rqd!~~77O#$ivQY4v`qr?s!!!Eg#fuwj-aau^0!HoN8feW;Z%j4mh;hTL=`wF)
z_wJOG?4?J_2YB0W78h4#6iuaEipPk@c4oaT?d<7^Q}E}|=OF}X$LV{!=rqY=B^23F
zx;)o5vuiSEKlk#esi<UbCWYX0xz4=2oE$yc=*w-ju2tOZmv>k<AJgW?SIyUw*+we^
z)+7hw7s|c#8ea3hs-WfvQua~a->Xdy(atj3lifx>>E9>X_w0#87TM#azbsk3^#BO=
zZGPv0MwM-+Uu8)Be9UBTUG$6R&(-hkkTbqrTwFXoJ-tbF9&X+BMEl&HsYpQ-U|NRr
zfX+Nt)}ueiG+FQ2eM>j^)~$!ksd1|k-#H179v#`Y?_9$Z+<0t5mIxiYKKZq5=@@iN
z_2$LSjCLbndL|eog<@T<b5!cn%e&D?Qz2nt>htH%iRzgd4k^~-S%eFA;H8<L?%9g)
z^OsLaNlI$j7JZMs^do0?V(2AHmTX-dUDx$GPVML&l5J;A6J_b3HVt>Sn2z4-m<Cy{
zTw`mFt>B15A@&VxPS{&}XN$aT!H}P^57U`Ss`W;d(3AEABbu)yvpL<{ab~n+?$X|K
zth|q9$ii8_(D8`aY1^Lk^5tQ&W<ZYqmdwNtA3lUio>ox_0MsV0u`N)C%4tbE<qrx#
zYX8NHznn5NOUTR1yPJYRC}LsQER;?J{|X|$CO@awegw+DJsrN?;*}Dw+X3?0;MW=8
z1_B#5ibzPPA(06qS5)*H9v%+uBQHa*e`~o&$<FRE+2I}UJW;n_EAe^7%v;29kk*oE
z6$fz^56H$WX#CP)s!O;e4{kGiX`a=Djoe#(3NZk|r4lBiy+J8Z3-zS13;Poy*@2?8
z0z%MacCx7@*D-E#v@cpe-7LPKzz`LM7W3@cK>))}WpBJgh9{ODxp$L4x6Ouer>v9i
zw`7gDMG<n$Yp7yXsH{g;{p;6b&>i>FjmkRho^{&YANLHY7Qhb`;gJ)w7^9f$0~~hR
zv^FP=RFszsZQQ5`Jg3D#JBecIzy0<bz<>~!pPT}u;wME#5h9?(h1Z~EICLNKLd8U3
z9SM{jdATh_QupM^7$VM7i)Rel@?DxH#|9drm6ne6*2#lrgzh*K`Uh`qGaMr6TKYC&
zFW^8wGc&U#$37a-5Oe-KOPl9NjY_EFx6g0$plO{T<0>jDEM}&^%Q$|z=HTEUY1_WL
zqO#KRIZdF(i6n`6ALb*+9MuGpC>Oeq)u-3)N$i-N8CUVv4GazrR`f`yH7NG>y7`Ih
z#K6hK@pE73kr()dL3xSb#e#MZx&BPJ2A*D2+(d{5NL6i1%LF7oiBsoIcvSQ6h$Wwy
z9-l)r#Q?5{<ieB`RaEo8gIrne%Y{9zUn4*c3Qfn`JyQQj8bK`C>=55<B$8EBKtKU1
zn&v$nvHXX<6Y8JLmho<|@k?92&+yGT1i)j{8uIN3uY7i6gH9D5?Hlm+1N*3vOPc^>
z4N<?c?xS?D+9a@Ak8l%Ig-d<7?3iGd5FIv4PQ!2b&F<1;g>yxNS?`hvQu&HrO#Zf4
z0t9g3!i8+?>^ksLAWY#f{G6TGm%ZU*<;&H~A)*f*KAiD04P7EOtY0rwi!pMNtesh9
z#KcrgE3ef-@3r+BNht<$TO9wH#-lI*R{r{H#L_(zCA&V^Nx}EP7VW_FiM)^j7Er8T
zA+O<~B!YK?ohj3hIuBPFyFUFusMX`;s2Wp@ViO3&D3BRZYJWpw2#17mYm77k%YE$g
zTTNkMVU)X@{HWE!+qP+9I;3PfJ3HS@2ed~xKh{2%lAPRt<!Lp4*-_k=ueHOb^$|)p
z@`kGMuf%`*_w9HtT?P<k;#jl2rP!w!o6a8}5wRC3STiH@4{=#)=0m)@jKf)z0;ST9
zpM=1vJidKacXG<eJ3E*FPbiV<D8tR(W?#hvqi3?l);qej)euSS;qQh70YkXAVGl+y
zn1hMt<VYohyAhcW&uRrHXY)sCIhbq;N=gIyj3P~-q)ow4rDE+aZgvI_W45~6;d!?M
zl?JOJ^gfCH#`Fe`z;u+Z=1txDw5aH40cvMwC)bt}7qe_TQvNf2pqKuUL|Ue`e#E`%
z=Qjn%nByHD#WQfMi2w0VrFrJmVLWT%dUBSXWH!i20OywdG%XX8RP44DK)(N7cIrQK
zlS2@1=bPFbMFMX0e62~pl3NuL8A_q&n#-3~{eh;`;uvXy19@+8kx9|CV7Q|ccZeZz
zL49YdIp(<|>2;Lw+>({09P=hkk_~e;OuWW!!7eH=B{s~@)N)@Hi1!~D%<gFX{5Jg*
zr(e5U50BT>V75*bsXYPMSf9zhI0J}f;>n=^SXSQ_Yk~Fncd>+GyUsuFmtT5el7ZhO
zkIz-ZlYKuOLtO)N(M7$9RSipc^F~XWn_VST@_D}=zJL51_E7MF>wb7pO^?*0B|RtX
z8X6k<5)xC={w}nW;|vBf4-5uLXCq0dk>o9;8lSYo0T3q{=l)=9-@3BR{M4!-KHW2C
zk}Lj7Wn_90F<0yBHB|0yt=Pp2vr3NHr^c}nf$rOe<=qvgigHsp^(CKcyJc{LJw`uf
zn$<El?8c2qVa-?}<SN+VB$MHjYFZVj|5tUFwf}Hn-u(FlM-9?WUO9jM{OEwmY|GDS
zS$9fV&CW30$J<;AgJALMKa7dd`&sID(Xm`~!EnzHUIlI=W{kMyB(H0p>bOK;e-k|z
zg-!VkX|M*szRH7ah$+0n)xf~Ut1CsF-j~b?;Mt{f>((uzfo8sC7so@osCW}bj0%@G
zmuJO_7`V2UFvL(;L$EfHB_$;rYj;^SZbE@}=R<<rbpQOMD#Ny8&m})UO9)dH?;lCY
z4xBmNQ7SGZba3CkeLsrYPK#GqJQ6!MH>%M1p>LL~+2Pvb{6nHy*IVJK5o3=lI09uo
zzV*UxCbngUQ=bN70<o>FtJ@Tl%jY$!UN~B>JUKJP5TN44*%mImQ9K*u<Gf|d77<ZV
ztDmDNoPz_Mnt7svt>MiZ0;26hB_(?czkj+HvYUaqui)iP+2b{74reB|)+N`_A=i)=
zgIcTnip2cwGLCsNyLKIY`}Qs2Gl-fF>k*sw0(}JQ>jdpYh{8A!`hgotyy^%?RG4m<
z6PO&pB7A%rAk+Qn(;-@YeSN6dDoiVhaIzZx2LgQ+t4y?Dzq0f0epKPVYT(;D+oq$o
zg%6)zU`+(b{bcQCUPfdK46%kqTTLq$7Y@G=y&4px3tMBO>V0!(LVYKG7t6AnD2>>6
zLDHPQma1WBcpo@pb7n2R(<{C05&>sDAyWRbB|TzaK?Si8Lw|<)jT<-Gt5OSwueN(I
zWU&pSqfTvJ<82(VB$C4r$%G4l6S;r5BIwkhtt-sWe*lg)am2^xYGWiMqrVt8JnOu^
z*n1??#0+zefJncQoBgf^Igs|ln5ucfG0%YZC2v5y5RwPLc!!8Io3^|@Vm?A3H`HM)
z5M&`bS=aYZFKl<}J%Tz!@B=FTMr9yRBC7OJo~5{M#|eK~x|$p_gP<uGfyJ}y>xg_M
z?lf%I-uSo&^&cwylLH9}zK3ed&d%NlnWTO6XyVteb`013`&qh4qYJNfREo@;cq>Ik
z#y+!C9mMXXP}ZbZ@bnUYxp=CLkrYTGLDU!UK)x&5|1ZWDzWt9nM3XE4MyHXQorx}t
zb{M9F2Z>T;?B4r=6t3jsOWk5x*$V~vCo&(2kca>z5zpf6m)3X`N3c7^quDkM|C#B*
zg36duG-*k#Un(~p!RI}(587nS4r@9ji~v<NP7~XsQ~3+#uv9%lA_|$0T1`A*Qpj>T
z=c9@6tYfxHZC}upE9z*PR0v~6??qSi`hP_O$lm@5>UyTv)Y-we={={z2mWby1Y`2V
zaI0OhY*`e7*9sf>Zytx>KRphJnLj)ZShu+L20yw6hoGjlPn?LFzmzu#jFYP2?Ol|X
zo}PSdAs)at4C7b%(~#>&Cz%F1JputHBq*pvYCsULtmN6xda}&X*FQoK2$~oT<SS8v
z%}7twl^AS((sh3bHT_ISQQP?g80|6FQ7L@<M|<|)*VS8P*Gpfwgm<?dc;XS_lGcPm
z1N@sydXv09aY^3;N=VNm--Xc#J9+!*SKNc`g-8F@j{n>KVu^~13Jwhwp|UMrtcAsd
zbFTGIk7Mf(bM1F{96Y;SkL)`8`3Njk&>!rC0QKs%Yg%q@xj!o1dr>0ym#E&Cid2tQ
zw16&!*<d;PKN}qII7sKl{;eZk`2BDf`>4V+v@I!lUvhhR1GuF%{>{vQm+`lmF`{pV
zZ7|%kqIc2e#qyB7<eg&=-~GT-`{D%&0DfX%0D5&4Nc09hxzxq71KP8gIcVGeHa$*5
zLxZ2{Gx<9F-`&tnsze}~zzr=ymu`L8=2ug#kUn1b=jTjw=~$=9uoy@;{^~0_=l@kG
z*mro13#0!4mEZm(AwdNCSY+M0eQ3`lq@-kRUWd1b<?r>|@{Nb-q1&yv0$^PWZ`dF&
zARs_`G2qJzgpj<qLYi_)*+qb{yfIIL*hRk6kS4Gl%@s(WTF;loYvI4{@D+0Fj~yNa
zcyO~<xrICaRr=UWDoF&MD&^@_{`FtXDEoixA9`EZzl|$;m9N9|n4b3hnYLo8*sveV
zU!G&%^`_Xz$BQ}EboI)WI&c%RTggWGI*<r(tit?F1qX4mz=@m%h*h~~=KO=!F7Swu
zIPtUJQ5bX+!$k<ZXD2iIzutMDTt;9R>RA_kGPn;(KWYTA+Ryv>&6k~ce}VMFh`|IW
zi*%y@;J8Zx=T9L83C;t5;IOHJ?U9bf4~qMKF4pacj!vw&QJyQ`pLTO8S<$a6(HU&M
zTis2=4?kd62~>?kaJca5RW@Kbys~6Ys|6}gC--S)mTwQ<F)&<Q3FwuW_0KqgF@OED
zqs13;`?iRT%;fzp2S>+S_wEf;n~SL$7=*UkWH)S~Uwz_7Uy;NoxzcJuPZ>v`GCkYY
zgj7$o??^ir;*RKzL{rp)64so#UJIw0%&``$R8nm%&CThc330IU|GYqm9%A17_&C=N
z>pKlw3<^ZO7A*=5+<r0k7&9VXPhN`e@^MbT1`^c1jpt?mY=ZWRwl?URwDk_1pSK8h
z-%={=6NoBY9}$T0>{bW=@z)$X&$Prq(tFgpK`anbz#Ivm9SxKP(Db})TajA$E%ZZI
zSGS#%(`av<*}YXYak)+P=rKn_w~I<k>w_wXdrx$-QtQKKr%!_~3t?>J968#ct$iyf
zW%z}IfTz!~wiZnwUuDpbof_+Vtpv8w1kN0X79DBFqhjU_RjwJE=nC?v-LhC7mG&_n
zjFSoy#OC2Jo?{@xCh1)df;txDGc6B)gm5NFxVw&lZk?A_Wq^~T0qSAdQB=tPeDEv|
zIOOT4FUDDofyQJ_YisNDSBpAwf-{|JWyj5HAIMHt$<0K9-1MdTFcbML`rpKhh>Po>
z7#n*IDYTt^cLsmCnJUqDwW*~@?(m(>`(Scy-7YbZ<-B#VJX%ZH^Vzu{!5~x&_nFBh
zn+Z-ztxwKwPXMbxO_UeaNur~7B0x>lP1Qkb7wiP>RO0sOxve(u6rxX1LD$<;!2pFX
zm(}%tCzS}-HV$`JH$=*_QZ?k|FA(_^V(cDY_k`fo+@-rC(A^b=$&0tCClOL!en92U
z=BFO5tbCpkGt4Q`?c24BCmT|oK9nz{)~oDni$SA%X4p_jJYQK~xy*0L%^waraZx@q
zs~9F?1{nw8#E<MLf|zQEQ{w~vxyh9uq@mXs4dFNTrR^x3!5S|&-5Al%skxXxVNb;T
zgZ%WUg3CXa-@#0&vxho5a_G?2@tYU@uU>V6WB#Jk?h$&agn(3aNuQo^X?IgXWc0V?
zleH7-n@t@o@|sG}O*MA!3fRfd)ceHe-n))iw#5?_5|Y-<o8i(JrYK*QSe+A<Y4N(U
z7k7@U;jG3rOqM=1ZIxPRR$fuD8iJ&<z0hMx+2(5jYj$$sXn`Vd#Qt;h6U7m9H(tBK
zGhU1@19^72SB0s@Z83Yv+4b`7Hboq4`Hhp){BuiYwYj<ZX7qL-(#~?CR###pEpbkv
z6vvpJRab9|RSR!7pJ{H^i;jx=kTz`;VG-&f7Dmonbm4%`);*KG(TZI#G$f(5175s1
z#Atq-YGq{=fCTOWORPk|oGvdf=j`wtQHeISvp+jI^WI4&0IzHO&1$#5vsD=U&(%Q~
z$1VIFd_nu3Fc&Vr<D4*i6ye8p4!M#E>|zwWo<N2aR<(Oi4OX5G_d5HE8%wGPorC;W
zN8^Um?>MZJHI;lbUs=mGFOE)0dDRQWIfVV@B#nU~f0OCbN2#ehL0zuBe}Aki%%$C;
z$p6wQAkd`&ScKzMS4nQoeD7H!I}2aL5`y?go6EIx9>W#uiq<=Jzt9vmVk2J2@zv~b
zZQejiKF$ohL1#?iz=219p6S6c5~r{4_nm3W(<Thdch1~k8Z8P{jMx+*J*9P@LmkDK
zi5-x`T&S=%g=3A1UDY8HO}LVa$vF-1(q^<5l0zh<w~WYa$#YI5d6#08e`uN9j-8R2
zK0dCq%<;mhz!0&cvDBLzg#kaXc0;<kW>*emb-R9j%G!|cVmmhYPM0=1RK-Vh0uklc
z*U96Ei-l#;oVoMOOD`^ThlT+2UE5a`9=ORchn`;qjjN45IrfAk$j$KZaB}RV4GUl6
z?e%y*V{SnBDTG+Sw+onZFXBZ5);U;03?b=o%5gv`90(>rh1Cm1*@Rf<f1b~p2Ec1i
zt=*H7ni_&J(QnH%uUpS&yxMcHC+T520xbnqh!G?|+W?KVhRNv(;Y`XnIT_~bOCcUR
zOuvxv6KmtUy?r>hp8Or{qV;`MXVw#lfJha8xBWpyUZxHRaJ$#I5u*!+T`z)=xQC-%
zyuB+$bx0BFZ%Pe^eMPpUx7QpC{LI4`aRfJQ?2V?HMfZzav=~5HP1r77xW%tu+!7`5
z;Dx{d0H&X={g#rQu9IxgeL>W3WeyethYCJ0-EkW_Tjryl(8R=1U``d!=3n*qcSb<>
z!`GF8O-v*?X2OsQlvBc43xeDXY2pSViSPBCI3}EKMWxM5HWl~7ook9x@P8Al(q*en
zx(}E6ydR<0rf3=&8OgR_0iCok&?BPCr6iJ#LhtxGfBW}{yM_xp`1G=qUqjxe<mJf}
z!lGz@($do67&w4#3RqC=>TPBRX=!O`8FrofA<SV&il7X^`<(3Jo9$&Tyu?4)QEY-x
zBB?-5t)XK=O0zq7Rv{xV&4pncw%bi(E-U++cecgJ;3;VNooUUfWA&N)b>5e~U4n(*
zTCB*%k>)g=egtP&!W|EH6nSm;nel`LW`U^6>-1YA{(Y>pv^3+`Qx>4XtEt<%5(Qt+
zlXDI@Ri!OaI7$oTmX0zmbhkkHZ||^Cs&)~!vKapO9LEy0k(&Zk^j6+%tKGeU>JHHT
zah4D)sC@*d(GqaDPXx!;2(!XM^D}c&48=oWUB?rFw|n=rpwH*PaK@pDnUcISVIZn#
zfQ~hr2K7A*G1h=<K$?SMFJ>+Aw8_-jNeLopP(p|=vjUt8dsq`LH(M~Gh7J$16LEo3
zLGYcNb-KX1R0PGy@3_#WO^3m%h&_vLVg#?OdoS=}BxO;Cx92iY7T=|$yh3X>AxM)0
y&u=`jh@6Ev#^UnF*-&z3^8decqj=?w({dbP&SR}ggTz&#P!&`Ur0hR-{(k@gYhE4z

literal 17440
zcmdsf2UJwqnr<nrbhm*vH!5I2zytyUlBEq;l2jA~2`WjlP~=pGR_rzx0+LaZWF><H
zBT7akNJhyyha%PcHhOR0nRV~XdvDEJGw-ohF9fPi?X&kEzVM&^<Cv=AvY*%gOrcPg
zQI8xvL7}j9QYZ^Xe_V{;BnYtY#<w46zfn*Ah%dJv4R7G*pDrHJqERRu=g5Bxp2$a;
zQYf1!)Po8qodSCsTpv{1UzVTi*tB$$)YW^%mzE#-bBo@{$~(>x!Xjcqt=HuxhEDWv
zueev(8rSJ^Nb-5(pS-%cp^;iX+G&5XAAhE?|4k8NZS3L8Uzcuw`D5XwJ$dp{CH(@W
zk;~m0Z{2+-IY4vD(;FxnxjGgdHBb~4mzkL<i8bw~+V8N(cglVRb_(UvZ^z0PP$;34
z)pq!m0A&#uzWw~|PYeG){!?i4+Fhsb3AuMEaqc<0Kg+(yhs)vfW<@UhlerT)1F06Y
z%2S$cxzqhS&mE;ShwU)8va%Ai|B-UGCD`ff9~+%p?{_uj(GTBV@7>+qy-O>Jm)2ff
zP$ZNW=&R?lotji4Ue~ImKXoEfB|xM01?7cjRc6Hk%B2a%HgC<yq$H7%_7bJCVA(Sz
zyU(P>9gL=2QhB+%pK|Q1P~ME~_-BK=6AwmsyFNv!F~Z~GcExHYmd2+y`X?nNJ@EE!
z8~voF#+d0+qZdl#4&CletF!8A%yD&@=nwKzlar*+wD2_1vRPA!timszJi+}9&rRNO
zITo+c)VUw69u*lE*DG5mKiPa$a<u%wr!d947SH?AB-<Yv49|>|rBA$LzjSwvX8D59
z7mJtvdaowgq$53V)^fa_R(9KpHr<!r5vxCICO6xq*3p|>rpoBe9x0_VoJZQ*{@@RO
z|Nj2g<d!?q?RV)j23A>JksEE^tZ@EvZC%n>)fs-Op8HrO*|_aHcNS&Xw0y*8&rIjB
zT(cbdvObL5UN5%ewI@fJShbE<dZWuui3U6M<&v%UWJZcbYfPhp#7p!SQ|vAVE3jW$
z7$9tL_>SCECvQQN{A|}om(EAbUK9C%t#@TabhGSCYx-w{73f&(rySP_4FjAE)zZ?^
z3is*Owo{+V2y|?8X3kCIg+zsTNKN-Yw`LeJhub*!UO46>q#u0f=31kuAcw-3Q?V!1
z)xF4`N;J9Aii9_;U!PQG^}LN|?hB9hZ8i#BK0SkN|79-wisOE=t`51AZ-v@3JA#Ye
z{IN>ZZE~>p>FedpHX%Bb-p|S_Vj{wZ^EX%<`szARaq=TL#q?s85INb@)YRQ=GE|qb
zZf|up+0@D*<Kl<+*X-12eW0#d(36<gm}HS!slh2_cY)h|^n=4>1!Zukb*&>sAcKuV
zX3s=DO|2?HH=rZfRdAxw)jHvMf0kR9l73-wnOwl;1J|<$a;@`b9r%<V-oxp)S8D2)
zuySpz+6ujQ)8{^g%sr^Bt)1x2D-wL3qD>!qyUvUE_4Q?Soy;9OcF2{Ruu<kbuPmpW
zxiLIEywR!ovgF_kw&b3K=k2#dOW!v*HU%jOX12OdNMaSD_T3e?I7RtzJJE4*%34}l
zhwp6kTed;mt$8WEu4BB;TFu^m%(Y$Gah$%J()@Zj?=j`+BB|y}Yz*GSoN;YNUT%TA
zfya1qh*8Cp<IHsWbh^(oG4m?Jw|}mQPuqYTEwYYbSr2ee*A;b4jdpouWo2=&NKOqd
zm7=Wv9ClaE|8`E<*kyL^BMg!hR0g-KOLmkUe}Tps%dz0uJfyYJ?9^C%{&fzXU=rCi
zCLxxbXL8+~`!ibZOju|hKYl+pA!|oZRlI8%b8e<Xo)r|bl;R{YpAMz66ZP4MW67a6
zt3;g!o8xq!i>t{D9rPIX3Na+{!f6tCpSOn2njIN@&CyY7o~V-NF)MPh`NdK$slD&G
zT`jEgX2(xzXc)i0x0%EXx5wm{QDkSOvtUN^m6=XO`NHFhT;8-Y+4MG(s;P;AkhUV<
zbe%_uPoBJ2k!X03h~TZkS&G=UzPG@ow6dLS=33<EwuUf^O@?O&==5Zk?1KG@%_dcp
z&=-+vA;HfFa;3=<^_<(h$#Q3>C(1r3@&w%!%o1>)8l^IJ>16CN32}F5FA4BXZ}N~1
zJAAt^#q?9s`L}=G;*u7hXrj}r;<UWV+}cGmzg*uSf~Y1b9q-}OhrKIt?g;Yb4{;Nz
zv(9rre3v^2``wCqpyclE9upIzo#x*A?49jcrKU+~prompS>bZ&-g^f)WcPLi(p1Nr
z=&X|2>K`avR}f&?iNXf?t)B8T=QQHAz2nnrV=8PYp$A5AviDjO<qO91a)$CZWX_J)
z(cW+N6zHHl<#bxiMLEK#sH$qevGvaT{^tY1*RNk!_7mB5`}S>4c@MYfumvyPk!$l(
za|f>pwVKFJYgy$o+fhM{uzN)hwyWEIKIiEzwVvX(i(Cv{?}aEDI^`JUXDNJJeaWO;
z-dNn?^SuZqU*kH<bQQ#ObPN8n?PGm?AZmgrqDd)AHKAec<tuZ-$#z{ID=I1;NC6HU
z`$frTSHT_W{?pyR9#mj`u5qZd^qLvV`+d!tHQQ9#D3k|B4Gav-%*?c<7Dhyc(u-&M
z(#HS%^Jmh02RqF#JC);O{)@WsXi4kPJDxI|cU-<6vr`f1Cld{h`mxOR_I4e&p|@Ri
z>DF0}1HnkA?8&b;<I}W|5UR3c6(>)g3^yqYmh6d3PBN>EDPam@NuLj{9xX8(>#oH4
zyNj5948L8f7PwDXqS48KPc1kwlhzTbc$b?~&dmvd^C_t$Ix5(|BZMv|YFzw2(PJ!D
zW^S@Ya=JHJqFNDGJG<(8U%hcR0L`u5bs$Ib%kQfgGh^M(NS@(%YxkW^bK~h=ygM}$
zluw;Hg?8$5f2(4<W8Spx^`)Hijpq{E@*@mo6jJ@kSL=ogcP57G=XubpGdFD5psuCm
z=e<)i^6}&IpWOBX<P6$f{c5u)hq89J_JgUZsfx<V)LW;JlAYy`NJ=sc$9rqKP#TmZ
z)fm_iyY7nZw2o5uCIbmEv505So>eXcY8tW=w5Krec3vu$6gz}i_sPVMqiEQbZXKmT
zpU|*<EX9pq{r*$t=treGNt<csFJ5)NQVY*K+)SUOMQ?q2ZAkVtdz<TZj>3thJP!sq
z8}D9ee0h6=grLE@JC<7378c4XDk?v%N==*7yz4DRK6P`(+Z;+vY^)dWvSrKS^U;j7
z+wo(vMNN{%$UydpDtcFAQc1pDU#%FTEPgHa_fJIyp-Y$9*dw#Dq(n_hEjZ0BEQ;$h
zY3gUr1n@@3$Cv!MW~T^!ZpO!w8-3a0^V7`-JeN20`T-AK(6_Q_QYP65Kof1`{^9W<
z=kXd7a%-alP4Wi$w+aZNP`J9-sQJZZl5kUFJ=MdJh^;2s5Y|k`NNS?{ZZCedQZzo>
zQppa(*7x`66C3Im$lZ}S%0{u+jzDTfcog1O<gQ7YOYeB42ry;(DYEU#62bS~m9aE@
zoC%6QFu(y3lq);bWup-3{*0}pe#lMNltrpXsd5URKWLRblwUHZg&V#=_h6qjh`f+z
zpV54+Ly7D1?qa<I;{x_Sr0vNynGAI+8&q(<Da~{8c;i>S6e}b+LB7BIKY)BkuQ-pC
zNFa+IF8y`SriSN^WiQ!z0)e|yizOr^I=<c%Qf<n0CwK=j-S&2!{;qT9u^+WpujuX3
zpY7HZ#evm(O!sL$J*yx$ky#dzkRX_tJ7vZ&{QBmPOw?9U=i%1kAZZC%m(llwUtT90
zm-w&b2R&#GJ667=wJ*J?w2VGCqKsVk<Z?LeP-EmfQ<OKOME+54_esUn%Bpm%T0exe
zmUldKxA8{TfqE1@Ex<9<=G+tS$hEeTYXvzr79dtW#h*#nXd!oE?J*LtOFKniY-}t#
zHrCuPHZHEX{LukX*YVyGu61B24&kBe-lL-j0ITSDFx$vQs}l7?cxF2lN!seU_N7)-
zSMyUz*Xh$db*fTi!V39gYGGlKid{G+^RvrKa`$AFqzNVHK5qvWE8ZRz$WRTEtCx@Y
zX%$~MDv(ihbhK&58@wlleNoO?heqd)Vt=unu^RD(jBQH8tzK%z4O#ZgR(^Ruz<yrp
zfddC(-Pe;;Y_zP;u+hr2H7ZInEMTb2apD@2&8tVRbbkm_oEq;7$eU>sa_)TOX*)CY
z&IBnFosdvk>?e9A?XR4C{-NdHWN#ELnytJR9rFQdHGwqFT@k*+K?835SJ^pA3Esex
z3Fgw&u{pGO+4=`ab>o@Ip|0vgebs$mem|l3w@o(&-s6<8(jiD-w7YWW*2A~ndu!_X
z3F>C?QVINs)60<}NI)a*^!)?NC3!(T=w$&-NRQ7+55y2+zEVE~sG(==*ts*=wyp5D
z`u^O1m0>x91%30WcTavMvtk8fb>fC*&UteA!#8p{-)09_#IY(3BJNP^en|*LN8G<a
z;NLV~bZ!b|QRsjBpY%5MR40~A4z-8^%GnJx%1RGp(**NohWAQHq+5If61XXpyU%@k
zTxFuaL0wbRcb&f5v)`jh3FY|N<+2W~L@cu_GJpUy`Ml|d0oxL>YLax!#);hOq!N$?
zTFY$-<Kfl<KWRrRm(kAqUtW8fl=zEnQ^pRwP^<)Toa7-e?vDK2<VIStFqIMHP`h1z
zZs_g^0+Nr}Q#k^fAh%A}?l^e~Z-HfcgAY=n-m%G@Q_|M3APL($_#xwp-X?Y)w~NN5
zfqtyH0hTI52C$u3V(sAt>?}o1y71vKvM4GcA;Bo5u&|K8)?|}1L8>thqSkX~=OSBB
zxyWbFgh})PQ+uP}cc9-ElP1p_otUV`IC<(+YHcDu+o-?6-`}6$p`^yopELgWNS1J%
zECGZ~jRr#6S}HwI{7Fqd=!|ua7-6ga_{Z1*9vZqP2!w?e&{3^0!5n>u6ol>#xOGTp
z4!&9<3iR8XTCK+{)XZg!l*!Al;!{4tIH{@0PbGQSS6yCSj_hyAzg<-5z4MIEHYFd6
zS`M7X3(v=>BKQn%)gW!w)*7?e_@r~Muk#kXxXh-ms~hAeYLZqgAuev*U!RHA*w88b
z+@a6^&p-dvSb{aYZQ&eUC<&fmO9(ujHCJPBS9Z8qw5%0Czr@cZ*pHXc+K0jhwr|&H
z8kUDDaPGT!HvcY92ryUM@2gY=Y+P^|xt|y1>?ity5z>^IFOhrSsK_S*c!e27pNb-I
zY<Hc`#g~h;a$M|7<)(Y%PCwh?xzY9rLGxA(S>Ahye#3K0i-xc015cTI#l?x50CcqL
zO<(^KJ*V7?6`x=AiW%NLp?G#l%U}3H?#W}fveohbE)oClRGv_jCG4k9uMbmqKF@C-
zI2v^6vF-nra&)OQV)=`Ab-PEN`w)1v!SjI>4$ETw|9(wyx!jxbNfo>AS|)<@;VtkN
zGbjH5x00^2NLD#~>(~3B1}me0%R^-aIv0spH)dy#*P8DxLhWSmg2C7_S|NUv%}z9h
z5YB|=SNJhPiT07o?=2!iA}}1cJSd;m<T0($qp9z1Pa4BU@HbNMBbEIas0pIr%}NaA
zASp3i!nM8F&+F=PsVA%NlS06s5W23SHr0G2jAvHBB-q&)=|yl?+?jMidwcucbxXNr
z_7Qg0Ta(;0vG3otE(Jg;cdcjV52?z}j`C8$t_q8aJ|Q(V>VQ8Za&zb0?~0y#^*gv_
z+Ac9Mu>cWq5+#f)ddJiwl?w>tLlC-$x7-50M5tJ3`RvmAsQMFGnOd-dP$)@B=gBXg
zpq>*cw@ztln)KDCk`IiH?F86Fr}Ahrum`TxeeTc@?u_?McIqftza}ja`JWOe&usPL
z6g5^R{7S{8k;qJHA#y}i(#|JD56o!WJpMeBaq-+*A-zB}aut6xr%NhoxDIQAHA{CM
z{h$bn#5ovtM{0+f9Mf=)8N!l;T2jL&*>Qc@1`(n!KqCs0W|b{l?(OYOh>3_p?`N5*
z&M-fDx<~20mRq2ME?1?|OxQz1)4veqibz%bRFA1oA;!47;#F!g_X&jj6sc0j+u??r
zPIeml8tB~aH{b{5q{kMN)%eYg6_QYIPUz}3bc;J-FVl2_u}lD2?cgSyy64DSArIjK
zeR?CQM7YjOeTSScdi3p5`<|-ZUVJKr;Z!g2f=uj`0GDQY-BZwiLsWR?SRy?L8`U}z
z!g6C6c72TK0s#8s6JA$-U!c=6qut+p_%+Ym9>t9|2UR889zbK0qcZg8hJAOZ=aMGr
zdhzD;^aa2eetBlWF1_r11o;CQJom}OOQliX?0=fS9LoKioPR(-GL)V0+}zw#U^Jv{
z07!2?PGVYW%c5C#LU67wURJolI{SFH0DkMXlq_SaJ9dP?jJ(-CYbCBlIH%>MqMere
z`cgg(SDmpo3UME&Iy*aqHa41hKswuY6uGruF${Y?)srY>{{yZh$G3V^90wLtu>wC$
zJ{(*+#ro3bt6zF(oi3Cu$!$eGyOkI+5F(xDm>o)iOYtdp4q}_*w33sm^<0fmhoPpq
z67I=cfY9`X))fAx3m5A=f=lZ-eYuwT_0OFmK(-HjebaAoY!KfCsc-wa`MX=cC}-J{
zwY%O!mYfV8l5?MG1w0i}*n+o9?ZDq;UB`>fp*sZ+As_jqd!HWoX(c&^<J&tee?`;G
zMXQd?&6P7UGOBrc?!erLJWG^NQCh8Ol(!T&F2i->O5=5+7oi=hB4r6zg5cH)*{_Iz
zb3oxI6+>jX`)8KBY(UP^igu^B0;-o@S+s=QmRq-aNMWeL)uUX?$QB<~U`eB4hDSz5
zR*=5e`J(=7<MkC=54ZZ8Cpec)5~wg-!0)?>Y$%sEJ{?F~hOpj7Q4eHK2Wb93+#x{g
zQD<kSYyNY)E=GQS{z_iZ$k4D82uEyh8wnuS{DG{Vgr%G!;XnLAKJZC4Z^E*ccX!LE
z{6P<V^*to#Yg%!jq$oifG;llLeMBhX`U4*hp;wS72?<%<F*A@~m4JVeAwhOPbuPt!
zdfg;->^A+3bAM=Xr$G4KM)mGr$*&;dvZRS1m6$zz6y>1atY$S4hi9N>84w+tXx2EJ
zLg=BDB2E7kV%c~Z;<(O^951gM#Fn;~AK)MWO!VxFOG`KI+r!4byjWGbk09c-X(sX~
zuT?aJshowa00!D4ls)*{Twl}NZpb>iMdGXYRf}&+G(Ifw;x|Fqbpf&yjS%W%iYqy4
zah~*y^e;qr2fH7E0=|RbJuR&oKG`OsR|D#x@TVH}(pfhz5nY8oJC@n`_$I;n{8S=5
z0~irAp$H1OgsB|4Bo~6POm5nuy(}c;<yCfTpOIOtvKZWp9I^@IlbUm$T*J4m-L~pa
z*H$xlZ?4^a6ty3`Lwg(m={rTTAYJ$DlJ>O|;y6d%#JssFW9cTIxj|1t%ZR4U4SC#c
zwnx8=$aMs4QBUaWhww(n#F(Y(qk$CFKS$H}OWpWhhiqRiTG#sQ9d{9G!?S{${V2gn
z=znc7u?dAya%0hf?&CGg+Ui^%Y#P1OK=t9C%_z#dm5>S2e5BRp;6hJgV~YYMX~Q*T
z9^$Ax1z?1}eQ9-q$Ucj)uW92S)2*9G_*H%eNo*X>pZq(kviv8k%ApK-@<OqLb5xuZ
z{eUqyB7G)Z{PIe}fHm913R=Lnn3vhJqy@T&C`v&uIx&@mJ|VXHtf^`Vmts)ANb~B1
z0P%e|1Rg)c@oNrwx8&df0su(bl8j~za<J|ZWs?a1q^a>zq5b9G;`G{?kg<gnOn^pF
zkD1BQo@ycfz>BYpOW+FbpFUB}S`HAaL0Vj@d3D<QCMa5UV>0aAPG|dbCK^ctMn}_H
zynN$B0@<U2UCfHa*XdkT&`8t^A_aV>MjQmrjSNUhQD@Sv`us`pEm1;QB$`T~O+Kgd
z;7fK1DX3#p4)Ps1p8x>Sarcg*!a`ec@wTKQLGuY3TE#Iby_MD1A*&JrolpxR^5L=_
zDBvY1I>SZ!b4lOp6KUo!PZ;_$U)M=Ok^hJePILz1hQSCB0k|N<BjV1FxJ&jgKf7!o
z>qSxJrIOw=!z`o23PLY*p6ZHH;{0CCtY0JYtA>+Jp!cbHOupuscOO~{yguR7S4Lrt
zp<8Wkx^Bn?2WiKfl^8m<1>2a=-RD)Fwzf7uaZ7E}*Q6{zLNHCfqtCqaksruo%`u5{
z(ZrgYASRsL&>w0K_7K9&b^FTP0*jonDx#Ag<-#YxQpLSNZV<JsO(BI0U8e}4tX(2M
z{C?{{^sRuV8DpX@8r+fObM4PG0YW><^%J?kQCvXRvB5iJZZwMdYF+NJ`9^7X(Ze`?
z8P<qR*ZuZUcNZX~o2IJLZC@{w2#48PaU5&vY_9)0C&(XUcOrI7k=vP<LvHdcarl6V
zA}Z@)DkMvWkvqJ+y>s9MN{ZR+n^kJNj<3X9<y2RXI{Hc3n?-NMPlLj(-&iySAA_%X
z21FEHD3pZFXcParT#xh&&??!R>x4?dTFG8R^7fV@7vy??cyF3%mzW14?h_5^e_vwb
zb$|c$d6x<0G=B80kC$Crr}vN){)jt7fMoALjd5cZ>6a1`1r<k<pbb>ikN@K^fa%1@
zEK#R$h}yf+o~``!lZ{n~MFFGKB)~GQ;KqurAHfv+37`5r130nk^wUlJuy$SZr{OJC
z>(kx60%^#s_W%i1N!YSHQ$K0e%m|G&JJEq+RZ3tA?mf`%BQ<~*P?i$l>-#445fpRh
zzSL^s`)IV3$&4sFz+_OdEiDa_j>^iKl5GJ1FDBl@9{wF~|27NI5#(t6`RO?-19@kR
zo?3>QX<M$qZpo0rg)O)fLetMdK4tbh<^@1#CKd)dGqa2}Q|2~Z*Hj97G>uxyL?4VK
zT$0FBXq{pR0aZqEh`c<MQzczp-J&!CRakSgL|Ri}%$F8LME5C8U2n9YsK`h@qIBKc
ze1JG3ZOvCX2w1D=@7MtIkt7ExACUo8g^v&2w1t$_ULw(?{10?w8eQ%i9jSp2g1Y;i
zIwNj?&^mn1SQS8-3|zEitzl)%DP~2KoHz7iOU4KnjLqG2NQOD0cAe#<@(@Zw)JzWh
zVn}nXAC-hG6K{~c*z)sK(<3aD@SNFfg5uzK`V-Tg6f5F{!a7aK`m?dT{OJdiDwZ*1
zIZkN@iy{|@S5!H|<t=G9T`$LLr|uwTJCWwV`{dyp7Z5qFeVw4r!?!orMWv;UyWjnQ
zDh?kcMH_Du;MRnCS=Th@0ADAAEK<o=uw<LB4g-GGj;KJ}qtM1^Raz!Sy`P_TRmN(f
za@zK1c0{J7RYCR!?nL75CkkAkZQ*ujh(1a1vpVh_SYCkRCjaKA68Dy1a_DnM)Ya7q
z2;rG&X5XfYlDglOsHb>^87&vB!EdkW2|gA76h8u6pJobl4h-wF#`NhEY1qaIK4qlr
zi+6uqnG+-iHxdD!RrUvFUI}UeaX|wl=YgCzUI1XR;~QkZj(S6;A->m}@K6?)9$+@%
z8R(goVsqd<&koK*#S6%pnJ*l40<a(AL{nS@b6_n|Wv-%x?0?>MvYT>f9NLm7L1IZz
zMPZ$GFgrrZTJFjRvn}7Z3AuGTx@ErdWuPh%KNe0v0Y9oa*?}M_`$EzP$?3yWED|k~
zQST-ZK;**obGU<fPmMQZ@(n-9LteW3{rXqemVEC%`D;eSdqE$ZaR7hwV1e?K!CFW^
zcY|rW%bp~c)Gw%}L^>lKjW@G5sfbX(c{n-Hm_w`9b@cuV@<N$QpK(S9FBqV+rg!P*
z%FN@(0CKvavQCS0mJxE_=r&TKVUXi8R!9mZ7M_d&4>W=s*s^+;o_J?WCNzo8!4dC3
z_uQ%QjgzEjSRX`)g0c}!k+_Rcn$7*r<axMr!+9s>((i<_cVg#ncVfG6rYAHse9)T4
zWJf?#h}EQYZ=OaSRcN}Q$m71Zq@;vX+~Tyy+(<~2W}@D18u}GEU-U!GSHF|Z-23y?
zVIpKu85nzrKwbq&IiKV~+fyN!<I(~);sp~|l}MKiUf14dilvyyfYjjgeOr`XTRL?|
zM<tDvN_T)v6b-i(CA%?)b)d=D0CAEz2jA~}X(uQoO6WX$POrCBl5(8u(XWJd;zg`M
zA|P<nXNFetgFc4BcBmVyP*+Rc(LwzPTal<@1ezRuIDd!vyZbGlz!?Af_GXi4E);h%
z#j@=P-nm~f0__Xn#<o<dCvFC{umCy*5#{Uct96;(ntA?UpAp}Y59yZl*=joRcEJ}?
z5Ngs?24wNp9(|TeLhe-ecc22EaZFDyc!O2OW+;yD$)pc2xxZw@weP-f;Nwhv-Nkl4
z6H*q7fwF>f#%p(C8#L1zU3%sP_F(6BKG_M(U<7>=?jwyXAoOy|y4VrcW0_XNM-u<Q
zf!BN`IO!Mdq!AOi27`Szc!)Uk=?xBbeAqGPk$If}sqPp%L}XCn)%Sy85W9lZB`~;@
z0$UvMQ>H-Yh7X@VKO{yd2EeRyna1lkY&hdHZ%8y>X8#za7NTO?e=1g^4f|cPJH4I(
zh@N2N2G77x#OP=OZeuC9_4Gg<JyJEuA6;?_z;z9;7Owu;P3CIu3Y1|o@PUX-!f;5g
zNiHe@N^%#U<V5F(Fd)Zw;qaS?OHnjA!yHM5lsRv`F8fRNqLFqt>O`b32*1Put>-2B
z$s`GPC$L+R70_$jFqPqqS(!=^O&kXRKXH5eJQ`6rcdE(H+F_JoBAB?b=r~bm#_{>F
zF)<)KtlW2jd(OS$oAw2$qonXGa^JdKerAxj0KM#E9-W1rYB@|3#2^>JN5n28eq+NO
zX{-7S%tE9Y<v{pG)9G_?!~1JJP0H)`fNPeVDD}!~e3c6qwioBq1#?x(u|H5N*|-e`
zlQQU{`l(ai=px2sjAh*X!i6K4f|=${M?H;x^5nuddkXU48_Xxv0EC2AekR5a+HuG8
zrdccO1(Zl%f?P>x%uj5PwvKtG?OT9v%kRKq`V3l?45=h-eGix`La9N-<Dd4Sad-V+
z099W298hSwqT|?L-&GN<PV`Yo4oi9DWaKAbavbGCkG$j%8p;aN1C1rri5MUlt9h8~
zK5aEjlLLB<uX>JwIjR?*(Z_IVz)U^1#6l09`3pH6l6So*1zPD=x@*x*@v=6|3klfo
z#ml6Aw6<xwfnZD3&O^Z>4n}(ce+ZOx884!}0@AQ_9uv37?jUMYF9G}LIv<i6Q3Q5j
zSk5Tt(WA|%V-Yuq=EUB%cH*`Zun19nl>@}*#kXY=_4hD$RiuU5386NK^b-hvJMERp
z0)D~KJ3z&m0i96vz_eTk=5^-DF+mYRIin|tw@C01x@^Q9gv);RuYfS5Yx_x9>ji@z
zkSQ1<3ynknXxm8ktK;XXX_RU9VVs&EAv|*+q_RXwChiC^SP1b02h;o3#Wp2`TE9up
z9g8QNj$C%GD`ajvG#6(QmNK)WpSr5z&wyvEwSkrH`fKR~hk;LS!s%&}xgFF7S?USo
z0!k7wbkv9*voea*(5=VIB_3lOLo4ONaRv(DwgKmMH%uQ9^soEp1+vG-u*b@v?igUw
zMpy6!8I~c7h~A3aD}U)stv?*-(j7fd%H+}SiJ;&z-mBGpguMGAZn_f+CPGekUfOX}
z&Uj03#eeDiXwH5zOV|tkjEp&Q=+GfEbMp~^)3oA5tz<q-5vrn7Cu8E$5W=nLg$n}T
ztnZsJn*h`dND(I>SWR~9B$%HM0Ek?fJoMLpsrl~{!!7hZn%DdZ)`RA+=WIk40ffj;
zuKB;K`D057gfK$ABI6eSQ2Q*%3>lfEAQs!X6ZpW|tFYLT&%V3@5gYEiSdDRjr7fTw
znHM9Y7wS0UU|H9++T56!9q`HrT(gn3d&R_3<6v7^CH3MOZcZpA#h~3jFmuI03`?r>
z=EjI@4#n%0r(EOjKdmApGYL9liUC2jwb-UT-)U?TYi`;E_8Xb)ndh(=8$q2SQzF`y
zP37hHpw3LlHez-&2^O%A5h5yJYMz>IhOLs0>WtMGfgq{*3hIiWNBCr3Te8l0=<8eD
zynE9u47J+fWJ~%Pe69d@8SFS*$`-Xf3Wge)Z=it)wG+*Wc;4tHf%a8e1iX{kNH}qU
z$UIzWY_TKZ<@@0Z!aZw)`BMav=7ZOh7@7CjD~>B}`_I9&Vc~=AmZgsWAla)EEGE9L
z?7!W0{y#LoT`rR(CNfhR00OWgKUAsTP<%g|1&|ra(asOVX-Gp%-TGB{3)F*cM}edE
zyB;rJVhE{`KoDYC0fQNm0SDq|fWf3&(4ftC%uNrJkzr4CH6y4Nl3n3mM<1@o8H|_D
zo6{H;uuheOz(PKFq@%11!le)3hB5Xx6M#w-lSRA*Xf?^;B4ij)vNQDR(G?p=-E4+O
z9<P%rLj8!zBHn_3l_Au_Bvy$Dor*9cGY|99ADSK?wY#V1@41ooEjKpEs)9a2mesDp
z7?^}EsdaLD!QafV6U;tqdH^wdwsMhZ|GU||eKTIwM};9qE0yHWo2pD4?hLzY46oXP
zGZE7oFz}LrhBfLW_U;YO%5n*Nw;m%E8VNczNf-b(^fD&X7O%0R(0(Mzj1rOjWFRvU
zD>#<<k1R8v^u=#1(;@6#FRl!-I4|I=zH5Rb!79ndjoj3q*#d#hvWV6O-bTi}{}Yax
zb`kqYlwT@i>z#d%D<`q*G>CHll3E5xZYOuGf@Fel&7uVWn~|mtdil|q;OQ6bk3-qY
zA2o>NnUlG-Oq(Oot0;puZQE?{zQ2!3{!f|K_-A?Ja(RC05a#{IMEyJfJ)*>$i~Vz7
z;n_h`biNB}*hT``Nhc>P`vdcjRs8caL(nfp9sBEXLR-nWI?QHq+I;yw({X~mnZ6NH
zh9tyk38`uT*8gZOqqY4O5sjs3i`<ty`-j<!I2igE!z&?Puwx!6@#5lX-vkJ}^7!Ah
zdoevS=aeo?<{!ZRiZH95Or?SK4}KBo<bvw)SCNjmL(cKAK<Akof7jQCV2S%=-tPK0
zp?4+s+OLS!@AA3Nv+vS5M!L5t6SU<zxOQTpR;t-4VEYK4F;K=f^a_<ky&M9e|AR@S
zm65~|Am;lpI7Zw1t#>W7cu?3GhM4sxD(}CTK~hV%YB;(E?CO#Onci5TxbZ%sbRtLy
zvlh@@ePC`{h(_Y3&(g>gI2qdgpB#zbSD3rVAL0$cpx_w0$ywtegKv<^4WM^zJG?n+
z{vECWIG((O)S&x^5PS~tdgnZfmRpEWph0JBH`kxmCeAl8f5@~Y;hpf)=kziEj96HK
zL%&Dp$J`8)4A*h?pdcU!M4$k!<O6|d4(rzGrp(L|tRJVJCkI1lo-z$jH?cgzrqR$N
zFcW9w4l+1zrKqeMWjj4R9uRTo#o9;#R*g~7mwz{9GP7F|(<y&=4NCm%!Y%>pKnKAx
zz?m2wy<K7bWT(fZT7)Xa&Mw#Fl1luACQ-siQPh;Uu41VV)ZD(VB4RcXXU(#K^eZbX
zH#C_1O~jibE)>*-8CG8dUAeN7=iSh{EqSdmn*`UiDjrci{>!36TYZapCJml+{;=%r
z!*dbgf?~GYA8uEFc|V}0m}f{)DsZz_*xpU{7v_FkW29B8yiNU+cY7I=t#smUZoNI5
z|N7=iS)MwR=M#Y=k{xcP4rR01105rh)kQ^0%a<>|U8^sX1w>R(U9E`;u$$VcW@8!8
z*_yY%kzorJ7Jetgc}qCUV>T;0Cx^#TbCCW+r`j(jCMMTBJ!8Ye4+5$BjB38nJ#peO
z;Fr{djg8Iyn>TN&_<MRPhE1$l76w-27QNuY$zoB;`c<@Bvn_?0?wj}1Uu84f1N>#2
zleTQxQW2%blLMHR<vbF}@bcnutVYmND69nwrdIK7-mKs&q_5%Pk|n5@{p9J>J(f+m
zIaRMnp{88oUxJqi04+A2sU5PrXk&By@ZrPp?^aPlAMGilTu?3+ei3G~CoY<Uy!wSf
zb*~Lj7Ww*xfFWeqx*y53ECdC%A3Di*%I5_O?LwDBL`{OSuDjE{IziX=;>B{{2K%pX
zembP67>jukTYLLcXU~Sh<xF&)=+`O|(#be>?%X+X`v*Qgg|A;907LpHUwxstsL1Q;
z)dNA664x;oY6i%QCr%hwuU>_19uTfeN%OZY)q3s_kH;c})YR0RU0lpCV<qJ@7+NBS
zgR!@Q{BT>8Lh-LaR~4hRD=%8KsGzd)6b4z;Po8`rH#2wuxJHSw`j=nM&i1BX`0(52
z&6~ART@3mq2S;11t1c<Fwzh_chw}oeD{*skOSz0bIdS5|NnPFe;^JbdiS4^~RiLJh
zLbcl|E8D;~04m@wYH|k)YFgdk%iECctPYHQ<!eieH-qufuy3K0(=xFRw00#_=yFIj
zLYkVIc-l+(!1bjcpw<ZS@bFMS0v}s+GrMCG_Zk(hz73thhGXUauxSXT<VgVz4h|z@
z;~l)b2L%NMW1c)Ihp5Mk_ubwod4iF0;lq}Q@bFI1#M&gIn^<gXOG_Bgk!*VyC0I{K
zft~v}b{<hu8K%g=oSrs^g>Y0^`J}0-X;7PlWo<aj!3&t)NWC%w)U=I{uNxo_hjljI
zoo(U56sMu{+qP{hIRET@Q<HuOeJ)c_KR51p*kSxeVDH|O=sDukeKpC7m|OA<2ncX<
zcNa6rUy6PB<S%dfkvAeL${a_uY|j?Tdi4V(k{cwfD)4N{>a}YlA+KG)X~ZNXXi%wt
zq$*F3_gSV{)}?KilBzZKII5(i4&9{k@u4M)7B4n?d1aBf=E$}#Oj#ZS&ui=Mj=^hp
z2nr@#R$TPrwZFKfzrTca0-EOAckg16oIfvJil!2F^7QE^pFf{AukqP^CI)sretl9~
z+rKu&H0<%?Edf)*BO{gZXCyK+G6Xhl+7I7#vU>csuy)49M>8`stzW-xv9`7@c_{qa
zpWac*RpKYAC?_Xp>*%OyXc!I}W2CNVWtD(eR$!TL-@Y}4yfZvFI5}{YZ9Xk7&p`LA
z&2@jShnYTYXCAk~%Qy&N*ToQ7s|xi~F5HJ0$f@Ed645wCv$b8lW=+Ia#XA?kV9^^G
z_}Xjp^&pp=+i3{lI9kS78WYp*rg%_cy|hEz9WLp}k^l*FJV<7r#?f3bqTe2vM<f|c
zPDF`Y)>injpbJDlefqR&rXtBmsja=e8)l$(rfmc)lGwn&jpiQb&OQ42_3K`@$%|Z)
zwg++4-Jhb$wr<>h>H09Fz8*lAD&YCH{(e0ZPIh)rZBvxu!m_ez0Q=V7UJDg@yUu%v
zb}^3`r&A|SUik9rn)cJPmu+dZp1uPf9=W!*wkgIX>w{#RPhecg@3-H6i({gb)TUd9
zAhqp!KMUh=2URq`N(>Neqp~mGSPon4fxmyu@USh4hdF5D-F;gwafwjO{80RC?CkdT
z_xFnj#Ky+nH%!pUJiogxQD460>(}$Ve(=rw5yR$shYw%3ac~f46oYy9Sj48JR7_4z
zJ}mMTI#b6jX7)gZlTN3LVrP4%m#h=~m1>OWs!Bn0a8lc09_0T@nEiREoZ7MR;NX6E
zHl{8-M9#hTh=}5yjm{6xEiB)UO2mwULD=xa;>A_|>&zfZ6+5}OR0Zg71uFU7)3a_=
zR9JZE@)RZ1w6`g*>4*q)q*Us}@NiP8G^PW5gxeEo9B2zI{ry%5T`OJypxTFM-r5tc
z81OD`etFfh+h5u-;iu)x6Vf>lntOtz9mO$h74Pi`?WLu&^Mc3hm}bBzrd#~+3;?IW
zEx$eIuHdPg3>A;xe!F7YkTu43i*lBgab<bMRw*fM@O$wyx;Mit8nd-LIe4m|A(8aJ
zDF)0Slw#lwfIY{I>2zSER}9G;ulrmB2{~Im_xtth+WJr{RT&l@HVITROT^m+QMa(9
z#QW{}XI8?(!pP$)xD!0*c`19ZMLwj>nSvyuZmTOOT-MRm?b(%#ut#~RDs0MjPJaCO
zaYe4XtJbcyB_Z<qU^P`UIPhf~_YEp=@Vqdv8*Ena_Vta$o@LQG?ufZ#9%qM;P&EvZ
zu=nqew6(Sh2nf9QXF;eE*s-H9GU3S+0TN@tErrF!DUuwZbEHt)O^;imeC*h^?GRc)
z7M+#L8~yn4Z+KaEeP(ir$8*PlDnw13CAwg3N)3oX{`&Qy1<#IbJjUVZ=olchnF=9S
zC*BsbbByrj=CgRGX?I1mKNOgaTwL2`?KfPov3a(bgF6n#(N~|T&Lw4cZ?xL769o5}
zb7O2)mNugz+nE-$&*nGeFk`fO;c+oDEBV-+@y09^!%o-$sm=+p*F!>@UN2pAiLZP$
zt+dI=#KauChw_%jTrfC*T1a}SWn*>?pF1A#sInJ5fBt+0j3FY~$hda5w~t?BqlA7S
z?e67E1%B1Qr_@9({jVV*b)-rP2pmyRP$2CA`L}w-ibu`O&Fp3a7`9M%bITsfNl#B7
zn~BWGI1N>eR8rUyDXtS%4H-5E?%cT}g$%b>jXn45G<Vy*%lyk1NR?j#H?YSLJcdN;
z+@2xqWAJ9u(O?y&&^lC|U0q{gKcH7VjWz)0Bxnc>*;IA~eKj1hk3%i*%!3BmoLppj
z1Jw=Ce6*7c|G*IEqZ>DFxb#5E&w2IwHR)FYqwqX5=H}*Hpb5p9Ldj3vzH0sY7zShI
zHhz9q_3IzX_h)Bk?{ykH3oG}ib@tHVoU%vNZ-LBy`spX6RWR-!+CXh;?0A2Jxw;j4
z2RV71T_wU_syh85jizp45qI9eK+Jv0u7`OH#T17OwOz1sI#C?hr4spN@o{{hiS^>$
zg3r<m?V4BMDUai-s;X6a{1`a?1RLhTy?bYAM$fN%dU`^kTD@vj7(9s!fI`FL;~G%n
z&Aq)&?wmar5fKq>-IRM&MWu3QV)hn97?-?<4t5<7(wp};{Y{vaV~ZjGUx2Q=@^i<G
zj3P=y<b|kJ`W5CHR*TQduQw~zO0)3$_~}!s=_kHLKm5=M=Zj&O<DHTZbe{qpOiI(x
z17_z*JfoPMEi=~NaOT>zYY2d++r94Je}p@ULH}4U?~&c#l$Qw5b^<V9nT!a9J0H(l
z0c=Fvxbfp|fU`%B9+_3e`GF!!&rY`7{h5O!(Xk1`@iBmJ!0GA-5B_cpTLn@?59AGX
zl$-Isx+p+bQy{0Ys`~)}FkE&VI(#?=^l%FDertO>AF3{#*X_c>$Ej4R_OtWXpdm-2
zDic_Y4_q&N{&zgu^a#035&-SH_L+fansrm-uLl;%aNnKY?DT@WGWA%1aDH!n=BKQ)
zMDKyc?)&WYuKb3Eh9B731B#241Fk&~5fM?o%PnisJqg4`9>7&TcI=)c_&Qr^-@t%1
znv%swUcW0>7WPbMt#y$Ge&m#KdYXC{i!g?qSv7;*|9m7S)yU9L=JWmq7k*KZ?3kRM
zmM`z{fFzp_LscioCACzpEJ(WA$=R9n=G*ftSFTh>NW>rV2c%FreflA3Kq#I-L${pY
ztXNon4@YSYdR~G7FY~~_z*}$63u72g7?8Cc3V!~>tTNeK!hC#u`LNlFr{*xdBlhy;
z%lv|Zf=MRRFD@>QZuvYk^u@u$hy5{R%8B5?RFA&8j*dl1p6f&cRZUG|!GZ<k^_ZS`
z`R46gKRiXKRRlWMh(7xvBQw(vNXRQFsAP+cli9vsM_4RFaAt4<Wq=s;xd>=)a{UPu
z7H4e6H?YtCcd%dom1ruIOI*$KPsLIw4wooI;-^rg_h0({r%n6G01V~D^1pF~zux{2
nT!FlV(o0!DSn&VbmtnDKtP$H_o_$Z2+zIuN>cNEH&R+Rn#&aY)

diff --git a/intervene_test/Enhancer_predications_updated.ipynb b/intervene_test/Enhancer_predications_updated.ipynb
new file mode 100644
index 0000000..31eb030
--- /dev/null
+++ b/intervene_test/Enhancer_predications_updated.ipynb
@@ -0,0 +1,175 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 9,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "([<matplotlib.axis.YTick at 0x10f6ed710>,\n",
+       "  <matplotlib.axis.YTick at 0x10fdad550>,\n",
+       "  <matplotlib.axis.YTick at 0x10fd9d8d0>,\n",
+       "  <matplotlib.axis.YTick at 0x10fee4350>,\n",
+       "  <matplotlib.axis.YTick at 0x10fee48d0>,\n",
+       "  <matplotlib.axis.YTick at 0x10fee4e50>,\n",
+       "  <matplotlib.axis.YTick at 0x10fee4ad0>,\n",
+       "  <matplotlib.axis.YTick at 0x10feee210>,\n",
+       "  <matplotlib.axis.YTick at 0x10feee790>,\n",
+       "  <matplotlib.axis.YTick at 0x10feeed10>,\n",
+       "  <matplotlib.axis.YTick at 0x10fef72d0>],\n",
+       " <a list of 11 Text yticklabel objects>)"
+      ]
+     },
+     "execution_count": 9,
+     "metadata": {},
+     "output_type": "execute_result"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmcAAAE/CAYAAAADh2QWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAHq5JREFUeJzt3X20XXV95/H3RyLlQSU8pDEQIhQiDLaK9RZB1CJgRRGJSCnoKtGhplq0OjpTaev4NO1aKFarrZXGYokzilJAQcaiTIRaGUACRh4VEKVCA4nyKFgE/M4fZ185ZPJwyL3n7H1z3q+17jpn/87e5/c9uXudfO7vtx9SVUiSJKkbntR2AZIkSXqM4UySJKlDDGeSJEkdYjiTJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mbnSQvTPJ/k9yb5K4klyT5rSSvT/LNtuuTpA2Z1XYBkjSdkjwNOB94M3AmsCXwIuChNuuSpEE5ciZpc/NMgKo6o6oeraqfVdXXgIeBU4EDkvw0yT0ASQ5P8u0k9yX5UZL39b9ZkuOT3JrkJ0n+e5IfJjm0ee1JSU5K8v3m9TOT7DDajytpc2M4k7S5uRF4NMmyJC9Psj1AVd0AvAm4tKqeUlWzm/UfAI4HZgOHA29OsgggyT7A3wGvA+YB2wG79PX1VmAR8NvAzsDdwCeG/PkkbeYMZ5I2K1V1H/BCoIBPAWuSnJdk7nrWv7iqrqmqX1TV1cAZ9MIWwNHAl6vqm1X1c+A9zftOehPw51V1W1U9BLwPODqJh4xI2mSGM0mbnaq6oapeX1XzgV+nN6r11+taN8nzk1yUZE2Se+kFrp2al3cGftT3vg8CP+nb/BnAF5Pc00yT3gA8CqwzCErSIAxnkjZrVfVd4HR6Ia3WscrngPOAXatqO3rHpaV5bRUwf3LFJFsDO/Zt+yPg5VU1u+9nq6q6ffo/iaRxYTiTtFlJsneSdyaZ3yzvChwHXAbcCcxPsmXfJk8F7qqq/0iyH/DavtfOAo5I8oJmm/fxWHCDXpD7yyTPaPqak+TIYX02SePBcCZpc3M/8Hzg8iQP0Atl1wLvBL4OXAfckeTHzfp/BHwgyf30jik7c/KNquo6egf9f57eKNpPgdU8dlmOj9Ebdftas/1lTd+StMlSta5RfknS2pI8BbgHWFhVP2i7HkmbJ0fOJGkDkhyRZJsk2wIfBq4BfthuVZI2Z4YzSdqwI4F/b34WAseWUw6ShshpTUmSpA5x5EySJKlDDGeSJEkdMqNvMbLTTjvVbrvt1nYZkiRJG3XllVf+uKrmbGy9GR3OdtttN1asWNF2GZIkSRuV5NZB1nNaU5IkqUMMZ5IkSR1iOJMkSeoQw5kkSVKHGM4kSZI6xHAmSZLUIYYzSZKkDhlaOEvy6SSrk1zb17ZDkguT3NQ8bt+0J8nHk9yc5OokvzmsuiRJkrpsmCNnpwOHrdV2ErC8qhYCy5tlgJcDC5ufJcAnh1iXJElSZw0tnFXVN4C71mo+EljWPF8GLOpr/0z1XAbMTjJvWLVJkiR11aiPOZtbVaua53cAc5vnuwA/6lvvtqZNkiRprLR2b82qqiT1RLdLsoTe1CcLFiyY9rrWdv2vbfT+pGNvn1vWtF1Ct30ubVfQfa99wl8F48V9aOPchzbMfWjjOrQPjTqc3ZlkXlWtaqYtVzfttwO79q03v2n7/1TVUmApwMTERHf+JaX1uP7dO7VdQuft89q2K+g296GNcx/S5mTU05rnAYub54uBc/vaj2/O2twfuLdv+lOSJGlsDG3kLMkZwEHATkluA94LnAycmeQE4FbgmGb1rwCvAG4GHgTeMKy6JEmSumxo4ayqjlvPS4esY90CThxWLZIkSTOFdwiQJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mSpA4xnEmSJHWI4UySJKlDDGeSJEkdYjiTJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mSpA6Z1XYBkiRpuK5/905tl9B5+7y27Qoe48iZJElShxjOJEmSOsRwJkmS1CGthLMkb0tybZLrkry9adshyYVJbmoet2+jNkmSpDaNPJwl+XXgjcB+wHOAVybZEzgJWF5VC4HlzbIkSdJYaWPk7D8Bl1fVg1X1CPAvwFHAkcCyZp1lwKIWapMkSWpVG+HsWuBFSXZMsg3wCmBXYG5VrWrWuQOY20JtkiRJrRr5dc6q6oYkHwS+BjwArAQeXWudSlLr2j7JEmAJwIIFC4ZcrSRJ0mi1ckJAVZ1WVc+rqhcDdwM3AncmmQfQPK5ez7ZLq2qiqibmzJkzuqIlSZJGoK2zNX+1eVxA73izzwHnAYubVRYD57ZRmyRJUpvaun3T2Ul2BB4GTqyqe5KcDJyZ5ATgVuCYlmqTJElqTSvhrKpetI62nwCHtFCOJElSZ3iHAEmSpA4xnEmSJHWI4UySJKlDDGeSJEkdYjiTJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mSpA5p6/ZNkiRpRE59/6Vtl9B5H2+7gD6OnEmSJHWI4UySJKlDDGeSJEkdYjiTJEnqEMOZJElShxjOJEmSOsRwJkmS1CGthLMk/yXJdUmuTXJGkq2S7J7k8iQ3J/lCki3bqE2SJKlNIw9nSXYB/hiYqKpfB7YAjgU+CHy0qvYE7gZOGHVtkiRJbWtrWnMWsHWSWcA2wCrgYOCs5vVlwKKWapMkSWrNyMNZVd0OfBj4N3qh7F7gSuCeqnqkWe02YJdR1yZJktS2NqY1tweOBHYHdga2BQ57AtsvSbIiyYo1a9YMqUpJkqR2tDGteSjwg6paU1UPA+cABwKzm2lOgPnA7evauKqWVtVEVU3MmTNnNBVLkiSNSBvh7N+A/ZNskyTAIcD1wEXA0c06i4FzW6hNkiSpVW0cc3Y5vQP/rwKuaWpYCrwLeEeSm4EdgdNGXZskSVLbZm18lelXVe8F3rtW8y3Afi2UI0mS1BneIUCSJKlDDGeSJEkdYjiTJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mSpA4xnEmSJHWI4UySJKlDDGeSJEkdYjiTJEnqkFburTmTPHr+xW2XIEmSxogjZ5IkSR1iOJMkSeoQw5kkSVKHGM4kSZI6xHAmSZLUISMPZ0n2SrKy7+e+JG9PskOSC5Pc1DxuP+raJEmS2jbycFZV36uqfatqX+B5wIPAF4GTgOVVtRBY3ixLkiSNlbanNQ8Bvl9VtwJHAsua9mXAotaqkiRJaknb4exY4Izm+dyqWtU8vwOY205JkiRJ7XlC4SzJ9kmePR0dJ9kSeBXwT2u/VlUF1Hq2W5JkRZIVa9asmY5SJEmSOmOj4SzJxUmelmQH4CrgU0k+Mg19vxy4qqrubJbvTDKv6XMesHpdG1XV0qqaqKqJOXPmTEMZkiRJ3THIyNl2VXUfcBTwmap6PnDoNPR9HI9NaQKcByxuni8Gzp2GPiRJkmaUQcLZrGYk6xjg/OnoNMm2wEuBc/qaTwZemuQmeuHv5OnoS5IkaSaZNcA67we+Cnyzqq5I8mvATVPptKoeAHZcq+0n9M7elCRJGlsbDGdJtgB2rapfngRQVbcArxl2YZIkSeNog9OaVfUovWPDJEmSNAKDTGtekuRvgS8AD0w2VtVVQ6tKkiRpTA0SzvZtHj/Q11bAwdNfjiRJ0njbaDirqpeMohBJkiQNdhHauUlOS/LPzfI+SU4YfmmSJEnjZ5DrnJ1O71IaOzfLNwJvH1ZBkiRJ42yQcLZTVZ0J/AKgqh4BHh1qVZIkSWNqkHD2QJIdaW5EnmR/4N6hViVJkjSmBjlb8x307nu5R5JLgDnA0UOtSpIkaUwNcrbmVUl+G9gLCPC9qnp46JVJkiSNoUHO1jwReEpVXVdV1wJPSfJHwy9NkiRp/AxyzNkbq+qeyYWquht44/BKkiRJGl+DhLMtkmRyobkZ+pbDK0mSJGl8DXJCwAXAF5L8fbP8h02bJEmSptkg4exd9ALZm5vlC4F/GFpFkiRJY2yQszV/AXyy+ZEkSdIQDXK25oFJLkxyY5JbkvwgyS1T6TTJ7CRnJflukhuSHJBkh6afm5rH7afShyRJ0kw0yAkBpwEfAV4I/BYw0TxOxceAC6pqb+A5wA3AScDyqloILG+WJUmSxsogx5zdW1X/PF0dJtkOeDHweoCq+jnw8yRHAgc1qy0DLqZ3vJskSdLYGCScXZTkFOAc4KHJxqq6ahP73B1YA/xjkucAVwJvA+ZW1apmnTuAuevaOMkSYAnAggULNrEESZKkbhoknD2/eZzoayvg4Cn0+ZvAW6vq8iQfY60pzKqqJLWujatqKbAUYGJiYp3rSJIkzVSDnK35kmnu8zbgtqq6vFk+i144uzPJvKpalWQesHqa+5UkSeq8QUbOSHI48Cxgq8m2qvrApnRYVXck+VGSvarqe8AhwPXNz2Lg5Obx3E15f0mSpJlso+EsyanANsBL6F189mjgW1Ps963AZ5NsCdwCvIHemaNnJjkBuBU4Zop9SJIkzTiDjJy9oKqeneTqqnp/kr8CpnT2ZlWt5PHHsE06ZCrvK0mSNNMNcp2znzWPDybZGXgYmDe8kiRJksbXICNn5yeZDZwCXEXvTE3vrSlJkjQEg5yt+T+ap2cnOR/YqqruHW5ZkiRJ42nQszVfAOw2uX4SquozQ6xLkiRpLA1ytub/BPYAVgKPNs0FGM4kSZKm2SAjZxPAPlXl1fglSZKGbJCzNa8Fnj7sQiRJkrSBkbMkX6Y3fflU4Pok3+LxNz5/1fDLkyRJGi8bmtb88MiqkCRJErCBcFZV/zLKQiRJkjTAMWdJjkpyU5J7k9yX5P4k942iOEmSpHEzyNmaHwKOqKobhl2MJEnSuBvkbM07DWaSJEmjMcjI2YokXwC+xOPP1jxnaFVJkiSNqUHC2dOAB4Hf6WsrwHAmSZI0zQa58fkbRlGIJEmSBjtb85lJlie5tll+dpJ3D780SZKk8TPICQGfAv4UeBigqq4Gjp1Kp0l+mOSaJCuTrGjadkhyYXPZjguTbD+VPiRJkmaiQcLZNlX1rbXaHpmGvl9SVftW1USzfBKwvKoWAsubZUmSpLEySDj7cZI96J0EQJKjgVVDqOVIYFnzfBmwaAh9SJIkddogZ2ueCCwF9k5yO/AD4HVT7LeAryUp4O+raikwt6omQ98dwNx1bZhkCbAEYMGCBVMsQ5IkqVsGOVvzFuDQJNsCT6qq+6eh3xdW1e1JfhW4MMl31+qzmuC2rnqW0guLTExMrHMdSZKkmWqQaU0AquqBaQpmVNXtzeNq4IvAfsCdSeYBNI+rp6MvSZKkmWTgcDZdkmyb5KmTz+ld3PZa4DxgcbPaYuDcUdcmSZLUtkGOOZtuc4EvJpns/3NVdUGSK4Azk5wA3Aoc00JtkiRJrVpvOEty1IY23NR7azbHsD1nHe0/AQ7ZlPeUJEnaXGxo5OyIDbzmvTUlSZKGYL3hzHtqSpIkjd5Ax5wlORx4FrDVZFtVfWBYRUmSJI2rQW58firwe8BbgQC/CzxjyHVJkiSNpUEupfGCqjoeuLuq3g8cADxzuGVJkiSNp0HC2c+axweT7Aw8DMwbXkmSJEnja5Bjzs5PMhs4BbiK3pma/zDUqiRJksbUIOHsQ1X1EHB2kvPpnRTwH8MtS5IkaTwNMq156eSTqnqoqu7tb5MkSdL02dAdAp4O7AJsneS59M7UBHgasM0IapMkSRo7G5rWfBnwemA+8JG+9vuBPxtiTZIkSWNrQ3cIWAYsS/Kaqjp7hDVJkiSNrUGOOVue5CNJVjQ/f5Vku6FXJkmSNIYGCWen0ZvKPKb5uQ/4x2EWJUmSNK4GuZTGHlX1mr7l9ydZOayCJEmSxtlAdwhI8sLJhSQH8thdAyRJkjSNBhk5exPwmb7jzO4GFk+14yRbACuA26vqlUl2Bz4P7AhcCfx+Vf18qv1IkiTNJIOMnN1XVc8Bng08u6qeS+8YtKl6G3BD3/IHgY9W1Z70AuAJ09CHJEnSjDJIODsboKruq6r7mrazptJpkvnA4TT36EwS4OC+910GLJpKH5IkSTPRhu4QsDfwLGC7JEf1vfQ0evfXnIq/Bv4EeGqzvCNwT1U90izfRu/uBJIkSWNlQ8ec7QW8EpgNHNHXfj/wxk3tMMkrgdVVdWWSgzZh+yXAEoAFCxZsahmSJEmdtKE7BJwLnJvkgKqazhudHwi8Kskr6I3APQ34GDA7yaxm9Gw+cPt66loKLAWYmJioaaxLkiSpdRs95myagxlV9adVNb+qdgOOBb5eVa8DLgKOblZbDJw7nf1KkiTNBIOcEDAq7wLekeRmesegndZyPZIkSSM3yHXOhqaqLgYubp7fAuzXZj2SJEltG3jkLMn+SS5IcnESL3MhSZI0BBu6lMbTq+qOvqZ3AK8GAlwOfGnItUmSJI2dDU1rnprkKuBDVfUfwD30Dtj/BXDfBraTJEnSJlrvtGZVLQK+DZyf5Hjg7cCv0DtY32lNSZKkIdjgMWdV9WXgZcB2wBeBG6vq41W1ZhTFSZIkjZv1hrMkr0pyEXABcC3we8CRST6fZI9RFShJkjRONnTM2V/Qu7TF1sBXq2o/4J1JFgJ/Se8CspIkSZpGGwpn9wJHAdsAqycbq+omDGaSJElDsaFjzl5N7+D/WcBrR1OOJEnSeNvQjc9/DPzNCGuRJEkae126t6YkSdLYM5xJkiR1iOFMkiSpQwxnkiRJHWI4kyRJ6hDDmSRJUocYziRJkjpk5OEsyVZJvpXkO0muS/L+pn33JJcnuTnJF5JsOeraJEmS2tbGyNlDwMFV9RxgX+CwJPsDHwQ+WlV7AncDJ7RQmyRJUqtGHs6q56fN4pObnwIOBs5q2pcBi0ZdmyRJUttaOeYsyRZJVtK7ofqFwPeBe6rqkWaV24Bd2qhNkiSpTa2Es6p6tKr2BeYD+wF7D7ptkiVJViRZsWbNmqHVKEmS1IZWz9asqnuAi4ADgNlJJm/EPh+4fT3bLK2qiaqamDNnzogqlSRJGo02ztack2R283xr4KXADfRC2tHNaouBc0ddmyRJUttmbXyVaTcPWJZkC3rh8MyqOj/J9cDnk/wF8G3gtBZqkyRJatXIw1lVXQ08dx3tt9A7/kySJGlseYcASZKkDjGcSZIkdYjhTJIkqUMMZ5IkSR1iOJMkSeoQw5kkSVKHGM4kSZI6xHAmSZLUIYYzSZKkDjGcSZIkdYjhTJIkqUMMZ5IkSR1iOJMkSeoQw5kkSVKHGM4kSZI6xHAmSZLUIYYzSZKkDhl5OEuya5KLklyf5Lokb2vad0hyYZKbmsftR12bJElS29oYOXsEeGdV7QPsD5yYZB/gJGB5VS0EljfLkiRJY2Xk4ayqVlXVVc3z+4EbgF2AI4FlzWrLgEWjrk2SJKltrR5zlmQ34LnA5cDcqlrVvHQHMLelsiRJklrTWjhL8hTgbODtVXVf/2tVVUCtZ7slSVYkWbFmzZoRVCpJkjQ6rYSzJE+mF8w+W1XnNM13JpnXvD4PWL2ubatqaVVNVNXEnDlzRlOwJEnSiLRxtmaA04AbquojfS+dByxuni8Gzh11bZIkSW2b1UKfBwK/D1yTZGXT9mfAycCZSU4AbgWOaaE2SZKkVo08nFXVN4Gs5+VDRlmLNAqPnn9x2yVohnMfksaLdwiQJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mSpA4xnEmSJHWI4UySJKlDDGeSJEkdYjiTJEnqEMOZJElShxjOJEmSOsRwJkmS1CGGM0mSpA4xnEmSJHWI4UySJKlDDGeSJEkdYjiTJEnqkFbCWZJPJ1md5Nq+th2SXJjkpuZx+zZqkyRJalNbI2enA4et1XYSsLyqFgLLm2VJkqSx0ko4q6pvAHet1XwksKx5vgxYNNKiJEmSOqBLx5zNrapVzfM7gLltFiNJktSGLoWzX6qqAmpdryVZkmRFkhVr1qwZcWWSJEnD1aVwdmeSeQDN4+p1rVRVS6tqoqom5syZM9ICJUmShq1L4ew8YHHzfDFwbou1SJIktaKtS2mcAVwK7JXktiQnACcDL01yE3BosyxJkjRWZrXRaVUdt56XDhlpIZIkSR3TpWlNSZKksWc4kyRJ6hDDmSRJUocYziRJkjrEcCZJktQhhjNJkqQOMZxJkiR1iOFMkiSpQwxnkiRJHWI4kyRJ6hDDmSRJUocYziRJkjrEcCZJktQhhjNJkqQOMZxJkiR1iOFMkiSpQwxnkiRJHdKpcJbksCTfS3JzkpParkeSJGnUOhPOkmwBfAJ4ObAPcFySfdqtSpIkabQ6E86A/YCbq+qWqvo58HngyJZrkiRJGqkuhbNdgB/1Ld/WtEmSJI2NWW0X8EQlWQIsaRZ/muR7bdbTkp2AH7ddhGY09yFNlfuQpqJz+8/fHD+Sbp4xyEpdCme3A7v2Lc9v2h6nqpYCS0dVVBclWVFVE23XoZnLfUhT5T6kqXD/2bAuTWteASxMsnuSLYFjgfNarkmSJGmkOjNyVlWPJHkL8FVgC+DTVXVdy2VJkiSNVGfCGUBVfQX4Stt1zABjPa2raeE+pKlyH9JUuP9sQKqq7RokSZLU6NIxZ5IkSWPPcDYESX661vLrk/xt8/xNSa5JsjLJNyfvgpDkoCTn923zF0kuSPIrfW0fX/u9p1jnW5pbZVWSnabrfTV9pntfSvLZ5hZp1yb5dJInN+v8t+Z9VjavPZpkh1F+1nExDt8P7rebpyRzk3wuyS1JrkxyaZJXN7+7e5vfw3eTfLhvmzT75s1Jrk7ym21+hpnCcDZ6n6uq36iqfYEPAR9Ze4Uk7wYOBF5dVQ81bRPA9tNcyyXAocCt0/y+Go1N2Zc+C+wN/AawNfAHAFV1SlXt27zXnwL/UlV3jehz6DEz6vshyelJDnqC7+t+OwMlCfAl4BtV9WtV9Tx6V1WY36zyr83v4bnAK5Mc2LS/HFjY/CwBPjnaymcmw9mIVdV9fYvbAo876C/JO+ntzEdU1c+ati2AU4A/WWvd05N8MsllzV8yBzV/Vd6Q5PS+9X6n+QvnqiT/lOQpTS3frqofDuNzavg2ZV+qqq9UA/gWj32x9jsOOKPvfb7U/JV8XXoXgZ5sP6zZp76TZPn0fbLxNQ7fD+63M9bBwM+r6tTJhqq6tar+pn+l5ne2ksfu8HMk8Jnm13cZMDvJvLXfPMnvNqOf30nyjaZtiySnJLmiGXX7w6Y9Sf62GU39P0m+kuTo4XzsdnTqbM3NyNZJVvYt70DfNduSnAi8A9iS3g4/6UBgL+B5VdU/LfAW4LyqWtX74+VxtgcOAF7V9HEgvb8qr0iyL73bYL0bOLSqHkjyrqbvD0z5U2oUpntfmtzuycDvA29bq30b4DB6+9yk/1xVdyXZmt5+dTa9P+w+Bby4qn7gVNITMg7fD+63m59nAVdtbKUk29MbJftG07S+WzOuWmvT9wAvq6rbk8xu2k4A7q2q30pvCv+SJF+jNzq3F7APMBe4Hvj0Jn2qjnLkbDh+NjnU3gzzvqf/xar6RFXtAbyL3hfjpJuBAC+dbEiyM/C7wOP+Ounz5eavyWuAO6vqmqr6BXAdsBuwP70d+JLmy3IxA94+Qp0wbfvSWv6O3vTEv67VfgRwyVpTQ3+c5DvAZfTu4rGQ3n71jar6QVOHU0mDm9HfD0le1hxbtJJe6PuHZvnyYXzGtbjfdkSSTzSjXFc0TS9q/r1vB75aVXc8wbe8BDg9yRvpXesU4HeA45t97XJgR3q/xxcDZ1TVo1X178DXp/p5usaRs3Z9nsfPv98JvA5YnuSuqrqI3l8IewI3N38Vb5Pk5qras9nmoebxF33PJ5dnAY8CF1bVccP7GOqAQfYlAJK8F5gD/OE63udYHj81dBC9444OqKoHk1wMbDXt1WtdOvn9UFVfpXexcJrp0dOr6uIn9tF+yf125rgOeM3kQlWdmN6JIiuapn+tqlcm2R24LMmZVbWS9dyaMclfAoc377VvVb0pyfObtiuTPI9eUH9rs8/9UpJXDOkzdoYjZyOWZGHf4uHATf2vV9WNwFHA/0qyb1X976p6elXtVlW7AQ/2ffEO4jLgwCR7Nv1vm+SZU/sU6oInui812/wB8DLguGYEpf/9tgN+Gzi3r3k74O7mP7i96Y08QG+/enHzRYzTQ9NjHL4f3G9nrK8DWyV5c1/bNmuv1IxKnkxvVBR609nHN8eJ7U9vmnJVVf1538gqSfaoqsur6j3AGnqB7qvAm/PY2bnPTLItvSnT32uOSZsHvGQ4H7k9jpyN3luSHAo8DNxNbxrhcarqiiRvAM5L8pKq+v6mdlZVa5K8Hjgjj512/27gxiR/TO8g4qcDVyf5SlX9wab2pZF7wvsScCq9s+8ubUZazqmqyeOLXg18raoe6HuLC4A3JbkB+B69/9wm96slwDlJngSsZv1TURrcOHw/uN/OQFVVSRYBH03yJ/QC1AM8FsL6nQr81yS70bvrzyvoTVk/CLxhPV2c0gT3AMuB7wBX05t+vyq9X/waYBHwRXrHKl4P/Btw6dQ/Ybd4hwBJkjRjNdPr51fVWW3XMl2c1pQkSeoQR84kSZI6xJEzSZKkDjGcSZIkdYjhTJIkqUMMZ5IkSR1iOJMkSeoQw5kkSVKH/D+/DMGAeurWSwAAAABJRU5ErkJggg==\n",
+      "text/plain": [
+       "<Figure size 720x360 with 1 Axes>"
+      ]
+     },
+     "metadata": {
+      "needs_background": "light"
+     },
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "import pandas as pd\n",
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "\n",
+    "\n",
+    "raw_data = {'Cell': ['H3K4me1', 'H3K27ac', 'H3K4me1+H3K27ac',  'GRO-seq'],\n",
+    "        '1 method': [float((0)),float((0)), float((0)),float((3754))],\n",
+    "        '2 method': [float((153327)),float((12510)), float((12510)),float((0))],\n",
+    "        '3 method': [float((54616)),float((54556)), float((54767)),float((362))],\n",
+    "        '4 method': [float((317)),float((317)),float((317)), float((317))]}\n",
+    "\n",
+    "df = pd.DataFrame(raw_data, columns = ['Cell', '1 method', '2 method', '3 method', '4 method'])\n",
+    "\n",
+    "# Create a figure with a single subplot\n",
+    "f, ax = plt.subplots(1, figsize=(10,5))\n",
+    "\n",
+    "# Set bar width at 1\n",
+    "bar_width = 0.5\n",
+    "\n",
+    "# positions of the left bar-boundaries\n",
+    "bar_l = [i for i in range(len(df['1 method']))]\n",
+    "\n",
+    "# positions of the x-axis ticks (center of the bars as bar labels)\n",
+    "tick_pos = [i+(bar_width/2) for i in bar_l]\n",
+    "\n",
+    "# Create the total enhancers\n",
+    "totals = [i+j+k+l for i,j,k,l in zip(df['1 method'], df['2 method'], df['3 method'], df['4 method'])]\n",
+    "\n",
+    "# Create the percentage of the total unmarked enhancers value for each cell was\n",
+    "meth_1 = [i / float(j) * 100 for  i,j in zip(df['1 method'], totals)]\n",
+    "\n",
+    "# Create the percentage of the total H3K4me1 alone enhancers value for each cell\n",
+    "meth_2 = [i / float(j) * 100 for  i,j in zip(df['2 method'], totals)]\n",
+    "\n",
+    "# Create the percentage of the total H3K4me1+, H3K27ac+ enhancers value for each cell\n",
+    "meth_3 = [i / float(j) * 100 for  i,j in zip(df['3 method'], totals)]\n",
+    "\n",
+    "# Create the percentage of the total H3K4me1+, H3K27ac+ enhancers value for each cell\n",
+    "meth_4 = [i / float(j) * 100 for  i,j in zip(df['4 method'], totals)]\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "N = 4\n",
+    "ind = np.arange(N)    # the x locations for the groups\n",
+    "width = 0.5       # the width of the bars: can also be len(x) sequence\n",
+    "\n",
+    "p1 = plt.bar(ind, meth_1, width, color='#6DA2DB')\n",
+    "p2 = plt.bar(ind, meth_2, width,\n",
+    "             bottom=meth_1,color='#D2D5D4')\n",
+    "p3 = plt.bar(ind, meth_3, width,\n",
+    "              bottom=[i+j for i,j in zip(meth_1, meth_2)],color='#D52114')\n",
+    "p4 = plt.bar(ind, meth_4, width,\n",
+    "              bottom=[i+j+l for i,j,l in zip(meth_1, meth_2, meth_3)],color='orange')\n",
+    "\n",
+    "plt.ylabel('% total enhancers')\n",
+    "plt.title('Stage')\n",
+    "plt.xticks(ind, ('H3K4me1', 'H3K27ac', 'H3K4me1+H3K27ac', 'GR0-seq'))\n",
+    "plt.yticks(np.arange(0, 110, 10))\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "<Figure size 432x288 with 0 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "N = 4\n",
+    "ind = np.arange(N)    # the x locations for the groups\n",
+    "width = 0.5       # the width of the bars: can also be len(x) sequence\n",
+    "\n",
+    "p1 = plt.bar(ind, meth_1, width, color='#6DA2DB')\n",
+    "p2 = plt.bar(ind, meth_2, width,\n",
+    "             bottom=meth_1,color='#D2D5D4')\n",
+    "p3 = plt.bar(ind, meth_3, width,\n",
+    "              bottom=[i+j for i,j in zip(meth_1, meth_2)],color='#D52114')\n",
+    "p4 = plt.bar(ind, meth_3, width,\n",
+    "              bottom=[i+j for i,j in zip(meth_1, meth_2)],color='#D52114')\n",
+    "\n",
+    "plt.ylabel('% total enhancers')\n",
+    "plt.title('Stage')\n",
+    "plt.xticks(ind, ('H3K4me1', 'H3K27ac', 'H3K4me1+H3K27ac' 'GR0-seq'))\n",
+    "plt.yticks(np.arange(0, 110, 10))\n",
+    "\n",
+    "plt.savefig('Enhancer_percentages.png')\n",
+    "plt.clf()\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 2",
+   "language": "python",
+   "name": "python2"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 2
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython2",
+   "version": "2.7.13"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
-- 
GitLab