From 09c573f2f2bdc05af30cf548870f47aadcd73134 Mon Sep 17 00:00:00 2001
From: Paul-Christian Volkmer <code@pcvolkmer.de>
Date: Mon, 11 Aug 2025 15:36:20 +0200
Subject: [PATCH] feat: implement info subcommand

---
 Cargo.toml               |   2 +
 README.md                |  16 +++-
 docs/info_subcommand.jpg | Bin 0 -> 16769 bytes
 src/fastq.rs             |  30 ++++++-
 src/main.rs              | 176 ++++++++++++++++++++++++++++++++++++++-
 5 files changed, 218 insertions(+), 6 deletions(-)
 create mode 100644 docs/info_subcommand.jpg

diff --git a/Cargo.toml b/Cargo.toml
index 7445e32..cf65ad9 100644
--- a/Cargo.toml
+++ b/Cargo.toml
@@ -10,3 +10,5 @@ readme = "README.md"
 [dependencies]
 regex = "1.11"
 clap = {  version = "4.5", features = ["color", "derive"]}
+console = "0.16"
+itertools = "0.14"
diff --git a/README.md b/README.md
index 0803fdf..5070713 100644
--- a/README.md
+++ b/README.md
@@ -6,10 +6,24 @@ Application to show information about and scramble FASTQ files to provide non-se
 
 This application provides the following subcommands
 
