From 95eb020548d6af355fd5d01f24e3880ecddf8eff Mon Sep 17 00:00:00 2001 From: FlorianTrigodet Date: Wed, 14 Aug 2024 13:46:46 +0200 Subject: [PATCH] documentation anvi-script-find-misassemblies --- help/main/artifacts/bam-file/index.md | 2 +- .../images/anvi-script-find-misassemblies.png | Bin 0 -> 43685 bytes help/main/index.md | 34 +- .../anvi-script-find-misassemblies/index.md | 115 ++ .../network.json | 30 + network/network.json | 1097 +++++++++-------- 6 files changed, 745 insertions(+), 533 deletions(-) create mode 100644 help/main/images/anvi-script-find-misassemblies.png create mode 100644 help/main/programs/anvi-script-find-misassemblies/index.md create mode 100644 help/main/programs/anvi-script-find-misassemblies/network.json diff --git a/help/main/artifacts/bam-file/index.md b/help/main/artifacts/bam-file/index.md index ad84258..99a80df 100644 --- a/help/main/artifacts/bam-file/index.md +++ b/help/main/artifacts/bam-file/index.md @@ -29,7 +29,7 @@ A BAM-type anvi'o artifact. This artifact is typically provided **by the user** ## Required or used by -

[anvi-get-short-reads-from-bam](../../programs/anvi-get-short-reads-from-bam) [anvi-get-short-reads-mapping-to-a-gene](../../programs/anvi-get-short-reads-mapping-to-a-gene) [anvi-get-tlen-dist-from-bam](../../programs/anvi-get-tlen-dist-from-bam) [anvi-profile](../../programs/anvi-profile) [anvi-profile-blitz](../../programs/anvi-profile-blitz) [anvi-report-linkmers](../../programs/anvi-report-linkmers) [anvi-script-get-coverage-from-bam](../../programs/anvi-script-get-coverage-from-bam) [anvi-script-reformat-bam](../../programs/anvi-script-reformat-bam)

+

[anvi-get-short-reads-from-bam](../../programs/anvi-get-short-reads-from-bam) [anvi-get-short-reads-mapping-to-a-gene](../../programs/anvi-get-short-reads-mapping-to-a-gene) [anvi-get-tlen-dist-from-bam](../../programs/anvi-get-tlen-dist-from-bam) [anvi-profile](../../programs/anvi-profile) [anvi-profile-blitz](../../programs/anvi-profile-blitz) [anvi-report-linkmers](../../programs/anvi-report-linkmers) [anvi-script-find-misassemblies](../../programs/anvi-script-find-misassemblies) [anvi-script-get-coverage-from-bam](../../programs/anvi-script-get-coverage-from-bam) [anvi-script-reformat-bam](../../programs/anvi-script-reformat-bam)

