From f0294a3e33257aee83aa7c507d2abddccebf0239 Mon Sep 17 00:00:00 2001 From: Tom Harrop Date: Fri, 5 Jan 2024 05:20:19 +1100 Subject: [PATCH] Check if cutadapt detects mispaired reads (#5696) For issue #5655 --- tools/cutadapt/cutadapt.xml | 25 ++++++++++++++++++ tools/cutadapt/macros.xml | 2 +- .../test-data/bwa-mem-fastq_mispaired2.fq.gz | Bin 0 -> 11426 bytes 3 files changed, 26 insertions(+), 1 deletion(-) create mode 100644 tools/cutadapt/test-data/bwa-mem-fastq_mispaired2.fq.gz diff --git a/tools/cutadapt/cutadapt.xml b/tools/cutadapt/cutadapt.xml index bc5fe47f3bd..7ccacb4c1c6 100644 --- a/tools/cutadapt/cutadapt.xml +++ b/tools/cutadapt/cutadapt.xml @@ -1312,6 +1312,31 @@ $read_mod_options.zero_cap + + + + + +
+ + + + + + +
+
+ + + + + + +
+ + + +
4.6 - 0 + 1 fastq.bz2,fastq.gz,fastq,fasta.bz2,fasta.gz,fasta diff --git a/tools/cutadapt/test-data/bwa-mem-fastq_mispaired2.fq.gz b/tools/cutadapt/test-data/bwa-mem-fastq_mispaired2.fq.gz new file mode 100644 index 0000000000000000000000000000000000000000..89b7139660b9de5a1634a1ab528fc175fa1e2999 GIT binary patch literal 11426 zcmV;TEM3zdiwFSjnUrM$1MPj;cH>H}<$L|cdF;htP;-#7shC+0eV(-X|39U%f#4+3 zDdnlURVUL|nU*zBw5-?$HV_-Y?f>d#f4a1nR*#p;?(^{f+NNz9O;?9wb7_q+mwNXP zZ|~mvZ0|YqEt~V*-}$D;_q(MV*7HGhNTyqNdf=Yv>rOAiA@{q-i!xvC8PEBQb6P$M zH!~i&^9ej(m($bh=r+#caT}+kSI|>*+#T;t-oqVVo?msRGtm*&@)LN!bb9(uw{f4| zujTiEcRqd5&!$J{SU#})|HuFNhsPz!b)GK6!zCxa=BMQL;{o(!TS zz($vF93YYgml;F!M2xsWN2MVELpT-P#!)tASl@tEx)XMviheR`i0JfMillk2&`m zD96L$)Hb`_0XX@)g5W=3l-Qtz5ZIigM}v}QYcYS@GB@Sv#Zx)u4V&|L)YJJHQuDjW zYZGEu_gwP5igx?SVvqD_d}gcCD;MZ@jltG$;N)^@+SB2m1XvTGZC(EIkJcMaSQ78gz)&KnCW%bwZaX91Aqw~dvRYg!^UnI~)!96_zG+gO;J3)sliyNJ^LbU*>(%S4z1G+3 zhbC|1ZBPw6F1Yg~0Pc#G!Gsx^cH)d3ddZ|bgA@Fuk)k7FJA z6nd;yowY0K5laYDZ zgfYqr;X{_c@j?dQh1vZj%^Gmak zZ3JY-4_hTZxxcghv3zs-;b17%35f+%^6Ih$l(b?YD+E2WkS%&rpo0I4AD5Q*R%FyC zav*WSjZfhzZl3H#`xHRisYHx<^K z7UpK1y{RY5CGc8iXXi|A*U_qfgdk6L6+N!vIpX;Z zO~q5t?6-pEB|&i<5wFQfUO16evTSk*vhYdPkvss0GqKMS({MU>O3=!QIIlbT8$M&0 zNH9)zwThfh*0JEqT2NGxT4EI8qh{1@ga_OC(gKHC$?+`A}iDCT3}D2Nv;9g)LyF#9gX-3}`h zX(f-5s6^AGcmzBMC`bB#Jf_?J==MY;>6x^dcbV#t5II2|=os%V;9P~*OL)o3-dA<)W5me7PkU%Sp8+}wn zM=Hg8OoJvc5e|dbwYN&=y6qJG>pHIuXh)_XhqVLP0;)C9RBZMbqgw3ewAj_@kjhOI z1zQ_bS5-GhB&O~?&%+Tx^e-vhLOG-MDs#{*^hVRK7AdA7OvRG8&AIc$ zTq%yBg)iD9ZLG0HjRfzcTV!T!M^Gty$-N|WjdSr0pzQXb$cor6vm!6^ z_kPOXGptJfeoWpkVxGmtF%J3|=Ui*ab&O#uVxE*F0$wWJICToe6I-Aohknz06O7iDRD62`8gX@iwtJ7%53wgg+9MGB?dfmyaCN#3u2kOB2O} zkS&Vo7H#8VbEB|z-ycGWWCz|P$bEnlxm~Z;iwC0zn^c4a7FU5BO`>bvxeTqOtyP;H z-KBuXaF|0>rzehLqjMT0Rze{S2@A{$A;p2as?I9KQb@8auq(E9WmE7qPEPuKXivw3 z5n0=}+17aokC@2X$)G=O@r*7aO#3lYMxG1>!6Ht!#)-Hjaa80R5E3#SbWPNwj#2JK zrm?Y*lYsmzD{MpWZ-khuPZZiVQ&@zbarN>pVy}pqK)a+oXl&|D+i4b5!vO&*hTK%& zLK)>)Rn1a@UiF@m)+?gvU`?luuh#d9KU*s_HAUb z+i`BiV~_C*9c)_e0jJ6Rgw^4B3~{CmTu#em-GZvbIY8o!2}q5YEg~&YQ%i1m5mc4e zXW2PK#-;MDaSMD$7&b@dgK-#AUw3x$S=dUJZDYq*@8HBxWV|~xhBpIwF?fA`VFmtQ z+8-31Y_P1U4q;=@c)7ddmJg>qY^!PWm*VR*GMYzC?U06d#MjF$CNF9h@8g$Z+O@RW zwg`|p(?<8ZU{m`BP_zSW08$j7#_Ph45 zmpY!m54rehChG5gG%4lXGRXOY#H?74^4&6|%UfX^$qZ#PgEGY&^-(408n`KCoY zo(lu*evWunn}NeZxHjO=cw1qe7lOtyMlK6JOxYLY>SN&|dBtx=2X$aRLL}AOzBoYo zh@MHV+OR2eWHf{oAQDnyewX&V2gkmIsAd?_WODYYm%7oKBb#hx z570`6e@|{G;`o7Bf}~w@ef703nV2BC*Q`ISZ$5Jj6_WU!*RA%g@AA;%qMf1mCSfap z_lOnw4@gT4Bn7yUtO?2@i)KFfmmWl$4DN0d*Cg zrPWa3Lu?|ueVsf)(Db|=I37yG#?{umhROMSJP}sTr*E%88O4t)x0BMWXp^nm$Iv6C z{g{YL0h+U6R5m)~=t>l_3{w_L%DJkC??>WCVHk?YgzQ!Dh+h*&-HASn5Tj3hDGXs_ zm06iUOHEe9L#5=kF`DcT_KahV3S{JB9Wt_sOESp3p)SL~2O+y(zt*4}>1y9VTX1H6 z@=b!^C%bbOmF(QOz&~pXz7nBq6p%h>67p7I>8y49DhPhq8Yoz!4@7dpbJE^+PV9-+ zs`?;rx>grXZWJ3qH0h~}eC@hnG*8;e43)_C!k02F0mJ5*Q z%^!kF$C(p&)=0D^x}vV8_=kI6yaAN`@l4ldcRX^b@vpk(`7eIsj8MYkFi#v08sG)EmavpkXUhf> z`8N8qq*Vv(+obGvpXrKPi<`e0SOdhnx5^M(VV_ui-nzf%+Yr4L zS)=P%+-SI*xC05{3_gOJ4i#QCXHVLylQIL7-T=z+h_bKdAQ8{^=#;aSo>_=cHd&0C}1YLJ7)4Z>ozFscXxVXf9?WVhH9 zrZ#nnULT~*#g+f3%oiFFsoX^d-4|(Wx;R~_S@Kwk&dk|)ol90D- zC)`tKlg;KQDLoV(TW`ycuW(~C1z(w)O{G4ESZXx}_3@#W6o45CI++v-%wXww>!;`+ z2{83*5eg*t$J3$3j^O@(d*+Ont=xTHmvzmsnKPU)&5;_eH8|qqwoQ&pyTqm` z+Y%`gX9_ShHig1i6H(QUSVu*1Zja()SKV&VstHPzs#@h(#3toCPP%lVV8gB2=JdpO zE29k&+&6%7X*6`Q+BA}D`!-P75LAQ4Zggu~rcu&j`G4GAX@W{F*|X>FOq)A=lh!Tp ztGw-c;w$8*Bn>NeB0M*56QGPHMM?LHmFIy5p?IJ1rLJ^ddIu=`J>5K=F8k76>Hb?u z3i3M%6+r4e<$a#fk|&ol4FCKmeJ!)LqKGWY);uh_n_z2jl-6=FD)0Oem8IXiXO)Yz z>WMz4-Vd3pjPZc)W0ZHJJ}nC$VyrCWvKm{){H7yO=Z@vGQMI1po9?D=W8scs@aPso zl~O($EHfAFUIS&PkLNSN<-q-wejSYXml(0EfPFjZMnLHo+53YZ%QDOQ`^bH3Xx%-` zJ!8c3$_3fQ<2~X2=yGiAC|?2;9OvAyuq?61fjIYp1u8u@i}z`dwo5(`I4s@BwT}VW z_1>`>N|TSdU2Pf~A6Rd#)YVs28)!8>NQNguiVl4a8+^t~WuT^Tslzs4Ydh|3A+0aQ z+6-kBA|E5W$8P|ofuMSKF~`}U+<#k{68pjpyl`GNge@++nb%13PHzxZlh~A)(SO1> z(Xy60jueLTDdn^r4(IX7*#Si|M3y&sfCa#)j*sG!mD&){!fs*q&h28K~U|CFb@)2|pM zTJl#arXQ7j{V-*^ije2ST}hgT{Md&nKQu0mGl3nYEW@x*>-`^CDsY!je zJO;{hweb=tX)R4mPg^@&=qFlhKSLo-`o<)U-KHnlYJ_FASjn%T7IZC3VVS}J17M|7 zao?e_QgZ!`90#+|=v?@ue}-G(x z>~>Hr?Jp>&1NKaT!sLo7$7 zUE-7);UU`4XRay}xzntz4kFIZ$gkmzAc&Au^WC(kO@R*s?b#(ahW%8vOn#P?Sa0HeYn0uJx_ju+e6=+xFInZ zrmb2OlP~)`%I%BGcr@KV`f(0yypApW)u__tA6y4}!*#!!>(0epHjv~&+ z=@1yVbQM|~gU`KZ+L}3>E~zu1!a3L+vOW|Kw_HgoqY|(kL`H=Oin+C7T2N5>21D7M zG?HzH%kK2A)dv1|-OW4tYhQYVN_x!}!ozKf@&i!9f}yOsl7lMvWimB!oLfR%9i<>S zjhK*9oRPM|?WFSb+&320fw-WVAX|x1(~lmq;I-;XF?+cB&;eca2Gk8cgA9o#BSL-$ z?bj2m2?Ad13bIs|O@>HS=p*SWG&+zTq0-F&>PXa*V)NO&0hIGepNY^Kqq#=o?@L5~ zLX6MN(d{V4zjNn4R^0QX4RCNfRvKsd1?;j$<00mmRf|h~^;6r}$7c$9{A$XqenuN)Cr(L;mM1!&bgmf6!%o{E)qWZc)9ND9cMjvn738bhDl^ zbL*K;X;K+1^|^W7*qEqq(X$jTYQyc5Vwhr6DDZ(r36(O?vQ?E)JQ4vN&8@9E)QmcC z6@bKXKR;H`IvJHvr#FCdIX5WPIA6Hx>s!bla{*8Jg*SS{ z50iWR_TU|=UOkG?wXb zHEK(sp?t&h>~56cvNP81=@_zV`~4d!$P@NYwj*w1nV4JaG@SF>k$?8alj{ zhWC)}T3UN;WCOx+97J4&xog+r9*C_9{bzgoXiarPWvO(}As0nA-T=zs(txP!_P<+e zzb}ca{ZYrZXNr3!zjkHAW=q~eriOK6_M+Q=52Pm!lL(xjcW_zLha z6!O%zP3XpisT%4npzJjm?Lq$*t^Iy3x0}~MSr(>;2XnIOwU$>+EYfy!u2uW`v>4kY z9;F{`4y+UTo7xj7h+`npR|m2vJq5eqz1rg#4J$>*pk-uM-YV3Uc8W!BDYsEM!?VF^ zC!(=P&gz}0W&cC#?Ek|#d@Hj``4wNu z?=5D+$;bM+^gYroX<3j3sfhHlNplZbDLsa(snL{+EClnPyLa(#|dbPBin;EvgCM>W+Rj2Yen3`^*n2e-3*aY=eUn7c_*jI zq*!GlhQB%cQ7OWZbh#)ckS#K9^bh+zDCOv0ZwT_Fm7^h`G?zwt`+dJXVMw;nB(;0CR?b}QFfYI zYeO1GAjp_fwk%e#!u5fz6Zt(_ioPK+Ac&myM@>QUiDlK_wRMgpW}o!;{;7^4Z_dB| z=OW2ASD)w2sl`7uEA#jtwslqA6BIBjLtr#<mUNz|xoSugBtV_RU_fc4D4y&c#<;6O4pUx3Zc`b{bTW5?J-=WNk zTQrYk;v#`^O@wvqCWvRo7Eh(8AHm@B7AEQ1grTleAz0qo=^!|`DP`K0?WOIiR<&(0 zEu~@g8$j9bk-e6Z`Z- zX+H}=!YI<046n_7pHb#6d44xbw532O3P{q;hsw>QmPL{Xr7Blhh0=M&+dUOpRJW=O zQMVQPxJkygmcNa6fO6EBTYG7C(lPXRD<*!sX#MLE310;+c~xnV9u6rl^LeE)Sl6r% z`QdUyA}y}lCH8a;WVqih;=0uM~~A?)QnsWdzTMres0_xURO7Iu~+n?^U}s_Q3t zCl{|kfTDl3zH0j4ekxU6uUA&9yuyj0tKC^6M&`aKp4Fwl1?Bij=-j_nqdZijZ&hrp z1LNUw^zqjv#zX(u#9#|BYun)^q){QdaTuv`*Yw!i;>o}`4CUq0@fx|32AmR4=jW$_j+>jvo*?8 zZ&w>5MJ!&nt5L_1HS!}oIRhsx=qB*c1HMBaVhEUc#+O_mlWk4CxV+ZO8Wd_Yx;^MJ zk94zf+amPIzia3-_Rkq5&dX`#`KfJN8FKIx&a7hGvgAyHSc8;iVgAZt&Aw>WGjeW1 z*{A%}DNy{DmyrrNkuM2j$G?$h^|$1e{Rp8_B3Y?EI~Bt?ezH>1?vMUB-oVLW2aSQz zJTd(@n=@(O+aqM8JSI_4tKmm&{bF%*GmkNoM}*W| zNj~uF7H4CS`CYh5XcE|>U@VUVUi7C0h=_CM>ot(&2ti-T-Q?VLEr!yLN!)-*g0+=? zw2+|id`X-}cicnPxEi$@U58f4b=}%TadU^31~tzVbgq3ZMj@aacY6X#bKaF~+b2`FdGaymEM;qT9-2l;dK zAbS^S_9#0xZJGAD^F7wk`jGTX;KvpbpvQIm5rceGw*{fN5yG`r=f!jSg z+#|zwqksj$a?dqgeXQ6Fc2kL$T+0NfssBA}ilKA5MH_qpA?V6$oLtT)g30A*zTP*^ z-=N2z4284qs}zcdx)HI;q7E<(!v@<(Q}5Y{q@=X4L-A>oFsJ3l&{Ew%pRII>B};Mp z;Ph)>aZ3K?#|~VhDGA%zl(tO2uEji&@fCbr{?U%;e6^5lX|!OQ-1=7OGo9SHUp z%+GpGp;#w0Hf-nG^0F2FbC)i8YLB)!bLu2LW!P2|T*eK@t;|W~u+!(>%(2CFCN?f% zqqzs<$1OY@4L9XF-F7egwgD)R7PPz5kwuWdb#2#uX1lp49~X+aiOFsW4D#68#IM?1 z#8QFJ(t!LrxGroA;z2ZLQ{K01(YQF- z%lK@3bKW2!<(3{hEnjPBQMa-g+u0ENeqHuU(nlPo>C z9ta7)lq{tu$x{AWvV`o$M$SJCGfp-HEUV0HSj<=6Fd*JlA9GJfM>NH@?MPtVM&W74 zDy25X{^^MHmPI)okET5z^}oyW{GvO9}Fs-sHrPbpXN&d&3K_pI> z@ekadm4rRYjpG}B=V8gCQp{F^syh+3p?;&b$=i?0m_-{M*SgZ0*tGB9o+$iPjT~IUM#tLYzJ=13+CFbcLrt1@Y<9%^|hPbs*A)V*E3tOiSTBFg>T}Fd-(LC zj&;4ISf{IQpK6pUk*TPp*(?&Vkis&DB?h9=*KZ1aa8GOA4CDc~SH{ED81)8DF1k4& z1985>!O zk7u&V_QY2)#nJ>H8+*?tVj`dJwQ=c@DZ%K3Do`Mt2sEvOj9PbEviA!AP1V*K4;unb zr_>wvWVhezJ>9$UxBIV$NacQ z_>_4s;+2G%qkNZpGnF%^adcs!Dk5qMLU_Fgo}BiB=CmNEOX}A4JQvVfMI1JI(Vb|X zW-JLg@I-EXa_er162Piw043r+Rh4t?3RAaC2H zr8i5%y5yMUka?*-;)-Ku!|I`BNAoU)XX}F$9H-&pYMw&6xm)!HP!7Axx!vhYbLQCi zf2?6>{<|5O``XBSSxNe3v(P->$WP5KTrn4UnkhG>Z^uUFW%`DjQgeA4t)2E<5XN$* zf|h!h^DXu=JvU;q(-lu)N?g8$LEWN{qep|aNz+~?SF&PSLezE?v?;XI;ezB2vnkDM zprCKt0T!h>{Rcs^3iX&p7UXk zy5l3~P=4bW_{Xj%+r_L&c%YeYJ{yurxa1-@9;kcQQ~^+5;8&ko*K zSMr0>6hVfXDqmkMcRcDY=0{1+T ziRoz(7AG+;F^yUL64wlK+AZE#CaGOd=xME@_rypb?8m+@Z)0XcjEs`wp2Uhf zZAEFmsMBL$*|AXgIDo z6VE9vA$>q%1;V9|s_K+NRNRfyMBhDPxuL z)6loPnD#YnGE20jLS-<*e=LMsGRSEHBd^sXg@yOM^7%ra%;z!s$iL$ zn4&+K7;x=ono}F@+f})`Q*Ywrbk>awIR36P61E2<%tUszj*kHl2eH{kzs`{4iSiCP z4{Mn|4^^4#m($$5YmvZ<<(e%?L=W1ri=Tg}>vD_wQ5$^fO&@vpCbHPkZS4mw?!1Bt z8eAL+Czn%m8MgwzUlA$Z#+16>rYCu!C>WIjPX=Nor51$vF!z9rGwDswy;9P;1ZXN= zkJ+Y)8nK~qIZX;npz9WR)@1b7eOQ@n9YgvcV{;2i%G7PuauKPib=zsBD(OAkVmDBm zmy>x2C%fi!HoVc({1evDjd} zg}6HDgO7GO5fcIpD^aSdZLtZ|(SNS$(5-q|#{-}okK7JLH>dw4DCh$DLQXxbx)H3q zg??3Knjg1!M)MH2wiG9`^w16^TK|wenNs(T$jdFF-J5LX(2>fyF-PyZ(IeHOot=f* zF{n=-a{3x5$MXRroQ&vXzde!t5+{#cQNkC^mh!3@#zq6Kb&1Zics~tK;{E(M58G`N zo=3Y>tQhxn)JC4p2J0$9&r6&fC@^l&1F{jIU@`OGUdZ^fb;McJC%NqRXK_t;fxWy* z$;fWjq%S#&6mR9Nxr@u$>U%zq)SrXrL^nI$vyXrou@EFN#oQ1Yx5CC zYBb=Ni7{m5ITxorDoQyT*471dSGrLI(>ixk7{)_fM7R8d)DC49ZcTPUyzPy+Xgm z*xC)db|IoM-Y;q=gwAdVG(>4@txPq6TeB*o}?(POrR!;XhZ%G3;IsXlTU5ae|+ zLk#12Qe}`=_;DsJPnqj0qNX-wRbd>QHXVzSbZhFCQ+&ORZEHHyLmUe(s~}ZknUbry z7@ae57G+;=0OfSrQ_OSR9WTEQQ1bskGiwfx-vG)}!_cz$|K!GP#~{&7KgC8dmjTD4 z^qvK5h;VoGq$|4nm}_X18;Q~opqY~J*r!U`Cg85D(sMq;n?{N&$BydLd%}aiWF%4%6rBB%}z?$ zKI>)I24&0hEXJY?DWzwY%H_!bM39OKvJ|&lSZ4lxsnx(O=X#%vi6Vo{eK3Amu@{+* w+UG_}<3&C0B3GvCN>&Z^rj8?_)vMcfB)AImzNT$uSO4Sx1GW`9oLJES00)qgpa1{> literal 0 HcmV?d00001