+### Info
+
+To show information about compressed FASTQ files use:
+
+```shell
+cat file_fastq.gz | gzip -d | fastq-tools info
+```
+
+This will result in output like
+
+![Info subcommand](docs/info_subcommand.jpg)
+
 ### Scramble
 
 To scramble compressed FASTQ files use:
 
 ```shell
 cat file_fastq.gz | gzip -d | fastq-tools scramble | gzip > scrambled_fastq.gz
-```
\ No newline at end of file
+```
+
+This will scramble headers and sequences and write the output into `scrambled_fastq.gz`.
\ No newline at end of file
diff --git a/docs/info_subcommand.jpg b/docs/info_subcommand.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..cee2ce6f06a4247b1a15436dad71db83571fbfd0
GIT binary patch
literal 16769
zcmb`u2{fDEw?7<3QB$?lR8whHO*Ky`RYg%0HJ7v$C5^P^F{)~wilQ|ZRccC&F{`<%
z#vmj}LQ!+fArWuB|M#x@-uJHi`@i?T>;EKst*nzLC(q}cea_kYv-dumKEnYn8tNJ7
z0cdDw08gnez!?Q_@g>yF-TB!|8+Bnhc@-6E3;Vyd*s2SwD$7wD*u0&+TwFZWh3`L-
zmOH}$bO5yH&i(VIUUbwqJtI9m9Uc7z28Q#DEEiZ<m@hChv$9>f$jZjS#>{;2%0-UL
zoLpR7EbQDoS2%etadL6~(+LeN^*40%O!V|joUF{Ood4tN><55@5%8PlA}!4|z&Q>Y
zS`M1C9)J)4KtoRj_z%W^yJ*f)pYc2cBhv+D>J9Z50q1CFY0uHo{`1t-yF;np1L!#D
zFNw%$pTBJOjN#gAPWkYZ0!Gn$)ty}CV+1h;C$G0m7p`#g@bX>1A%0WhmZFlfimIBr
z&V5}y{Raky7LP5htZi)VoS(b6y1j7s@b>Zb^A89Nig*_p6&({Bmzwq=JtOnur>w%F
zFU4O=O23uW)YjEEAR3#RySjUN`}%(l{2E72OirPuXVCbC#Xn2SD}Ptlc6Rsn4-Sur
z$0z^bq5;tU-?082*?+>tLB(~Bj*ga&;U8Qy=lrQc%RxskB6t3hwi&~-*O#x!hcj~C
zODU-CWD->{CvZ7=ja|4Rrij13^AEKDh3vlv?Ct+AWd9EA|G|X@u+Y*_Cy$l`pameR
zbUCt)6(H;4Hv`>To2ER0i5kX;%wH_D@8k_hHf39?G}=IOBD^z8X(+1R|8sF~;#RCy
zdZ~Nu;EK!~2{@fomfF%N@&-7$&9s)~7E^ZyV2-WlH&(GH{2>)3Kgt|<yUVWdn|T@c
z?x~mB)#coUN(*DcQc2vGWV5=7`Vo#_odw9{JzVXZf#uusI0mOzs;hCUS;?~46Z?Y>
zXyeXASukX78({<jNqVirC_#eOj;jSB41}v%%&p=JX8_(1(<VYiyC})F;X8tW#NxE&
z5@S=`99sn+`)<b(noqPSeEE(9_8&sgf?Ndd_${Bxdy~E*RkoebH^9Nf(WPAF2D6Xk
zcZ)zAD&Iao<yJ6)XnU*y_DtTn(`cVPk2#?r^m76%4*iG`i&z#1Vp_K&f74_gkD#Sd
zeLzuQkF)I3{b50KhDYJm<F)&yZ()i%oOu_MO~=|tDbnspT_KN1&537DHOZg)hQ4s@
zZH~M|wTJK$2MK$Igq5yb=N`=)&13cOXn7AAp9lLhZV>wszZBp%>puxg4Tz8WhS_-b
zzL97vKJIP-?;b&srm^~t=uEzf!Giv;gKCqygLNq}#h_H+6|{fKcy+E_wO38EqbW2P
zq){($1y#ER&yY`?+J>$dw;<p;y+9R`LWt20Gzt=;{UzYpc3pYt)^$XN#^9?xMxKfS
zQF#Ve^UB>wsk=YvyzTR8J6ue8L)=@&sut_A4GC2Reut1M`F`1qUtema&fgE@yh!7M
zb4%+DV(6Ala%3H?7n+%e=-VoP!-Fdw9bSZHxsinP{9aJj1V|r}n_2{R?8E3`av=)U
z4Mgp4#Q46;#NbB+>E*9!?%hUtIEQ|3(PqOMJw{7@!jHGJuUL%@>r>M+tum5Z>-X9+
zv7zC8e5D;kttBXtPhpW1(;6S^G`yK$t-2oJuo9tp`3%sdP1x%dRGzFi-H@i$XNEu|
z5PSKLI$-8gO0yl^8Lk7<MXGAOuMnwH2=C^ca*d0jdf>I~nf1d$f_PK`G4L}H-p$SD
z72qKI+wofT<M&;u!>P;GwMplmJMlUXFWl|?D`!9XVNv(jg<VQMncX^A!$w*N46O|#
zLe~*_B?A$aLH=($+~6@`93gUAF(DfG_*|n0vAJ~XQV)>wERC@O5x?%PiLyGniLgp?
zj8<Q3|3t`4)D#TJ#zVVXWH|y>tO-_?xIm?Iko{~4$F)<p%ptvoT+IjhS`zZkQ)Ly!
zJdPSCj+k|q^1l@=*a9XQu-5tta$--G5jFz4NDjsuC0}}>$<s8U81OJ+#P*urhdu@2
zIOV6in*x54PE|zzLuX?8!q8>o$$dDR^kmeR*y|Pf=qsh;k}ZAF&W61HK+gzam+uTx
zG>5Y9!gyhV?r8C@Q?TQX(pE}E2r$z;jl@mTw2xJK222!#JsI2;8H7NHt638S_Du%d
z9dT{KOO+GsjHeXD$EP`d!yEBobluzGp}eT@X5?h!uy=S`@3UL!flcNlgO~e`hwZ`Q
zOTl*xm*V36XH?d<>L`f0mjRlQBn7h3j!JapL7PI??wmV{>9v@Y%xBx{yS%9`Hd%K!
z)-2oawj{~ZT%G);$w~~`##x>LlI!w$jY~a(9D0kZiVFsd>Lj6JON1-l64=D(_Ynm6
zm;EpvqEI&<Lnla(7}2x0`re%zCNlBzWg)USXYCpMIfk$13{d8rPQtMW7IERQj3E}d
zIS@X-N0yEV@4*tHl$%1GOVJ+Yt9y<*AMG5*YPy)u-H%}4!A-p`xBBvY(=}@q?FV)^
zo%5Mgux7j$Jm=BD;C6&@$R#Jqc4^T7gBLqk2H<}k1v1V6Mi#;bbUUKM&b;q+_*Dbr
zlG^PWBwGk>8_>nDt7OqIe)qz$FR5!cH}+R$D;D=~E?(nx$?C^0hT%)iQEESN)?tLB
z#gcdAI|O8x@<rl!e>WNx0ebYHIJ=8u&ELAqAnSokZPsH=ZH<QSuBMLQVak<&w9a-B
zqRJvWEfL;Zj>EK^x*}~FCd83$t#M{oy^8v*=%q?l7{r{ad~elj>Voai8klU5gHp|;
zWaY-x>jaBrm0ZdSwC@+`&Y@L^TjUwQ?F=A9M0W432gycXZ$}oon>~4QFInThf^d}b
zt3Al}QdDo?Dvv6V?*xg@?9XR;k$<6;+x-l{<eM4Us*O7XAR;dl4W+HNwtuG=$8L(J
zFuuvf*eHm@m3Do0>-GjDuhv#Oy-qTyEKqAsO~{CckSTWCp4oD$y?Uo%@HLWVbOBNS
zFn}0I9a#MZR`%}7;IGEZr@4u<9E09VXMl>`XMuF7o1MeGo_e>0m_lfGj*?ncZnPB;
z{lJ@Y56hBQbg!m;?0bH(XG_qj2%b?Hceovny`m|!ie;^1743%!k1GkNJOp=WXit(A
zq<hKCqF4IzgV|5o{=D_C!PibWRfH%JkMR(3JTwA=?6u0{K!sd`9_stzReiyj6sbxn
zWALlL_oeqBklzsrdo~m0Rx%!vkUCfxY*2UhRf1i780t8vwl{={q#R;E{8V5M7UU<i
zqNB2XUrl{CLEY9_{@D@yB4r}oapU~hXo_u+ey0{67)dnY7vLXB>oo-ish^f{ku2`@
zUccHnQy(s*IeR!1s!ft5c613UPqj15m@cU8Z(kwW9hpavIBhiI@vp`W#cZQEIM3<O
zj$EaEa?}-;5XP<u;m22%FR*hwSyjh&mnT}Hbd8m5F_zqy6jP4P;RAxq5RIR!hL&Y|
z=2n=2&28|CUr__97L6a^=KgiCO5DbHp5$G|gWhOr?I5CEQHmfoqu{T7<AsfRVu8P%
z<r$sb*yYSngI2xv1k}dNj%c6emAaiPiTINCai~aMu)FydV6Q^T>QTxk+FA*V-N@O0
z&L77~0_N}yVK|r2(mAFQSJ5iEAj?#L%rP5aB4D@mXHe3J-|mSz2NQJ*`W?}slusC0
z411?2t#cTZZ!wkFNT~0tlxctsDujM8RZw)h5DmlW&(0+gOVrcTQsdUYCfC&;&-jiX
z1|{}oJsH4;<sXNl_kqEkpxTih)0jM-{2om~Lho?#m!5IqPcE^1R@#|1Dbv7WeKQ`X
zTeAOO>QtI&&$UJo-0A3sbCFdETn!zIm7va`lD8!j*+0i$x&)UW<U!T=(7wk`A>=Ou
z$#z&GWakXf&F9v&&2TTo9*j4=NHm83v>&;kQ;u~Atf{NMQ`uc`5iwMckY^*0sn?sQ
z#Zl*?su>!pXji_b`Cv{G3+f=ZIVKnE8vL{B_MwcF5uuB%{`eTOb^*~B?laJet*vf(
ztL7Qfagg^=Ryft6QEj^NXK&!2|G3=XWfX=G8)8^Vi1MgnvH&sS!GzD8Torkf!J6<J
z!A+1#u;&Zi&{{i6Wttz18ZnETuPeFW`N;0{43Ls)=5184pm~YNu(&M{;AtGJL)bu}
zojzF#4acd8%O(62iJb-Z>Vu^1x#j{C4pnwKl3*MiiFIQy#tnQq{epZlr(Pcu<eS{>
z^m9XMitp~u>0JHuO0ty*-r4Iyn|5K(V@l#GfudE!``C@7^{!WID8w~eJDaR2{Ai5*
zG8*%Vfp@u;gm;St@-GNQz-|+?;>dRGtSuc=5T$s?y-?-gpT!*gw}38B&40mOdFe?q
zO;oz?mfKB@wW-w~Anm2io9Bnzc{)|+j>x&KM-46iY<BtoY<99Sd5Q$W)IO9=X*q9G
zz6<*jkz}8CC1m8rjrTyzEA8FMX6x+tkKgMtuCa!K2u2Ht?o9>!RAf=WN{^l9w_<lM
zhghFG_TBO$p_d-?MjwX$0sS+Cd%S)PK?&Bs{8{-;MPlT8Y(B7f-5ufm_dQSOfO5aF
zKVv=k3;-gW2Jz=zCi;6X7|CujXHP38-u^AU*w|?0pJzel(0fyhxDzQGM-g}u;MZHv
zZ}e9wQCY2hL~&HI?{N1F;AJy=UZ->2Au1Y7-2(#{%)aY4*Nu9%wi%-7<bk~TVOeI{
z*KLE}z_uiOgKZV1B^TJrL;2Cu;Iuk*Nb6tc;E0v&OLR{eEgik|gZ5qMijF85r0o*2
z+Jcm>z6GoSF+5MT(usU|+79jU1nWMSNM`5G>WK<ksJG8$OdnkUs{i=$OMYPA@7tm)
zn7o;)Yweq-`96~Q7m2}#{yPrw7z6wnfNroth{LgM;zjdQ_dH$ssHb{M!bxgFfd^~B
zHzkhOS6$BlikZI;#=V;DF{a8P5>m!=aQ?)q1k30G0kHR7LgwLNBz!FlW2cvUn!@58
znnYRu)z~(R_p|@X{;^vl=I4*)_Icv<Mb>+)Du=HKbFI`3tJO%k15#!?a8BLvuIKYE
zO`OTRcd{cm<20Ylx&!V5)mryp*=O9dq6%;%(m9~ki*``D^v`xv6*u7(#^c%7ofX$(
zJmevE^P$<MyQ6;(UmwOvtopQ@LWy?}no<JplP97isp)_cBLb>3Sc|>aot<jiI9Lhs
z^ub1&WHa<HGWmy4o<3?(-m?)pR~L9NB#9G`+nAr(XQ!Q8Ph8MIH6?Nzvnu?Idi0xd
z_6tY9q7i4|(%1OPlVJ3ThNQo-)-^8WMOa%)o6i9fq<}+9ZjXekm$_rlwTcJhv1}s?
z#<va1wBkNT*i=`4cV1Q=`^JRxce%Yet#ebL2_C-1E4W~l=%;6s{$zy~$TBrrKB9IH
z6V|!9i)?1n)^*TWKMm{LfczD_Oqm}&I5j1KA0RtwM3SJ@*ttM!?v@;xL7&b&9=F%L
zJXfb2LL*3O0pbhT#4~`)c1$zZ*W6MMO^yklUdoSpkd&^{Xq%4mu=X(5UY@EFI)oFY
zl+V7JUz-_mUE)f5r2*XiWGYutu6ed__9u(^B_`pk1U7{s7YmM!rJLm{Xo~Urk^Qn^
z+^dyf2Q1+Xpe0k%Y1st9<$COOBHAUR24UeY{sQlgU>>Y#KP!G_GkdBr_6uK>8-OD^
zTz?)oxZ#$g3(@v#?E+&N5JQd1F#U|Kb)%=hXVmNK)*B%_mIF(#^gWg}X~??7rUk8;
z!Jhf6#NaB2HNk76_3;ljUT!2@d*HL178uF(;}<t(Eyz6Yda>r?HtPAO;);c??$?f)
z*QKUywYzJ;KqdC@B)HO?(Iezb#r4X2Yb&P6?syG@61TW%4bL!XrY=uU6TrZqNveHh
zM4VX=dW)seWcTnpC;8D+$>F6(P-xE!bBc(ReL+fG>cuTJ+IIG5Nh>(3eDT=2<Tvwj
z$nkL@wgBaI`?Q|Rp=FB35{_;ea>A~3=Mt26h8~wU>JC-P9CW1llg2?>?P`ihSbTe>
zJ<2%wT`7DlIf5eA(Lhk?ZUKY8+AG(NUDbc1zg@rn@i{PgA@&n{S+&5BIT7N^9!=qM
zhvMO^MpdiICsB1kGk0+L^tAip<)0i%^jQq9AaU2<^5IEp6)W+tmXG|r?9TwQu)%}T
zF>~WD`MVutCY6%V<PpQ!27*a_AP$4&>UQb4Tho%%72Kh8e&CPQbK^+F)pySC8F~Jo
z&{*0K+c5Vc*bs#m5rFER7Y+48J^r4`nIiL<bpK=Gpu}~(f7T@fPOHZ>K2P*CH^dA-
z5}h@)go>$?b{*H?XV!aX3@!PP{TnS*2F0We3_aetWR@LA1X-VMS%<BU6zBA1c^P!9
z4Z#O@9X#753yYiG?UtQ2ANCG}9L(i~*^x?cvT^m2vh0*0N1h~`;fMv<%|T<UUk}X2
zcZKep0sKhe*8_(YP+9Yp`X(0)9NM(bp8<3v#pVt{Q&byV22o1%evFu?%O+LInjs*g
z^R^F1hdBD=tG--(^#FHN_Of@n_R?<?gFD8sDA<4_e`5nd|E=^j?yQl>1I&=#9Y=n(
zmrs$K7@2tbHK4Il_jf?){6~7fmGhNkAry3<IzZ~*8=78hOby=eQKe&7acqh*DIylS
z>a11hd)pXYcHb?o9LIJjQ!<X{m#gxno?~W|DaNv>KQM1p0o7R5kxf>lx|#8hrXWBC
zeCT;4m7wAU#q8yfW2+s#*489#+$(@sKBxx9@1R57jH9ybiQsPY^=9Mtf=YHyJG-?*
zWxjx(T{>fJ*YKIlftgyC&(2BU6(Kg8p@rlb+CHlwrIy**RZZ!JNhsv$V5?OkL<}rK
zn_55=Z}7G_$W%4xwVUloNNwxjJPh6dRWD^3AAq-vDbluwO*qTH>atxDA*O_KtDUIC
z3N&c$Z*q&(sGq!p9p=h~%v+{|V!=uLP|3akrGxy9A4-s5Ty(Tz{0k)zyF32@eG6(o
z_N7c6V%8>Ova-aX`^AqZdDe7(aK_V|{lP@>ZXw=QU9U)>7W+H6VqeECD_$A5RHdlv
zcbMgPX-6%GTO7wF^|``8vbmL53ss%$l=dntEI;ZDAZ>cCuuTJhL>AW~l)}WQEQ6-h
z834~Q<3e?17;^svHtw0GEo~oD4!vty`<juhIl*D+zblKo;%=5*G}Gk-45?VPjaBLP
z<p;iYes!L%VZ08KTq-=i7SwiF?|qnf?{MKH#CfPG{maXs=cyjWoLh(D>7XohB&YQx
zws!litS|c;^nB)}=|GqY_^*SfFLQt!v>LNryIYvjle*RH)X`$qHxK5k=ep&Uy?P4M
zd!GzkWZO8_i#s9V^2!twMRC~Zg()$l&l!NjK>RQ%THs?UtM~7T>xgtud*Npyx|l{~
zxjy7jU`IEU7duFBTNsISHUc{cqIB&ui~_ixH)WDcCM1W%=e^2$A7757!#z*+_I3D^
zirM!IBWT?#Vy?@cD?Zvh127x_`|;5|cl7c4$%<iDRxjGkdZ<77YU9+V_rfM!a*y6c
z>{;?X$($?;;E3o`g=*)d9MotSp#2&c%^rOz_!zREAH~O=e_DM`Hyl26+bJXUyqKvZ
z%Pj)ki(+D*Ny*00&gCAoXRa@<v{ip0^xR5MMhzyw0p>r>k@NG0Urck9)Z;!enl3rl
zTZRdjvm<g+CuR4L(0psiOAZ6{)<9d>$WoZB#j>5`notCmu2m<HP&Kh?=U18*@b^h?
z@Pd)CvVAs0@*zUDS^WJo%atst+fuh<I<UPVZFYKvAC~BmCue~0wfx~Hc5!LcQ_;~%
znJj^;g!Ua}Dx+`YN|Q6b-<j+dx$yp8N05O)%ww}3gB%02-?;J}GmCy5YSjwG%ujrr
zC-o{Tqh984M@m<T-G?6fFIBG<jkoDj9pX;tMNr7CI*7Q*<R#G6${znb(R!pyW5f9V
zpx07i;CU{`z{2ix-Hh?4<j0T8zU?)=tZ0XbA8U1PS3gSDl5WNlAh7&6JI|043iH$=
zAqOT-8BVWbl3^{uEQ%}|yczEAyv}qr@rN4&Z)b(wfob?|D_Xy2PE$|1WCU;Y>1)hs
zx0rd)j+|h@PQ;de&~Edj!Cqtg!08bv>v?g&USK^@{3Mdo%<92EYW*L@3Ka=-uOjlB
zrINZ&>gQ8@rDUOE88sjFu04y=6W}#1_*r2mp_o<6T7LnZ-<P#H5-HuPaHG<1Ua%pz
zst^er7q-ogxY619@^#wJO`Ez$<H1p#x||k(!nm9N$hmZ`bxn>WY?qbhMI-Qo8&e_r
zub8p3_iM6kLa<U&+ToA-Q*9#lO9CpT$qtZ3{@@m)fP-YVUn6V<bViq-+`GOzbDfdP
zh~+I0jlE%R2Gvtxr7T7|!8N>Au`7@k(C=Nt)t3}eVzOx~p>49^%d`A$9?F1A6wuzu
zHCq)Et?68S{m4xf2$0Nz*R=DJFOF_gJpwFUKpqmeoV~;Fcp{SIVU425q#4KsTC=1?
zCd}&bFEMPv33~%o>XV{|S9j4Ahx}masvwT|>?H@&9H7?DVBh?eF_-bWHWx_6%yKGl
z2GiOwbhGE{yj4Oge>6xUQ%NgA^G-nj4?&?!-7&W@Jo2k~+I3WP9jWPtIsR}!Y*eFp
z>%v`EVJ@CNdVjeWl=C~Fjn0M2pphSh7kLtdt%4Fe0Y_<*HD{mNx|YHFTh5Lb*{)>b
zGdF4>i`&crj$224Q_bu#Ax}_7_#xi?V>kP5np^WP*M8_^((CC2O_`^Qm5lJ)e%Wu^
zU$%EV(KuaBG@A3D-$q*HZa0*B$B%ZrFy#*j?zU<dtZy1y9@})y2#|)8HA);A^Q`7X
ztg_~&pD%XY1>i1x!*3^PDih$*_B!JS1&|U7>DT#jhAhkw&j)P<kyOBp8b{5_jW+S_
zuo~8t*j6QQIVK=hyQvN-BZ5#?lf6FqC!?l^{nHsh>fytlGr+|~=Bc9k(~smkkRpsD
zF?s=(kN-Nk4d2*jZ*U*yfm%&jTs8t09AYF@e+-jdV@=+A#T%6cL7^iQp<24K7pTld
z_~1#?qT`jFDg>_03S%}_&-=1!^C<8kz}MlIWxF1|iH=}Vo0Ij}8Q?;g)28!6*xR~9
zh?0GaeeUZmi-$?b`+0{WM^cZn%NvIWpALVK5UEWdNZw+l?SzS(egs;ScjFp;lDUa)
zkB}A5r-~-3T)-V_#e2MkZo9}+8=+g1hIQBG>fOIkJO{U+j^*Erm<+L6ouC~flY~7G
zEm=3SOWn1;XwE&Oc`Lj^uiU@HA@kAohaXE4N23fb);>*QbfvrKQ4&dZuOup1Hca?H
z?%7MWgw9(41gyQ%p()`O@B2%)9ZQ^x0>pc@YOt}sgV6IjUx>BcnP^Bh*lbNjz8?2+
zur2FKm^AxORZDVtPx<{t2@`3LI#|$+j`OW2KH`=<E!*&~4P;aSY8A_IeGEubN_2aO
zguGmye8OsKMW402O#gJ%Q{L~(vk!7;`5qaQ^OEo3g)kO-vRVuA^5Pt{M~JnuRl6#{
z7UCfLWl9z@G=lbIc`C0jY&k8L!1j8`Pm3AG6C#fS3J}1$Q>|A}LD*!B`CBN5{ThP9
z)OF{H;%jxC#-3?@ZTi6eN|#$@*_qATKFuD-_DILUpYDzl2p9O^$VrQ)ZejadV(?Ok
zIFWe~#6N9|b|8L*P<x-jI{c^&jk?T|YAj~>BV6<ix1ap^*HeWA{S5?rZ2RTaJmra2
z<KVB_-%iVal7ynGl$ed6B5`Mc;@SJy##!zDpR?YN*ZY7&8<p^EYs?tD|NCi9or44Q
zpo0!+t>)3|2BSLdOWTFXw}`Pg2iM{8HPrnvnO37BYz*g<nc1$fZOiwpeQKPz@=L08
zPOY`w^L(TAri=k;k0*d;Y-#1IgG|VM|GGdSa~a~R8upAk5qaJ@hvh4EIUpmgG~}M?
z-gUiIr8#(3Ejddq?q*@V@Xf+Fz`v4pviRbVGWH>Ie=`WZn)$GiWRhFzzQl6)+eXeD
zzT12SfZH$Lu}=1xshczG)<(3<m75+W>YGm}Gi|YLcAgXe72mi*Tlwm-RT89JA^q96
zZ+m2N3nl1ez0um68fl9?0}#?f<OxFW!fIiIX8^|UZ9*^ij!sK9FEhcG$7((0CTmei
z>D|_)L<?LpS?mbJ`iz1msgQ4&&-S)#v{3q&44p@?g08Rg@%?SyBk5vGV-)4E;^E%V
zml}AXzB@uc+PQYLB125^Ld1ta6}SQ9wRD|b_UfLlB7L0vuTpZz#)b6#<ra1t*kn-H
zTOWsW+yvJ4l5(Oscm^n2I#^~HVQL>GW&WaTVG%4RoSOaz-4PyE(|+}@<}H$c_UX4$
z$LY|tU;XKazTt@Nx#4ZCujiEy)E0w(<e&etgXlKpCRTWqyU4_pL45Y0t!wm;2UOV`
zU+ymOkFe$}Ce^1u82{j69_))9r+NaJn|~e$!g*v`VCq4RCy)uM=>lqSKy?o>bn-pV
z<2JUYaLmX$+S;7~21(EyheEnQ<Jw=Tde4q(DML?*7KaIAFGDyAU*t&<{Y!{bd)?!6
zgGCPWSOd=6*-`WBKpMm00?Jb+kYObb`bJ$mo~(gayRj2FmAzofDLsR%+barec=KZ;
zzpd$V)Lp+7{}1<)Olg&}@nI1`f(3Zv%N9HAiJu8WTOHSXf*xL);e})mKq9w1fxQ<+
zH6LGuBgJTJe<@!5f<{$Qg6*qF;9UtrKf@EyNe5Yg70N~#MOz}e=fu17#0XM?V$GxN
zsMC}HC4&llq?^-5tv1go!TS1#^ACMiZKy^tO#TIhfzS~LlwC+>q}s>N@NGdx;8=*$
zM80X3#Ahk93`_mRm8Xy4B}(7$>~HL};zCsLDq-A1&I`$$r{8?q<mz${lL$Y@-lhIA
zZZb;8#l-OPsxE&AY7jC_QJOr&w4UvSmKXJ*Hzpjb!s%;`MRg88Rt_4=EZM}mYg^I`
zTrjx6cmrE=lo@)oeT=n2IoM=symkx<IRlhSu8#46SlSi-3bKp=LHJ7lPU8g-BQa`$
z8NTsnGH*E8Ei$J!c1n2`;KgR0$Ii8k+U>k^(3A0<-`*`H>n5M+?Y?*9GJ9IT^XFeF
zIK~qeaZ}NY7~?DBru?1G3zdEb&Pc@osrK=*pAXD@Z@%Xb-C)C}5WBCQV|5@x2$@}Z
zN`#F1^0hqK$vWOTuz%ih)SJvJk8)C2U}5N#Pk&Tt?mz#NS75-oWRm?McOM&Onms;u
z!sKl@<^amv&%}9z$P$-u@N;D89fo&G&}eQ=;d~HBRZ(gG0(ae|m|}UI202;?`oS>&
za_>vx<-?M(*2}mNx*3)3rAeR}ZY;nj)ML@cCrU!*y1PnSO`po<>}}@5p^!btL>!3a
zuW~+ntIi^b18rQ^3HD}7FL6&Xd1hNxY;biz5<XO;n2=lk_V+LO^X`@<=Zhu8E1TOW
zTo9l6IbfIcoOa+DAiKlYB`l?yfL*YTmeCX-*8Xr*CHbg(miuJ&MwOTJ)fl*y>Lq+`
zeV#bh+<=To=@i6(=+6M|i5j@&l#<i@&LeL1!@ND9|6vyS&=2LZhuHKhtEi<PK7Zjk
zfGr4>MFjs4!e}z@2*t^6chu}x-qfd){d1fc94K-+*UUPY{6VAEyWYA$>;j4ajq(xV
zCu{B~Qqj0944of42-vpr8lA^T77jA^!4^M<uKfz62^V0!^}fU%3z0nv+g0vJHJvaq
z&78zcc1%*F$qgZwcSO5E>=f>R#+k~c?PRhgR3=SWPRcLHH@>+pp+qfE3{%GfDC%(T
z2?*=SqFPhzm)x_WBk_~<p}NobVBe3t^w945gXsw&HPEXe(|sB!IbF!c6uhm2wF9m)
z`((t_o{T*p(sP>GUF0;=1^g>9FKos#ku{oju4f)b4sG!1eRRFDGy6nRV?=_&SajDG
zlh>8m8^%@z)Nn#47Vq=*Wnaa*Wxe293=RxylJcY#!5x&pO0|B`jX`?<Y41QeB<Zz(
zBpQ?SsLPY~D_L~x^G;O3mlr`HN0keb{LboKRt7U^zjb*sZ|cvn{sk`w9VF9Jew_g>
z!-R<@{XJ;GaLDUI5HEfGP+zq(Wcj66^1BQdPacl(rt+Jru$`H(A8rrrP3LAu$5CCH
zN<>-w9D~`;VLIN}+Gxs`ze)N)ga0_>XNl3%gi0v%+Ck(l6GJi8A7zEzCW`*xlMN>;
z_`5u-v>jCbj9#wkIuPvi+yjee^3_On{Qa4{M}6vGT83+_{kBB<#avoEI$-(eK{!xm
z`)U{@-n>3NJ?DLRp~1uW)RMBK^p(l0?@NohB&*7o5HQ3U;Elr)zlQC8n9A-}av6n(
z>_FuWqk%NIMCI*;r{72SqdFXdyO`aAWe<PLoBS^2{<X^=bk$YgylI5an&Io35PK_x
zaMX3hQEcb_ynrTO@SjHHwVDiY_FtFG;RO2!rDIV8H$%ohJ8P*J?$Au-nogHDIn12|
z9tXBj?yamswBYN;Ad6M@XqYNNGx;-F3gs(wr5jutAh%tTeb=7<o6J^jiTkw&v$}uY
z%seuBAeYv;`edbx2PniI4)n~z_81FwpCrfBbv|Q_oQx`v{4p3vhxzgXlB=8^SzRAH
zz)PDGknWGp{E^S^F2wqD=U_OEk&x3Tj2K(HGJaeB(nUohkrabY(l2%M52Boybogou
z6l&Ea4=I<n>hS;Dfgj$+djyDfYtpo}HxT2CUcl~5@PQRHi%ZsIPjg+p-eEpMzF0{8
zGI{dej19oj6Az-J(r)a4Zs*Fu#g<dWJ_Hjyz_M~B^^VJ{1+jRiKwaJ^l1$TF8ZpCn
z1^vH;Fzqy?5YX6ul=?%_E^nT<Sp8D_2S0v3xH~2lr7O;`eRK7R16z)VcNKU)*nJ6+
zH-dU){@TnB>AQ03zdLIsTrK?6tvCqr1|Q$o!4P7ABY|7&Lt$8p2!~okORfla^%wbG
z#n-j(zv{9jajQM$=M`n(HYf5fI*8(11@APCKLcHO=*%1{SO5{F6%nDmU&G6pBh*?{
z{K~ST(p>&+FpO34;JDa#Iz6Cq7)!8d4vQQqgt3LZB;1I1!zIf0@VQoh9mMa12~SAN
zA~(!U)qet4D|>4nY}WM)U1H-1VnOi}#(7%#38cQfTN7~UXUUf1Uq|a&MAVgfCw5z5
z&WxS2)qUZyet=7lSyEpwaXwZz*i$-0kfBF&i3U}mtrqwYOhh+`yBqXQp15~!hdZX8
z)fL&T$vx_9;a}z`YS8XBk`Cnr5B~fk$N)HJaIl9)eNt0B0KpGVpJZ!*KASm46q2uH
zQXQiU^-!J{=%M1%`500?gu(gxWM%c!Z(d)>HTsSChTT@T2*=QSsgSBtsm)h~akt+6
ztKBCk>HkEGXt(lzLVSVU9~(oeIuhZ>A<1gbo|m*|cH|{z@2z7-%!@;9%ZdG+^H)wY
zHCe@DDT0^)bnv?=-8$qN%1RH(9fG^}L(_m^nX=s$QljhixMcHEe6R+EH3V@-bufMQ
zL<1)qQ6EM%Uu#jlASR;Y2+6+3drjT=kKGgV!XYuay3RQrj`uB@s;eoEOC<H_b#+iT
zHu{R4kDZSu1EO&fPfQo+pJ@=4{I2cavHT-JW){YEd<L+GaT2VW7K2-|A?zj@6eD~8
z=EmnD=8V?c2BOs$;#7H^)P#MnRaffg8Y7+`_&XF0R<+oIriM?rne5Yj?d?PDg97LN
zlH8rC;Tkk7ZDBh)jC~<}LA2++PRXTHt2=^+3x6rEKlXaVE=;xhFMy&1`6ns7|L~7f
z>QN<?4gt7zP3>B#y^Q!xQhS)}01!e}gGVsf&4tl}-SIccW~aH5f7xF_G_kzwd=Jn2
zX)@h>@W+!@Rdtxv4JhyRAqU}Ki{1`CYR~dF$mDA<lnBMlh=cuV$n1@pv>p;)$c%(U
zRCKF`cud_i_OYfUqIT+Kh0m0xGsNNUle=3#pJBdlJ)!AWTF9pb3YC!`5ovdvyV^Oc
zDaakVlf3lc^XV+5Hv%|H*naUOcKXY_UXxrtk<^c#@;n1@bufARV+t`QU&=D-n?XCq
zp(m%O?8TZ)WGxigmmt2d6232er=Fc#ro4FM^nL2JVRP|Zm!~`lkz%w-pSbomPZLgD
z4tjGZT=}yD`FNc@OV?igIAv7IQJ*YB-}dlI3&w04$Ra6BtW}o$kl-L??D`o1a3YD*
z`bQD+Yq7mT7XqA+|LB3oKYH-e{MwM%O;@%x>OV*|^*@X|xqSGLkW^Gel<ZW0vg5dx
zb)_<}6tmvYP&2Ll%>AaU(9CxFLSU7chV~<YR~I&iSz!TVHdCM1zLl>^f=tHDw*TbA
z?O#x}PrtdX9n#SOYALbUHhrrpKB2fNdH9-8ajBo`K%D_%#84FnM!X{WMvbsBVUBj3
zfmzPtoxl<@OS>^Q&Mi`$a1<@R!}f_YtAu>rY7~4$S(z@D-(&4sfWS;hZVRJJpXzcu
zV_{){rPB3KOjts0UUP9+XAoa7C}`FFILPni5%#~u%{y+BZnN4u5=HPs*U+ri-Dtg2
zZO?<#Yxc~Q(nnXW*Mr?ruh6%-#nU?l$`?Qv0+0>_<)Z6H(FfqRYil-+^9Ir`63h3y
zSrXLc^L0{e>V03R$GvV(ttDm3j)32n+4de-lkLOGj|E92wbe{M(<n`afVia#0WcvZ
zNB2b7;|lpFBVu3{!=Nip@(iD=bpCsEdt(#Xy=s<dj_H>7Ktb^9?&vjXg`sQopBm*~
zUcJ)K3rw1fw^N5SsZ3r!(pcYvBDE!o!T!+g8swk%X1ZKWAO9nV`5!*h7{=ZFe6Fn@
z4@#`(&ksKXaFS)pA&@Vazh-mlm>=J+|L!-hJ7<u;t!C*|#?<xpmArwbfMX<XB-v4p
z7=}OgF}?&56i|1Okaq|OHfag`&=^`{pvbn>^7cJ%485LQxb5Tk96nv`-C_9+B=E4!
z{y~Nj`%H9CvPYjv`;Kv9Jw}`!PxYWU_)l}s0LphB2U%wrD?Bk4gCl~go$K!<q~!r{
zq4d6a48%m)cZzLLy%waj?;NddqSdxZi_B-kPjaUG1uUr|uANse{$><Lmgan`Lm&z6
zlnlnYuldGdoG5&7-Xg&pmmeeB&WTn|`??k@o7DT;KToXge-X}gE&rKtw$3KArZ10c
zio5U1K{UgdW``7wk62QSE_d|=m52XB7vox`7agFTLG=wqMIQ1q4j9kfY??`b39Dls
zE;W{*qyu<Kh;HV6V%WC;Dv=~jC6aPjO|9wC{{@jGaSw_P3$zXN>|hnDoZl!NkTDha
zUnvVqtX%bsS^QfVx2nMDTW<N-xBT({i--OXbNFAi)2!zGpJ9=01TXI_XG6cbpZK*w
zd<MTh2`w?1P~43Qc9r7!W_|3y40sh}PAu<leu?%he)btNw>4Qp(kuGr!B2O@;D-yp
zyretp9X{3UcC&FCyN{Sdb^P9$yp@MOoFW~Ac9+JngcW)agv>z=aSAO8MdY&`Kuq!7
z9)5P?YvMhcN_>B&@~6joH7dL|+<`@r0vEgy6R4i_Eg0izNq*;v{+s<zK;`D7OgxFz
z+GwKE@yKay3NdWEsMVrrHGdlhdn~AKog0&|Vlr{mLBYsO1!LS+Mh|i0sa}|qV+dBK
zkmyF@S_54MN|b<=gC;k^Fy0RRMcI(^3tMAT+|Swly`mwde&*aFZ6@@tWUkXxH{u|^
z+7}F;N3ZP}zlPFRc0^ZIRV)NVz4i@gOk(YM8cLbz1B-TJ>ty2H2uB8lsWzM{S^q}F
zVf~orpULHh2bcV0tOAN?N6>%1cZ-)8{%z-ip_GN#eL&|rONp0B-WG(AOs$!t6tFz{
z{@C!wu;9gn3S4z}8j{x{t}D#OFr-7e_22-#x;i~h7&|#4wQ1NpF2gt0NpeI`Pf&h$
zrBT12@VjKr`MKi!$lkie&G+_O2E0|y&z6*Tx)yWW^uWtTKc^&{#`L%x-Uvo0?Q?fw
z(Ss(0E4jG($q*X$n0gSS!`#r9H7#_FJx@WS)41N)I~Tv%ezu@6AJmBdrQ@SZgNZHE
z4GPE`^>4pzimf0_FdgHl4ttEYQHpLUHMm1Pf@lvWj*_FoY9<lMt+M#|7A(KrUrmu6
zZwJ|EA2;NmTqp}&d0?%?xvSrA$@5R;;$&x%b$3kPkp;@Ypo;~B-ZMZPT-IluSz^B-
zOYSvo*_%Bkz+<OwDgD8eZkG4hsIdC$h&R-pgYB_&81TFZA>IKIt&}+-7XI@0pW49h
zw^Ed^W(dijSmzVr3uyr*RJw>PeBz%<5o+#e22=AX{L#jPMI5y+W)N;iN!}RFtcsQE
z!IH<2s(HWqZ4`uF(8MzI41l{shJ)?xJ%d3TX8@86(ajbG6PhEg_gHO(7?nEexZ^({
za^j~U@%^Iwr5^*`f^v(5ZVTvn<aY~A4_Xql@d66=(H3nr+4dQ$ilvY`qxiwQB2J%L
znspkc6(lUhR37GB_(e1#S-E4`_Mx<$WyA|Ymps~88=FGK%_)<q6aB26(qa%EJuArX
z2hK>+*V=s{G`B&-4?gSYR$bhSRs-*X?2e82zDd{mB(B`}nh{Rtqs}K(<Uu|aU#qz@
zCHFB9KZnCM^5$r1b7&ao4!KPUeQL~3L-~!CjRQ*JPEBoz5$sV=vPEEU$@GkAaot!f
z&i9evGvNyfHZC#8UaFoCvhg_3J1h;*eGzP=Aql@?T-8q<c+5w0-|sQqVrc^xpIQe4
zIF<jQoLezH*vndYQKBslOU)Q?J<X?ZVB*)48;(isV~3e~eM3dAb03Ena}KMLsR2|i
z7I?QQX9#fTb30;Oi9hDRLBj&Mc*(#e?2<g|PO6ID8pec8pN^V0sU)C=b*S52wk9ox
zs&B4o_nK&FW%fdwb~{A>M>O|CqNFc(k5D9t1IS!4gLbg{n{^!(>l3_lE<Ri1qurgg
z8k)xbYdKs_`hWxr7&8UQ=ifObYD&_I*pl(xyt5pzy{UA(u~{{^0!43H+PX65uyFna
zTs^1&(UBhzP~B&M$axuMaA5Ev)%dd?{PL5<NFXoQX{%CHV%X6@d6O1!b9=fhuIMnn
zxNQ!xy7D!#;g+XlN*fY|glC`wXBU8X@RiJu3DL1i#yY(hSJm5dzUOM`2kWOCT-%Oj
z7Cv^KR$EaU@mQ8Dp}P06iqf)Of?`M|cINJKkCV7A#nWe7e-i7)lQSiGj-St~9#<#$
zS2;CL9_Jq*w@DN`UY`p$3rn?iIUmvd3!%Stv+&J@!x0+T$ca!wofiEu<#sEw)o5W=
z<{DCQe+)=^wDZeJ5b`5O{+ylxhw|9Ks-rl_C`943KsrPgmrNfZ1LDN4<_WYS4<IT%
zx6z*_A$$w^yPmg;R-XS-8ED*`ZQHT2@_uxrLZvKDgQiQ-Ya~JUf)L^FlkI%Oq8=||
z&D0%dy#|uaD5MhY{~~h2Z)1G0$x5p3b1BUrfXaZ14@qKR7wrO+2Y=+LMA#COi(dFs
zgpucS>2rIdiX?jm0_i&E3LAXUgG_x)TDN~*UF>z$jdSM^vo*5`zxI}v<$xEkw{60@
zr+5aqmh}YTfSoHLUmgR};}g##@kcSOkOd>PdK+n(XX8F|@)Ex#m1e(*9;Y^wF`^F!
z(2Y{Ri-cWkLKLXuPjc&oiR~Cu<E)2{>0e@Sk-s&a9p@i}`FloqEiSQb%x+r#H8|It
z4AACIY-c5Ab~d2=8s>!&lhOz48sMgJj$l3=DGPmmwF}Q5{l$A{&{U`14hBVl_*)qY
zZ3Bdis0MMxYFa?><}Rz=@Tq36*j`1%aMrXvFbVU}Y)2T+{dvgE+&n)2R?GN&hH*#x
zI;z>~>Lrh>^ZulLvh;5IP@XTYgXM*$Dlu>&K(sb%u(z4MXr05cVyU3L!CE(iJK?nd
zE@|?^F6}vWzP-4E{HRC0)Fg^_r5(rOKrl;CRl0IS@eI>gmFG(O;@(v~wxuLR7NhTg
z)X&e;l5bR69HH%PLr<vbkv~b?!E>i*7~<co70YiT|3bh1E4^dT#%o3UYl_uSMw{9q
z9w1q|<~8CY-2#`iIyql3Yj{sn>$32Z{n0P&&Fb1CY>>{`8+*p*TA@(cqw9%3o&yEz
zr(gKhfYID<eEI2u%Z<;k4^C8m))DRebaY(IdC4-fkzs2g)4D$G-gvw48e*)iGO9^K
z=PPD7N~KW~wj-<uh_qZ0^RN1(eUmzs#qFfkGXPtlB!}4qwz6^$Q8|vHoG*EA{f%gS
z2qI=LLeJX~%XJ7NhZ*sHX?I$>!yqfe*PLozzWI6mB^QuRJV|npnwj_}vt1}2-~Z#I
z%NUgjPRi)-QGs!_<Q*(D*#F()P33iZq!N9H4gdD@)2+qfp2z98LU!h<8uc>Y7t^%S
z>1`l>>Vqw*{6r9y(7W}coc1ldnE~jt7M)V;28=Plky=_onCp*uwpVWpbOE!$olB-a
zJ*?baGXk>ph7`8k{&ik;xlZdj>p5QnceqHntWwwo;@fRA2gmteU`8L3xzRQ4AULzJ
z*G7S1;_4e~&%g@txK|dC5Y6E2BXNk3g{f;Ds>39HFBGMN%cl=eu%$+ybG5qLg=#6I
zQ_$Id2xR#S%o$*lID9XOW=;3G0l%)OD2se5<G#{~07<eW#)&Y}`JC{+M)0zOwIzrv
z{4gXQE=L(!h8`bn%+x|v;Ul3s6s@hW_3|(UyDr2X1D!(B5qM|L-^wXn)24xX_WDv&
z>{2iupmQORFCj`1|2txhHI$!lxwE;yhg##bKrI)Mb327E7}N{g*p-%-;BZalK|kFp
zi23O-nB(&k{)))67hp~GIWHSh)83otbDAq)ay9X`)z#IPq!P|~+<WodRnThgojNt$
z?_0UU_)1%D7Q*<dMa!<q8p-jeS+sg&^Y*%b&7K{mq2z~ir5(IQlb(opR5uppI>Leh
zy5q5z4B&Qa8SEDxWIhTl*x}DGH%wR&Y(?piB=r{<$qIzmz-tXUT}6*4_2U5~P2EvH
zlFFx!6>P2RQ0v6VCES-|KL~VY0E~gCOn%XFTRF8i^2ZDsn?_W&UaoWuP(JE?H3J=$
zs#^S>w>_j?6Ug=H@-s-k2+yB$=hYNFC`?mzLacw)PY*IvzN{DeA*+BBP~bhvll+@b
zgHemiw0CMG*Wuw>926O1{$aq5bu>i?Y>mhh8&fKOTvJ{D;{K@nMp|IDPLc&|HqQk7
znaID$>Nu!-n$s#XejlwBE96)|^>P`Jr0p++7W<{<BJOFRO2;yuckT6uR3p7gxITmI
z5Ol5iiD!;8KcaqnHw~0&(z05;3;G{0e<!rsNbw!^jo3v8tq0*-TXj%wjE8>%;oOKo
z0V4JKhw$N1@@0(DR(>sDvVQxm;MH;I&yyi{zCb6gY=g$He{UYgzIL*=YRL6xEu!%X
ze8Pb)EA}@<n-(@?LpisJn58wr5u4sE44L1GEW4-qa16#_v=A0)he%L_xQ<h{?}wrr
zBC1o@rl!9?ZZe12iq5{m1Uf2o!ZqQc|D8NLYY_s8ayg&)SnE<4Ug$~~A!dVGcJ(e|
zco#SuSTeACM5a>xYBP2uy<(XGaWu*gvA*{pot5jT1>~E*4oz8BG}HCrrG`)L{wrZY
z`hw^9{ycF%%DkcicxnWZ)SEte%PLb+js;t5y_k^olr8B4Vm{gF^?x<EQMQm+_iy=j
zvAQzoBcaGVr3p_MwfKwBTv^|c1})I{S*N>@`sBB-<Xy2(l)ZGTTc2}oiqndv<@{OX
zFM3erezUz&wnL};1svYav{qy=Mqwqwq+N|wOm{d_!(RM-JrgU80uJ4?57Zre&LVcR
z#s7StDmoCZN`x#*9>Z8eYzY<a5O5LtZ@Je1^jt;eC2D*OO&|E2G;C>8_V(519F|lQ
znex_Ddsxy4h2ihxXMSnHb?v7eI7-Pow)z%4h#m}6X=!R0Q(T4qsYe8tN8RHRj)d)}
zYFsQW3eONtWiXdy=Yd;~Uk?_ViSL1WJcNcW#Pfq2OqjbeP?xT$kTiCu56o8WAA8>R
zas<U29o&%2GQ6`*yg^NE{RM`QxhhNE747IB`oKTm>8vA#JO{ep=A9;r6^N!Zvfvg~
zmG)&jN8aWMw-YU1piOei#M(8u<>muYQN|ey7{l?+!Hn2HJeGmW@hr?N*JO`&9Ki%-
z>+A=_+>)QYb;o3JZ^bc{S$SP4Q)Z5}s`Knk0_zf%jwc{D>s{<F7&RD_038Uc+dOK4
zA@&cDzoP99gd(w5WhoqJ6SQnJP(!oPQNLk_{`5;K)|NZ@*{<nFQ664VFJ)%FUHnv~
z)p%Ua6PH>M4aLz4YVz&XYBjN=FNJXho{#{GzO3!=dS2O9t#XG|Zu73gT?4R$sr(I|
z=g-DIoLr6TEsV3WZVcqNAR--_2*rq~=aW^K?kQybgq>T~eP)GM&!<(TljH*h%CE%T
zG?0En&C7WDuNd`z+8+y)e!l7y2ap_X<=t876JjT8$;1{K7j>6BEWTj`XLmm`mCs~v
zFt`yuM~#8d$^lCIZcC21dt2Q6CTZYsZfs5f72U2Hfa=N;%QDHVkgog`_&omf!rW55
zyPlZxZk@~K5nH>3j7Qk{kf5F9H$W|!_`(rt@4x@8cj$l7d+nxy0WzX(e4@&2HiTLQ
z1pIe~X_}2;)Y=vo7E?pie>k-tTq%!o7*&dy^xR*sS)==w#<|mcap2mOmmHleKi;x>
zosaXBsIvZEwLLU}As`ZcJLuOM&`F7<gD(EimS0b4%eJU>#kX&H@f--OTbmH$Wn!sz
zqB&vs(^tEhp>kEeyej7Uq2Cue$7i2!KIp#bT;~$=B!|UXNS94wD_pt$UQ)eB8CCwb
NHu`_Ip%~6){x5rUsyzSz

literal 0
HcmV?d00001

diff --git a/src/fastq.rs b/src/fastq.rs
index deb2ab9..4d25d62 100644
--- a/src/fastq.rs
+++ b/src/fastq.rs
@@ -32,6 +32,34 @@ pub struct IlluminaHeader {
 }
 
 impl Header {
+    pub fn instrument_name(&self) -> String {
+        match self {
+            Header::Casava18(h) => h.instrument_name.clone(),
+            Header::Illumina(h) => h.instrument_name.clone(),
+        }
+    }
+
+    pub fn flowcell_id(&self) -> Option<String> {
+        match self {
+            Header::Casava18(h) => Some(h.flowcell_id.clone()),
+            Header::Illumina(_) => None,
+        }
+    }
+
+    pub fn flowcell_lane(&self) -> u32 {
+        match self {
+            Header::Casava18(h) => h.flowcell_lane,
+            Header::Illumina(h) => h.flowcell_lane,
+        }
+    }
+
+    pub fn pair_member(&self) -> Pair {
+        match self {
+            Header::Casava18(h) => h.pair_member.clone(),
+            Header::Illumina(h) => h.pair_member.clone(),
+        }
+    }
+
     pub fn scramble(self) -> Self {
         fn number(value: u32) -> u32 {
             value % 3 + value % 17 + value % 271 + value % 911
@@ -213,7 +241,7 @@ impl FromStr for Header {
     }
 }
 
-#[derive(Debug, PartialEq)]
+#[derive(Clone, Debug, PartialEq)]
 pub enum Pair {
     PairedEnd = 1,
     MatePair = 2,
diff --git a/src/main.rs b/src/main.rs
index 769560a..4c98b95 100644
--- a/src/main.rs
+++ b/src/main.rs
@@ -2,8 +2,10 @@ mod cli;
 mod fastq;
 
 use crate::cli::{Args, Command};
-use crate::fastq::Header;
+use crate::fastq::{Header, Pair};
 use clap::Parser;
+use console::Style;
+use itertools::Itertools;
 use regex::Regex;
 
 fn scramble_sequence(value: &str, seed: u32) -> String {
@@ -53,9 +55,7 @@ fn main() {
     let args = Args::parse();
 
     match &args.command {
-        Command::Info => {
-            println!("Not implemented yet");
-        }
+        Command::Info => info(),
         Command::Scramble => scramble(),
     }
 
@@ -87,3 +87,171 @@ fn scramble() {
         buf.clear();
     }
 }
+
+fn info() {
+    let stdin = std::io::stdin();
+    let mut buf = String::new();
+
+    let mut headers = vec![];
+    let mut read_lens = vec![];
+    let mut quality_lens = vec![];
+
+    let headline_style = Style::new().bold();
+    let info_style = Style::new().bold().blue();
+    let error_style = Style::new().bold().red();
+
+    let mut line = 1;
+    while let Ok(n) = stdin.read_line(&mut buf) {
+        if n == 0 {
+            break;
+        }
+
+        if buf.starts_with("@") {
+            if let Ok(header) = buf.parse::<Header>() {
+                headers.push(header)
+            } else {
+                println!(
+                    "{}",
+                    error_style.apply_to(format!("🔥 Invalid header at line {}", line))
+                );
+            }
+        } else if buf.starts_with("+") {
+            // ignore optional description
+        } else if line % 4 == 0 {
+            // check if quality values differs from sequence values
+            if Some(&buf.trim().len()) != read_lens.last() {
+                println!(
+                    "{}",
+                    error_style
+                        .apply_to(format!("🔥 Invalid quality string length at line {}", line))
+                );
+                return;
+            }
+            quality_lens.push(buf.trim().len());
+        } else if line % 4 == 2 {
+            read_lens.push(buf.trim().len());
+        }
+
+        line += 1;
+        buf.clear();
+    }
+
+    if line % 4 != 1 {
+        println!(
+            "{}",
+            error_style.apply_to("🔥 File contains invalid or incomplete sequences")
+        );
+        return;
+    }
+
+    println!(
+        "{} {}",
+        info_style.apply_to("🛈 "),
+        headline_style.apply_to(format!("Found {} complete sequence sets", headers.len()))
+    );
+
+    // Instruments
+
+    println!(
+        "{} {}",
+        info_style.apply_to("🛈 "),
+        headline_style.apply_to("Unique instrument name(s):")
+    );
+    println!(
+        "{}",
+        headers
+            .iter()
+            .map(|header| header.instrument_name())
+            .sorted()
+            .chunk_by(|value| value.clone())
+            .into_iter()
+            .map(|g| format!("   {} ({})", g.0, g.1.count()))
+            .collect::<Vec<String>>()
+            .join("\n")
+    );
+
+    // Flowcell IDs
+
+    println!(
+        "{} {}",
+        info_style.apply_to("🛈 "),
+        headline_style.apply_to("Flowcell ID(s):")
+    );
+    println!(
+        "{}",
+        headers
+            .iter()
+            .filter_map(|header| header.flowcell_id())
+            .sorted()
+            .chunk_by(|value| value.clone())
+            .into_iter()
+            .map(|g| format!("   {} ({})", g.0, g.1.count()))
+            .collect::<Vec<String>>()
+            .join("\n")
+    );
+
+    // Flowcell Lanes
+
+    println!(
+        "{} {}",
+        info_style.apply_to("🛈 "),
+        headline_style.apply_to("Flowcell lane(s):")
+    );
+
+    println!(
+        "{}",
+        headers
+            .iter()
+            .map(|header| header.flowcell_lane())
+            .sorted()
+            .chunk_by(|value| value.to_string())
+            .into_iter()
+            .map(|g| format!("   {} ({})", g.0, g.1.count()))
+            .collect::<Vec<String>>()
+            .join("\n")
+    );
+
+    // Read Orders
+
+    println!(
+        "{} {}",
+        info_style.apply_to("🛈 "),
+        headline_style.apply_to("Read order(s):")
+    );
+
+    println!(
+        "{}",
+        headers
+            .iter()
+            .map(|header| match header.pair_member() {
+                Pair::PairedEnd => "R1",
+                Pair::MatePair => "R2",
+            })
+            .sorted()
+            .chunk_by(|value| value.to_string())
+            .into_iter()
+            .map(|g| format!("   {} ({})", g.0, g.1.count()))
+            .collect::<Vec<String>>()
+            .join("\n")
+    );
+
+    // Read Lengths
+
+    println!(
+        "{} {}",
+        info_style.apply_to("🛈 "),
+        headline_style.apply_to("Read length(s):")
+    );
+
+    println!(
+        "{}",
+        read_lens
+            .iter()
+            .sorted()
+            .chunk_by(|value| value.to_string())
+            .into_iter()
+            .map(|g| format!("   {} ({})", g.0, g.1.count()))
+            .collect::<Vec<String>>()
+            .join("\n")
+    )
+}