## Description diff --git a/help/main/images/anvi-script-find-misassemblies.png b/help/main/images/anvi-script-find-misassemblies.png new file mode 100644 index 0000000000000000000000000000000000000000..aaba8ae515c8b16aed54bd0473f435712e2669dc GIT binary patch literal 43685 zcmeFZXIPZkwl0cot8GBK+kgZ?lpr9W1j(Qxl9LEX7LhDD7s+lFxXoI!GyoP(m| zoO2KmD00ql$E;Ge6KS4x9L@XhG z_dXHP(Ly4kpZb0|3ZJ~3Mvk`r0qiC*gW3OYS zM`Ukr&uV05Vy&xVp~q@&We_?qKu$z-iAdt^?FSBFvxANf4<-f=mgLr5PUwr`=xqgEBhmD(ypHp^wzrwhL_eMLyp$q;OQYhY%&p&k;r}`7 zjo^znw#z$L%7fW?4Y0>OeLrl>6j^EV4rec0`aU~+1>iZia zqLzm*G=99Z;QzmG{_h?T5&iG>VCb9m><;+(`57(@mFAhUv9axaAuYEU{R>%A#on)0 zhtU7!QnRw|_xASYnVOrsavp|{o>e_Ayo0X6zX~qiVg3GZqW}Ab=zn4l%vDH;HVdUB zCGV@M22YpiH-=1mP`h3qRMFB(7j#CH?&g>e z3;Fx`_25kVaujTAa@yM4{ey#zQmQ|FO0=EnVq#*tucj7~Ha|I;V3JX)eBM)(-)`1$ zS@2y~l81qreHk0`_xF$UWsvK`>FDU-TB1ct1Z`_?8e);lco83*za9-;wWgkyZ#rgmi|G>b6D`MU~xIg~*1J@Sk$;8Z@YtqeH zY_c}ntEi==wLR)VZL`0-#oPv4LbbfUZf4`NyX$;W$a(YC>R(=7O84($qdoj#pNi}k zSS9-MOhX6D%m#~UhljO~ou-txw9HH!N>*a#GU<{|kc;e^Zci{>OH90)Yuc}D&>YdZ zzqb>PYt1!LIr8(*oJmw~-@cWUl}*pDUZ3k%QdHEPLOulX*$=#lH6T3Pl=OsV4eZmX zBdpr)koeuZFKCmJlk3u1Yu7#NW@n!kYm+n3(J@}TCjIi|OS{RRj~$c3N=iw|I6Y8N z;Vu%pefy6=-rKit#|wU}th~7X9xa)ck%2)BuFk_~XoS%w6%_Q$YnYpx=b8FMd& ztyNUq>hGCZFVYZ9EkB!{A!cTlt{L=r&T@TLRfI&s3)fm`ZLqbq#XGC8Vfj$CGNowP zAR^)*oG8rSUypQvmVrSEyS%aj^VLgHFE?)r=QprtF~B_0*XQfJA|Rl#w6x?-)n(0Q z;Dg2fo`BKV+-ZX(3FxA z0d+xo`g_Jixb1ZQb15 zxXgz%V3Q=}}7}z3=#wA8Z2{29Y20w7Gr)%Kf zj66IK9UUFvZC0m;ahD)Cd-nnz< z?3pt$ZEdjGdc2)6di4Qd-*!bFN{=2zefjd`0Yr=T1i5}(L_~yfXR<%-G&p~9Qc~ue zN$b0p;QD^xc<^g=K_(_i@I9t$*B-dKmPO-NI@0zVivJai;5E4?^;7w5rlS7(>lcVI zC=|W#jCuF2qr}k~7TnuOmjr=gd2z7?oMG2;0#mawZ3*V&vV6KNF$%+Wvc|I z7cX7}ts4%PmDXeK-n+*)tJ4W1z*45c>6utr^E#c@XOkd4=skHdHPsrEhXLP&Fr%cQ z@zlOCD(Vs)J$*9lXA8uf#nDPTnil^l%c&NcVk$w*s=S58dv*>E+wq#;im7hgxN+L) z;kEN3KAW2kj!sVegOxD8tfi$T6AMdjaeRJ$KFfGtv3&+?s;TyqC$X>@HnTl9@PTDv zpfLO6$B%x#zA-L4>tqq%MS zi+p}$nm?{EtD2jyU2EcBn2J2uO{m;+b#*1b`8236$0!C?8O~=bNpbV3JS=^BY1r+c z)ODkHxYJ4I$Im%g%+6U|6n51h5!#HJ7J`)27{;sUBe))P`O1~KbtM&*ndc-z^2Wx= z)5|ZwgSuCz+aZM|FteNe7(+wLELUT6pwJqfc%qY|qfu{`9%g$wf%s6M?E%=++FG|Z z^Ua%G@5|%k-Gl#uc+EbU>RZ9#-CHPG@t!Ud-kX-&o?vt99pUQBHA(2!b{kyn@e$e< zhx?i@l0)iI%Qp+}%Fxb$T}HukXlUpaq?2mzi?%C`+$Jqiw;}LYJ9x!O1=c;z?CFG9 z%jLYeh;cF>Dp7@(U2ag&1!8i}no(X466&s3=a_T!>d%80CPhWv%FN6(?#+rbE1ISs zblsUdaq3iR{sNq|EizR#&+{xD8QQ`Uxiq6Nw_Jy%$0&7!hhn}48%R^jiHg3wX+D(N zEVLuNvp&~e>&u|@;6W@Hzo?iPGsHc}oh=oA9Og1?x#Wo9i>%z=9P`;8c}|j}TYWlQ z$VCm7A}KE3^4D)y{IJ+}Q&Wbp^Cqy3?$&?|jLF}0UEG1Jor763j*E+9WM;1K8VfN= zZ;h4k?I?A!PpCc#rxBC{L85&lHYwMCws<>odjeutNkhSi-nsZ z;Cq(F7rO`-W>rrV#1}gYFKyeeZ4~rT$tliS(pAsRH?P#Uj(yZ^TvL$@js4xq2xtz1(OP64)Gai%LxrkVNLwT^Z<;WXqolTG0qvj_V- z3%Ra)n8r|Ux#P!=V|JISxomL{AasyTJ{C25U zR#vmPZ2d-d;zI9>{PC*U`bM=_I(G(nfHbjv8Pi~dTn0^)f~zeerw_-yJV}cCOsbue zg3b-soe>FG$hKS4e%h0*?0tV^M1)zsrlMjD04Eu@y%}v=TnA&jkZrC}`@NN~jqPdq zLuw0_kcTm=tzv6k>gHrrGjQy=u+@88&2C+P9y!VR_|q?gmem*CO`Lb;@^OR3`O|AA z#>T6#Jd^d=UOva=Se^3RYT^AQlC`;hm4of>@`u~M!-xyY03JDJ7Ot8XEXftaV2Ti4OTfG)=aK~TbEhhG+E4~4(k=*LloJd! zT-Z%(raPlNmddn0PbrX9^ORufPe1*Hc0={r`82iMIS>n3V zLS-X`(k}Ab_26h|Xvo>0{8d+950H(omeT+n2su*&Hg zM)pl3AD*8hCr^e%)HmOCNwe5a38v^55Dr?8E>aRirHr;uTKcu!{D9F$aVFSa3b-<#>uv$qv2Kdk zy?Z!!3Nv}w0HM8A#bOgLDnT^}a%$zSrFC_6d;R515df;o#;XXcQ9;85+3*fnOM;UF24f5*R$!bM#G8(M^I?pW^Kvaia&j{JY7eR^X#L$W4~Lmg zk+FZ$VSVOaC{v@fA#9z(s2E=}JY25Oje)Zwrg=^4!4cBd*2bydaG{vHW{np>B9+6) zp|ppYcRWeil;N&}`^V4FUw}+BOyor^JhU2u%$;Tg`-HsbY0vxHV^v-pnthq+*#^x_ zTwF>iDj8kLN-1Z5{k3I;R~e>z|6qSlraXW&IW@K8)8A)g4CFf}0d0|U=)WFv7$Y9E zA9jha246I+c|#M=_LLOtvc14^?38d0XKYD{#>Sw1tSwH)y? z^wNe8*uH1-5rO<)3Ss3rjw|}@E=|`IAF4zOxpXARQ73rI$jNm>0-M`^+@=aap}Xeo zl{E;uJx46TfI453am?<&lT6-AxgbWPGFTq;*7%?;W`5ZLz(-4>Ld>advzmlR*JS?r z0+ISaHtuj2Wfd72FNhZDM-Q_~gz#F&S3YB=rM(N04ItmmJ@W33j;RmFs5{`t;KT+! z+&Uj$z@;h4$xJd40=-GPbhQAPb~l$u&%GXW*%$;=kFqDZ?GIMN!h3~Vz$zM*s504k z%WG?;@9IB7-~w}P;&hgPY#X`1%D|#k^3c3|FE3s$G9A)&rDcDfX}h`W3PrIRL<)#S z#yoZazUWC#J|U7ee=U7=)g~@B_QFrd5M{#ol+f=3ST&PLwH$6k-U%PYpwpYl zC^zdXPiHINZ8PgVA_i4}873p!Zmw?{lB5z;Ag3?##X{+b;HilFsT-8)fH~=;0xsYH zY_MOWyWGjzeBEq-hee}6p4w&g?f(8od5U_z8P?dqhIgixA(E50pvZP69`;pDDv-4k z!cE38YUh?kb8Q8M=!31ugByEgN4}qJU=xeJm1e>G{LA0ZoS7IKv;Mrzs#){^R?E?0 zUbdz0vN`HuH}hl&5~rk$OeY-gT$j0jEPyagfMAJGV4I}K1YQ5|UMSmEaHr_D+~8VW z8qQtAZGR<_m_`-?X|7J?(c&R2r^~j@9!{m4J#}V#XZS$Ne#kN5j)(~H9;=N`B;YS} zSh(Kogm>@WLE6_n*$SyZpL_sWBgFdL6x-Pz1=xW5KoSs}kWvd`I?RgWQsw?x7Ij#ZT!BEQ0`O&j$;a&~9i8Zvd#_u|ezIS~_w1^v#ZDsM_;OU3Tn>Dcsa#B83KjwUEF<8KmX1rbNfUz`s@y;kNHs zOra&2%TW-4!7;W^P&%#Q;)eIv^WC?h__)~A^rPq?ZDzT{s#4R~$QAzC_F!)#^r(h1 z&{CYGB<+A{=kkEx;IbMQue8+P7-CQ94F9|aMs^{9@hMO^MkFccS=6TDL(SMbNZGVz z7sF1{(b5_MLIuaOaA?pz*jBm7YxQYfw>d)4{ZflDfD&1jrO7EOQx5$G3kPptYUaxi z_K#+;1iLj-LuB+A3y6gvdchdDB3YK19~?P(HA?{VYj5wHX{6`RYh^}2PJH|M@r{mS z#N-_#g~(B^i8^d9>cDX`2?-4;hhEqDQzxKnX{nlmwA?lK0NT=xJ=AU_1tGXgY-j4O=K9Ew7BHF5F3 z`tr@S7aVN^#W}43Gy?J_!yO}`db&76Mnz?|Sq3Fff;%xq*CXxk_5>01bRob`B?Sc< zZSCRuAr20X_MC{Ipq`OJZ|~Zgn)fFQ@A*)rLEr#1)^{-}EPo;rc2BqXYv%zG5X4BxTkZKKAH&b;} znQ3UEKdgjsn-7iLfULouCcHD7O-$18smx`E6({gxqCx7{RC==Xq}q4Pq1H8qQ}Zet z0>GSDQzKqkEyY&)gr8nmczAWd{m;m4?<9%_ptNQH9dl9IkIu4Y?^t$fCd*qn4`+Eo zq0=#Pqv?g@ebYxlP`m&YG-%x{7$WReW;fe2dVK^a2E)Y={f$EQVijdg1qB6Ssg$T`M$vj4lGWGOf1X`P$zp40*Qjn*C>k9djVDKpjg51xCM1Dx zNJ>p@nBN50lmaNueI}?;p8Uq+7w6fuG!`CSJWfoUW?!C>k)g<|Qd9X}r4Gp#K)Q94LMhh89gB=n`x7E7atH{h$!Hb$j0p(4UtOo(e)o z-O+4%C%t%~OxlG+BC!(xnvfq2EZDiW!%Ra1DO0`A(0pYe6h=g z#F$W7&P)SP7G?Pks9lo*HrMv|@AhZJ#>Rdyv7hN`l+S*64FmCKZLyMsUtj^k_tV8b zs9tzQ2~SHelUMwqCkaMFIQKHrty$6Rx=m zVQMPaAmYX#x8w4Y;yK`zt}-#zSx~2a_9Uh7E|o`Y{O_;02P~dLtts6c3=xB7eEW+h z+bOAEF%AwqrwZLqk-dZL{UP0Mq~Z{u`vy*p?7GIrPAL9njiJ(E*Ze&*^E8!!A!Z@9 zqvL_cAr)oifTl>{;gMWeBNOysD(^p&6WWSJTLiF8+1>j5wTDr}6gN2?<^WbeR96q> zHqT4RfJ$X!$jQ8P(iKpq`y15Nw9{+P}~7f>TzL#Y>)zUDevG9+Los}?M^ztfzp*G{^eCp>h%yV6U49} zI(Bbo1L4cgYy&xpi;58yop4oq!U(`R7}*n}tm5 zsnfvz5FT?o*cm47HUh@d=#f5PotS>}^1)G8X0hyo*gicHx1IiBa@+yHskIdCJxN51 zX+r#;Bj+;q8*w$ zl_h-X896y0KnY#+HN=F|d`O@@!g~Nvj(e{P;xl-yVznG6&K2u%NJ&X?x$Qd>U$|Z| z?{o(XgpU4p=?{d)>ZSe?paZ)tTiNEj$|yhuyk+~_CVP{@2dx0^4ana-6a?Dn^x=4v ziMJ*rxjX9X>Nu^xO2a0lKRmD%uEUt{!+kJuL-5Q2OLW^Pr8#r#-VRG zQt^lT_VPr%?Obj*_?jfD?16$!FW<+cspmIY%=*d+Irh%;cg9JPii(Po;h2?EWnI0n{h>e_u~9Q zdxDQz-c!{DE=VqG&`(MRlpMdava({kdfP%}{ks-mN+F!zzRBjmk5Q2mnrTO?ao3;? z)84CxMdX~|(x)>ae74E@O<_r}$G!S}!7K(3A3ik1gb(k2<;IM@Bt6A3-yj{v(`pkR zeg!Jb4(5!K$9(FkMMY{1auLvFOOQ7B3v!v!*7C#&Vq)7O>hG(yh){Rez;e3oIv_GJ z=AM`FY9Zj@mbN%46z6IyD`P&4MLNy%vOIhrrN9w}_4QqYfE_Z_U;7$bF;!j^lNN3w z*#)W4HCXLb(VhlZu)t+XEpiM>aJ5L6^43 zs(h#yR&Ak@Ku3uxWVX^(?9Hb=QDI?}m6iaQ5pjsvdFUjs?`Ft2XXcgv@mUF((=(z5b-;*^y+h<16i4iI;lZHhJ6VT6EGA0DywzO z8CQjcwE=ViDf|G+gTZoJ1Ua_WjsyYe)A8Z?u{B&uhA~i-s1FG((>ZV8WfaO?_iUO4 zHw>4@YvZ$9zaM(V#8?+iU?)&ATn>wmpdN*rO9uxB*EKaM0~m$I$vNl_K`WE->eVUi z_KTMyrWQ+p<0UtSV#NT@))AM@h?{X{89!|*CC)fZ2-(r z;Gld&RgFyZ7NCqZAWO+{(AO)O9s5xZsdOoE`JNPJ19q}wGSbaV1-OMaV1tusMBEL> zYu>(s9z3KvWMYsKFU%jaYvk%ikX(2)@w?$yOx&tDlFwz_akVk%h6u7q%|W0g01WjN zSn8HeR##Wg#P|q((Ng<9(E^$P<`)PHM1Wgfp{9h?Lxuxe?Vlq1QMyDKT^|3a#6UxH8+<993TdtIuYx}xrvLi= z4Uw70n;(1rpMFyv6CWR6)%l#rwF%n6N?KZGMQk^3DuYAC?*OZp3ftS02MMDyRh0wu z17uA`(C9>QAF?a3dg_9Fd}>g_)ipP(fLWwbfjPf_e;o=aCa6%9H&#}F(rgf2y9(@i zq4ne`JSq9-pFbr~c_@LmB9r*qKnCci2ep5}i|&J;HAM)Tq#QeXl!=`k^P z$}TFv_=vo6_5z+8wSw#FB%yBW$5ByH8N-Tb8LwOsIq&&7F9RAHz(>JGL;!{~Pw{y9 zk`cHsIyyR_G)pcELXZQ>w=WhWejgZgBNAq@W*=b0UwU}V>~5_9v6tQH;_Q5uf+7WW z2kPKH%@oI#ugj~e4*-<{e|Jjo#L1J;=X^Y+u;FOta{6Yx7@Rjy_hIuoE|eezAtNJ; z7eu{KRaK*@(b?`g+~bJnqobpUa}W?tZTakWMy!Z~Y1Vq;PLYB($&7W2f=^MiwOEEs7nPS|1V6^B zG^kof&t=bgbvJQ`GN}{Kd9TJeY{?4Iau;mG>}IZdR6% zFrXB4eCmS$bAgW!q($vA(*0lbTn2&4fE5I%nYlzn3OXqRUJ*7TnK4n+#XO9hv~_x@ zT^;U1Nm^FbHK)_r;v_G&K+>ocntBk(gMoQab+pcCyBuiq!kfJ6&n4d6c^>{83)UA) z%{bQu6Q4K_7W7bgYlX6Ywd7pq%Db@;mquCIY`a=yIl5BW+f8pmVxS$$I=M68>9}Vy zWm|$SVC~)!F^G8HR91DfBc`l7KE9(8cYo+PAI8*bop-7MW;Q&j)i4tovoYhXb%5JZ zgPWVo!XCOdk0pwUELdxHR<`b}Rh6l{xB9v_8dui5$LOTdu;3WB=c42yNLd<~XSKTD z5*ppO`N3}7a6hh+T@(3kB{b;@i& zjU~mbl3s!p`5t66ZiCW51sVb=9;11`6cTa!LnjV$S^?(5H;(z-H~I*C~(kGfg+(8zqSgwhM8*TBVOvvN)M#cDr05N#6HqiA;Cwk5ZXJe6`xU) z38E(NfE!RrAPVY|Z-#OYNTiIB8o1TSi8K|`tGJ1rKCjyF*h!B9yIW`iEq5~sDiltd zYO~`ijtfD3UiKRj^(2MgzyCh{FVoS7N$8X5KW5lDeVuGe%4WTY8iGWeq6Cp=%clqT zc5qX#m2)XXC)!qt?{hIXqTQ@``EA#ZGfh!Mm#BX+<}Ry!qct?qcL;^|TW3US@CAcR z$KN~Z$(I=u?R4k!+nSzspT-ktA-@vzGKTN2xkm5_WP2Gq3&b&ba+;PtxIc>1!4KgJ zDGN527Uq2esIBInp3)q*Ovvk&QUIV#JuR#F3*+F)p^dqc`zOjP1N{>t%8be`nUa zgc+E-)_KBHhXojw@$G7Y}3Q#f-JZ)`t)dU?D$ zTTIyr~zMG8e|__&|+Tn9b;OMy%P5bgJ{{1_h@G z@*z&kh{m0*#B^l4z(yL~1Ft`W|KB<@7O0t?XWcfgVDxskBvv$E zoaIkV!E5Lm!~5SYW=TLLZd{n^jlHG+IL-7+n)?|%-))}q+nF*Yn{=!PUpc{JS!Tz- zM%q?zI{o4qGzzD}H@mz`TV|)B(|^L4+su@?5jXeU?!M094iOdb0KJ52n6>Nh8;jr< zX?$rtO1-%hn!2%p9-XWzcq6z?%Zcf5h5J!HGgbKCnLQ)yhpVJ^yp)`);?9;|IyO*cZW25=cP#ex?@6NXXR7er zmf82Mdv(3Mj|45?*7I-YE(!^QmV`}(1lh|o&&|!vn#$*)OIp_mTVU}P+T3+tDoJwS z>xyP2s5}Vzp8^&+kZ=BY8GgGS|H1#?`q)jzxrQtkoPv?%=q)sXic*3#R>&o#_N^y* z2x#!&A(@3a;nWrtjn5nItk;g5~Y{^@;Wfr7h7<;LXhwjs?FO^q>51aRmRt z|CevZklnzmq8W)bjh(mpMd01YWL>*&$xbbDcV;Rt` z*KFeg{E51xz+V|ca|Tdjk#vl;IS_4akW1%`p;~Xlb*5-?UPgCDFN>sK3ezuhadGuQ zTL2pDNd}QF%7|s+dQ?FKj;JCIk~-N4foYr$Ts^{mT);tjPE!zX`W9DSXs9GXLIcK7 znu6PmcB)7cN<#p{1P3aRhN@Ci5Fc19n9(x?>0yA;kSK-Eem*S41D!DbWuXcUEn)%C zn}HUJqXQ^Cc6N3o7l7&_8Bn-|w7!147Br^>x2Ix+Kg;6}2mR9Z_ZJ#wBZzAQ+0=|6 zG>OYvp}!15U4x$dC$+y5+y`B0p`=szG3D-XnakR^uUyD{KX{a=W_E1s_^Sa@K!)N_ zPzJ7rz~lq>LP)_K?}Pkx`YA8^EqX1sUl}iYH2?)n{90BPJIGARKCI-M_QxYa1?nr` zpAH{}lJ&B%urN?acR=CMilvhRS>#aS4Iqt7Ko*Q@ePBZ-=H}8C7cGEXVFanb-+%vI z1T@fEkFJjq>Sq+K1H)zmR5sFHz-}nY%fErHBjQ^?Fp{UxT_o@5SO{cVte@Z6T=QWq z6x>KiRG{L8ngPjJK}t#Ou;@@Y{LeduWWDyAUxRvf!58RHK3V{g9Z!I8oqjUVs5xx8%SV|Fkz5eA~yv1Y>9dLV`g70 z<@hgA$O#%z6V| z*#RpxFAzh`D#W3Rlou56AleV9K#(xQ1=Mj!=>&a)N81h!KNTZcTS`V5dxxN)AQv>% zU4YksZddoluv-S09uqC?r`10FfNQFcA1b$k5(+vl-N0CgvKh|9=uX%=N9CreIGPxCn}M#+Zp=wRldf0qEKK!&4s zUi)0s;5Dhv3fOo1z$Y#juy zDL}f~ae*uiQGQ5o2h%stUjTB--6ZIxI1(2X*(8}2&)rmK0qy}hW1S%R=qPt{Ln0~* zElp7KU++MgA}~miQtfhAq_0f?4y&(I_+WSC=nh(t_hwk(2uY3ixk9r`==})&FNFgY z>gJ`LQ`Sh$2BbXdNyhI)Lc3gK4j9C|lLf3u*)r&~PzvU-ucxA&WB69H2GsvCPzF#m zf!e5`GP5iE99=DRk`~H#WSutVxhT2y>dqEFo4JD~3e6axrg87{l@#7C1^xnJe6wwxxfaYivhk7-{oP$wWDTeQy$c6c@idUc>~|+Vq%@a2GN6z@s8z zUSCMW5FYMTX-p@PWLquG_2-ceJCG!{8p*LXkdBxRwIQ5PmDm*LOI3UWle4Z+BpLm>Jh}2_)$eD`$15X4lxkY zWU&iy8rEMM6rMg!1s#uBJ{t&RbDFEE-JBK1py|SqWsD(vPP!l{uE$hkOw9`m&i=?Lzp|>_S zIeGBK1lm0d$6u`q{{V-32qens4G^t7fF@r?+3vDbWf?<+!yWVv_?M@-V2t&AWgrvj zHw9^xn*E?{AP5Gtk2bFXdjpLa;8SM=({pk<
    8%x9=yxqLZ{Zyw?U0;>2X3S)w`16ktHz~eSK0=g~$2>ipso~q+4A#EdK>D1g6=l;e!yRs>QG*WEtXL&OE4~_%!2z~)2fLHP@Bq{cO+$xhFnkRw zHe@E~)p8wVSfT*ho~kv29HhStdWL>)-b6v`y|)uIw#jfCo)T9O$Dv*ZaML};sC@w< zLnu?CAWtD`(qhOoe~=zj4j$doQc~7E?MStGJ*|ATP3m-|pacf`!XR$ZoC3BGxb*7{ zSNn5}eAA2P`6|K2moU)vkT&?l)$xPL)u3tsLP5AV zaIqRVJZjQIH^T%pYp6H>_Zfy9{YG*xO5V0#FE~0lOf6MY=bI8J-L@8!4}TunhEIss zfaaM7r89@iHK2&x`ye$MV`iZ(huqk7xJh+EE=eH24Nf&*w`#qy5?U*{|}fI3^`hb$P~~F;h*I2gzb8nMHlC z&JSkyzsjr5E6}We#oCEceZi!ed^wlt*H|UByeg~gTjH1g__cOVC^%SJf8D27CV==C ziHO>bEm^G$YpbS(jRs$8YU|UJn~U&oN(=igaM@*V`;xI|8{Fi?!_@>Q5)nivN(P?sr3d9YZ8v28Ty&+-QD&T zJ45Sv`4bRDd8b$=Bb>}YGr653teyE-6hOF5#S8$n1`to2-i|y*bju`Uv!W?%nSaA8 z%z2C(Yfb8R5fW$1mw{TYs%sA|MvjfLsOOE7mT!JOlo1Ic!!3vvyE(&wYxQBE6DY(b zr<``!I}^Ei7j$426*KVYWg)2hSNt=a>EKC#Ip3eUZ&XuSSL{IHvF_(1v|C@kTc5Bu zpP9f8Pm7cwonX~QxJ{eHR$q77x}V#5h`Y7+Wh>#%pw({JDBCDaIHe&KaIJPyXrs`e zl;BD{0$y#04bws2DL4$3ABn@>G2zV*ib9VUKy1?VPSkVx4-)4~5cIn_w8~Rbn5lP1 zj;(mOZF-P#n++Ti7KUE4H9RWtwGjlU970U*i;MXqcgGnx;Bf&Bn^sY|ZB=++O*!@9 zTX3l*iHznDF6R6}JN2QRp^dNV!`oeKrTd`yvj!oRz(ly7+=Q38v#wd`T9<}ZEuGw! z!{>yyp`i%w@~y_LkECU9fZ7*I8E(G@0zHQ=`^`f{EqC%rVojaguumN3$74;$C-S1# zB<7j(VR2FHmsLCsImsQHorYHYg!zN`(}tD<@}EEiJXW$_|7W|D_IBGf^}K7)bk;4~ z8C>`(CcLfCm9AO5;2j4tzj1(q&7fI`1c@0Bz`&h2iY8EW4S-@~3u0v?g?Y&hNRkU> z3%MaTbni6@tTdowbcDy1!T}?us^`x$9^`|1X8l=>KP@2Qt(o%unE~gO#w|c^5fD(m z{>>0k16>k0?l~Csf}tawke_cq?4;iu5e~`*F60$MtMS}=;8*3V&AK@``S}T@gN;B0 zb{2k}A1B$@8Fy&h!KQn^9h+>DPIqdoSt84h+}?n){Jqb@_!lEcg;&DE!>4gLzF8B% zNi3kecBm^Z5799+yb66Mxt&sdkCcoIHb8()5S&7(0#2{-1A(5{5y|3&L(EOZZA%7H zITPq6%M^DR3fEf>wzhHb6!7ZV z#34AryCA9Nc{K-*y{x`S=Gv%#32&KFwMzLwZ^0&rZp!oOZ+NPU5(~Xr9_W3=8K8hJ z#`yLm(K8Wt9-bf%D$G@Iz_(uT?w-kWB)$-xz;{b9#%wa*O{au9Y=tv%%8 z&gp*l6Z7zGVi^A0-28lls@G(32DeO%XAdsiq_jtmUQ#llWyGr&Sm+cL7fUri>stMy zyK(Wf&&pev{D}P4T^L!hF;Z9wsWb(Kr7#fk77cD7M) zFLM4Fe9G{AsYmg80QHOzJY7@|`D!E*o*_RHSH@aFZ;t!r! z_X@&yVu$djY0(%Ow!1o4o9Twh6bXtZm^jt@4=VWk}kxSIo%-wf>Sd}n$6rqS-$B!k@YVWPZW?z|ea?jpWtZy$YD_>=T#p4J@I!9McrxNB% zYxI$Ofq30#*3!4M^5^>Jq&FikF3p5#oE~6 z_fVeBhZ&hWHgAo=<#}yL_}>ZVV0=n;badzrI0w`jy+s4S$jY_Cc~retI+i||9U-)`d^t1n#8|ZSUZGGh-8L9Du$i1 z-!Arm8|4^tJ_)HDm4g~}L_DM~*oF(MORqm)QYc{{HZ2O5D=n)Wh5hi3?eC`DdtQq( ztsTx_ECt0*lp0~yR&(XF+-a1&u8r@;2af1%z;7NjjGaKgc_EfZcQ9RjrMgMQV3kgJvnk`}(uCV;J-*zjFPJHbN{HL88II z!8n5l_4q;W^(P$MP|ahN$l7{$q}M`H1t0y^&fq$W@ zcM+DhJ538AYO_t0EUfKnDZqk94R{L-nQlJ_x~3jT>-c4uGW4Q>G>ZVVXdXTN9jzv0 zg1U=7@?eryV&Ti*j=QESguCJ)(lh;@j=+)?bKlZRGV@)w2+>^Bt>p4xioI|p3g#mggurDn`f=PdL$(jQuKApIY57=tkI=(^iHzpgWLaGl zu(2zoE~>l7>f7@S5($x20v}o9hg~d}x6k;}`pyqV98pPQB1AUAij!Ds{U_-S8YQp| z+aIv#4UxFs@8DA4t=^~2Sq0wvG zW9O!|0Gpw>c0+_<0L=>$LI-uREJnjzL1)<9u6xZrk5fTMe;jYhh8Qs6Hc~nXKDXAu z&u;GpOr6~E**P&1A6WZbH$K-=$Vy_Ri|n1wYvd-+!k^cAb^~E2q6?;`5<)vBwbR+Y zwArrk;97Zc`Yb(sUB@K}z6BuGh8qi$)3SA0{*>NOj#1*1osk$;<&a{7-9KUUPhQAt z_js9flA(MlEhk)v-?~}&O7H+afZy`R!-Em=B1^~@egnCzThnLax?%JOQfLl5_XC7N zcDDwvRI>&7mN;wS6W+O%f;%vcD?07-mv>@vhE&7eg9~ygs<*Hhb1sC zFRzUw`I|S6(w3KUs^&9GEPIj&bNehj3?ERm0|HXM3tCAb^Yim?=>#4z#mv=j&=0%O zO0+o*Svu~j(4RxuTV5EQ*CS5B0^Y4pa#&iCYT26QjiztLjW@J=^%X<%$n^uCg=>AE_Vyxv?wzjHhYtN82)_YING2%UAmhj_BQy+Tm5`KxSo$*@7 z-NtJCU9aa6EHux^ewo@##@#o1-aB=7>fI^sKFJJpM2ZA|C6`@}j;k%wl)!dUJIHQ( zQ)t_PETQ58Epi{ucI^S*-I2GnwO<1pYvQ#-l?d)k__F?asBis(IQG=4&DT~E7AC@+ z2)3*9$b8wmT!yfq<<`$+VIBVr87R|BD{Hq7I8-y~B(m07th2`)74gRP-6PlRedspK zF7TR}*8>U)L_o*l>^G;(po(-}{d}cbq*8S?X@V@@@ zPrB#6JcfQLga0LmtsSj)=_S?cbQ0_b^*(xi;p&7rf5UyDH-rk;xoXFoo*42pZy}HC zR3`+3Z*g<7Nqm2|Q$kFg&u*dd%P_cHniJtLIPizjYBludxGtLn)+TMC`HLKjTz^C~ z0URDBn9$z(N$kIwIswV#3G*`W?iA40HiU9NX(BekW76-=%Eyn*(0X$%2umA$Fj37u z*TjO~la4e6^ygYC6m+UK$4y9SNr&Fd_W8ba%VK9v7Z7x54t|Is0}B>cS681tu7po3 z-hPS(blMI|%iardE+0C=Y`T(acNy`K;{iUX@nBu|cx$!8B1;~-vL|J)l7m+f0KZtD z$$5V*xSxFKdK;W5dW}57!z}pM{}N8XcmFPd{~wWw;SrywQbT)tZ?!i}ui=yQ1pZ!W zkF2$q!_QqQ@0w%!r;pJT5@0(l47Yc7vP}ZZM^6%6SLVz!Z}w}u0$(y+wz?e$&(*z! zhm8E8pt|KN?d+U@M>oIXLlIo_4EKKoNB_muutb?IEbjY$h_(_zbHJA_InS0%w$~#E z@-qjZo_QwuulF=ObeMyDNA9Z(D=L$s@YVUPlRf;71N*+JgA?vB%!50~Rp9j>ksEyU z-`Z;YP{6=F6KRLq6COGV#BY#|j2W{aI8YSRWzq-6b$*P&?0S-@80M@`C$4?%M#W3j ztzZYy@zxQ%11K-^qb4t;pC74UXU)2Yr3hEIHkI3Dy$^yV*&x{uf&AX=q>(uk6%G93t0jT3_}AKV2> z2h^L)K&(69XI}e>m|qlCRL@87K_6Z8zXHPe(UFC~>%X*$Ll^f%cI50%`%%C-$sflj z{3A9(t>p(#K6m|-CBV1@C>(WUwRw8L&FahQKYg3?EPCnUjvT!=^FJ&BP33vBEdk_^TP%Lb@i20bT9G-|fEU zrRl4H1fLmi%KoP-55sdJl~>AP9EW@O_{{t;k2To3vnE4bB~S~(OUp7cD)jnC6%&4M z`cOr5FZAW$uNi4$j#5Uo8a{#UAFx)5O`BZ(>=~h6y{F{i*i>_@!SFq z$8KRr6LXa4fx#E=@Hhq-^ucF51c!h9_gx0f8@>L+iM0I0m$%UCmZf7~NCV&c6ZYxW zWWa|#2A_qi|L!!^D(1JwHX_4rruH?VexRA3z+)_N1{H_;`uf=6@9m(SO31WU#O~lp!bJC)BmTJAS!n|U1A_m*ySF?) z=|kBZdHK23+PAmovk^-<^Lm;T&w)hc;9XC*yNr&Yxu;#$$vfwRX~; zG|NxJ@`+gF@aW)(o2%S%=MFY9O=^Fkfb>{21*a2&>P}%A=*x~@%|Ck{f-46Z^w*I_V z+ZSHzVJ=m+JQK$N!6%0nU(*}X`0Vv{tVY-m)YQr_`<9oK&t8II&f}p-PSC_ccxbdP z!|O0-AluHV*PVt~V0asE>Yjk7-cFR4UNSP5!fen64uqYWWf6o8h>F7}0~W){l?44g zii2j|lM7mi10($U2I1G3*)NXMf{x2IwXoQISDD`1Xov3NcSxbqD-RD8ZS84ADr}{FYTdyvt%{{ zt-ornY}Np_s#ReV=i?^~O+z@q;6)4b^RD`O?3`iJIO+7#>%pCx7^q2ZN#QGLHF4eg zTj#a+ANAQU$n^DN5BBTE{etnU%b~>IKs0B8tO-HP$?g5Oq(Vz=qE4RI7=vQ0v)gDjBL5ND&}zTK?p(>Q;IT-OFYN7c8`jd%WA zDHn#tz3uab;IOmt&Zs3zdyTv`1=^RA zcfZOUZ(7Y2|0awBo@2@daapQbZZc|Uggp6jnBQ7*`ZpN+95jRAfAIPbQ5!z_?+4O< z<=|R=DvFpIc3W-<_?L9Pz8}5EP`v&xs9<)^tly2jF#DgvKFgy=>zjKFp7H+wfDbUa z|9Nl{5fc*|)q|^SEzqgoI#hnJQ(n$yL07A3Kj_cq zGA(t46kP&NeAdrNQtL0XF!nwaUh zV?U?U_P*B<5!$MiJO2>z1k{8s)@}4BWLA(e&d6=Zkkblev=bdb{Qa2%cfjF;6D7^z z!-sk2a@v+ur~G6!+!fShnxFuS$a^DMK2RRA$Lc(LjZg%w)=t zu`*;RV*?qY$W&+$$xt#CWzI}Q$<$}gEQttVpZBZZT5GSp|5@u;``GXCJ379;-{7c4NU6r1$+l%Te%a z{{s&E-@_OFzyE3^Y;G7h@-HpGe<~;7YX9R1`@iuG|NUDL&`Uzl;P>w}s||h^x%FMR zaDlcF|L$}A1!weAlN^G4#96j4Hs<`7AIHD{vj5_f6($6$C#yAt4v<}*2qqi+{(&n$ zcN$pP*zAK$uoeU-Nz$S1sE(2teL&{rjl5Ba#TRJ zG@g_{leEGu0B!LHHTAEU5-|;KYD8X!)c2{syk1e@1LUNSeu>P^%-BHz-3KvM72E{U zoJQN+WUPc)B_Qy$A1h(*6Z04s6T0{{%gpVcP#(?RZng#Ii!iR3n3@`OnIi%jY|F4W zIp6qdXLJ@C(RdK~CU*p+KP1A^5DjTU+v@(Yk!u9L#kKF+e!U$2nmsO`BA}bK54ud= z49npxxV74sZZ;d1GRU-4hmu^`23X-^BEg!@%OQi17_SzB4tk11;&*ui1Q=5}eFw9l zKzulsi>FOQQ$K4Ua@yPuRrFM^LyC(Y?1wtdTt^AX3^iPFyP~H?8^UL9#GU&D; zJ~#IuG@dRh!!zZ_rGt51W_4{&`Q`CKH+#QXw1dX$OYIK>0|UQ0wqCdbYn=IH509iQ z2z|z{Gw1ApZ?r^*8M~%Efj|Z6-v}4dv03-Zx6s)ojfEs)MM8D6@S)A^n4ClI==*%{ zKp#2<)eW0gTdK**SvScI@3!O>3~Ms>4Tv!E^<_OJ^gMOA$DT&(E+b2EN?OXe5&d%+N(TnuHud@%`YNF1`eIL&C-;`zf`S~^2fIWs;(R1Y| z*PRdU+O-sYCxvkL67CJ~@$rP+Msk6vxw#LPPz)sP=AWN$ShHrn-Q)W$Z|@igdYmWv zgGzJYz{>&|Gp7w4J1FRIIC1)P8HBoq7cQ)&{Dv^!7WP8X-w!J(t;3bQdwNVL1rWzh zD=}qJI$*WONXw5S)fH<++;Nx%y+XbaHsy8Jp=}+Xxi_!MVX`hO8+%M&e>=o;egx}? zm%%xQo`r?81vUf)@S@nTZCe2JK)b}nS)i&efhv5(#*G^{h{1U`GyOV;d=ZrTV?Vzs zjd&) zt^{o?o&Eis{z>XdOz0^h_IhJ2uv5Zhf>oKu^DuoV9)m@uUlg&5aejoG_G{K2+C_aVPQ5-P67;1l{t%tAXZ(ucCEV0!P-YMW$;Ln z3$lTw>q7jb27r(2r=vTl-)ZKwi_-$L0`(!k`3hiFoSK;_BK^B#=?gc23m>TNm??K3 zJ`BWgg%}vOcDbHEzvue(>#{EWAD_HRORJCgR9CkFb{oa$b))HFZVyq2v##Z*IFEJF z!sHZaHAzEVsvGtqN?T=M9^1qsAG{BWS@;Ajp_B%0^FMrJxy9DLq_C){+wj2Iv3KtZ zQBl!Xsi|Cc>)VcDlVmxKZck25R(^VHZKiE7$4}H4yyc7+Ho(yvoe*|- zwlKc0n;%ZEh3AD?Zdw1x2%q#N$V1P_eoZ%94_(_eSX3QOP%}cyl~(%sQ$h1kwwIt~ zF%>q6h~h_bUMh2bgi8s$TbDzP=!s{l)D7>CLNr8?ICUmg%@pfXDO&OzRBq;8XYq;bM#NfmFw1(lAniZ`k4-HLrY6E)U98` zD!dXAsd0M70p2z@aQElXL*jG$Hg8&TLUy)nhE0h;r zR$vc6TflBgS09Eb=*%L{D3q3##?&I%z8Igy3-EJ?+gNyXYtEH)Yu@$S4}9e6n;!2o zokPEiDn~Cmfx28FApP9?K-ePKXBH93KPumG~D>|iYmijL+ZmPsm{n>p$w~~gpqy@H1J$i&p zQjj4sSKJQ@+8F>l@4nh+$g1kGW+rA0g|#v*LrvghGd?>7eLET^xJ5P5>E?gnDBO8D ze-4Y{^gGdxIHZVBqto7wT^-+kmh1)DQP!Fs>SMf03s0>e++)kHp~lASO9(m}n3$Fy zlD1#zX)PLzPVdrhLuz=Ql`UT5OYr8_tX+$7A-7-}m6@$(wpv0`va%*FKHeAs#^`l3 zpjBOlg=%!EmIf@m1aU?T!*y`-ZD=SXtFSr6*w|QBPOjqnVP$1{^wasnr;a!xr8wkx z5#whqi~%T(fY}=|^%d}tW=ug2(vdsWb^|+>Y`X5!wLRay(NOpe@^0hNcOj{?adhMt z?L%uBJw5#nmk%{H>#%lTWn_>*gG+nCh?g>5wToXrhvn1*;X)jZ{vU8vAGE8n!KkK7 zWY?};S#RGK;Rgi_@``cJ1U-1*hmPHo_xIDnS?KfCJ23VP7-@Y|rEhO%$KaB8(D}!q z$!YIfx8CEpYD_UW9;XnR(c5T(9UmcR>^+>R5XxOxQLze<5nXsC_{pHCDBc>_FES&V zl@ULx-t~wu0I`16Yj~_t+0aTg@){PwPD5B&3eTu4_c=)*Vj2pz(^0hP{6NFuV|ent z)#qN^1EZ%jo2kcDo8#{L6%O!G(SC%TIP}LIe=^TmTU(p-0g3Oek&u!aa@Yo&o|sXz zclZY+Xr?^;qHP3co)F53J$v@3!0V}H3|&tr6yj7;1jD=_9HI`KsL1f=bp`)lCXXn_SAoSEs6ntqqAwnx%s zGI=-$m<|vnlYP_!Va5W^cxg@8`c{8ZA#{x4x?BCFu6FVIZIWsWtM`(ibX9S)Cv|aG z>J|oHq(#JW_>Qp*+=1%hS9qc{A4PeMGi^7x23yW$FlVKBq6?f%De2OM(;Ha>|VQeZ7OB}kb^!4ro56^_+bNCmw)}}Jo-6NeKi~ruiv_LJXH{h z17jILB*TMEC%kG*81Lv)L{P zqWBVYT8i{01IwbTVZq`yjVC&YZLI-tfRHj%aeDklJHD(GT&Z~Dv>BNBb?M9SV*p9H zrN7~6P|*U+*vPVgFBK6HK@dN>-n@@I;biFU>)ZPf%_LSvSVm~1UpPIB^rSK@|7)u8 z3Os1xz#EGTws3L5*6}=k=@YxUVaW}Zgb+shO74`m2M|rQvh8+Lc1TEUM*G-d^r?3C z_Il$w?A+XRlppB6dx9-Q@r0T2Z5F{L*kTk<=*^2-ThI3B<~R!??o;0o0z)a*^(%NU zRzEzr<>kwlO#nE)De~$@M(pTs-rIspABFuw#-xyl2rTGbcecO?uL)+M_ira~fR|!c zOHhexvahSF6Gsc4IXtvA+KhyOxTqPQM1y1T7^1TO@CS7*7)0E?d)E_-958ALf}>ZT^ZD$aTDcvvjb zuHsv_Zk4S|j)e~z4^MKx%1zFFh6R*}$Vi-XODQ%uDb}z_=csjat=hYHF9swW9hpY5 z;MdKqp7R8ms<@#Dwc4=bhYd)FZG zTtj*F>Qw>GB-!uSUS~Se@HZ$)N~#m#^IG`n(?&q$L&)FJC&3k|r^bMTow*Fj=qB6C z=N8d!6udT1ylhOsiYwtiquGH_)1?I*-c{@?*#{qV=CVXD#ZwFn4D{Q#2TD18?||SV z%W;GirU81f^kN!9~r@+9OvQYS}1Z0tV}Vu%hjf zCL4Z@FNV_A)J1Axf^HK(y9E#3^Ln-0+^2#RGo7DRF__VLEx~J2l{(8QD7XP&j5$3X zPRhl><`meeE2TPpL}>Eo-6oD{1;3}KJIM`Q9$}Uv_lK1!M82mu(_fL5H1yc?@^KtbdW_KK6&z_?@wl4 zy;T%o)W)O~{xK9uyzkxHNc^drr#kR@R<2x0Px}g8UQhAVSpSb3H4ZA7*NQMHt1-9K+A7Nqbg#Wjiq8I~NNV#jX_EHb!Iv!=`({b00 z?(oAXvlt5|)U;2n6^W%}pBi^)+V7$(ibj||5wmhX!0T;{NpaM77=4aoigCc;)lD6e ze1%#bKKDF9wi*PGex1T)OiT!;Y-b%cBZ#N;HGP(pNz^Jkez6QQHrXQnh zWa;lRgOIrw(^!0AWod|T3CbXpF{}E+2Lig${IU+BlB?goOim}YHP+wdNT`0en;MM7 ze1#SX8GRb;dhN3gY}+numvK)Py<<6x=T}})QD>2z27p;H*~1;BX@LBp0B0nHLZUxj zbx;JpneA|IuZH#ZhPY;2Q$gezV6JaB>_>BFlWL^7lMY5r#9f%YkIo_Sm#Ai-$@CpV z?rASNMw1cju!{3u>d9IgU>tiCrzS~Hv$C=ZkdZ{m*mQW7{sH)Mp8q33nWM|0`~Zm+ zZ44CVi?Lg4U^=@Qo5}`Omr>Ggtc8VzS+3J*!}n!8<`3YUy@5Qfpsh_mMUb&6%XSy? zdI3Kj#BSsxOW!(<&|wde-&9w3rkaK4z_`Me3y<;x0s@jTlZGk!yeZT1a3l7`NE5h7 z{KRM$w6w4yV=~eRvn)Iw3F(%f@hLfC){n#qE(*(a?3vsb>NjNrxZi;HG8CNPc;F$e zeJ#P0Uv5`}`^+XWzo$A`<<-^2fa9i$^p+d9Z|4u(wsot5NLX-v1{?DZgUZ^EAL#)Z zxPQ6eP3UTL{2Uv*mY*+^{NjaTR5~)}F85=Y+S8!AZa$T!BS zf6!va2DXG>^tE9M|BieXjjatafv=a`5e4a=HVF{u@O>5x*-dY6KS%VDDvE?@mcYkk z>v6bZ98^pc2lv+1R9Q)fGTld=j{4VZQG;{u?7q<_6=7G_V9XGH#;N_lg|GkoV_LPpEz9=l$?=WUJ#MA-H>s7AbR6ak1T0GV zWakRclF&>{H{!5rOwv3(xx%GkiI-$B$2l;XF`z-`7;4pE+F1fjPeVaV{%^EXHFaJ0oLHlzmswLIqZ_D_$cC=3R$){s937J{h^zl>kB7vHZ7sp-saR!tflGzd6K z8OTt1czU1PkG365#_)g)gbp8^X5p5+DSmNry9qa0r~k_En$qkYbeb{)LM*$RU;PBL zPQ=5(Y6PmAwI_`CRBXhVe{52xkinf6b2yhF6|U2VZ;Mh%9QeWV$ylE1rY5b|pL=n? z8NJJT(w8Z|&3QEF7R({9v58N~(NPpfJg>`(A8y>f;CB4G(L>&WV?$^5fcqa7*vZ^6 zpMvSt=WZ&GdUeZR!hOOo>)pGj2)8jKqoa3FBw*YU2`&*Zl^s5|IijYOAx7apB609r zPv^sVsVYod&0*@h-pPo4hciA@SDSe4&_gO~R2?DMb!URGE@?8#=Kzy8l$9`cfszN| z1DaL8GUgQ-pVMsp5?Z!R_qWA!3kz=ovJB;$ox%*o@T&vGnidEr`-Fs+QVzKOrbihX zM5k}dZ^T`g0LW*|VnwvJDmt{Mfz9qhHf?DSSE?8!qbE+BxSOl*EOEr^$Tg(7{*Jo3 zTQNH~5C>qwhg?kQ!7j9eD5Fy;MZQ z%v3^vz?90rx~JZ(`IrN|G`Z=@kskIpr|$)U@?2D!i)$4owSW-aG!GK-8;;zg+6CiG#SyTqsfcK+# z0(~;xuSf^eF8Dqf1v;?E4OY4(sc+ttVWxpx5G+JjV>HhpkSVVKw2&^YKpp{ni#ZJ& zKjQQRFXcNAB>ZKdrHw;VQLlYAYx?qY8>L@pF-SQs`xd|2*zGEDDXHjyJ1G~KsLjn;XWwP96 zMd7(-gZYHm_Dk>{LPl2gBj21^gj?erWMQz^b4o~PFWKRa{QsHe0{Zejg5f33d+XUF z_(3S(JvBsNO3PWCeRUFhDHcCE3RPg9+}4Qm5DoV8@wb&VTl-3Zma0V!0EJ-r)b zxCL||T@Tfgwa$QtZgkOofceCVSO{1Dd4xKUBB3!AgJG(><=G4 z#Nh>2>iY)<^6&4(4weVhxT=?ru$cj`oH>~tjL1Oir{@TR426&)%$s@bE+REyHv>p5 z8A9$@e}xzScn_!zryZexfs|DAWEMe)_rHKzDyDO-#84vH+^fGbN+&^x!bIi>@RLaC z<}F){VGaKqE@H|rbtM8)0amnLTW>~9cHPpYbMvcgP_c7G#hUUnT=FYy%f3vh{(_y5 z5jp`ujUeSOK;#$1n?%rEjynbkRe(r!CSGYHFx;f(N4)KO}t zCW1=f6G2cZZ?W}`z%L@RcbmA)4zfH{c`3pdzp~(#i{Wc<*-z69YdZp=_sb}2- zTW{k#R-I1|jl*W?1OV#CN&2M;$Hk6)`{;qo-F6%AB#Z8<>ZV||{m#8K<~h#3jZ(&}$9O3WMk4E2f?6vB~!*nH-+ow5?87My^#*o_>! zcNZd|2mA6KGwe_;){BXH3lAJQL&F_}oJJ;!f<3~t zM8IG>;@NRPT-?nHK$e`7Q`-JuhOuC5rhl3q@R;9<{1f^R^QKojabofFxVSi_+@R={A=gKwyJ7uK1)$VZFSSdL-8M5ZDFBc}!nqxY&j$fw zK(AvDbqz58J>q3c;FLr^tIoBFj!XeLxtf*=kYr`!{*qo~-;qzQeU?7W;y*6_nQQ+S zPLk{46XUe*FS`1Bg^8MTixovGS&2W|BzW!q;}z0+0#%!Zxtb9spLpOW`Q(g4TVq(uys@3b0uK{(|CVCg5oy` z=d6z5VcohMo4mgJKiAyLb z7yKZV0{J#DGVD+^P&$T&{Gt2_mDlIaMFj|xSbQ8&C3<^#KqWB>2ZdKsQj$3g^`Bs= z{oD8Iu=65cPaHS&YRB=f-W#PnUhvh+i@My;6>tM=BnO3|~-})X%b@JWwj9FV6&&_F!&qYWgunCB5eM+MZt2ApG01vDrfQYIZb9;E;%5}n=Al~8qKi5pzj?+;=M9B?16Jrq=2ns-|QyDIB zhmbPBM*GCXJYE=IX#mJ0ixq=#N%{El-atXKGT zUqn>Y`^Jr<_iq!40EUAS1pw0OCy-8D04m{5Z+3LA=S~@K*REjqlSB<%s#E*4Jx+jDkvBkuEjmxk*94 z?$7iahbXbE*ucyy9V934T1DE;a7$p|CVbwE-?C<0P5 zkup1f-W$iJ5x_BQk)!pyKR4n)2*9uH1Zh`ZQ^Rl9R$RP{96gZNbwEWR75?zyI!wB# z)0pGE!zL35@agY}p#1~o`5_Te&mz4z5H?In@$#N@ zM9FOvauY4W>x*4X7A-?b&lXiSO2Kh9x0?jhZ<)d>pgyl$7cs(Y^ti%BjF9mlh1090 z4e$-dblr(^2vQ@0KXLSRAnTBOs4XQWMI;4zd3n$EbKlmC-h+~bhzc-|GYg~kg2qvl zA+iLd(_RoG426tCDT8$NC;|oCko}^Z9377VYmrj}R4Q{MI1AO`@Z}V6N!JiWJ^z^Q z5zu?PjDpiE&gTs1zCIKaeZ(mcFI^j+BFcl#o}Sn-9du&_fQ4U+Uwjm!H-@7CgQ1@i zJXf^rG8(mvu?r^+p-Y0s@P$qxPAm%He+Weio31oxXuUG5n;8&S8XrWGz;S|SF@CHi z8+F0fJg$;zf1w=)`P(q2KWEt)QV^qhWam-Onb43op!ytr|(KZi9eLcibvDj}z z!g+wi@(^@brZ~U{kNpd3h;QBJ-u85KJXZaVZ>50pi_WMS=9Us?M|sff9|oY50#~7A zRfF)RfD@tq!D&=b$>dqsNAGcQ1C-QhKN<;=KmpJ|@}H#Zk3iCw+hML>|JOG#X1#i+ z4k0c6D~{dpUy;yDi1&4O?}B3O1V2B&&UnOw2MQ#^ZW$vrmExJ18Ky1W$VVV@QtF%C z57%8%>Kk=Y)xtuMs{wPQfHy3e>j7>s7M65Jip0>nHtGRh=hxo|On zMesjD-Iqz+rRpT4*ma%Tu`L>0QY=->#%xi zw4HejoQf>4B|u~*|9EIn&^71<{o){Wy4}&)Nod&E`1oDCylZf_F^&OJ-vaytqWa8# zCj#OqHa0e-X*xLw^4_4`-xn6r;MJniatyJ}^Uvl3$fk=C?@M@ti>}|e(TU;-EH}q= zI%nxoRY9rXMDR@{+l*50;)H+x{P_{3$7Tm$pseid!n(Q*$QVe(g)6$i%e53DDngd` z7>Y4d$%mSlLb(EL!G_u#$Tpy$RnE3EJpdpGoe^WY3sf!pu(qH++X7DM(`+upaX0Y* zk0hx1H4YCCf5gBkBxV~xUlADtLUi#oSnSQ%$R&tG3RmJ5YH{dOm5tcXkian;i#vVa zO4Xl10K9sng5MA z=C420^ZTf^ zC0*^>9+kQKaQwz(?7jZ#+@UIPHdjMdraerCK@W>0)!;uMgRkO%0hoPd?0auRF|td0 zm{=>feS)pIk2r#l1r7yNH-dThApXH1{T@zMS!+Bw|J+eA-pnPfVL5a0zI!BZFg|dl=3yEE27mz5s zm$x^Q??I3GoSx3ko&F>6^Gg+)Rx`8HK6_RXwL(ZWLzwt=wR+x#gop2iR9a>7;0Rz( zP!tuVzBO7*7xdqqd3vnEPXiN>>sZ(zb(?+f3>wwHI46<4@lCsXlC3q|Tvd=cnRuuP z@H^Scy)053WDE-qW+%S+U;6NS<+`xU_W0Mmj%~%_HSzMqqGI&D4>b*MA2honP}*0A z5TM^H$B&aesH&{&xK++zx;j;v37(6|77JgXfo1h9V2s>%T+Q;fI-ecQUGNfdh;c~) zwNSI63dz@nJ|q8tPij4W{{EkIn$A$wMt_jS371-_W|c#tyZfJP`{Txsuxu2N8phlq znh4v(#9^R>GG>j1X^M@s=_%DYZISP9yf?N zl);SI!gX=wnHWrJj`b_r7=P}aRf3l3tBl*|;1ig2nM3P(wX0ixd#Z)kHs;Z!X zf_fW7Km2hK4=EcN#LT3&!IBiJ->P8NOVlX;~hvjyZz{-%Y`u}+{uTD6t z%YIfZ8C)pja zS04G`j;^klPvZJ%b=x=KUH%iT{Ye9;XM0Fi*|g=B!l zWX@+FfKiHVA!|VCNp==sJto2!;(8qR{UPbP*>>iwuV3(Zb%&s$s%dwJ14!o3_VzrYUm6-10jJ_B^^`qT^(2*oIhw~$r z4i0iC4LW=u8XTlUiFhp8T!$U1M3wl>M$`kx)N7G`y81#bJnm{kdaAQ8 zWxZ6>>kC^*BMZJHfGqPCGjsFj4hiLNLI>NdQXW7+1XM~e9P|ARj6kHQW2Fj2rl$*G zw!{tK)~8r8rjqO2^34)!jZwY>T!EzDHzPy1*?NatzOmQp+g{45&&JPFa_teKg3FjeU zDyH$ILoJ8B{kWJo!v{Kn__$M@+6I`IXc~6*iE{bRq~03HAb7$h6D3yrSx_b(*lIlH zUHBLc>f(!N#&Fp7F{8PvpMI6;i@9Ze_O78;?G?n3M%y)N*4#R1N` zX$v#jcF+n$3;lKxc2d>l6d||#)q`Q#JUL&Z^J)z{`*N~IM{zd!Af_x^w#*O#_&K&_ zPy`eQ2Dub7|~ zgN#qy_4h0156#V+NqGrxyYBJfrQ#A26?C2`n)BS7nim-R`SUtDB|tp};zaENB0+9a ziWh5IrPbThV+g+vs7eZfX3p*z)q5~{|HzX=;BMeH!mS?6cQ4^mTal1}KaU~JbbdbS zjDD!=vIBP`#gIcIAwv2MXs^r9Ye<8;LL7e&VuzH#8=Wg3436GjRGqZ3HMFm%zTr!( z|2lxwp{E0@zmx(#Y~U;W(3;s^@_|<}GE(X7=0R@sOw$!Dl*unF?&o7$hiPAgz}_|r zFB#{-KM9jTRC6YKp%Jooqu)O|iQdlh z&uAeE8)NwkbEm62bPxrJDZ!qqfHu?^fw3AT^nblGG&bJlmyToo)tfiRAz!1+F^}E{mE3L|*NE=Qn7+f#$$27WEJUPXXw~}lw>ab_ zM{h+yONS@#=hzviO(hrFIQTVBkf)InimHov-ql};ExDm|Xn7$lPLw_6xtmE8_t6y>B_>>0}=B{r=v^7vIYGByHEqhMbf?<^@m*957g8) zxD{_o`PSRpMP%UnncN;S%LkjZdlj{(0&w|uBPr!pN=!=Q)x`yiPWkPiTvTSX1|%-f@$Tq2p)c$0?M;(m&!>Y31TBt|CHT7$B+oH(Xy300+rV1`UEbyGWjHwG z-r>-De)i2vb1}=$E3sRMP#{HBp90Qpo35@dPWS6Xz;_kD-gI`7sBk3v2>iPo*wBk7 z-IWExP;#XGGU~H(xR#3#McxAAFB_@VkOpW-Y=U0MYckQtfBTG z3S!oxy#o1vV(dxZhrcn>AQWs~o-8{?R)!kuVtt~{ip&FPe#nnk^7-wks*Vn;`=7Dn zb>ACP*pN5qH?yq%gF3z>-`OlwL1m8}^P5CxL&l^w`ZJW$>qIZL-NvC`ZkNQ5$@Atu zXv^@O#5|JA4?ZxWvjX%c#S>>YL3lOTk8-bX>K;0DD5%+_q&dT4Bbme|qNJIaZUD+I zqqH58MkGQTMy}7@&!`i5$)RuhJy?Y=*!$*9pWXAqT&@dMMhQoD)q7*7yjMkn3-o<~CI#phdMiGw>Wvm$c0nu#`fhHFFK#e!mibH6OrTG51a=hrCgC zc`y>?<(Nm!*8-^_kt1hJw&uCLvkq$|$-(4d?0wV^|0NGeQO~yfw$HxzWL+S*a>$Uz zkXjjatE#A=-O&(2Q%bo&x$De8e2Qr@2=)FiFCvr6)gnOAk{Qq~cmsqV%_IlQVa&?` zi;yhZ3ofJ(p+pXfr^*Oem;_~EO#2gv=_sDap2?($eWZ)o_w{+SL( z0D7`M8HEh;W`U0nQ2d&^3KwQ_(Z>EFojYlB-|4D8d?>-~$+d#2N1==>eQfIiXL`(x zReN{YNI`*)XjO68kF@}{77G~^0RG`vi!(=`APdN-kPj=>Lc!O9(5^f$gH6YT!gX-T zT9$o>wyvP1?ZlX|PtFHH#H4LDfXWI0YvEB_V(oLHG_O}Q!g-N4IzbpU7MACt?x+m2 zh**Z0YeHRt2X-Z)bg?G(B_n0(Fm)-HLgGz|fT@k2@GE~{=bi%=Alni^AP6#UBG}1gX~*dZ!LWQl!?-zklq+lyX7mA69H1Fq(ALq`d55HN%O)hShkIP+8zN zEW<7AiTq}W^-aL`5idGcd-xRY5)olUE@|GLcaSs^fOmlYcJm?k&P#{@P5B=A5xGx* zL68cbeejAj8o~(d4(T<57z;OV-j*Xleh!k(0{k56pP+jQT)>{#hjx-~SObJCT8+ zQgR)ws{mDm4Mej*Cy)xd8<5PUT~EOJ1ldX`dF>}ph|*9TNeBS+F#r}yk3hh4ovhE{ z=l~K&Pa+M{nwu!;VGyve`Ht51!BZ$Bm5?P6>;P@y5m0huFJKV)Dga(DU*A=vPYWGG ze=rmxZF2yct*=q*JaHXYRzbIjDhFvmCiFO>Dd8Z|N^%pw+>c%Cu1xwMrW10Bwxe~? zhLq?fFRwf_Kb(i9XN1A+XqY?<{`xNHzm5(=BCbno=R+df1o=+T_0OmUjJ;nkvxC}} zfa&Sy$Rff4jXPL|^a)j{U815)WS1cE4+;)mvLl~U24!i6e^CGSB~k(MtafWq2-tj29S^o0*OvfS&5_#;mWa1m4w)MBg_ zMF^(U)Hp_LhqSc5C7LXk9qN08g1jqlg}PkcY>@|U+qLMf90x-qly zonhgU4MX3*^UvXY2}TaY^#)AW1S+YPx$#g_1_R2~7(<60LS|z7j1b)u6(#*~Cpv0p3E{`a<`>g;i zQgnfzn8m1J=j7BlJ*S6;D^wLNsvkXk$d9Azbjpw&)TqdT;zav!&hJJAF4qpe017x( zc!PORbY)`Rv`KB)rT5!6dAxFN(Kisvlz4l2Rn(vcoRAxaprE}UL1bDI;*B_H9Fjh~ zf%t2Km9gtsAIiy7xAkma>$e*%qZO*5JOBH6_qT6eAi<1*#U=X)=>TS)sKf$nk@6N0#;uAek(@p^Z{8#$ zo>;rlGfYYd03&NLTk9jsJ2YM`qosv_L)};xRXEalOWG}l#2ri_mU;YmN6y0hENeH; z+kKbX4n#>gtwU))Y((H$Ad)8_A}^e%VGTB3pdY%hZjM7}UUk4p8=|i=0E<)*WYLrU zY1oVapG41s7LV{|ab7+?TWGRe(M^exsqz+5CM4`O2DOQW%YKmEJw~Zv&J zY8j>(egyUVA{Cg1bUR3+k!t_8rQ~6s(LQrRh08I8szK)_eT^j2p&cRs9?i!6rT6af0z6}ic;wGz%FgfcQ*uo_g5NEe>GEpcf zZ4dzq9&{&mC+eoR5OSlUqw6pM5cAD7heY^wv*?j5F1fp~VjORDX?gpWK9tgk{0OyW z6dh1fW+ZKwxThkxW5kO%*+SU%%45MDr$7gAczIRTjevl)B#d0TbSYIhMw)~y^jLL4 zY6caJERpB-^xVX9egZZK+<_;Dd|K|iuV26V4M(dnpukgtH-;vK3mkj*zOSxcfl7HN0g|6I_#H+B3k0CfDq(!=gx5eCPSDx0Omf%HK3z{BApb70Jpl3FQ6{%9Tc>I zWTB`mKgG$!I_d;L0pb>a3n~#hs{^2!^zrwOn~#sL zN?5G!T`Vdf`$R>TQ_#GZ(^*ITfx34mcghDt2}H#@1n-bf?E1eM-r`@;#r=P%hu$wX z@{!#~KL39tO#inpF|*K;zLYx1{qTapn_VA?bsRoY!La7S$8+0XH{OCkVl{=2;c;Tf zh0m*9aE~FMrXXU1hR^Vczs9SM9B>*Vpa0%vmfZ?I`@OAcDwteCUQS8& {:.notice} -The help contents were last updated on **02 Aug 24 22:54:33** for anvi'o version **8-dev (marie)**. +The help contents were last updated on **14 Aug 24 13:42:12** for anvi'o version **8-dev (marie)**. {% include _project-anvio-version.html %} @@ -132,7 +132,7 @@ Listed below **a total of 134 artifacts**. Anvi'o programs perform atomic tasks that can be weaved together to implement complete 'omics workflows. Please note that there may be programs that are not listed on this page. You can type 'anvi-' in your terminal, and press the TAB key twice to see the full list of programs available to you on your system, and type `anvi-program-name --help` to read the full list of command line options. -Listed below **a total of 151 programs**. +Listed below **a total of 152 programs**.
    @@ -4029,6 +4029,36 @@ Listed below **a total of 151 programs**.
    +
    + + + + + + + + + + + + + + + +
    + 🔥 **[anvi-script-find-misassemblies](programs/anvi-script-find-misassemblies)**. This script report errors in long read assembly using read-recruitment information. The input file should be a BAM file of long reads mapped to an assembly made from these reads.. +
    + 🧀 + [bam-file](artifacts/bam-file) +
    + 🍕 + n/a +
    + 🧠
    + +
    +
    +
    diff --git a/help/main/programs/anvi-script-find-misassemblies/index.md b/help/main/programs/anvi-script-find-misassemblies/index.md new file mode 100644 index 0000000..6e5b668 --- /dev/null +++ b/help/main/programs/anvi-script-find-misassemblies/index.md @@ -0,0 +1,115 @@ +--- +layout: program +title: anvi-script-find-misassemblies +excerpt: An anvi'o program. This script report errors in long read assembly using read-recruitment information. +categories: [anvio] +comments: false +redirect_from: /m/anvi-script-find-misassemblies +image: + featurerelative: ../../../images/header.png + display: true +--- + +This script report errors in long read assembly using read-recruitment information. The input file should be a BAM file of long reads mapped to an assembly made from these reads.. + +🔙 **[To the main page](../../)** of anvi'o programs and artifacts. + + +{% include _toc.html %} +
    +{% capture network_path %}{{ "network.json" }}{% endcapture %} +{% capture network_height %}{{ 300 }}{% endcapture %} +{% include _project-anvio-graph.html %} + + +## Authors + + + + + +## Can consume + + +

    [bam-file](../../artifacts/bam-file)

    + + +## Can provide + + +This program does not seem to provide any artifacts. Such programs usually print out some information for you to see or alter some anvi'o artifacts without producing any immediate outputs. + + +## Usage + + +The aim of this script is to find potential assembly errors from long read assemblers. + +### Principle + +This script searches for potential assembly errors in a [bam-file](/help/main/artifacts/bam-file) generated from the mapping of long reads onto an assembly made using the same reads. The basic principle is that (1) every single nucleotide, and (2) consecutive pairs/triplets of nucleotides in an assembly should be covered by at least one of the long reads that was used to generate the assembly. + +To find out if every nucleotide is covered by at least one read, one simply needs to find regions with 0x coverage, and that is one output of this script. + +To find where two/three consecutive nucleotides are not covered by at least one read, we can leverage the clipping information reported by long read mapping software like [minimap2](https://github.com/lh3/minimap2). Clipping happens when the left- or right-most part of a read does not align to the reference. If 100% of reads clip at the same nucleotide position, then it means that not a single read is covering at least two consecutive nucleotides. AND THAT'S NOT GOOD. + +Here is an example visualized in [IGV](https://igv.org/). All reads are clipped at the same position. There is no support in the long reads that the left and right sides of the contig should be joined, suggesting that there is a misassembly here. + +![clipping_example](../../images/anvi-script-find-misassemblies.png) + + +### Basic usage + +The only input file to this script is a simple [bam-file](/help/main/artifacts/bam-file). But not any BAM file. It has to made by mapping long reads onto an assembly generated using the same reads. You also need to provide an output prefix. + +
    +anvi-script-find-misassemblies -b sample01.bam -o result +
    + + +### Outputs + +The first output is a table reporting regions in your assemblies with zero coverage. This output includes the contig in which the region occurs, the contig's length, the region's start and stop positions, and the length of the region with no coverage. + +|**contig**|**length**|**range**|**range_size**| +|:--|:--|:--|:--| +|contig_001|1665603|0-498|498| +|contig_001|1665603|100500-101000|500| +|contig_001|1665603|1665106-1665603|497| + +The second output is a table reporting the positions with high relative abundance of clipped reads. It includes the contig in which the position occurs, the contig's length, the position in the contig with high levels of clipping, its relative position in the contig, the total coverage at that position, the coverage of reads clipping at that position, and the ratio between these two coverage values. The 'relative position' is the position divided by the contig's length (a value from 0 to 1), which tells you whether the position is close to the contig ends or somewhere in the middle. + +|**contig**|**length**|**pos**|**relative_pos**|**cov**|**clipping**|**clipping_ratio**| +|:--|:--|:--|:--|:--|:--|:--| +|contig_001|1665603|498|0.0002989908159387321|48|48|1.0| +|contig_001|1665603|500999|0.30079136504917436|120|120|1.0| +|contig_001|1665603|501000|0.30079196543233894|79|79|1.0| +|contig_001|1665603|1665105|0.9997010091840612|45|45|1.0| + + +### Additional parameters + +By default, the script will report clipping positions if the ratio of clipping reads to the total number of reads is over 80%. You can change that threshold with the flag `--clipping-ratio`: + +
    +anvi-script-find-misassemblies -b sample01.bam -o result --clipping-ratio 0.6 +
    + + +Another default behavior of the script is to skip the first and last 100bp of a contig (only valid for the table reporting positions of high clipping). This is because contig ends often have high proportions of clipping which have nothing to do with misassemblies. You can change that parameter with the flag `--min-dist-to-end`: + +
    +anvi-script-find-misassemblies -b sample01.bam -o result --min-dist-to-end 0 +
    + + +{:.notice} +Edit [this file](https://github.com/merenlab/anvio/tree/master/anvio/docs/programs/anvi-script-find-misassemblies.md) to update this information. + + +## Additional Resources + + + +{:.notice} +Are you aware of resources that may help users better understand the utility of this program? Please feel free to edit [this file](https://github.com/merenlab/anvio/tree/master/bin/anvi-script-find-misassemblies) on GitHub. If you are not sure how to do that, find the `__resources__` tag in [this file](https://github.com/merenlab/anvio/blob/master/bin/anvi-interactive) to see an example. diff --git a/help/main/programs/anvi-script-find-misassemblies/network.json b/help/main/programs/anvi-script-find-misassemblies/network.json new file mode 100644 index 0000000..78925d1 --- /dev/null +++ b/help/main/programs/anvi-script-find-misassemblies/network.json @@ -0,0 +1,30 @@ +{ + "graph": [], + "nodes": [ + { + "size": 1, + "score": 1, + "color": "#AA0000", + "id": "bam-file", + "name": "bam-file", + "provided_by_anvio": false, + "type": "BAM" + }, + { + "size": 1, + "score": 0.1, + "color": "#AAAA00", + "id": "anvi-script-find-misassemblies", + "name": "anvi-script-find-misassemblies", + "type": "PROGRAM" + } + ], + "links": [ + { + "target": 1, + "source": 0 + } + ], + "directed": false, + "multigraph": false +} \ No newline at end of file diff --git a/network/network.json b/network/network.json index bc16507..33f6551 100644 --- a/network/network.json +++ b/network/network.json @@ -65,7 +65,7 @@ "type": "DB" }, { - "size": 75, + "size": 76, "score": 0.5, "color": "#00AA00", "id": "contigs-db", @@ -677,7 +677,7 @@ "type": "FASTA" }, { - "size": 10, + "size": 11, "score": 1, "color": "#AA0000", "id": "bam-file", @@ -1126,6 +1126,15 @@ "provided_by_anvio": false, "type": "TXT" }, + { + "size": 1, + "score": 0.5, + "color": "#00AA00", + "id": "genes-stats", + "name": "GENE STATS", + "provided_by_anvio": true, + "type": "STATS" + }, { "size": 1, "score": 1, @@ -2186,6 +2195,14 @@ "name": "anvi-script-filter-hmm-hits-table", "type": "PROGRAM" }, + { + "size": 1, + "score": 0.1, + "color": "#AAAA00", + "id": "anvi-script-find-misassemblies", + "name": "anvi-script-find-misassemblies", + "type": "PROGRAM" + }, { "size": 2, "score": 0.1, @@ -2266,6 +2283,14 @@ "name": "anvi-script-gen-user-module-file", "type": "PROGRAM" }, + { + "size": 3, + "score": 0.1, + "color": "#AAAA00", + "id": "anvi-script-gen_stats_for_single_copy_genes.py", + "name": "anvi-script-gen_stats_for_single_copy_genes.py", + "type": "PROGRAM" + }, { "size": 4, "score": 0.1, @@ -2373,15 +2398,11 @@ ], "links": [ { - "source": 129, + "source": 130, "target": 0 }, { - "target": 129, - "source": 1 - }, - { - "target": 132, + "target": 130, "source": 1 }, { @@ -2393,35 +2414,35 @@ "source": 1 }, { - "target": 137, + "target": 135, "source": 1 }, { - "target": 145, + "target": 138, "source": 1 }, { - "target": 147, + "target": 146, "source": 1 }, { - "target": 150, + "target": 148, "source": 1 }, { - "source": 170, - "target": 1 + "target": 151, + "source": 1 }, { - "target": 180, - "source": 1 + "source": 171, + "target": 1 }, { "target": 181, "source": 1 }, { - "target": 203, + "target": 182, "source": 1 }, { @@ -2429,15 +2450,15 @@ "source": 1 }, { - "target": 206, + "target": 205, "source": 1 }, { - "target": 225, + "target": 207, "source": 1 }, { - "target": 241, + "target": 226, "source": 1 }, { @@ -2445,27 +2466,27 @@ "source": 1 }, { - "target": 246, + "target": 243, "source": 1 }, { - "target": 253, + "target": 247, "source": 1 }, { - "target": 260, + "target": 254, "source": 1 }, { - "target": 261, + "target": 262, "source": 1 }, { - "target": 129, - "source": 2 + "target": 263, + "source": 1 }, { - "target": 132, + "target": 130, "source": 2 }, { @@ -2473,31 +2494,31 @@ "source": 2 }, { - "target": 136, + "target": 134, "source": 2 }, { - "target": 139, + "target": 137, "source": 2 }, { - "target": 145, + "target": 140, "source": 2 }, { - "target": 156, + "target": 146, "source": 2 }, { - "source": 188, - "target": 2 + "target": 157, + "source": 2 }, { - "source": 215, + "source": 189, "target": 2 }, { - "source": 218, + "source": 216, "target": 2 }, { @@ -2509,23 +2530,23 @@ "target": 2 }, { - "target": 232, - "source": 2 + "source": 221, + "target": 2 }, { - "target": 260, + "target": 233, "source": 2 }, { - "target": 261, + "target": 262, "source": 2 }, { - "target": 129, - "source": 3 + "target": 263, + "source": 2 }, { - "target": 133, + "target": 130, "source": 3 }, { @@ -2537,11 +2558,11 @@ "source": 3 }, { - "target": 137, + "target": 136, "source": 3 }, { - "target": 141, + "target": 138, "source": 3 }, { @@ -2549,51 +2570,51 @@ "source": 3 }, { - "target": 147, + "target": 143, "source": 3 }, { - "target": 150, + "target": 148, "source": 3 }, { - "target": 159, + "target": 151, "source": 3 }, { - "target": 161, + "target": 160, "source": 3 }, { - "target": 164, + "target": 162, "source": 3 }, { - "target": 181, + "target": 165, "source": 3 }, { - "target": 187, + "target": 182, "source": 3 }, { - "target": 189, + "target": 188, "source": 3 }, { - "target": 191, + "target": 190, "source": 3 }, { - "target": 193, + "target": 192, "source": 3 }, { - "target": 201, + "target": 194, "source": 3 }, { - "target": 203, + "target": 202, "source": 3 }, { @@ -2601,12 +2622,12 @@ "source": 3 }, { - "source": 206, - "target": 3 + "target": 205, + "source": 3 }, { - "target": 239, - "source": 3 + "source": 207, + "target": 3 }, { "target": 240, @@ -2621,41 +2642,41 @@ "source": 3 }, { - "target": 246, + "target": 243, "source": 3 }, { - "target": 248, + "target": 247, "source": 3 }, { - "target": 253, + "target": 249, "source": 3 }, { - "source": 130, + "target": 254, + "source": 3 + }, + { + "source": 131, "target": 4 }, { - "target": 130, + "target": 131, "source": 4 }, { - "target": 138, + "target": 139, "source": 4 }, { - "target": 146, + "target": 147, "source": 4 }, { - "source": 147, + "source": 148, "target": 4 }, - { - "target": 149, - "source": 4 - }, { "target": 150, "source": 4 @@ -2669,15 +2690,15 @@ "source": 4 }, { - "target": 154, + "target": 153, "source": 4 }, { - "target": 169, + "target": 155, "source": 4 }, { - "target": 175, + "target": 170, "source": 4 }, { @@ -2689,65 +2710,65 @@ "source": 4 }, { - "target": 185, + "target": 178, + "source": 4 + }, + { + "target": 186, "source": 4 }, { - "source": 187, + "source": 188, "target": 4 }, { - "source": 198, + "source": 199, "target": 4 }, { - "target": 201, + "target": 202, "source": 4 }, { - "source": 212, + "source": 213, "target": 4 }, { - "target": 212, + "target": 213, "source": 4 }, { - "target": 241, + "target": 242, "source": 4 }, { - "target": 242, + "target": 243, "source": 4 }, { - "source": 248, + "source": 249, "target": 4 }, { - "target": 259, + "target": 261, "source": 4 }, { - "target": 262, + "target": 264, "source": 4 }, { - "source": 130, + "source": 131, "target": 5 }, { - "target": 146, + "target": 147, "source": 5 }, { - "source": 147, + "source": 148, "target": 5 }, - { - "target": 150, - "source": 5 - }, { "target": 151, "source": 5 @@ -2757,19 +2778,19 @@ "source": 5 }, { - "target": 167, + "target": 153, "source": 5 }, { - "target": 169, + "target": 168, "source": 5 }, { - "target": 174, + "target": 170, "source": 5 }, { - "target": 176, + "target": 175, "source": 5 }, { @@ -2777,63 +2798,63 @@ "source": 5 }, { - "target": 183, + "target": 178, "source": 5 }, { - "target": 185, + "target": 184, "source": 5 }, { - "target": 197, + "target": 186, "source": 5 }, { - "source": 198, + "target": 198, + "source": 5 + }, + { + "source": 199, "target": 5 }, { - "target": 198, + "target": 199, "source": 5 }, { - "target": 201, + "target": 202, "source": 5 }, { - "source": 211, + "source": 212, "target": 5 }, { - "target": 211, + "target": 212, "source": 5 }, { - "source": 212, + "source": 213, "target": 5 }, { - "target": 212, + "target": 213, "source": 5 }, { - "source": 248, + "source": 249, "target": 5 }, { - "source": 253, + "source": 254, "target": 5 }, { - "target": 259, + "target": 261, "source": 5 }, { - "target": 130, - "source": 6 - }, - { - "target": 137, + "target": 131, "source": 6 }, { @@ -2841,7 +2862,7 @@ "source": 6 }, { - "target": 141, + "target": 139, "source": 6 }, { @@ -2849,11 +2870,11 @@ "source": 6 }, { - "target": 146, + "target": 143, "source": 6 }, { - "target": 148, + "target": 147, "source": 6 }, { @@ -2873,11 +2894,11 @@ "source": 6 }, { - "target": 154, + "target": 153, "source": 6 }, { - "target": 158, + "target": 155, "source": 6 }, { @@ -2885,7 +2906,7 @@ "source": 6 }, { - "target": 161, + "target": 160, "source": 6 }, { @@ -2893,11 +2914,11 @@ "source": 6 }, { - "target": 164, + "target": 163, "source": 6 }, { - "target": 167, + "target": 165, "source": 6 }, { @@ -2909,7 +2930,7 @@ "source": 6 }, { - "target": 174, + "target": 170, "source": 6 }, { @@ -2925,7 +2946,7 @@ "source": 6 }, { - "target": 182, + "target": 178, "source": 6 }, { @@ -2933,15 +2954,15 @@ "source": 6 }, { - "target": 185, + "target": 184, "source": 6 }, { - "target": 187, + "target": 186, "source": 6 }, { - "target": 189, + "target": 188, "source": 6 }, { @@ -2961,7 +2982,7 @@ "source": 6 }, { - "target": 197, + "target": 194, "source": 6 }, { @@ -2969,19 +2990,19 @@ "source": 6 }, { - "source": 200, - "target": 6 + "target": 199, + "source": 6 }, { - "target": 201, - "source": 6 + "source": 201, + "target": 6 }, { - "target": 204, + "target": 202, "source": 6 }, { - "target": 211, + "target": 205, "source": 6 }, { @@ -2989,11 +3010,11 @@ "source": 6 }, { - "target": 228, + "target": 213, "source": 6 }, { - "target": 239, + "target": 229, "source": 6 }, { @@ -3009,39 +3030,39 @@ "source": 6 }, { - "target": 246, + "target": 243, "source": 6 }, { - "target": 248, + "target": 247, "source": 6 }, { - "target": 259, + "target": 249, "source": 6 }, { - "target": 262, + "target": 261, "source": 6 }, { - "target": 271, + "target": 264, "source": 6 }, { - "target": 130, - "source": 7 + "target": 274, + "source": 6 }, { "target": 131, "source": 7 }, { - "target": 137, + "target": 132, "source": 7 }, { - "target": 139, + "target": 138, "source": 7 }, { @@ -3049,15 +3070,15 @@ "source": 7 }, { - "target": 144, + "target": 141, "source": 7 }, { - "target": 146, + "target": 145, "source": 7 }, { - "target": 148, + "target": 147, "source": 7 }, { @@ -3077,7 +3098,7 @@ "source": 7 }, { - "target": 155, + "target": 153, "source": 7 }, { @@ -3093,7 +3114,7 @@ "source": 7 }, { - "target": 160, + "target": 159, "source": 7 }, { @@ -3109,12 +3130,12 @@ "source": 7 }, { - "source": 166, - "target": 7 + "target": 164, + "source": 7 }, { - "target": 167, - "source": 7 + "source": 167, + "target": 7 }, { "target": 168, @@ -3125,11 +3146,11 @@ "source": 7 }, { - "target": 172, + "target": 170, "source": 7 }, { - "target": 174, + "target": 173, "source": 7 }, { @@ -3157,7 +3178,7 @@ "source": 7 }, { - "target": 182, + "target": 181, "source": 7 }, { @@ -3173,7 +3194,7 @@ "source": 7 }, { - "target": 187, + "target": 186, "source": 7 }, { @@ -3181,15 +3202,15 @@ "source": 7 }, { - "target": 191, + "target": 189, "source": 7 }, { - "target": 194, + "target": 192, "source": 7 }, { - "target": 197, + "target": 195, "source": 7 }, { @@ -3197,15 +3218,15 @@ "source": 7 }, { - "target": 200, + "target": 199, "source": 7 }, { - "target": 204, + "target": 201, "source": 7 }, { - "target": 208, + "target": 205, "source": 7 }, { @@ -3225,7 +3246,7 @@ "source": 7 }, { - "target": 215, + "target": 213, "source": 7 }, { @@ -3257,7 +3278,7 @@ "source": 7 }, { - "target": 224, + "target": 223, "source": 7 }, { @@ -3269,11 +3290,11 @@ "source": 7 }, { - "target": 228, + "target": 227, "source": 7 }, { - "target": 240, + "target": 229, "source": 7 }, { @@ -3289,7 +3310,7 @@ "source": 7 }, { - "target": 246, + "target": 244, "source": 7 }, { @@ -3301,27 +3322,27 @@ "source": 7 }, { - "target": 254, + "target": 249, "source": 7 }, { - "target": 256, + "target": 255, "source": 7 }, { - "target": 259, + "target": 257, "source": 7 }, { - "target": 262, + "target": 261, "source": 7 }, { - "target": 268, + "target": 264, "source": 7 }, { - "target": 270, + "target": 269, "source": 7 }, { @@ -3329,19 +3350,23 @@ "source": 7 }, { - "target": 131, - "source": 8 + "target": 273, + "source": 7 }, { - "target": 148, + "target": 274, + "source": 7 + }, + { + "target": 132, "source": 8 }, { - "target": 167, + "target": 149, "source": 8 }, { - "target": 174, + "target": 168, "source": 8 }, { @@ -3349,87 +3374,87 @@ "source": 8 }, { - "target": 131, + "target": 176, + "source": 8 + }, + { + "target": 132, "source": 9 }, { - "target": 140, + "target": 141, "source": 9 }, { - "target": 182, + "target": 183, "source": 9 }, { - "target": 216, + "target": 217, "source": 9 }, { - "target": 256, + "target": 257, "source": 9 }, { - "source": 258, + "source": 260, "target": 9 }, { - "target": 263, + "target": 265, "source": 9 }, { - "target": 268, + "target": 271, "source": 9 }, { - "source": 269, + "source": 272, "target": 9 }, { - "source": 272, + "source": 275, "target": 9 }, { - "source": 132, + "source": 133, "target": 10 }, { - "source": 133, + "source": 134, "target": 10 }, { - "source": 136, + "source": 137, "target": 10 }, { - "source": 145, + "source": 146, "target": 10 }, { - "source": 260, + "source": 262, "target": 10 }, { - "target": 132, + "target": 133, "source": 11 }, { - "target": 136, + "target": 137, "source": 11 }, { - "target": 145, + "target": 146, "source": 11 }, { - "target": 260, + "target": 262, "source": 11 }, { - "target": 132, - "source": 12 - }, - { - "target": 135, + "target": 133, "source": 12 }, { @@ -3437,15 +3462,15 @@ "source": 12 }, { - "target": 143, + "target": 137, "source": 12 }, { - "target": 145, + "target": 144, "source": 12 }, { - "target": 149, + "target": 146, "source": 12 }, { @@ -3453,11 +3478,11 @@ "source": 12 }, { - "target": 170, + "target": 151, "source": 12 }, { - "target": 176, + "target": 171, "source": 12 }, { @@ -3465,333 +3490,333 @@ "source": 12 }, { - "target": 182, + "target": 178, "source": 12 }, { - "target": 260, + "target": 183, "source": 12 }, { - "target": 261, + "target": 262, "source": 12 }, { - "source": 262, + "target": 263, + "source": 12 + }, + { + "source": 264, "target": 12 }, { - "target": 263, + "target": 265, "source": 12 }, { - "target": 132, + "target": 133, "source": 13 }, { - "target": 135, + "target": 136, "source": 13 }, { - "target": 136, + "target": 137, "source": 13 }, { - "target": 143, + "target": 144, "source": 13 }, { - "target": 145, + "target": 146, "source": 13 }, { - "target": 150, + "target": 151, "source": 13 }, { - "target": 170, + "target": 171, "source": 13 }, { - "target": 176, + "target": 177, "source": 13 }, { - "target": 177, + "target": 178, "source": 13 }, { - "target": 182, + "target": 183, "source": 13 }, { - "target": 203, + "target": 204, "source": 13 }, { - "target": 260, + "target": 262, "source": 13 }, { - "target": 261, + "target": 263, "source": 13 }, { - "source": 262, + "source": 264, "target": 13 }, { - "target": 263, + "target": 265, "source": 13 }, { - "target": 133, + "target": 134, "source": 14 }, { - "target": 141, + "target": 142, "source": 14 }, { - "target": 161, + "target": 162, "source": 14 }, { - "source": 191, + "source": 192, "target": 14 }, { - "source": 228, + "source": 229, "target": 14 }, { - "source": 135, + "source": 136, "target": 15 }, { - "target": 143, + "target": 144, "source": 15 }, { - "source": 252, + "source": 253, "target": 15 }, { - "target": 136, + "target": 137, "source": 16 }, { - "source": 150, + "source": 151, "target": 16 }, { - "target": 136, + "target": 137, "source": 17 }, { - "source": 150, + "source": 151, "target": 17 }, { - "target": 137, + "target": 138, "source": 18 }, { - "target": 148, + "target": 149, "source": 18 }, { - "target": 165, + "target": 166, "source": 18 }, { - "target": 167, + "target": 168, "source": 18 }, { - "source": 172, + "source": 173, "target": 18 }, { - "target": 174, + "target": 175, "source": 18 }, { - "target": 204, + "target": 205, "source": 18 }, { - "target": 247, + "target": 248, "source": 18 }, { - "target": 137, + "target": 138, "source": 19 }, { - "source": 169, + "source": 170, "target": 19 }, { - "target": 198, + "target": 199, "source": 19 }, { - "target": 204, + "target": 205, "source": 19 }, { - "target": 228, + "target": 229, "source": 19 }, { - "target": 140, + "target": 141, "source": 20 }, { - "target": 182, + "target": 183, "source": 20 }, { - "source": 216, + "source": 217, "target": 20 }, { - "target": 221, + "target": 222, "source": 20 }, { - "source": 224, + "source": 225, "target": 20 }, { - "source": 256, + "source": 257, "target": 20 }, { - "target": 256, + "target": 257, "source": 20 }, { - "target": 263, + "target": 265, "source": 20 }, { - "target": 268, + "target": 271, "source": 20 }, { - "target": 141, + "target": 142, "source": 21 }, { - "target": 161, + "target": 162, "source": 21 }, { - "source": 181, + "source": 182, "target": 21 }, { - "source": 191, + "source": 192, "target": 21 }, { - "source": 228, + "source": 229, "target": 21 }, { - "target": 141, + "target": 142, "source": 22 }, { - "target": 161, + "target": 162, "source": 22 }, { - "source": 191, + "source": 192, "target": 22 }, { - "target": 141, + "target": 142, "source": 23 }, { - "target": 161, + "target": 162, "source": 23 }, { - "source": 191, + "source": 192, "target": 23 }, { - "target": 141, + "target": 142, "source": 24 }, { - "target": 161, + "target": 162, "source": 24 }, { - "source": 191, + "source": 192, "target": 24 }, { - "source": 217, + "source": 218, "target": 24 }, { - "target": 142, + "target": 143, "source": 25 }, { - "target": 164, + "target": 165, "source": 25 }, { - "source": 193, + "source": 194, "target": 25 }, { - "source": 143, + "source": 144, "target": 26 }, { - "target": 143, + "target": 144, "source": 26 }, { - "target": 219, + "target": 220, "source": 26 }, { - "target": 226, + "target": 227, "source": 26 }, { - "target": 252, + "target": 253, "source": 26 }, { - "source": 257, + "source": 259, "target": 26 }, { - "target": 257, + "target": 259, "source": 26 }, { - "target": 276, + "target": 279, "source": 26 }, { - "source": 144, + "source": 145, "target": 27 }, - { - "source": 144, - "target": 28 - }, { "source": 145, "target": 28 @@ -3809,7 +3834,7 @@ "target": 28 }, { - "source": 197, + "source": 149, "target": 28 }, { @@ -3817,936 +3842,948 @@ "target": 28 }, { - "source": 251, + "source": 199, "target": 28 }, { - "source": 261, + "source": 252, "target": 28 }, { - "source": 277, + "source": 263, "target": 28 }, { - "source": 144, + "source": 280, + "target": 28 + }, + { + "source": 145, "target": 29 }, { - "source": 147, + "source": 148, "target": 29 }, { - "source": 198, + "source": 199, "target": 29 }, { - "target": 146, + "target": 147, "source": 30 }, { - "target": 150, + "target": 151, "source": 30 }, { - "target": 210, + "target": 211, "source": 30 }, { - "target": 218, + "target": 219, "source": 30 }, { - "source": 231, + "source": 232, "target": 30 }, { - "target": 146, + "target": 147, "source": 31 }, { - "target": 150, + "target": 151, "source": 31 }, { - "target": 210, + "target": 211, "source": 31 }, { - "source": 218, + "source": 219, "target": 31 }, { - "source": 147, + "source": 148, "target": 32 }, { - "target": 148, + "target": 149, "source": 33 }, { - "target": 167, + "target": 168, "source": 33 }, { - "target": 173, + "target": 174, "source": 33 }, { - "source": 174, + "source": 175, "target": 33 }, { - "target": 179, + "target": 180, "source": 33 }, { - "target": 277, + "target": 280, "source": 33 }, { - "target": 279, + "target": 282, "source": 33 }, { - "source": 149, + "source": 150, "target": 34 }, { - "target": 150, + "target": 151, "source": 35 }, { - "target": 151, + "target": 152, "source": 35 }, { - "target": 152, + "target": 153, "source": 35 }, { - "target": 150, + "target": 151, "source": 36 }, { - "source": 238, + "source": 239, "target": 36 }, { - "target": 238, + "target": 239, "source": 36 }, { - "source": 266, + "source": 268, "target": 36 }, { - "target": 150, + "target": 151, "source": 37 }, { - "source": 151, + "source": 152, "target": 38 }, { - "source": 152, + "source": 153, "target": 38 }, { - "source": 151, + "source": 152, "target": 39 }, { - "source": 152, + "source": 153, "target": 39 }, { - "target": 151, + "target": 152, "source": 40 }, { - "source": 221, + "source": 222, "target": 40 }, { - "target": 152, + "target": 153, "source": 41 }, { - "source": 222, + "source": 223, "target": 41 }, { - "target": 152, + "target": 153, "source": 42 }, { - "target": 226, + "target": 227, "source": 42 }, { - "source": 153, + "source": 154, "target": 43 }, { - "source": 159, + "source": 160, "target": 43 }, { - "target": 189, + "target": 190, "source": 43 }, { - "target": 191, + "target": 192, "source": 43 }, { - "target": 198, + "target": 199, "source": 43 }, { - "source": 199, + "source": 200, "target": 43 }, { - "target": 153, + "target": 154, "source": 44 }, { - "source": 154, + "source": 155, "target": 45 }, { - "target": 187, + "target": 188, "source": 45 }, { - "target": 267, + "target": 270, "source": 45 }, { - "target": 270, + "target": 273, "source": 45 }, { - "source": 155, + "source": 156, "target": 46 }, { - "source": 162, + "source": 163, "target": 46 }, { - "target": 166, + "target": 167, "source": 46 }, { - "source": 255, + "source": 256, "target": 46 }, { - "target": 255, + "target": 256, "source": 46 }, { - "source": 271, + "source": 274, "target": 46 }, { - "source": 273, + "source": 276, "target": 46 }, { - "source": 274, + "source": 277, "target": 46 }, { - "source": 276, + "source": 279, "target": 46 }, { - "source": 156, + "source": 157, "target": 47 }, { - "target": 188, + "target": 189, "source": 47 }, { - "source": 219, + "source": 220, "target": 47 }, { - "source": 225, + "source": 226, "target": 47 }, { - "source": 268, + "source": 271, "target": 47 }, { - "source": 273, + "source": 276, "target": 47 }, { - "source": 274, + "source": 277, "target": 47 }, { - "source": 278, + "source": 281, "target": 47 }, { - "target": 278, + "target": 281, "source": 47 }, { - "source": 157, + "source": 158, "target": 48 }, { - "source": 278, + "source": 281, "target": 48 }, { - "target": 278, + "target": 281, "source": 48 }, { - "source": 158, + "source": 159, "target": 49 }, { - "source": 162, + "source": 163, "target": 49 }, { - "source": 185, + "source": 186, "target": 49 }, { - "source": 267, + "source": 270, "target": 49 }, { - "source": 158, + "source": 159, "target": 50 }, { - "source": 159, + "source": 160, "target": 51 }, { - "target": 189, + "target": 190, "source": 51 }, { - "source": 159, + "source": 160, "target": 52 }, { - "source": 171, + "source": 172, "target": 52 }, { - "target": 189, + "target": 190, "source": 52 }, { - "target": 191, + "target": 192, "source": 52 }, { - "target": 198, + "target": 199, "source": 52 }, { - "target": 251, + "target": 252, "source": 52 }, { - "source": 160, + "source": 161, "target": 53 }, { - "source": 161, + "source": 162, "target": 54 }, { - "target": 191, + "target": 192, "source": 54 }, { - "source": 259, + "source": 261, "target": 54 }, { - "source": 278, + "source": 281, "target": 54 }, { - "target": 278, + "target": 281, "source": 54 }, { - "source": 161, + "source": 162, "target": 55 }, { - "target": 191, + "target": 192, "source": 55 }, { - "source": 278, + "source": 281, "target": 55 }, { - "target": 278, + "target": 281, "source": 55 }, { - "source": 161, + "source": 162, "target": 56 }, { - "target": 191, + "target": 192, "source": 56 }, { - "source": 161, + "source": 162, "target": 57 }, { - "target": 191, + "target": 192, "source": 57 }, { - "source": 161, + "source": 162, "target": 58 }, { - "target": 191, + "target": 192, "source": 58 }, { - "source": 163, + "source": 164, "target": 59 }, { - "source": 164, + "source": 165, "target": 60 }, { - "target": 193, + "target": 194, "source": 60 }, { - "source": 165, + "source": 166, "target": 61 }, { - "target": 166, + "target": 167, "source": 62 }, { - "source": 180, + "source": 181, "target": 62 }, { - "source": 250, + "source": 251, "target": 62 }, { - "source": 273, + "source": 276, "target": 62 }, { - "source": 274, + "source": 277, "target": 62 }, { - "source": 167, + "source": 168, "target": 63 }, { - "source": 168, + "source": 169, "target": 64 }, { - "source": 180, + "source": 181, "target": 64 }, { - "source": 181, + "source": 182, "target": 64 }, { - "source": 182, + "source": 183, "target": 64 }, { - "target": 171, + "target": 172, "source": 65 }, { - "source": 181, + "source": 182, "target": 65 }, { - "source": 182, + "source": 183, "target": 65 }, { - "target": 172, + "target": 173, "source": 66 }, { - "source": 234, + "source": 235, "target": 66 }, { - "source": 173, + "source": 174, "target": 67 }, { - "target": 174, + "target": 175, "source": 68 }, { - "source": 208, + "source": 209, "target": 68 }, { - "source": 175, + "source": 176, "target": 69 }, { - "source": 176, + "source": 177, "target": 69 }, { - "source": 176, + "source": 177, "target": 70 }, { - "source": 178, + "source": 179, "target": 71 }, { - "target": 178, + "target": 179, "source": 72 }, { - "source": 210, + "source": 211, "target": 72 }, { - "source": 179, + "source": 180, "target": 73 }, { - "source": 183, + "source": 184, "target": 74 }, { - "source": 184, + "source": 185, "target": 74 }, { - "source": 227, + "source": 228, "target": 74 }, { - "target": 264, + "target": 266, "source": 74 }, { - "source": 265, + "source": 267, "target": 74 }, { - "target": 183, + "target": 184, "source": 75 }, { - "target": 184, + "target": 185, "source": 75 }, { - "target": 186, + "target": 187, "source": 75 }, { - "source": 196, + "source": 197, "target": 75 }, { - "target": 208, + "target": 209, "source": 75 }, { - "target": 209, + "target": 210, "source": 75 }, { - "target": 214, + "target": 215, "source": 75 }, { - "target": 267, + "target": 258, "source": 75 }, { - "source": 275, + "target": 270, + "source": 75 + }, + { + "source": 278, "target": 75 }, { - "target": 275, + "target": 278, "source": 75 }, { - "source": 189, + "source": 190, "target": 76 }, { - "source": 200, + "source": 201, "target": 76 }, { - "source": 206, + "source": 207, "target": 76 }, { - "source": 208, + "source": 209, "target": 76 }, { - "source": 194, + "source": 195, "target": 77 }, { - "target": 194, + "target": 195, "source": 78 }, { - "source": 195, + "source": 196, "target": 79 }, { - "target": 195, + "target": 196, "source": 80 }, { - "target": 198, + "target": 199, "source": 80 }, { - "target": 200, + "target": 201, "source": 80 }, { - "source": 208, + "source": 209, "target": 80 }, { - "target": 243, + "target": 244, "source": 80 }, { - "target": 195, + "target": 196, "source": 81 }, { - "target": 196, + "target": 197, "source": 82 }, { - "source": 197, + "source": 198, "target": 83 }, { - "source": 198, + "source": 199, "target": 83 }, { - "target": 198, + "target": 199, "source": 84 }, { - "target": 199, + "target": 200, "source": 84 }, { - "source": 259, + "source": 261, "target": 84 }, { - "source": 278, + "source": 281, "target": 84 }, { - "target": 278, + "target": 281, "source": 84 }, { - "source": 202, + "source": 203, "target": 85 }, { - "target": 207, + "target": 208, "source": 85 }, { - "target": 244, + "target": 245, "source": 85 }, { - "source": 202, + "source": 203, "target": 86 }, { - "target": 244, + "target": 245, "source": 86 }, { - "target": 202, + "target": 203, "source": 87 }, { - "source": 245, + "source": 246, "target": 87 }, { - "source": 203, + "source": 204, "target": 88 }, { - "target": 204, + "target": 205, "source": 89 }, { - "source": 231, + "source": 232, "target": 89 }, { - "source": 238, + "source": 239, "target": 89 }, { - "target": 204, + "target": 205, "source": 90 }, { - "target": 223, + "target": 224, "source": 90 }, { - "source": 205, + "source": 206, "target": 91 }, { - "target": 205, + "target": 206, "source": 92 }, { - "source": 214, + "source": 215, "target": 92 }, { - "source": 278, + "source": 281, "target": 92 }, { - "target": 278, + "target": 281, "source": 92 }, { - "source": 207, + "source": 208, "target": 93 }, { - "target": 207, + "target": 208, "source": 94 }, { - "source": 244, + "source": 245, "target": 94 }, { - "target": 207, + "target": 208, "source": 95 }, { - "source": 244, + "source": 245, "target": 95 }, { - "source": 209, + "source": 210, "target": 96 }, { - "target": 210, + "target": 211, "source": 97 }, { - "source": 232, + "source": 233, "target": 97 }, { - "source": 213, + "source": 214, "target": 98 }, { - "target": 213, + "target": 214, "source": 99 }, { - "target": 215, + "target": 216, "source": 100 }, { - "source": 229, + "source": 230, "target": 100 }, { - "source": 217, + "source": 218, "target": 101 }, { - "target": 217, + "target": 218, "source": 102 }, { - "source": 230, + "source": 231, "target": 102 }, { - "target": 219, + "target": 220, "source": 103 }, { - "source": 233, + "source": 234, "target": 103 }, { - "target": 220, + "target": 221, "source": 104 }, { - "source": 235, + "source": 236, "target": 104 }, { - "target": 221, + "target": 222, "source": 105 }, { - "source": 236, + "source": 237, "target": 105 }, { - "target": 222, + "target": 223, "source": 106 }, { - "source": 237, + "source": 238, "target": 106 }, { - "source": 223, + "source": 224, "target": 107 }, { - "source": 226, + "source": 227, "target": 108 }, { - "target": 227, + "target": 228, "source": 109 }, { - "target": 227, + "target": 228, "source": 110 }, { - "source": 241, + "source": 242, "target": 111 }, { - "source": 242, + "source": 243, "target": 112 }, { - "source": 243, + "source": 244, "target": 113 }, { - "target": 245, + "target": 246, "source": 114 }, { - "source": 249, + "source": 250, "target": 115 }, { - "target": 250, + "target": 251, "source": 116 }, { - "source": 254, + "source": 255, "target": 117 }, { - "target": 255, + "target": 256, "source": 118 }, { - "target": 258, + "target": 260, "source": 119 }, { - "target": 269, + "target": 272, "source": 119 }, { - "source": 260, + "source": 262, "target": 120 }, { - "source": 263, + "source": 265, "target": 121 }, { - "source": 264, + "source": 266, "target": 122 }, { - "target": 265, + "target": 267, "source": 123 }, { - "target": 266, + "target": 268, "source": 124 }, { - "target": 272, - "source": 125 + "source": 269, + "target": 125 }, { - "target": 273, + "target": 275, "source": 126 }, { - "target": 275, + "target": 276, "source": 127 }, { - "source": 276, - "target": 127 + "target": 278, + "source": 128 }, { "source": 279, "target": 128 + }, + { + "source": 282, + "target": 129 } ], "directed": false,