From 368545f795cb313fb2cd895a8f3ab4cd9d101380 Mon Sep 17 00:00:00 2001 From: Redmar van den Berg <RedmarvandenBerg@lumc.nl> Date: Tue, 27 Aug 2019 16:02:01 +0200 Subject: [PATCH] Add gatk test file and rename test script --- tests/cases/gatk.vcf.gz | Bin 0 -> 7761 bytes tests/cases/gatk.vcf.gz.tbi | Bin 0 -> 172 bytes tests/test_dummy.py | 16 ---------------- tests/test_evaluate.py | 31 +++++++++++++++++++++++++++++++ 4 files changed, 31 insertions(+), 16 deletions(-) create mode 100644 tests/cases/gatk.vcf.gz create mode 100644 tests/cases/gatk.vcf.gz.tbi delete mode 100644 tests/test_dummy.py create mode 100644 tests/test_evaluate.py diff --git a/tests/cases/gatk.vcf.gz b/tests/cases/gatk.vcf.gz new file mode 100644 index 0000000000000000000000000000000000000000..df6722dd81a3e3c9425a9fadf40e965ff288fd60 GIT binary patch literal 7761 zcmV-X9<JdZiwFb&00000{{{d;LjnLa9<^NCavMjIeWt$xVL$9l1k|eg1zOAoElRRt zDe6MnbBhTCNHj@=L4XE8+nTT6lXdBCpj&Fd6Z;?%L?OGn>tvohnN`qt-`&m^<?U+y zFx`%SdinZSRq$`V`|kDoPw!5z#{d2HY<zk)ym)%JF4yB?c=h?Q9KT&|%e!(tJS#Uh z>-pn$zFLmIJzXry#r$S4D<8M_n*lC3*xr|eb$L6OE@y-3VvB$HNVyrzmtknNo|Wst zV!qjyvv2>hJLlPlJSYF!oR5#|)$eoMIZZU3#mNVMqSu?jZ*#nU_i-$=#$55wYkY`4 zz8Or{<p58<E!XAj&->HA`M5Xz>&0ri-JSeRxm;~AHu(57UCg(iTi3YSyT;{qJzw5+ zUu56-!+AD-^OfhnAk-G~f0wi2;KMy0<k{wCx+rH@)Sm&y{c1IfNZh4q@YowH((QbS z4<#6Y`C#){-pp_3H&eQ9xR<vVuir=XU*<bG?L)sD+^n9KTe@&H;V{f!PH*l9r{`D6 z1%`R<&GaGr0dB0}3ib4__onYaB2E5!UH<Y^E^j{n$8%mB;(S-D?Q}6%h7S%_w@vT9 zq1Ur|5q2;C{d7}){4`zueEIay?`Ho1rYrK*_3B~p*ZJaR_4{f$Si^j<d3qSYmA5qe zDdOv=O)>a{c)P{~gI|JI&i5xj`#w*}@|4p(pPXG^UVIq5E9q04`}yPE1^$D`d@<Op z9^xAtgyHAG{q)yza9x(mDirr7yX?K&zrqR(^O+VQ!TU%s%vQfG<LlY~>|L~6Z2H&z z^?G_kujsm;!VSM%Ew|J85-~iue4H-B3tX$DPwlPp>&yNLb^`YC*{{>ZQ#p9r&{bd0 zH}~cGzcz#Ncicq8c)P+HZcF^R*{+FEujkXfS%3KSv^)fy^XdAIE=dnGzk6CeZ3aKA z9#)U5#peh3GOj`o{nK<kpCX!5biG{;UV%X|el}m1H*KQ$VfwfT?tV%B+W-8ReM7;A zfUU5b!(9STPq%}I6+QKDtHGwcTbFmBCalbTo@DRQ-!8A$WjVt~Uard9+y2>W1CSOh z<~Q>a4+}=X$^^fi#gL=>%7MJYtB>XSS#+?$Z};%vo7KZ3k_uAY9lV6@9hN^#*R#*h z{#wqLxSc!T38LM0u<W}763(Z0mxFiH$49#7X3l9`gdI*?dsiC#qg=0|qwinzXsUYM zPgTDrstVljkZ$D*l=bnff9mY4aH`w+3_evp&kuI}8Bq%I*}W(|>4(+k|AU5ip6zi3 zrzRSy!u7``GHp4_nSVY0cDkG{J|h+;^!Ipw90}#|X*(g6HvVUz&f<Fd@UL+B$@SCi zZMmLopaI6qr^RCMc)wbflY7wY_~QNK)m8jhWY_UO|AmX;1-NE2{^!m8T8NPnOfV_< zGbOaJUdw;MOK%oWv+`~A(QxH#IoXuk$@&R*^W(+a_ZRj1;Ep@MI6gf;C)o#H1`POg zNXq5-`E0tKKCcFzKizC@p5NbbA6^yos<@qQw$s6SHT!<PoYCzrSIaZ7-ueN&i@3Z2 zUbEH1B_elxi#sc8n#as>{l5SE<nr|0hx1pB6foJYCbx;u<J0Vp-K>5s*VDTaT9F}r z{xW`fdj9gq^V6$Wlhez|#jBUEE-z0%{hfy1OdsDrmdg(%K-+#ED^@U~f*_#rSm&#| zIKk7Ve7~M!L8gn3*{}0Gm$%cW#Wv`GYUmjsU=ioDKn@dd{QPFSEH|6+<=gXD7gvKP z-0p1x8KPsdy$`-H7U8d($EWI@>Lc@2`#rc$aHYvTRuXwKjsz#UTd$rTC)bPV&Ce5( zMaiK+lHlKx(8#Z&l5M{8Q+4Z6>Rx?%`!Xo6)o+u^<t9Y!lP4t3S^0Z8qkBm!vRK`b z>r58qujOJ)K}v50)X42Y{y!IZV{l(C9+OGt<Vk~*Cl9L`>`A^nPT)Flw)EDkzn$WK z?#Y!_|NeYOLHT(M=olY7-mE9(dc9gt)?ocHFW9SRDqo&%imos|dwcov{ZFsp1FNSk z#bbE))w{P>lb=pMy*<6Sx}3bfIR87cGZq%CKY`1_uQDA4BsTY}1%eS%kwT;;l;l2# znUoLn?Y6R^Z1c1&FVEkfr3Fd)bNr`IF?@(om(aRnh_0vg4TbCX?|-~F`}2UNA&oeR zg!_3CsBO$GnE3a}EJ>O=U|7NR-~k>xpT;6D%}Ae~60t0?aFaa04%a7)7n4vAguAcT z6Ca3!duZ@NXj=I&h$su{Wcd`FDn}cqN6#PF2}7L+Y8gVD=qQ}=_UtW<rTY$_J1eiB z?(!RVODq7~<}^J+iqq@O^0B?1;qAPr`>Ynk^Mp(BJba8^MQZ!?=JwCzJLA>*uDD*W zrZW-=^DRIuHiT;tynMdGUr*PY0%h4>*U&@jOZ+(im;dnU^z7BAKW0wCYBO0a7oY!_ z{frhX#$q6ahm!bd(|8AZcol^W(L%}~r#o-PBm^hlKwjP4Q)mS*1U|tN)3Sey9C;Iz zX-L$_Qhc~5mv`Iyu}~faZ3G`qHyhL^cpO^WzqBn_HsMlhsjRZSCcI(7n9<zyn@sK& zj3IxD$vPyad<n%klfCAmo>4l@o!7l4vNGYgLGP^VHBprbvr>7jSg(n$O@x+QdeLiQ zvWaI>E35mT-ewEOxz~dC-<F-+Z0|HO*JF*9am<PSMX(tuxrdar${%5__RJlZ2i<rw z*L&vND(#M#-t7AsM+WOPclqPiF>8hDpU~wA?f!(q^!tL(pVpoW=Uu-q_&lH6o6ig5 z{Nd()e|>QWrr$oAbI<q9nNhO;?z8zm3m^`Z?`08Dvzz4JOWQl4sMyU&1gI9h9w@4S zQH~3L0Mm%fZJg~rm9WP3`hcjON^{ASHoePGGn;ozU@dwdswOrMHdMH;UJD`+ze6C} zGx&7BVS;#{K0w;G-$5($efq$itA78i>DV#PwJ^O)V`_roSQ}<~)1R44A7?u&ql^&3 z_fG|p_d6L@h%2fhZ7(gE$?W=?3}rndvSxq4m<cMaC0VN4(n(=@0|%Bjw`8d6aF9~= zCjc|~<$6Jerg5)uRMl^|*`SJnWtHxrVVv#<#5OrhI@AAjHhWns%OFGcFKd-0;F8iP z#9*&Qo^U@cNe90^*m^(5z;|pfKl%D@x!@Ai*=v#KYvmZmT>pI8x|L!`h4nscm7(_8 zlsgN-bg=n;hLVQ&r}KC9z3TqB{BfC?zoD|s0SsOn*22~qr+e?OHms9?Q9||-V{KUH z0%<KjV#b<Vt<b6$h}y70v7pIb2dE8g!;q^Zp0N&4a8G|hQybQ?Vm<FCyUMWE8JH1C z`EbM9b5MQx!*5mWc#b*TSA}fkNTm<QQK-h^7~e~>wPEEs2o&Bg+G@kvb5Lg<fky2( z!pVNhsSRt#5fJ3P8ND{F9mgXF9c(8nzfsco!$FjFp>L7#4!5hN2kQ32-Bi}fLwj|^ zeb%~z>Hfr1+f^<jl{x}IS*s2x&?D?>&7obt;;8K^hmoj94reG?`HPfJ9DbjQHJ|Tr zKU1|<1(6?)yQ=*~YH`GTwZBMtguf`&SWj{IdaBxAkbsYvuksh^P#7I<XY2V?zx=E1 z8vbIBh+9?ri_`t-<!_CoA}rJWWzgU2Di39Kr-Z3Gw`1yMj-Ee1c$?~*?_U1!>HWKt z5AQEe-kzO&diDC`^!)1N<BzB3C$Hb0U%mQtLK_(;@l?;r_vfb<|9f(A3i7rt{`Nfm z7WCUA?$U1`e=GJ)D5XzmYYy*<ll6wEAKAv8yg4~N5nl7c>yzT-^yOH*IDJj6dL>_+ zUO*PVXr42{fVBF>+4nEbK2Y8<FU~Gs&_Rf>3G(9gW%$ir&=G<$e(mh#BX7t;Yi?eI z(^RtX0<Rc8_pWh$D7_aS&&Im&?8WZxDh5$zpS-ymot}-(K8)Uc9DO)HVb9sf4DoNo zUBT3d4WV$5OEJJlsCb5)*^$s}<oyUG1|PDaWXi{36wS6CW%x8kwUpti<X9{ZvY}yY zsI3%f3>R*UA*&ix;c-l}9Xp6a?Y!?A)7B4B#Hl7ghD!j%xztYKX-{6Ba8haolMx^j z0aC(u0mMxKPpA;$u8ENBMM#MP+CCzL*UGdI5(VRo5b^~+w};R(URcT75D^kXHR5(; zL}AnjRVW(@Z{iRD;{q@!L#raFb_a@<rzHbL(oYqa3^{g)w3keQr_<RIHH@e00>N}) zTt*5Ei*R8b(8s4-9A*Z?x*Ep9(=e*eWQ}2AFf7Axyg1j=Hqo-0;^rlAlBV!3Lc}9P z%<?XXo^e&65RU)}+XRUB0YqVgc@+VKlDCi%AkhHm8E1uy5YdiX8zNg^F}uPDk`g{` zM_>cy$U8Y4l*J(^OcmTdqCjFoA$-ChI0jx!@e&77Rt({uA}%9b+|ZyNS&e(gPYN>Q zQjp=Z0V0!Zt$?l6X;$Fv!-Wm*Nrz~{eW*MyB8vrE3AeNuhK2({H4&1G#dv~*w9F{* z-J27O@G5j-JQvZ@-DsJD36tXq3Hxh|$6N8l5D?BsxE$_Lwc_a+A`Z9%EmindK)C`% zVP`x6gAo@a5GH<;*K2#iG2pR;5?mOa78!8}6+C3oNXwzu87QfBj20aS6(LuV6NfgT zVMIp?H5Eh6N3OIijM~hkZ0KuS(j0VYmghA?$LVfp)QWInYt#-*))E{5N2>xHa!4Qn zG$T_S@uN}^=5R_zT#U-N@LgOWBMBD(<U5ECbxkKQ7DOk2gD5O0izua;P#Q@?L(m!T zba!GvLh;5`$wo92ksey(Y%>w%(9=#1ZCv&c7-W~n4C31=D?9DH8iLAX6A{6DtrXbF z6=fo=wJAW52Luj?O1oMr_@j6Z^=07%i-|}e_H8?G5HBid1rVO(ZLwfvCqQxtstI1I z$qx~tA#Si#l5&uxA*G2(z<J4kC548~^v9vUiNm<EDN_h!3ZHY7Fe(P5&iPHJ{AOcJ ziV!+lG>2WjILO_^;WRhO1rA-41J9elpTfBvtRQ%R6=O={6+wlY(S;U=h&rsFcJ>QU z^{YAtPuE6es5EOsh13WV4IJlMu>@KKcXi|=l3%o{1v8X#)KKKqjO+sjuc^ufM&Vf% zU<JdHrD0A$!kA@2VTZB~Cxv=bW6{%?ZiWnr4V31lD-6RJa#5EaM8qRZjc-DvLO5{* zyaT1a0Z7InJwVYOAnOng0hg49AOMbnNF$oOu>m3mKWGbYRVXW>@HQ>Y;0tK7PO`&7 z(mI1oM>=Ra(~)Sn_BJRw@2W*s^<qiV60AxjZs<fK2?f#RTn+{cWNY9dz``)oRi5!R z?lQ)q6e)wDX0`x{aT@h1Xhr5k6C}D)<~*ShMpq}?Na!qiU_z3%A-=y_9*KRFEQA>X zImQuWhwjkc&IFJQ#7mp-NZoXV5vTHqLnHZuh4qGjwz?||kOH>cB~}9kLN!&SXKdgl z--w2WAT$C=37Wk-VIjF-bvvsg$&r+V=}u;pAa`YgS3%7g5AG$g03ix#$k`zoNfm~P z8lB@}NV)>x#1O)3(^hg246%0Qh#|QTkAxiP%tv4m_;p-j8LCqm5|SK}uuTgE?h=sE zHYX&wa>ORcX32d>DNZL(hi`jVB^+={sB3fG6O5~2Dd@78HX94dW-=&t8dk#X4l7t! z$xhOZ(bP&eU{YY5Tozf4)jCLO3le_nkk&$5Xj@WgjOW6&qqTq%8ns&zbgkaPGnBKC zVte>R9OTH1G>W?s@`}j_Q03_|JeI{Z>KZdLBy8q7gtE~gGl|4~4)Nggn6fm6G?!#0 z#9KvtS_Zsr>T*DfTc6{`vLI_@Q_CqOiPtaz%5w=)J3Y~OyATN%8F4`|_61r)iB|=V zFkh2$ZD=qLZ9HdaSi7fyEGT5`DfpH=khr$%3LK;Y5?fKPg_eQu)Y8z{s6I^A%e$Hw zwUM<_chVMWC)(-iyry$Q>$plG!BVFn#+wdFOX9<cNO>%*Qhx)LS{SxX1*brCKnRov zNg<<zqs`00@Ixe=ED)1_RKaOl;*f4YZ43o#SK<+o+VHBCc}T!@G<_-oXJb1-Kz3=q zn_J>~%9yF8&kKq_A=;{8B<VQ3Q}fjL3RI@nb>U>8o}H69z1SUmy&gJSr5?X~O+R!- zH#?__QqFoKokwA0gET-efi({6n>JF9^7_aKQiWDY8E{Jitu5aqYFGjYI&L2ylux`P z^)OHy{*7D_9vKguXb5CxbtZVlLRZYV?CS^$m3wB=#yUR-5Fl>kZw0Q?EH8fSF2xYR z7%)`Msj7k09HX|6Vmghf-7#oKWd8;*8s0?Y#x$Z_#g3^$@mKg1Iyyolwe3~~Ua2~C z+_9<v`O-&W4L^@*$LxfTj?FmaAKOu-O4gyQ*@(oI3$>4qwU1=J0s#N2*dR$vY~phv zqQi2z++!fkYg6h5FLCCpN>ha%jubKvo7gN*4No1KIH+4(-4B(lmJKZG!7TZMY*6*E z4Iz!3*o4a3diQ}I9Er$kB#2z7s+5SF=`2*iG{Svh`=Pg|;5@i1ls?ecxRP&0kuFdj zeJKwiTiTXDpeILN3V|R{B;n~pxuR^vDXiIJfrN8auG#GcLBJt#*Yde3t+hbv=Xsf_ zD_xBWV`6P9lT45*R>oQS#i~rXR3J&)ESyCQZsB*D37YD(3M3r^$`c12M0JKbV+W?9 z7S&1AJxV2v9922(9U`b)!bl>%_c20UBszD9+L>sttr-`@yWq5%r&xmz+e(me78*CS ziP{MvkTJ@t&M{;qqg-co1ZGFN$t}V?nGq_doe~u!$7^&{5DF4<q(B@VfQ~AxZqz{R zIJPs03l`Dg6r~+P0X#+Mpk#J!bc7L|6)HV6DI}X(J-{bBG-@JO?Jd9!gV2^f=?IiA z*5y)#+LS>jsLtv(EVOgA4s4{A0FX9ygJ`Eu>qskn0a)AsvwjC4(Evz>P7iL0Rh^zb zfXsgY5Vf+M{3<|Rs0fg20i;pEc6Q<_fK&<~<Rfl70!T0!YFMLv4Un(lW0b6vsQLuE zi)o-OP-J&E%c;M>L}vscanZ56?p<4ImhTip+G|Qe7P5gP9^wbp5QHFBCIVzrm5OEV zs#2*GG&RBL<63b_6~fz~7~rp63<u3`SZtP>7DTo{J=ImF3Y4ju3RPLxs6l1d``a?K z;VkRi!VaBJk%`Zkp{v_$@Ix4SscOVc&5A)q@TwV(KV7{O#F&qxIbfmH>pLQGM=f~Q zY#tce-;2!yr?FfrZw9f5i~T~CN>Yq>fzrRLs%YO@>bM;wpy@kyRSW5T?7g{G-hqki zL<P<XR3vS!=1n*U;xVFXtFe%wc~zfEl6Q6Tm(~rXbEZ<}JnT{gzO$h$0#>4G7>PWF z>+qofQY?&V!nX=i?(tb_orW=$>|%9F+kkYGswvxy3W!kK0Iq-v<B(>kfE$YdhY-h| z69F26B^Cp3s`lOwwUFGG_7kXP?mUAakV+t!6B7-o%1iZ!Bf*9DQjLK%*V39aW{2Tt z8pG1**inv&q9O1M#Fs(3PIfxj>12n?XwCD!wit$=SvgRzQ`=;R+7~}Xu9E`4i5)ah zHl+lEIu)S2#QOpX#vihkKnDPvt`G4@(F8~I5*S2)b3URYI}h>50@YoVXrjGyN`4_b zHxm#TI^iS2Q5+?5=^@fdjevrU2Shw`b=66dm*))~Ii#;@HbO?`f|b4|M_vEaXz+BP zLsj~Lg&lG3^3jDbEzhB0<933`+X+=Ln%JQYQW}E12^H!u>gpIq-=B$115RHts&=@) z-Vp^<R<jx^dsG}b+P8>@4}=7NwOiIugO107r;U^?R5wsm=>e$J>`G-`OhUvf6}{8$ z*pR+~l}`ci(3MWy??Qi5#I1<5WmIPdYr^)7;(6nTO4mBA*w}TBiY=w2dgF(JhW4VA zsf>NK!OE=)w!mZuDHBy2lpfqVs5YgccXI#j>!QGitW@Fj)h0r#=<%*_ipTyN!dSH> zfiJL(F`^KmEu4(&l+_6*hqa}%2$U+IhkQl)6?F+}b<#QCyx5}IenU15yWLfZ=~5RV z*lNM!esoIk5YoOHq^&(tgN)Z0*Rmoce`$O^f8v4zq+4=vFT~4Tkb-j4h#=Yvrpgov z5~%Ql&c?<1KBJ_L)7D7QK*6#ou61GBX~`oOinPPff{Bi0b;T1oHQ@v}DBNW{!D!NI zYa+V31dSV>!pX+jfJ%zd3B|A#qZ*^O9*z4DOw`d>Z#Kx#VU>B42W3qSln4hUq^f51 z5#kU|GBzN*_q#w*uN1<Y(-EZ5=TRt54%<ehgPe^sj7Q?+d?Brh+YQ_ou6-QA6x^HG zEVZf~O}cv<m7>WJN+ioEYf~K)76fS@4iQ#MQ;j=PHXQp&0-_O3g1YM@wxP|0bFFcO zKzBi9kY0<05FJUciIHr>S%y5wDqfUyfQ%#~VmT`f3My1hVcep22t9F&alfPQkbVd* z1)-;*HN(E0^o3!Eoo9q3bYF_AeL<mG>?e_#a#O)k;I#FGjZFnk=N(&31yW;Q!wZ5% zg>?b@r53lMQvr0+bR@wR)W?qdcX7>#;CMsV2@kAU+I)yDv_@NjS~e=T?x-8<D}heD z$dLh+tHAQMZaPTqhRS5!rxcw<sz-&ztO%c7oWkxmbIf~#4E5%OQ@aiec-T=!PRHf` zSXe=`7_B4wAuXzQ_Y~;put|hMthEP$4Agy)v4Sad4}vN;u1ZlV(~PDsbJ9_rT-utQ z!-~!X{?dXshGo@vurj6^n?xbjx<RTaLuu4;x0AxE8kKUG;x>&^gS~OMNch<B5N-!3 z$10#w2h|v`Q7)an4wQwI=QU7A9&H57%K3a8OSFL!g^1g&JZ$#uAq8yJ7iAwZ1&@hc z2c%R$ij#V%mPFJ43m|JXvEOqRNik-qRC6m(h=`<;tbz)<yw(iyEQ88G2Y%zVE{uWf zOa^*K$sm&kYbl%iM<_>BsBR%9#DSd0vBowPx5qVrG&LkN&&6%WULLNkd7``4IRyY} zJC9-!(cn2pmeaQ&!sbiv^iT;a@H21Ap~5uGr8A<`ucM!-gr$C+Oqqwiffs%@$8yPR zX<H6Knafl&)T%lUPFo^~n_MzO&d|3}vK_1KH0Si?h+Q}=)v**ah8>)O1RaY{=M#1S z@x5CPR2re)&bm&uo2GPP-?((7UbG~`GmaD{x}{xQ-H}01+OUR0mrr|u(CE~>425ZM zOb7~@RuBTfP%w<ZyNb2|Ef$!wykBaDRCxk5#pz%QnbL4xW9nhRBT^5Enr=$x{Zv%I zk}jRTPeC9nZ!@YkSPGpjEuC}R={``?gs2T1N1aGgc1qA`+e5nB$nU1^or)t*n$q64 zlc~KOj?1@q@Y3!Ol~nE6QMy~&g_bdDctjV~Nz{U2s+nmja>JSjzUP*8r31mJU4-17 zw`2E?bvr{R+VT<Xo#O}qVQ9GPJc<&IiW|i`uZdWPZfk~+RW=o~bh08I-A0W^rxna@ zFO^o4U=sWBgutnHA=q&lhy;n#GPbx4s-CVdJ>ZrD$Jx#k!5wF@iswa4HbctlIoD8A zCY+Q?JsCmY@2L-ti1<|^mHS9KejzA<Wx0v{9@BAUIKYtp%PLll)fF)!je@n&K&=kD zYKmP+d*e+&=qmsL5L|`xaB>ljMS{TWwijGC9&l$bzko;Y`4LA2PoOpWc6#n^@AgyE z>0eZ^EkOX0U|R>UG8%$<H2KJYq775UtNjNq)Az4*Qg<vi-|{)cd~7Oyd&`E+SMF`| z*p#a@oO7W$XRCuKB3Z|-jte@BPPjFOZEYP-NJ7muVX?l_)2IxnJ10RAh1uImHMSEJ zIyO}un`c%X1&W0tb|yA`n5U>}mpzQD2ui!`0papaP{PSLL1zqUmw_JLgxrEi%pFlV zL~2gmXJn!Hbpt7?^bmg5o7&pBlvhmH1~k<HcuToGXROdB_pU)$Jbe#^j%A0Rxj*ZN zh^}V$&AM<_Cu}faF`V2vxC6K9T7Q(CLZ+P0jw@!2bcnooTol5*a2>}*A%j)S_^d@x zM4YV?ReANPcOAbg!l`maiU^Xl+IEgoSQkDBf}zsMicaOmLY~w2<RnFdR%sC4Y;uc$ zA<_ua^mPZ4a6EoNM1whTX%W73RNYLvEHsi0U&N$x<(vNlck3~dv0(rJABzYC00000 X0RIL6LPG)o8vp|U0000000000I<x8V literal 0 HcmV?d00001 diff --git a/tests/cases/gatk.vcf.gz.tbi b/tests/cases/gatk.vcf.gz.tbi new file mode 100644 index 0000000000000000000000000000000000000000..26f3289b70db11b8b83bdb80424688968852fba5 GIT binary patch literal 172 zcmb2|=3rp}f&Xj_PR>jW{S0qUZ{$5>Aix?hyQxW1z&C->JF-D<<z`1IE-kl2hP0!7 za}_o^%=*4`xBC6L%YO?bb3B}oy>&^pT=CidaIe$PD%XY2kB)Bt7xe4LI{W|6Z*P8D k_3Q2g70*ey(YvR!^f%9G-aadkfk7V4J<<%!VCR7d0H+&2hX4Qo literal 0 HcmV?d00001 diff --git a/tests/test_dummy.py b/tests/test_dummy.py deleted file mode 100644 index 11b4eaa..0000000 --- a/tests/test_dummy.py +++ /dev/null @@ -1,16 +0,0 @@ -import pytest - -from vtools.evaluate import site_concordancy - -from cyvcf2 import VCF - - -@pytest.fixture -def known(): - bla =VCF('tests/cases/known1.vcf.gz', gts012=True) - return bla - -def test_evaluate(known): - d, disc = site_concordancy(known, known, ["known1"], ["known1"]) - assert d['total_sites'] == 3 - assert not disc \ No newline at end of file diff --git a/tests/test_evaluate.py b/tests/test_evaluate.py new file mode 100644 index 0000000..df922e7 --- /dev/null +++ b/tests/test_evaluate.py @@ -0,0 +1,31 @@ +import pytest + +from vtools.evaluate import site_concordancy + +from cyvcf2 import VCF + + +@pytest.fixture +def known_concordant(): + known = '/home/rrvandenberg/devel/vtools/tests/cases/gatk.vcf.gz' + d, disc = site_concordancy(VCF(known, gts012=True), VCF(known, + gts012=True), + ['NA12878'], ['NA12878'], min_gq=0, + min_dp=0) + return d + + +def test_total_sites(known_concordant): + assert known_concordant['total_sites'] == 37 + + +def test_sites_considered(known_concordant): + assert known_concordant['sites_considered'] == 37 + + +def test_alleles_considered(known_concordant): + assert known_concordant['alleles_considered'] == 74 + + +def test_alleles_het_concordant(known_concordant): + assert known_concordant['alleles_het_concordant'] == 42 -- GitLab