From 46c217fb27de33165a788397b0d1faef4f6d2039 Mon Sep 17 00:00:00 2001
From: nkindlon <nek3d@virginia.edu>
Date: Sat, 2 Nov 2013 14:57:54 -0400
Subject: [PATCH] new unit tests for intersect
---
test/intersect/a.bam | Bin 0 -> 1780 bytes
test/intersect/aVSb.bam | Bin 0 -> 1778 bytes
test/intersect/a_bgzipped.bed.gz | Bin 0 -> 85 bytes
test/intersect/a_gzipped.bed.gz | Bin 0 -> 63 bytes
test/intersect/new_test-intersect.sh | 21 +++++++++++++++++++++
5 files changed, 21 insertions(+)
create mode 100644 test/intersect/a.bam
create mode 100644 test/intersect/aVSb.bam
create mode 100644 test/intersect/a_bgzipped.bed.gz
create mode 100644 test/intersect/a_gzipped.bed.gz
create mode 100755 test/intersect/new_test-intersect.sh
diff --git a/test/intersect/a.bam b/test/intersect/a.bam
new file mode 100644
index 0000000000000000000000000000000000000000..bf9f1fc23f810b8d8904396b991ca91c9dcc602e
GIT binary patch
literal 1780
zcmV<Q1`GKgiwFb&00000{{{d;LjnNT2BleBh?P|oUUNom&KPP@l28#5L<lqcewh?z
z)HJEFVr)=JPBU}HnWnkrTo8*u2&2d>!3-3eZr*x|dI>A^%_Mv1A-Yg+5+nn?lqI3<
zgJ=Ke|L0#%>mN8Ud+p7~+Iz2Wt#7Y?T~~Kg(>(y4k8UXUZt37{tlaZN$ISTD#O~?A
zYHR0{o64IvbgX-1!`6w3(Wy%RV0G)ny1ud6+TK~w#%&vG?d*A~+_R-)V0bq#clC6%
z4^NB@wvP`!FRSg<zUjX9Fj1L$Ztvh!`|!+I-+0^b5ch3EgX3XY`SC3sLV2N?5j+jV
zC`eptEtOR^jYRS(ATee%HyNbL&?sX}Fj3jvH(s3>qaZ8BnM|V!nt#%1?z~PTNfHtd
zyY$YwG?F4AF)NkVDiaWzh9tD)(u*{bp&)rCwX!;M;Fbd9IQLrc%-(EeV0uVbhWjRI
zF=iZd!etTb5he>nn6%2{sHy;}<8tBd>!`W_Dz{2I+*wnEWSp@~I&l;f>zK7dW!!Vo
zJzK;*VKSlZDInT&;k?U)wvXPq0^WIHoJUc8k(&<wnhfksRPK2ZRA!Wv7;2Fm9a4t!
z>C{ik(Lzb?y|iif1gWkW84{uvY5ECL%2AFBpRF1qRyEEFoS3l2r9)bfPE2#jls4&;
zr`0s?m<Y$8J}_EI^PxPdu(>o)E$e%7*4Lf|`e*R!5nCkZ2Pt4Pq)kwZ<ovLGmEkO{
znmgvX7AAcgq+-r357)s=THcvT-(;118OjJDe8#cLz{H*!9w~s8^;%k<!K?Q7*YHSz
ztV|G#nM`CVVcB>MkQC6$I4O+G0uJn|;gRlEIVojEZJ5dn{i8KJ(iMjd2lA>6-qh6e
zL=BLn+rrzz=}g3sU~1J%I+Hm$8!1*u00jk@xnPn~2{e^Pg`AfLitZH=QqqMKfPof~
zKY<}OU4lsRhbu0)3?()V6um1Xe;9Kb11+R}YELYcM-{9{)#c8H_YH<xL{mw_Gn%)X
zwqc#ggHMk~nL&IjnXE2o;L+aC&N%H*JeuM#N=T35(S(PzmgD}>YOobrX%vrGibLI0
zC>|}p+Hg?#`Oy+$s3Vy+Py>&aZLH@Wy$_nA5B>#-isI4up|$krerZa)#t(z?SJ3$3
zg0M&9Ptf>bz3~#oBk?1oC2ny%8b5TXA}Ab<9>PhB(oN9lAw17KCmVP)diWrZFgQtK
zhmyvl!6#`AW*ej%3P%(9%xILllBV&sSNLXW3ZH3-!;y#~mBIs0()4|J=HhTPb+4U4
zvz{c;LVAxEaY<r@lun?VRWw$FvqMRyXm8&jCD1TdG*(D0aPMfWkP>HwBC#UetUzy8
z(O41WDDEAJ71D8y2MUQ5;VCHan<cWs;^|Y-SmAU=`D`QV9buWy^pv*K+Hum#+S6hu
z04xr_u%5%(^}j4zQXeolde<E-al}EDZ*5r_M+?6Q!pUO*ajk%>BiC+;BMH4beq9nr
z3LaRpD;a6%+;?{}GSKyT|Kd2ZaPiL{7OjY*#y2PU^5a#oGO;iCas+!1L4G!N1L64E
zru=MdWWax?OZnN@M}UXs8hdomx%1?rxHk871Wv4Ito3kY&Hm=NHutRpZX0R%l>6oY
z=f_HLXMJ|zrvo_uUQ=!|xsy{nTH@6TI57WmyjsCv^3(Ipal~P0=(p0!I9k{ihx<p9
ztq1|LKQ3#ogG#vf@TH~ofr9FV>eA))L8Dg;J4)r~Vy@*2J1>`VlgVApUnxcZHLQ8#
zLTi0s;Ng;NT~;3?jZpah)s2u?V>=VBtlr$P#-=1Zvwb5hkJhBM2sdUnH5SJPA<Ujj
zLBuv6yj4j-B+WXo>apmxg`xE?mFxS<wasAiVzNI6v%O!m)K}&@SWtQwVBusEYWD4f
z!5*%DJrZGZT_Lz!Nunj!-GM_#N(&|%O&U1*bAB=jLUkZ<pc~G1LOhjJ?*mS}Udo#V
z{P7cFmT+(%Vpj0?D#Waz=>x=UV9EQ4*}}{mVs`N23}W_h{cPSG_|v)vF>_d4LCgZ?
zkzEQ2zn??Q3g*8;%o;Ac=DgX!Tu1)Ag@xI_n(}4`&;Em$J$%!h2ayauIfIxvTpK~m
z0=_+tm?g}7hL{!fV5>Ep*^ihFeD@(@wy+w*cF=qbF?+c0mAqLoxP)LU4yQgs%mSK`
zGf@({Pa|dp{~krm8de`f%m!KyBW4R5-$cv~Mlfs-yWUCIxDDpkc9rXc?Qs3-Yt0YB
z&A-<0Yre7h_NMyyM*Xk$e;<cX=^!qxCQDc3mkQ_%wc&q4AX5fTB>(^)iwFb&00000
W{{{d;LjnLB00RI30000000022NMNl1
literal 0
HcmV?d00001
diff --git a/test/intersect/aVSb.bam b/test/intersect/aVSb.bam
new file mode 100644
index 0000000000000000000000000000000000000000..04eb46945a6a1e4a5a1147390ba7f9f31aafaea6
GIT binary patch
literal 1778
zcmV<O1`YWiiwFb&00000{{{d;LjnNR2CbP}Xk0}YfdA~aF`G)QqQ!!Wh#*o(=6;DS
zByDXinyRFRwpi9pHcgC4HZ~VZ6`>#%#ab=Y3N1C>dh<#2MXJ)i)na|{LA;=EMX&^X
z(Q1*_IXh?0ImwwP=L8lq{|q1h%>1|c&$^!8mX^B#x*yq4>EF_YJ5Z(X@vfPPsmAW<
zp?Z6F-;<TTEnR~nyK&{QEnS3Z!jT}j($m+~Ino#(>YNyQj@CQtwdq=?Syi2ScJI(s
z=g7=>ZK7jj80(JVp$W4r3qlbIQ_d-qOootfOhQO-EP_bY;W2~|L8Q97Hc@YkM`VUD
zq9Igb^`nYorFjU6Hc8m*QY)!KNUTXhk}}O13L$YzBAjBXNeD?qB#kI%k_R2Qj6e!w
z%?S?Lo2(8_5A*6sZPMD9P)L!WbC=G-XvV^*WV#qC%TN`jp}UW%JVV8jaaD}6BI{%r
zp`dW0B8pT<O2UHFQ%*gZRgZ`u)7BA$YfO|@L8h&pb}FNtCPL{#s5*PoMqd$xwN$Ke
zmWl+UK}X=I*%)m?2Gb#-pV+)cvDQ?Es7IQ4O;8gMISFw;(hNCFVWPvb260&pqs)mB
zDO89_ljezWOcCQE<m7Em(+UxD{2>QsQ<^sEk(teffTF0clc>)%GVUM5i!C-A^Gzur
zgPb;@W@EnDJ`3_$o6Qx{m=h6l8^R_FY_0=|kn&DcYm@clGMFHQ=pe<agN;2AuW7(c
zYEETH<`LEUz(B;KRzFiL21O!KHFFaYkXfLjM1({n1038H@wg2rl+vJTLsXw17>jr!
z-mr0A7Q~yHnr=kEX5D7q7K#U1jC9$HDiL&Y(z)1z0?1T=K^NStRE(Qa!A#5x0Xg-u
zg%njL1QY?43kp;$LsF2oTgC)alVU?aPP+x=4?-$PI|p4*{X|c!kb_KWDzw0Oq<P;I
zLCw}wRN$bRx7Ti%2gSjsCR~{jzm-H-U6AtZ`&kLaLpGW6Y{elMp&^+v<=Mi6N{Nfn
z&*orBI190qlozLBlQ%=`IpxWihgmV_X9J?iBSCGTlxL%j)Ho>3B|KZv8~sd)S_sej
z2j^55vR{(;hbY8<35{z&#tX4O!s7c7IwTpVyrzGcFoT00YAMh92RB)y5S(=nqG(Y1
zNO{&hh{i$B$&_c^L!0ub2pnyC2ctq4Vtmv#n58M*3c=YT9|>NF>rz|ebIpp`tgY}7
z52}e$T+=ZqW5pOAwe`Juwibi4)jd~1wS9_fx&_s`m?Tb{UO}lMg>04~cXow18MAL+
zQyLUv#>{#J<)j$xtXEK4j8`xg-z+KQW|{R0Q;rs+ov^Pg%!@G;X|-cPFToTSdxb0}
zd@}15iU(Cb+v0kMNW_DBO51Haj!OE!o)%ARs%+lSweH~!TN{nBsp`N`eQRS~Z9Mwm
z**{A<u<StXI{{#+`56#~wd;Rg-Wm@O9KG}QHaEh+DmS;Sbfd+e4dLW50JoOF)zND=
zxe<l_9lx}?5rg|%cO@eZ-TUrJMgn?18(8W_5-$Gn{gM@Klzy0kFFslYD--#m&r7iP
z5aefF%MgyQZOPBNZU_8(x|E-F%?5aIF5RPo?wu!>xV5=YC~#s$x>m!HHT%nMZSI2x
z+&Y@-l>0CN=f_KMM?Aau`vRPQw<R~3+{vjOZSHCb9GHK}T`fT{`N{dR8(|n8{<XBy
zjTZOL;l8nCD}=!856jCjl)^oSFD;7)4C)u^%WjDWX+s%ylqyceTwfS=UM}S(le?O~
zQgZrpSo8XY_IM!Rp%QIh9uJakDSY?JMo6Tt=LuI<Z%(DIa|usx-w3xjY0`g$8#9~I
z#jcA8v*$7h*9(L<s~JSn(F4gIby`apUjJex-e0b_29p<){VAC3|GX_;nQL=F=^cQ@
zlZmOxw-OK5aP_NEhsm{$;BqywmRt)64jm~inoK)6aP+7AWMYKuK%zh|ob3j8Dwzfe
zoOrF27ZdpXM~|4o!F?VvgTGdJ#2i}Q_lO0wzUL83n3?m46}&Lx5o@@9HZL~*)4s<e
z#;~;N5fhmA?2=IU?VLx<VE#*wn8RgN&Wi=ib>+{SNSOVzB`;R+%-<fdhOc|`5Q^aA
zGafO9Yoi`9fp3m`#1v*e^@tht`Brl{v)>~Y@a+d4v4qt=wu168k66RKFXzP!!6gry
zVL0`nM@*pXITJ>q_q0dM;Gd%&F^AO$Jz@dvhdp8m8{hDV6^#1W8g{*%c;hyhTia8K
z2ixKL)mO_8fc>kRU-`!7+ge~nb9^KIi~g?>=r#fWKjS3py(1+603VA81ONa400936
U0763o02=@U00000000000JiK&M*si-
literal 0
HcmV?d00001
diff --git a/test/intersect/a_bgzipped.bed.gz b/test/intersect/a_bgzipped.bed.gz
new file mode 100644
index 0000000000000000000000000000000000000000..b4832259598cd51ab13e40080c991bc2febb7da0
GIT binary patch
literal 85
zcmb2|=3rp}f&Xj_PR>jW77X6!G>o2@m^?5sdf>xiWU+PUBt<0!#RW42mkTb}6KClC
Uqq9j3s9qjTqcj6E*c=c60GBfocmMzZ
literal 0
HcmV?d00001
diff --git a/test/intersect/a_gzipped.bed.gz b/test/intersect/a_gzipped.bed.gz
new file mode 100644
index 0000000000000000000000000000000000000000..745284b5bfd5040c8028876367749c2d949cfa1c
GIT binary patch
literal 63
zcmb2|=HNJPUmC=~oEV>8m03`bnxdDKn!@0HPQ&PliOB;KqX#}LMiyIlPEu4-P+TxW
TaJk@eJ#mKKKRTP#7#J7;uND<F
literal 0
HcmV?d00001
diff --git a/test/intersect/new_test-intersect.sh b/test/intersect/new_test-intersect.sh
new file mode 100755
index 00000000..6f21ef16
--- /dev/null
+++ b/test/intersect/new_test-intersect.sh
@@ -0,0 +1,21 @@
+BT=${BT-../../bin/bedtools}
+
+check()
+{
+ if diff $1 $2; then
+ echo ok
+ else
+ echo fail
+ fi
+}
+
+###########################################################
+# Test intersection of a as bed from file vs b as bed from file
+############################################################
+echo " intersect.t01...\c"
+echo \
+"chr1 100 101 a2 2 -
+chr1 100 110 a2 2 -" > exp
+$BT intersect -a a.bed -b b.bed > obs
+check obs exp
+rm obs exp
--
GitLab