Repository 'krakentools_combine_kreports'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/krakentools_combine_kreports

Changeset 2:e0d3f5969458 (2023-03-15)
Previous changeset 1:c26130696a4e (2023-03-06) Next changeset 3:073c867c5244 (2023-06-20)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/krakentools commit 92c58a65005708fbd56b27b86969c01a1258ddb9
added:
test-data/extract_kraken_reads/R1.fq.gz
test-data/extract_kraken_reads/R2.fq.gz
test-data/extract_kraken_reads/kraken1.report
test-data/extract_kraken_reads/kraken1.results
test-data/extract_kraken_reads/kraken2.report
test-data/extract_kraken_reads/kraken2.results
test-data/extract_kraken_reads/krakenU.report
test-data/extract_kraken_reads/krona.k2.int.tsv
test-data/extract_kraken_reads/krona.k2.noint.tsv
test-data/extract_kraken_reads/out1.k1.11176.children.fa
test-data/extract_kraken_reads/out1.k1.11176.fa
test-data/extract_kraken_reads/out1.k1.e10386.children.fa
test-data/extract_kraken_reads/out1.k2.11176.fa
test-data/extract_kraken_reads/out1.k2.11176.fq
test-data/extract_kraken_reads/out1.k2.11176.max2.fa
test-data/extract_kraken_reads/out1.k2.11176.parents.fa
test-data/extract_kraken_reads/out1.k2.exclude_both.fa
test-data/extract_kraken_reads/out2.k1.11176.children.fa
test-data/extract_kraken_reads/out2.k1.11176.fa
test-data/extract_kraken_reads/out2.k1.e10386.children.fa
test-data/extract_kraken_reads/out2.k2.11176.fa
test-data/extract_kraken_reads/out2.k2.11176.max2.fa
test-data/extract_kraken_reads/out2.k2.11176.parents.fa
test-data/extract_kraken_reads/out2.k2.exclude_both.fa
test-data/extract_kraken_reads/read1.fastq
test-data/extract_kraken_reads/read2.fastq
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/R1.fq.gz
b
Binary file test-data/extract_kraken_reads/R1.fq.gz has changed
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/R2.fq.gz
b
Binary file test-data/extract_kraken_reads/R2.fq.gz has changed
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/kraken1.report
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/kraken1.report Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,18 @@
+  6.79 38 38 U 0 unclassified
+ 93.21 522 0 - 1 root
+ 93.21 522 0 D 10239   Viruses
+ 88.75 497 0 - 35237     dsDNA viruses, no RNA stage
+ 88.75 497 0 O 548681       Herpesvirales
+ 88.75 497 0 F 10292         Herpesviridae
+ 88.75 497 0 - 10293           Alphaherpesvirinae
+ 88.75 497 0 G 180255             Iltovirus
+ 88.75 497 497 S 10386               Gallid herpesvirus 1
+  4.46 25 0 - 439488     ssRNA viruses
+  4.46 25 0 - 35301       ssRNA negative-strand viruses
+  4.46 25 0 O 11157         Mononegavirales
+  4.46 25 0 F 11158           Paramyxoviridae
+  4.46 25 0 G 260963             Avulavirus
+  3.57 20 0 - 336096               unclassified Avulavirus
+  3.57 20 20 S 204987                 Goose paramyxovirus SF02
+  0.89 5 0 S 11176               Newcastle disease virus
+  0.89 5 5 - 139270                 Newcastle disease virus B1
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/kraken1.results
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/kraken1.results Wed Mar 15 18:18:07 2023 +0000
b
b'@@ -0,0 +1,560 @@\n+U\tC_II_XVIIIb-101\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-99\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-97\t204987\t301\t0:120 A:31 0:72 204987:2 0:46\n+C\tC_II_XVIIIb-95\t204987\t301\t0:87 204987:24 0:9 A:31 0:120\n+U\tC_II_XVIIIb-93\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-91\t204987\t301\t0:80 204987:2 0:38 A:31 0:120\n+C\tC_II_XVIIIb-89\t204987\t301\t0:120 A:31 0:60 204987:2 0:58\n+C\tC_II_XVIIIb-87\t139270\t301\t0:58 139270:9 0:53 A:31 0:120\n+U\tC_II_XVIIIb-85\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-83\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-81\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-79\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-77\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-75\t204987\t301\t0:120 A:31 0:54 204987:18 0:48\n+U\tC_II_XVIIIb-73\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-71\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-69\t204987\t301\t0:26 204987:2 0:92 A:31 0:120\n+U\tC_II_XVIIIb-67\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-65\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-63\t204987\t301\t0:43 139270:3 0:74 A:31 0:58 204987:14 0:48\n+C\tC_II_XVIIIb-61\t204987\t301\t0:71 204987:22 0:27 A:31 0:108 204987:12\n+U\tC_II_XVIIIb-59\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-57\t204987\t301\t0:34 204987:14 0:72 A:31 0:120\n+U\tC_II_XVIIIb-55\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-53\t139270\t301\t0:120 A:31 0:91 139270:2 0:27\n+C\tC_II_XVIIIb-51\t204987\t301\t0:9 204987:16 0:95 A:31 0:72 139270:3 0:45\n+C\tC_II_XVIIIb-49\t204987\t301\t0:41 204987:2 0:77 A:31 0:120\n+U\tC_II_XVIIIb-47\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-45\t204987\t301\t0:65 204987:12 260963:2 0:41 A:31 0:120\n+C\tC_II_XVIIIb-43\t204987\t301\t0:120 A:31 0:23 204987:3 0:94\n+U\tC_II_XVIIIb-41\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-39\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-37\t139270\t301\t0:120 A:31 0:116 139270:4\n+C\tC_II_XVIIIb-35\t204987\t301\t0:120 A:31 0:45 204987:11 0:64\n+C\tC_II_XVIIIb-33\t204987\t301\t0:32 204987:3 0:85 A:31 0:46 204987:18 0:56\n+U\tC_II_XVIIIb-31\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-29\t204987\t301\t0:120 A:31 0:82 204987:9 0:29\n+U\tC_II_XVIIIb-27\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-25\t139270\t301\t0:120 A:31 0:88 139270:1 0:31\n+C\tC_II_XVIIIb-23\t139270\t301\t0:44 139270:2 0:74 A:31 0:120\n+U\tC_II_XVIIIb-21\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-19\t204987\t301\t0:34 204987:9 0:77 A:31 0:23 204987:13 0:84\n+U\tC_II_XVIIIb-17\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-15\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-13\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-11\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-9\t0\t301\t0:120 A:31 0:120\n+U\tC_II_XVIIIb-7\t0\t301\t0:120 A:31 0:120\n+C\tC_II_XVIIIb-5\t204987\t301\t0:120 A:31 0:77 204987:28 0:15\n+C\tC_II_XVIIIb-3\t204987\t301\t0:99 204987:2 0:19 A:31 0:105 204987:2 0:13\n+C\tC_II_XVIIIb-1\t204987\t301\t0:109 204987:5 0:6 A:31 0:103 204987:5 0:12\n+C\tHQ630064-1017\t10386\t301\t0:32 10386:88 A:31 0:32 10386:23 0:31 10386:34\n+C\tHQ630064-1015\t10386\t301\t0:48 10386:72 A:31 0:22 10386:98\n+C\tHQ630064-1013\t10386\t301\t0:19 10386:101 A:31 0:22 10386:98\n+C\tHQ630064-1011\t10386\t301\t0:31 10386:89 A:31 0:50 10386:8 0:31 10386:31\n+C\tHQ630064-1009\t10386\t301\t0:31 10386:89 A:31 0:39 10386:65 0:16\n+C\tHQ630064-1007\t10386\t301\t0:1 10386:2 0:31 10386:86 A:31 0:28 10386:92\n+C\tHQ630064-1005\t10386\t301\t0:22 10386:16 0:31 10386:51 A:31 0:18 10386:2 0:40 10386:12 0:48\n+C\tHQ630064-1003\t10386\t301\t0:31 10386:89 A:31 0:99 10386:21\n+C\tHQ630064-1001\t10386\t301\t0:30 10386:17 0:31 10386:42 A:31 0:34 10386:77 0:9\n+C\tHQ630064-999\t10386\t301\t0:1 10386:1 0:55 10386:63 A:31 0:37 10386:83\n+C\tHQ630064-997\t10386\t301\t0:37 10386:61 0:22 A:31 0:19 10386:31 0:52 10386:18\n+C\tHQ630064-995\t10386\t301\t0:32 10386:1 0:31 10386:56 A:31 0:33 10386:53 0:31 10386:3\n+C\tHQ630064-993\t10386\t301\t0:30 10386:90 A:31 0:39 10386:81\n+C\tHQ630064-991\t10386\t301\t0:44 10386:46 0:25 10386:5 A:31 0:55 10386:18 0:47\n+C\tHQ630064-989\t10386\t301\t0:30 10386:67 0:23 A:31 0:91 10386:29\n+C\tHQ630064-987\t10386\t301\t0:67 10386:8 0:31 10386:1 0:13 A:31 0:20 10386:41 0:31 10386:24 0:4\n+C\tHQ630064-985\t10386\t301\t0:29 10386:91 A:31 0:35 10386:30 0:31 10386:8 0:16\n+C\tHQ630064-983\t10386\t301\t0:24 10386:35 0:44 10386:17 A:31 0:24 10386:13 0:38 10386:45\n+C\tHQ630064-9'..b' 0:86 10386:34\n+C\tHQ630064-119\t10386\t301\t0:19 10386:101 A:31 0:43 10386:77\n+C\tHQ630064-117\t10386\t301\t0:20 10386:100 A:31 0:25 10386:39 0:56\n+C\tHQ630064-115\t10386\t301\t0:18 10386:32 0:32 10386:8 0:30 A:31 0:7 10386:14 0:31 10386:1 0:31 10386:8 0:28\n+C\tHQ630064-113\t10386\t301\t0:35 10386:85 A:31 0:42 10386:78\n+U\tHQ630064-111\t0\t301\t0:120 A:31 0:120\n+C\tHQ630064-109\t10386\t301\t0:27 10386:93 A:31 0:43 10386:77\n+C\tHQ630064-107\t10386\t301\t0:23 10386:97 A:31 0:11 10386:17 0:31 10386:61\n+U\tHQ630064-105\t0\t301\t0:120 A:31 0:120\n+C\tHQ630064-103\t10386\t301\t0:33 10386:87 A:31 0:31 10386:23 0:31 10386:35\n+C\tHQ630064-101\t10386\t301\t0:42 10386:78 A:31 0:77 10386:43\n+C\tHQ630064-99\t10386\t301\t0:15 10386:105 A:31 0:30 10386:47 0:31 10386:12\n+C\tHQ630064-97\t10386\t301\t0:32 10386:88 A:31 0:28 10386:18 0:31 10386:43\n+C\tHQ630064-95\t10386\t301\t0:26 10386:94 A:31 0:45 10386:5 0:31 10386:39\n+C\tHQ630064-93\t10386\t301\t0:35 10386:52 0:33 A:31 0:20 10386:10 0:31 10386:59\n+C\tHQ630064-91\t10386\t301\t0:18 10386:4 0:31 10386:67 A:31 0:35 10386:85\n+C\tHQ630064-89\t10386\t301\t0:12 10386:108 A:31 0:102 10386:18\n+C\tHQ630064-87\t10386\t301\t0:40 10386:80 A:31 0:62 10386:8 0:35 10386:15\n+C\tHQ630064-85\t10386\t301\t0:33 10386:87 A:31 0:56 10386:52 0:12\n+C\tHQ630064-83\t10386\t301\t0:35 10386:85 A:31 10386:120\n+C\tHQ630064-81\t10386\t301\t0:26 10386:94 A:31 0:68 10386:52\n+C\tHQ630064-79\t10386\t301\t0:35 10386:85 A:31 0:59 10386:61\n+C\tHQ630064-77\t10386\t301\t0:35 10386:85 A:31 0:43 10386:77\n+C\tHQ630064-75\t10386\t301\t0:30 10386:90 A:31 0:32 10386:17 0:45 10386:26\n+C\tHQ630064-73\t10386\t301\t0:16 10386:22 0:31 10386:51 A:31 10386:54 0:31 10386:35\n+C\tHQ630064-71\t10386\t301\t0:46 10386:74 A:31 0:19 10386:22 0:31 10386:44 0:4\n+C\tHQ630064-69\t10386\t301\t0:34 10386:36 0:31 10386:19 A:31 0:11 10386:109\n+C\tHQ630064-67\t10386\t301\t0:29 10386:91 A:31 0:31 10386:89\n+C\tHQ630064-65\t10386\t301\t0:18 10386:102 A:31 0:35 10386:85\n+C\tHQ630064-63\t10386\t301\t0:23 10386:97 A:31 0:62 10386:21 0:31 10386:6\n+C\tHQ630064-61\t10386\t301\t0:13 10386:107 A:31 0:45 10386:75\n+C\tHQ630064-59\t10386\t301\t0:27 10386:93 A:31 0:27 10386:93\n+C\tHQ630064-57\t10386\t301\t0:32 10386:76 0:12 A:31 0:59 10386:29 0:32\n+C\tHQ630064-55\t10386\t301\t0:30 10386:90 A:31 0:27 10386:93\n+C\tHQ630064-53\t10386\t301\t0:27 10386:93 A:31 0:22 10386:4 0:31 10386:22 0:31 10386:10\n+C\tHQ630064-51\t10386\t301\t0:26 10386:14 0:31 10386:49 A:31 0:31 10386:10 0:73 10386:6\n+C\tHQ630064-49\t10386\t301\t0:34 10386:86 A:31 0:120\n+C\tHQ630064-47\t10386\t301\t0:31 10386:27 0:31 10386:31 A:31 0:35 10386:46 0:31 10386:8\n+C\tHQ630064-45\t10386\t301\t0:21 10386:20 0:31 10386:14 0:34 A:31 0:62 10386:35 0:23\n+C\tHQ630064-43\t10386\t301\t0:33 10386:28 0:48 10386:11 A:31 0:28 10386:19 0:43 10386:30\n+C\tHQ630064-41\t10386\t301\t0:26 10386:15 0:29 10386:50 A:31 0:34 10386:80 0:6\n+C\tHQ630064-39\t10386\t301\t0:44 10386:76 A:31 0:15 10386:105\n+C\tHQ630064-37\t10386\t301\t0:33 10386:87 A:31 0:46 10386:48 0:26\n+C\tHQ630064-35\t10386\t301\t0:35 10386:85 A:31 0:11 10386:109\n+C\tHQ630064-33\t10386\t301\t0:22 10386:98 A:31 0:55 10386:65\n+C\tHQ630064-31\t10386\t301\t0:25 10386:13 0:31 10386:15 0:31 10386:5 A:31 0:24 10386:96\n+C\tHQ630064-29\t10386\t301\t0:27 10386:46 0:32 10386:15 A:31 0:24 10386:53 0:41 10386:2\n+C\tHQ630064-27\t10386\t301\t0:77 10386:43 A:31 0:21 10386:47 0:31 10386:14 0:7\n+C\tHQ630064-25\t10386\t301\t0:35 10386:85 A:31 0:63 10386:57\n+C\tHQ630064-23\t10386\t301\t0:26 10386:35 0:31 10386:28 A:31 0:23 10386:97\n+C\tHQ630064-21\t10386\t301\t0:25 10386:95 A:31 0:28 10386:40 0:31 10386:21\n+C\tHQ630064-19\t10386\t301\t0:32 10386:88 A:31 0:22 10386:98\n+C\tHQ630064-17\t10386\t301\t0:32 10386:88 A:31 0:26 10386:16 0:31 10386:47\n+C\tHQ630064-15\t10386\t301\t0:59 10386:61 A:31 0:21 10386:21 0:39 10386:39\n+C\tHQ630064-13\t10386\t301\t0:33 10386:87 A:31 0:35 10386:85\n+U\tHQ630064-11\t0\t301\t0:120 A:31 0:120\n+C\tHQ630064-9\t10386\t301\t0:22 10386:98 A:31 0:33 10386:1 0:31 10386:55\n+C\tHQ630064-7\t10386\t301\t0:30 10386:90 A:31 0:63 10386:57\n+C\tHQ630064-5\t10386\t301\t0:15 10386:47 0:58 A:31 0:14 10386:28 0:31 10386:4 0:43\n+U\tHQ630064-3\t0\t301\t0:120 A:31 0:120\n+C\tHQ630064-1\t10386\t301\t0:45 10386:31 0:44 A:31 0:120\n'
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/kraken2.report
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/kraken2.report Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,73 @@
+  4.29 24 24 U 0 unclassified
+ 95.71 536 0 R 1 root
+ 95.00 532 0 D 10239   Viruses
+ 88.75 497 0 D1 2731341     Duplodnaviria
+ 88.75 497 0 D2 2731360       Heunggongvirae
+ 88.75 497 0 P 2731361         Peploviricota
+ 88.75 497 0 C 2731363           Herviviricetes
+ 88.75 497 0 O 548681             Herpesvirales
+ 88.75 497 0 F 10292               Herpesviridae
+ 88.75 497 0 F1 10293                 Alphaherpesvirinae
+ 88.75 497 0 G 180255                   Iltovirus
+ 88.75 497 497 S 10386                     Gallid alphaherpesvirus 1
+  6.25 35 0 D1 2559587     Riboviria
+  6.25 35 0 K 2732396       Orthornavirae
+  6.25 35 0 P 2497569         Negarnaviricota
+  6.25 35 0 P1 2497570           Haploviricotina
+  6.25 35 0 C 2497574             Monjiviricetes
+  6.25 35 0 O 11157               Mononegavirales
+  6.25 35 0 F 11158                 Paramyxoviridae
+  6.25 35 12 G 260963                   Avulavirus
+  2.68 15 0 G1 336096                     unclassified Avulavirus
+  2.68 15 15 S 204987                       Goose paramyxovirus SF02
+  1.25 7 7 S 11176                     Avian avulavirus 1
+  0.18 1 1 S 1305716                     Avian avulavirus 12
+  0.71 4 0 R1 131567   cellular organisms
+  0.36 2 0 D 2     Bacteria
+  0.36 2 0 P 1224       Proteobacteria
+  0.36 2 0 C 1236         Gammaproteobacteria
+  0.36 2 0 O 72274           Pseudomonadales
+  0.36 2 0 F 468             Moraxellaceae
+  0.36 2 0 G 469               Acinetobacter
+  0.36 2 0 G1 909768                 Acinetobacter calcoaceticus/baumannii complex
+  0.36 2 2 S 470                   Acinetobacter baumannii
+  0.36 2 0 D 2759     Eukaryota
+  0.18 1 0 D1 33154       Opisthokonta
+  0.18 1 0 K 33208         Metazoa
+  0.18 1 0 K1 6072           Eumetazoa
+  0.18 1 0 K2 33213             Bilateria
+  0.18 1 0 K3 33511               Deuterostomia
+  0.18 1 0 P 7711                 Chordata
+  0.18 1 0 P1 89593                   Craniata
+  0.18 1 0 P2 7742                     Vertebrata
+  0.18 1 0 P3 7776                       Gnathostomata
+  0.18 1 0 P4 117570                         Teleostomi
+  0.18 1 0 P5 117571                           Euteleostomi
+  0.18 1 0 P6 8287                             Sarcopterygii
+  0.18 1 0 P7 1338369                               Dipnotetrapodomorpha
+  0.18 1 0 P8 32523                                 Tetrapoda
+  0.18 1 0 P9 32524                                   Amniota
+  0.18 1 0 C 40674                                     Mammalia
+  0.18 1 0 C1 32525                                       Theria
+  0.18 1 0 C2 9347                                         Eutheria
+  0.18 1 0 C3 1437010                                           Boreoeutheria
+  0.18 1 0 C4 314146                                             Euarchontoglires
+  0.18 1 0 O 9443                                               Primates
+  0.18 1 0 O1 376913                                                 Haplorrhini
+  0.18 1 0 O2 314293                                                   Simiiformes
+  0.18 1 0 O3 9526                                                     Catarrhini
+  0.18 1 0 O4 314295                                                       Hominoidea
+  0.18 1 0 F 9604                                                         Hominidae
+  0.18 1 0 F1 207598                                                           Homininae
+  0.18 1 0 G 9605                                                             Homo
+  0.18 1 1 S 9606                                                               Homo sapiens
+  0.18 1 0 D1 2698737       Sar
+  0.18 1 0 D2 33630         Alveolata
+  0.18 1 0 P 5794           Apicomplexa
+  0.18 1 0 C 422676             Aconoidasida
+  0.18 1 0 O 5819               Haemosporida
+  0.18 1 0 F 1639119                 Plasmodiidae
+  0.18 1 0 G 5820                   Plasmodium
+  0.18 1 0 G1 418103                     Plasmodium (Plasmodium)
+  0.18 1 0 S 5850                       Plasmodium knowlesi
+  0.18 1 1 S1 5851                         Plasmodium knowlesi strain H
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/kraken2.results
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/kraken2.results Wed Mar 15 18:18:07 2023 +0000
b
b'@@ -0,0 +1,560 @@\n+C\tC_II_XVIIIb-101\t5851\t150|150\t0:116 |:| 0:17 5851:5 0:94\n+C\tC_II_XVIIIb-99\t260963\t150|150\t0:40 260963:5 0:71 |:| 0:116\n+C\tC_II_XVIIIb-97\t260963\t150|150\t0:116 |:| 0:75 260963:1 0:40\n+C\tC_II_XVIIIb-95\t204987\t150|150\t0:63 260963:3 0:21 204987:20 0:9 |:| 0:116\n+C\tC_II_XVIIIb-93\t260963\t150|150\t0:116 |:| 0:82 260963:3 0:31\n+C\tC_II_XVIIIb-91\t204987\t150|150\t0:65 260963:2 0:49 |:| 0:18 204987:2 0:96\n+C\tC_II_XVIIIb-89\t204987\t150|150\t0:11 204987:5 0:38 260963:3 0:59 |:| 0:56 260963:5 0:55\n+C\tC_II_XVIIIb-87\t260963\t150|150\t0:23 260963:4 0:58 260963:2 0:1 260963:5 0:23 |:| 0:116\n+C\tC_II_XVIIIb-85\t260963\t150|150\t0:42 260963:3 0:3 260963:5 0:1 260963:5 0:57 |:| 0:116\n+U\tC_II_XVIIIb-83\t0\t150|150\t0:116 |:| 0:116\n+U\tC_II_XVIIIb-81\t0\t150|150\t0:116 |:| 0:116\n+U\tC_II_XVIIIb-79\t0\t150|150\t0:116 |:| 0:116\n+U\tC_II_XVIIIb-77\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-75\t260963\t150|150\t0:59 260963:1 0:7 260963:3 0:7 260963:3 0:36 |:| 0:52 260963:16 0:48\n+C\tC_II_XVIIIb-73\t1305716\t150|150\t0:116 |:| 0:53 1305716:4 0:59\n+U\tC_II_XVIIIb-71\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-69\t204987\t150|150\t0:11 260963:2 0:67 204987:3 0:33 |:| 0:116\n+U\tC_II_XVIIIb-67\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-65\t11176\t150|150\t0:103 11176:3 0:1 11176:2 0:7 |:| 0:31 260963:6 11176:5 0:5 11176:4 260963:1 11176:5 0:8 11176:2 0:1 11176:3 0:45\n+C\tC_II_XVIIIb-63\t204987\t150|150\t0:39 260963:5 0:72 |:| 0:58 204987:12 0:1 204987:5 0:23 260963:3 0:14\n+C\tC_II_XVIIIb-61\t204987\t150|150\t0:45 260963:3 0:23 204987:18 0:27 |:| 0:108 204987:8\n+U\tC_II_XVIIIb-59\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-57\t11176\t150|150\t0:31 260963:13 0:1 260963:6 0:16 11176:1 0:48 |:| 0:102 11176:1 0:13\n+C\tC_II_XVIIIb-55\t204987\t150|150\t0:67 204987:6 0:43 |:| 0:113 204987:3\n+C\tC_II_XVIIIb-53\t204987\t150|150\t0:55 260963:5 0:10 260963:2 204987:5 0:7 260963:5 0:27 |:| 0:28 260963:1 0:65 260963:5 0:17\n+C\tC_II_XVIIIb-51\t204987\t150|150\t0:8 204987:13 0:95 |:| 0:21 260963:4 0:33 260963:3 0:7 260963:5 0:43\n+C\tC_II_XVIIIb-49\t11176\t150|150\t0:35 260963:1 0:80 |:| 0:30 11176:5 0:5 11176:3 0:73\n+U\tC_II_XVIIIb-47\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-45\t11176\t150|150\t0:28 11176:1 0:35 260963:13 0:39 |:| 0:98 11176:7 0:1 260963:5 0:5\n+C\tC_II_XVIIIb-43\t260963\t150|150\t0:99 260963:12 0:5 |:| 0:20 260963:4 0:51 260963:2 0:21 260963:1 0:17\n+C\tC_II_XVIIIb-41\t204987\t150|150\t0:116 |:| 0:77 204987:6 0:33\n+U\tC_II_XVIIIb-39\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-37\t204987\t150|150\t0:110 204987:5 0:1 |:| 0:116\n+C\tC_II_XVIIIb-35\t204987\t150|150\t0:48 204987:5 0:5 204987:1 0:42 260963:4 0:11 |:| 0:19 260963:7 0:17 260963:10 0:56 260963:4 0:3\n+C\tC_II_XVIIIb-33\t260963\t150|150\t0:32 260963:7 0:77 |:| 0:46 260963:16 0:54\n+C\tC_II_XVIIIb-31\t204987\t150|150\t0:113 204987:3 |:| 0:116\n+C\tC_II_XVIIIb-29\t260963\t150|150\t0:116 |:| 0:74 260963:16 0:6 260963:1 0:11 260963:5 0:3\n+U\tC_II_XVIIIb-27\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-25\t204987\t150|150\t0:106 260963:5 0:5 |:| 0:61 204987:7 0:48\n+C\tC_II_XVIIIb-23\t204987\t150|150\t0:47 260963:5 0:31 260963:5 0:7 204987:5 260963:2 0:10 260963:4 |:| 0:57 260963:12 0:9 260963:1 0:37\n+U\tC_II_XVIIIb-21\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-19\t260963\t150|150\t0:24 260963:1 0:6 260963:16 0:69 |:| 0:19 260963:15 0:82\n+C\tC_II_XVIIIb-17\t11176\t150|150\t0:116 |:| 0:79 11176:2 0:7 11176:12 260963:1 11176:6 0:9\n+U\tC_II_XVIIIb-15\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-13\t470\t150|150\t0:116 |:| 0:22 470:2 0:5 470:2 0:37 470:5 0:42 470:1\n+C\tC_II_XVIIIb-11\t470\t150|150\t0:49 470:5 0:42 470:6 0:10 260963:1 0:3 |:| 0:116\n+C\tC_II_XVIIIb-9\t11176\t150|150\t0:116 |:| 0:25 11176:3 0:4 11176:8 0:76\n+U\tC_II_XVIIIb-7\t0\t150|150\t0:116 |:| 0:116\n+C\tC_II_XVIIIb-5\t260963\t150|150\t0:116 |:| 0:73 260963:32 0:11\n+C\tC_II_XVIIIb-3\t260963\t150|150\t0:33 260963:5 0:55 260963:1 0:22 |:| 0:108 260963:1 0:7\n+C\tC_II_XVIIIb-1\t11176\t150|150\t0:36 260963:5 0:11 204987:1 0:54 260963:4 0:5 |:| 0:38 11176:3 0:5 11176:5 0:51 260963:4 0:10\n+C\tHQ630064-1017\t10386\t150|150\t0:29 10386:87 |:| 0:26 10386:33 0:1 10386:5 0:9 10386:3 0:8 10386:31\n+C\tHQ630064'..b'86\t150|150\t0:32 10386:84 |:| 0:45 10386:5 0:1 10386:3 0:3 10386:59\n+C\tHQ630064-77\t10386\t150|150\t0:33 10386:83 |:| 10386:1 0:5 10386:2 0:27 10386:7 0:1 10386:73\n+C\tHQ630064-75\t10386\t150|150\t0:28 10386:88 |:| 0:28 10386:21 0:31 10386:5 0:5 10386:26\n+C\tHQ630064-73\t10386\t150|150\t0:10 10386:1 0:5 10386:20 0:23 10386:1 0:8 10386:48 |:| 10386:50 0:21 10386:11 0:3 10386:31\n+C\tHQ630064-71\t10386\t150|150\t0:28 10386:1 0:17 10386:70 |:| 0:18 10386:19 0:23 10386:3 0:7 10386:46\n+C\tHQ630064-69\t10386\t150|150\t0:31 10386:36 0:4 10386:2 0:28 10386:15 |:| 0:47 10386:69\n+C\tHQ630064-67\t10386\t150|150\t0:23 10386:93 |:| 0:3 10386:1 0:15 10386:7 0:4 10386:86\n+C\tHQ630064-65\t10386\t150|150\t0:16 10386:100 |:| 0:19 10386:5 0:9 10386:83\n+C\tHQ630064-63\t10386\t150|150\t0:23 10386:93 |:| 0:62 10386:20 0:2 10386:5 0:1 10386:5 0:17 10386:4\n+C\tHQ630064-61\t10386\t150|150\t0:7 10386:1 0:4 10386:104 |:| 0:23 10386:1 0:7 10386:1 0:12 10386:72\n+C\tHQ630064-59\t10386\t150|150\t0:19 10386:1 0:5 10386:91 |:| 0:27 10386:89\n+C\tHQ630064-57\t10386\t150|150\t0:31 10386:76 0:6 10386:3 |:| 0:4 10386:1 0:53 10386:26 0:3 10386:5 0:24\n+C\tHQ630064-55\t10386\t150|150\t0:29 10386:87 |:| 0:15 10386:3 0:5 10386:3 0:1 10386:89\n+C\tHQ630064-53\t10386\t150|150\t0:9 10386:1 0:16 10386:90 |:| 0:14 10386:1 0:3 10386:8 0:3 10386:1 0:1 10386:5 0:13 10386:5 0:3 10386:28 0:9 10386:5 0:9 10386:8\n+C\tHQ630064-51\t10386\t150|150\t0:25 10386:12 0:10 10386:1 0:13 10386:2 0:8 10386:45 |:| 0:27 10386:12 0:71 10386:6\n+C\tHQ630064-49\t10386\t150|150\t0:16 10386:5 0:12 10386:83 |:| 0:116\n+C\tHQ630064-47\t10386\t150|150\t0:16 10386:4 0:11 10386:27 0:3 10386:9 0:19 10386:27 |:| 0:19 10386:5 0:7 10386:48 0:31 10386:6\n+C\tHQ630064-45\t10386\t150|150\t0:20 10386:25 0:17 10386:22 0:32 |:| 0:26 10386:5 0:15 10386:1 0:11 10386:37 0:21\n+C\tHQ630064-43\t10386\t150|150\t0:27 10386:1 0:5 10386:28 0:8 10386:2 0:25 10386:2 0:11 10386:7 |:| 0:25 10386:30 0:20 10386:1 0:10 10386:30\n+C\tHQ630064-41\t10386\t150|150\t0:26 10386:15 0:21 10386:5 0:1 10386:48 |:| 0:31 10386:81 0:1 10386:3\n+C\tHQ630064-39\t10386\t150|150\t0:40 10386:76 |:| 0:14 10386:102\n+C\tHQ630064-37\t10386\t150|150\t0:20 10386:6 0:3 10386:87 |:| 0:46 10386:45 0:25\n+C\tHQ630064-35\t10386\t150|150\t0:33 10386:83 |:| 0:7 10386:109\n+C\tHQ630064-33\t10386\t150|150\t0:10 10386:5 0:4 10386:97 |:| 0:13 10386:5 0:33 10386:65\n+C\tHQ630064-31\t10386\t150|150\t0:13 10386:5 0:4 10386:13 0:12 10386:11 0:10 10386:12 0:7 10386:9 0:13 10386:1 0:4 10386:2 |:| 0:18 10386:1 0:4 10386:93\n+C\tHQ630064-29\t10386\t150|150\t0:23 10386:46 0:35 10386:12 |:| 0:12 10386:2 0:10 10386:53 0:1 10386:5 0:33\n+C\tHQ630064-27\t10386\t150|150\t0:24 10386:2 0:6 10386:4 0:17 10386:7 0:1 10386:1 0:1 10386:1 0:10 10386:42 |:| 0:19 10386:46 0:2 10386:5 0:12 10386:4 0:11 10386:13 0:4\n+C\tHQ630064-25\t10386\t150|150\t0:20 10386:2 0:8 10386:86 |:| 0:51 10386:1 0:5 10386:1 0:5 10386:53\n+C\tHQ630064-23\t10386\t150|150\t0:18 10386:5 0:1 10386:34 0:12 10386:1 0:5 10386:5 0:5 10386:5 0:1 10386:24 |:| 0:23 10386:93\n+C\tHQ630064-21\t10386\t150|150\t0:7 10386:1 0:12 10386:2 0:1 10386:93 |:| 0:20 10386:44 0:9 10386:1 0:7 10386:7 0:10 10386:18\n+C\tHQ630064-19\t10386\t150|150\t0:32 10386:84 |:| 0:19 10386:97\n+C\tHQ630064-17\t10386\t150|150\t0:28 10386:88 |:| 0:23 10386:15 0:5 10386:3 0:13 10386:5 0:5 10386:47\n+C\tHQ630064-15\t10386\t150|150\t0:43 10386:1 0:12 10386:60 |:| 0:6 10386:3 0:12 10386:21 0:7 10386:1 0:15 10386:1 0:3 10386:6 0:5 10386:36\n+C\tHQ630064-13\t10386\t150|150\t0:25 10386:1 0:1 10386:89 |:| 0:17 10386:1 0:1 10386:5 0:1 10386:1 0:1 10386:1 0:3 10386:85\n+U\tHQ630064-11\t0\t150|150\t0:116 |:| 0:116\n+C\tHQ630064-9\t10386\t150|150\t0:19 10386:97 |:| 0:15 10386:3 0:5 10386:1 0:1 10386:2 0:16 10386:11 0:9 10386:53\n+C\tHQ630064-7\t10386\t150|150\t0:15 10386:10 0:2 10386:89 |:| 0:35 10386:3 0:23 10386:55\n+C\tHQ630064-5\t10386\t150|150\t0:1 10386:1 0:5 10386:1 0:6 10386:48 0:29 10386:7 0:1 10386:7 0:7 10386:3 |:| 0:13 10386:25 0:5 10386:5 0:7 10386:7 0:1 10386:7 0:3 10386:2 0:28 10386:9 0:3 10386:1\n+U\tHQ630064-3\t0\t150|150\t0:116 |:| 0:116\n+C\tHQ630064-1\t10386\t150|150\t0:35 10386:5 0:5 10386:28 0:4 10386:3 0:36 |:| 0:116\n'
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/krakenU.report
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/krakenU.report Wed Mar 15 18:18:07 2023 +0000
b
b'@@ -0,0 +1,1305 @@\n+# KrakenUniq v0.5.8 DATE:2021-09-21T05:28:21Z DB:./db1 ./db2 ./db3 ./db4 ./db5 ./db6 ./db7 ./db8 ./db9 DB_SIZE:43487608 WD:/tmp/nxf.LS2KtRM0fv\n+# CL:perl /opt/conda/envs/nf/bin/krakenuniq --db db1 --db db2 --db db3 --db db4 --db db5 --db db6 --db db7 --db db8 --db db9 --threads 64 --report-file lowfilt.report.pre --output lowfilt.results --paired lowfilt.qfilt.R1_val_1.rmspikein.fq lowfilt.qfilt.R2_val_2.rmspikein.fq\n+\n+%\treads\ttaxReads\tkmers\tdup\tcov\ttaxID\trank\ttaxName\n+2.18\t1681\t1681\t12617084\t1.2\tNA\t0\tno rank\tunclassified\n+97.82\t75426\t277\t1886797\t1.23\t7.482e-05\t1\tno rank\troot\n+95.62\t73730\t29\t1742828\t1.21\t7.188e-05\t131567\tno rank\t  cellular organisms\n+93.26\t71908\t380\t1704354\t1.21\t9.69e-05\t2\tsuperkingdom\t    Bacteria\n+44.5\t34310\t59\t797434\t1.21\t0.0001228\t1783272\tclade\t      Terrabacteria group\n+29.93\t23076\t82\t557233\t1.2\t0.0002134\t1239\tphylum\t        Firmicutes\n+20.77\t16015\t73\t379200\t1.22\t0.0003627\t91061\tclass\t          Bacilli\n+12.81\t9880\t45\t253369\t1.19\t0.0003098\t1385\torder\t            Bacillales\n+4.103\t3164\t10\t102496\t1.16\t0.00027\t186817\tfamily\t              Bacillaceae\n+4.088\t3152\t963\t100677\t1.16\t0.001024\t1386\tgenus\t                Bacillus\n+2.521\t1944\t1944\t46683\t1.13\t0.07898\t1963032\tspecies\t                  Bacillus intestinalis\n+0.284\t219\t150\t4091\t1.26\t0.0002188\t86661\tspecies group\t                  Bacillus cereus group\n+0.05447\t42\t42\t1083\t1.29\t0.001994\t2026186\tspecies\t                    Bacillus paranthracis\n+0.009078\t7\t7\t226\t1.29\t0.0004816\t2026188\tspecies\t                    Bacillus tropicus\n+0.007781\t6\t6\t208\t1.21\t0.00036\t2026190\tspecies\t                    Bacillus mobilis\n+0.005188\t4\t4\t170\t1.21\t4.36e-05\t1396\tspecies\t                    Bacillus cereus\n+0.003891\t3\t3\t75\t1.29\t0.0001177\t2026189\tspecies\t                    Bacillus albus\n+0.002594\t2\t2\t56\t1.16\t0.0001006\t2026193\tspecies\t                    Bacillus nitratireducens\n+0.002594\t2\t2\t62\t1.27\t0.0002833\t1392\tspecies\t                    Bacillus anthracis\n+0.001297\t1\t1\t44\t1.27\t0.0001063\t2026191\tspecies\t                    Bacillus luti\n+0.001297\t1\t1\t66\t1.26\t8.275e-05\t1890302\tspecies\t                    Bacillus wiedmannii\n+0.001297\t1\t1\t70\t1.19\t5.253e-05\t1428\tspecies\t                    Bacillus thuringiensis\n+0.01816\t14\t0\t1824\t1.22\t5.778e-05\t185979\tno rank\t                  unclassified Bacillus (in: Bacteria)\n+0.01167\t9\t9\t1141\t1.25\t0.002064\t2597270\tspecies\t                    Bacillus sp. P16(2019)\n+0.006484\t5\t5\t681\t1.19\t4.507e-05\t1409\tspecies\t                    Bacillus sp. (in: Bacteria)\n+0.01038\t8\t6\t3457\t1.41\t0.0003944\t653685\tspecies group\t                  Bacillus subtilis group\n+0.002594\t2\t2\t1721\t1.15\t0.002671\t483913\tspecies\t                    Bacillus inaquosorum\n+0.002594\t2\t2\t119\t1.18\t0.0002753\t230954\tspecies\t                  Bacillus thermozeamaize\n+0.001297\t1\t1\t179\t1.23\t0.0004264\t408580\tspecies\t                  Bacillus coahuilensis\n+0.001297\t1\t1\t61\t1.36\t0.000156\t2039284\tspecies\t                  Bacillus fungorum\n+0.001297\t1\t0\t42\t1.17\t4.106e-06\t84406\tgenus\t                Virgibacillus\n+0.001297\t1\t1\t37\t1.19\t6.655e-05\t1638645\tspecies\t                  Virgibacillus kapii\n+0.001297\t1\t0\t5\t1\t1.141e-06\t704093\tgenus\t                Anaerobacillus\n+0.001297\t1\t0\t2\t1\t8.059e-07\t2625356\tno rank\t                  unclassified Anaerobacillus\n+0.001297\t1\t1\t2\t1\t8.059e-07\t1872506\tspecies\t                    Anaerobacillus sp.\n+3.84\t2961\t27\t58221\t1.28\t0.001726\t90964\tfamily\t              Staphylococcaceae\n+3.803\t2932\t321\t56780\t1.27\t0.002658\t1279\tgenus\t                Staphylococcus\n+1.935\t1492\t1492\t25704\t1.26\t0.07286\t1280\tspecies\t                  Staphylococcus aureus\n+1.414\t1090\t1090\t20177\t1.19\t0.06293\t1282\tspecies\t                  Staphylococcus epidermidis\n+0.01038\t8\t8\t231\t1.31\t0.0002607\t1283\tspecies\t                  Staphylococcus haemolyticus\n+0.007781\t6\t6\t92\t1.54\t0.0002367\t29388\tspecies\t                  Staphylococcus capitis\n+0.005188\t4\t4\t121\t1.34\t0.0001435\t246432\tspecies\t                  Staphylococcus equorum\n+0.003891\t3\t3\t197\t1.42\t0.'..b'illevirus\n+0.001297\t1\t0\t8\t1\t3.411e-05\t2315196\tno rank\t                    unclassified Bastillevirus\n+0.001297\t1\t1\t8\t1\t0.0009028\t2282408\tspecies\t                      Bacillus phage Maceta\n+0.01167\t9\t0\t944\t1.61\t2.861e-05\t2731361\tphylum\t        Peploviricota\n+0.01167\t9\t0\t944\t1.61\t2.861e-05\t2731363\tclass\t          Herviviricetes\n+0.01167\t9\t0\t944\t1.61\t2.861e-05\t548681\torder\t            Herpesvirales\n+0.01167\t9\t0\t944\t1.61\t3.27e-05\t10292\tfamily\t              Herpesviridae\n+0.01167\t9\t0\t935\t1.62\t8.081e-05\t10357\tsubfamily\t                Betaherpesvirinae\n+0.01167\t9\t0\t935\t1.62\t0.0001562\t10358\tgenus\t                  Cytomegalovirus\n+0.01038\t8\t8\t783\t1.65\t0.0003119\t10359\tspecies\t                    Human betaherpesvirus 5\n+0.001297\t1\t1\t152\t1.47\t0.0003854\t47929\tspecies\t                    Macacine betaherpesvirus 3\n+0.1945\t150\t53\t13859\t1.21\t0.0002668\t12333\tno rank\t    unclassified bacterial viruses\n+0.05058\t39\t39\t3891\t1.28\t0.186\t2716355\tspecies\t      Escherichia phage 503046\n+0.02334\t18\t18\t1490\t1.29\t0.05535\t2783517\tspecies\t      Escherichia phage 204576\n+0.02334\t18\t18\t3311\t1.13\t0.1106\t129862\tspecies\t      Phage Gifsy-2\n+0.02205\t17\t17\t1867\t1.24\t0.05953\t129861\tspecies\t      Phage Gifsy-1\n+0.002594\t2\t2\t269\t1.07\t0.01331\t2859565\tspecies\t      Enterococcus phage VEsP-2\n+0.001297\t1\t1\t22\t1\t0.01401\t356350\tspecies\t      Phage 258-320\n+0.001297\t1\t1\t191\t1\t0.1322\t2558708\tspecies\t      Streptococcus satellite phage Javan449\n+0.001297\t1\t1\t67\t1\t0.001649\t128975\tspecies\t      Salmonella phage Fels-1\n+0.006484\t5\t0\t444\t1.01\t1.139e-05\t2731342\tclade\t    Monodnaviria\n+0.006484\t5\t0\t431\t1\t0.0004848\t2732090\tkingdom\t      Loebvirae\n+0.006484\t5\t0\t431\t1\t0.0004848\t2732410\tphylum\t        Hofneiviricota\n+0.006484\t5\t0\t431\t1\t0.0004848\t2732411\tclass\t          Faserviricetes\n+0.006484\t5\t0\t431\t1\t0.0004848\t2732094\torder\t            Tubulavirales\n+0.006484\t5\t3\t431\t1\t0.0005104\t10860\tfamily\t              Inoviridae\n+0.001297\t1\t0\t83\t1\t0.0001515\t456491\tno rank\t                unclassified Inoviridae\n+0.001297\t1\t1\t83\t1\t0.01283\t2686285\tspecies\t                  Pseudomonas phage pf8_ST274-AUS411\n+0.001297\t1\t0\t77\t1\t0.00601\t2732875\tgenus\t                Primolicivirus\n+0.001297\t1\t1\t77\t1\t0.00601\t2011081\tspecies\t                  Pseudomonas virus Pf1\n+0.005188\t4\t0\t210\t2.23\t4.643e-06\t12429\tno rank\t    unclassified viruses\n+0.005188\t4\t0\t210\t2.23\t5.953e-06\t2204151\tno rank\t      unclassified DNA viruses\n+0.005188\t4\t0\t210\t2.23\t0.0009295\t79205\tno rank\t        unclassified dsDNA phages\n+0.005188\t4\t4\t210\t2.23\t0.01276\t320831\tspecies\t          Staphylococcus phage PT1028\n+0.003891\t3\t0\t150\t1.22\t5.186e-07\t2559587\tclade\t    Riboviria\n+0.003891\t3\t0\t144\t1.23\t6.207e-07\t2732396\tkingdom\t      Orthornavirae\n+0.002594\t2\t0\t100\t1.12\t4.316e-06\t2732405\tphylum\t        Duplornaviricota\n+0.002594\t2\t0\t83\t1.1\t2.868e-05\t2732458\tclass\t          Chrymotiviricetes\n+0.002594\t2\t0\t83\t1.1\t2.868e-05\t2732540\torder\t            Ghabrivirales\n+0.002594\t2\t0\t83\t1.1\t3.981e-05\t11006\tfamily\t              Totiviridae\n+0.002594\t2\t0\t83\t1.1\t0.0002877\t11007\tgenus\t                Totivirus\n+0.002594\t2\t2\t83\t1.1\t0.002982\t11008\tspecies\t                  Saccharomyces cerevisiae virus L-A\n+0.001297\t1\t0\t24\t1.71\t2.937e-07\t2732408\tphylum\t        Pisuviricota\n+0.001297\t1\t0\t20\t1.85\t3.197e-07\t2732506\tclass\t          Pisoniviricetes\n+0.001297\t1\t0\t20\t1.85\t5.285e-07\t464095\torder\t            Picornavirales\n+0.001297\t1\t0\t20\t1.85\t9.652e-07\t12058\tfamily\t              Picornaviridae\n+0.001297\t1\t0\t20\t1.85\t1.591e-06\t12059\tgenus\t                Enterovirus\n+0.001297\t1\t1\t20\t1.85\t1.107e-05\t138950\tspecies\t                  Enterovirus C\n+0.001297\t1\t0\t63\t1.35\t3.295e-06\t186616\tno rank\t    environmental samples\n+0.001297\t1\t1\t16\t1.69\t1.073e-05\t340016\tspecies\t      uncultured virus\n+0.6575\t507\t0\t45700\t1.63\t0.006306\t2787854\tno rank\t  other entries\n+0.6575\t507\t0\t45700\t1.63\t0.006306\t28384\tno rank\t    other sequences\n+0.6575\t507\t0\t45700\t1.63\t0.006306\t81077\tno rank\t      artificial sequences\n+0.6575\t507\t507\t45700\t1.63\t0.006306\t32630\tspecies\t        synthetic construct\n'
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/krona.k2.int.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/krona.k2.int.tsv Wed Mar 15 18:18:07 2023 +0000
b
b'@@ -0,0 +1,73 @@\n+24\tUnclassified\n+0\tx__root\n+0\tx__root\tk__Viruses\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\tc__Herviviricetes\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\tc__Herviviricetes\to__Herpesvirales\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\tc__Herviviricetes\to__Herpesvirales\tf__Herpesviridae\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\tc__Herviviricetes\to__Herpesvirales\tf__Herpesviridae\tx__Alphaherpesvirinae\n+0\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\tc__Herviviricetes\to__Herpesvirales\tf__Herpesviridae\tx__Alphaherpesvirinae\tg__Iltovirus\n+497\tx__root\tk__Viruses\tx__Duplodnaviria\tx__Heunggongvirae\tp__Peploviricota\tc__Herviviricetes\to__Herpesvirales\tf__Herpesviridae\tx__Alphaherpesvirinae\tg__Iltovirus\ts__Gallid_alphaherpesvirus_1\n+0\tx__root\tk__Viruses\tx__Riboviria\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\tf__Paramyxoviridae\n+12\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\tf__Paramyxoviridae\tg__Avulavirus\n+0\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\tf__Paramyxoviridae\tg__Avulavirus\tx__unclassified_Avulavirus\n+15\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\tf__Paramyxoviridae\tg__Avulavirus\tx__unclassified_Avulavirus\ts__Goose_paramyxovirus_SF02\n+7\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\tf__Paramyxoviridae\tg__Avulavirus\ts__Avian_avulavirus_1\n+1\tx__root\tk__Viruses\tx__Riboviria\tx__Orthornavirae\tp__Negarnaviricota\tx__Haploviricotina\tc__Monjiviricetes\to__Mononegavirales\tf__Paramyxoviridae\tg__Avulavirus\ts__Avian_avulavirus_12\n+0\tx__root\tx__cellular_organisms\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\tc__Gammaproteobacteria\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\tc__Gammaproteobacteria\to__Pseudomonadales\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\tc__Gammaproteobacteria\to__Pseudomonadales\tf__Moraxellaceae\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\tc__Gammaproteobacteria\to__Pseudomonadales\tf__Moraxellaceae\tg__Acinetobacter\n+0\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\tc__Gammaproteobacteria\to__Pseudomonadales\tf__Moraxellaceae\tg__Acinetobacter\tx__Acinetobacter_calcoaceticus/baumannii_complex\n+2\tx__root\tx__cellular_organisms\tk__Bacteria\tp__Proteobacteria\tc__Gammaproteobacteria\to__Pseudomonadales\tf__Moraxellaceae\tg__Acinetobacter\tx__Acinetobacter_calcoaceticus/baumannii_complex\ts__Acinetobacter_baumannii\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\n+0\tx__root\tx__'..b'x__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\tx__Deuterostomia\tp__Chordata\tx__Craniata\tx__Vertebrata\tx__Gnathostomata\tx__Teleostomi\tx__Euteleostomi\tx__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\tx__Catarrhini\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\tx__Deuterostomia\tp__Chordata\tx__Craniata\tx__Vertebrata\tx__Gnathostomata\tx__Teleostomi\tx__Euteleostomi\tx__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\tx__Catarrhini\tx__Hominoidea\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\tx__Deuterostomia\tp__Chordata\tx__Craniata\tx__Vertebrata\tx__Gnathostomata\tx__Teleostomi\tx__Euteleostomi\tx__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\tx__Catarrhini\tx__Hominoidea\tf__Hominidae\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\tx__Deuterostomia\tp__Chordata\tx__Craniata\tx__Vertebrata\tx__Gnathostomata\tx__Teleostomi\tx__Euteleostomi\tx__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\tx__Catarrhini\tx__Hominoidea\tf__Hominidae\tx__Homininae\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\tx__Deuterostomia\tp__Chordata\tx__Craniata\tx__Vertebrata\tx__Gnathostomata\tx__Teleostomi\tx__Euteleostomi\tx__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\tx__Catarrhini\tx__Hominoidea\tf__Hominidae\tx__Homininae\tg__Homo\n+1\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Opisthokonta\tx__Metazoa\tx__Eumetazoa\tx__Bilateria\tx__Deuterostomia\tp__Chordata\tx__Craniata\tx__Vertebrata\tx__Gnathostomata\tx__Teleostomi\tx__Euteleostomi\tx__Sarcopterygii\tx__Dipnotetrapodomorpha\tx__Tetrapoda\tx__Amniota\tc__Mammalia\tx__Theria\tx__Eutheria\tx__Boreoeutheria\tx__Euarchontoglires\to__Primates\tx__Haplorrhini\tx__Simiiformes\tx__Catarrhini\tx__Hominoidea\tf__Hominidae\tx__Homininae\tg__Homo\ts__Homo_sapiens\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\to__Haemosporida\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\to__Haemosporida\tf__Plasmodiidae\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\to__Haemosporida\tf__Plasmodiidae\tg__Plasmodium\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\to__Haemosporida\tf__Plasmodiidae\tg__Plasmodium\tx__Plasmodium_(Plasmodium)\n+0\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\to__Haemosporida\tf__Plasmodiidae\tg__Plasmodium\tx__Plasmodium_(Plasmodium)\ts__Plasmodium_knowlesi\n+1\tx__root\tx__cellular_organisms\tk__Eukaryota\tx__Sar\tx__Alveolata\tp__Apicomplexa\tc__Aconoidasida\to__Haemosporida\tf__Plasmodiidae\tg__Plasmodium\tx__Plasmodium_(Plasmodium)\ts__Plasmodium_knowlesi\tx__Plasmodium_knowlesi_strain_H\n'
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/krona.k2.noint.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/krona.k2.noint.tsv Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,36 @@
+24 Unclassified
+0 k__Viruses
+0 k__Viruses p__Peploviricota
+0 k__Viruses p__Peploviricota c__Herviviricetes
+0 k__Viruses p__Peploviricota c__Herviviricetes o__Herpesvirales
+0 k__Viruses p__Peploviricota c__Herviviricetes o__Herpesvirales f__Herpesviridae
+0 k__Viruses p__Peploviricota c__Herviviricetes o__Herpesvirales f__Herpesviridae g__Iltovirus
+497 k__Viruses p__Peploviricota c__Herviviricetes o__Herpesvirales f__Herpesviridae g__Iltovirus s__Gallid_alphaherpesvirus_1
+0 k__Viruses p__Negarnaviricota
+0 k__Viruses p__Negarnaviricota c__Monjiviricetes
+0 k__Viruses p__Negarnaviricota c__Monjiviricetes o__Mononegavirales
+0 k__Viruses p__Negarnaviricota c__Monjiviricetes o__Mononegavirales f__Paramyxoviridae
+12 k__Viruses p__Negarnaviricota c__Monjiviricetes o__Mononegavirales f__Paramyxoviridae g__Avulavirus
+15 k__Viruses p__Negarnaviricota c__Monjiviricetes o__Mononegavirales f__Paramyxoviridae g__Avulavirus s__Goose_paramyxovirus_SF02
+7 k__Viruses p__Negarnaviricota c__Monjiviricetes o__Mononegavirales f__Paramyxoviridae g__Avulavirus s__Avian_avulavirus_1
+1 k__Viruses p__Negarnaviricota c__Monjiviricetes o__Mononegavirales f__Paramyxoviridae g__Avulavirus s__Avian_avulavirus_12
+0 k__Bacteria
+0 k__Bacteria p__Proteobacteria
+0 k__Bacteria p__Proteobacteria c__Gammaproteobacteria
+0 k__Bacteria p__Proteobacteria c__Gammaproteobacteria o__Pseudomonadales
+0 k__Bacteria p__Proteobacteria c__Gammaproteobacteria o__Pseudomonadales f__Moraxellaceae
+0 k__Bacteria p__Proteobacteria c__Gammaproteobacteria o__Pseudomonadales f__Moraxellaceae g__Acinetobacter
+2 k__Bacteria p__Proteobacteria c__Gammaproteobacteria o__Pseudomonadales f__Moraxellaceae g__Acinetobacter s__Acinetobacter_baumannii
+0 k__Eukaryota
+0 k__Eukaryota p__Chordata
+0 k__Eukaryota p__Chordata c__Mammalia
+0 k__Eukaryota p__Chordata c__Mammalia o__Primates
+0 k__Eukaryota p__Chordata c__Mammalia o__Primates f__Hominidae
+0 k__Eukaryota p__Chordata c__Mammalia o__Primates f__Hominidae g__Homo
+1 k__Eukaryota p__Chordata c__Mammalia o__Primates f__Hominidae g__Homo s__Homo_sapiens
+0 k__Eukaryota p__Apicomplexa
+0 k__Eukaryota p__Apicomplexa c__Aconoidasida
+0 k__Eukaryota p__Apicomplexa c__Aconoidasida o__Haemosporida
+0 k__Eukaryota p__Apicomplexa c__Aconoidasida o__Haemosporida f__Plasmodiidae
+0 k__Eukaryota p__Apicomplexa c__Aconoidasida o__Haemosporida f__Plasmodiidae g__Plasmodium
+1 k__Eukaryota p__Apicomplexa c__Aconoidasida o__Haemosporida f__Plasmodiidae g__Plasmodium s__Plasmodium_knowlesi
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k1.11176.children.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k1.11176.children.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,20 @@
+>C_II_XVIIIb-87/1
+TTTACCAAGAGTCCCAAATCATGATATATATATTGTCAGTGCTAGGGGGGGCATTGAGGG
+ATTATGCCAGAAGCTATGGACAATGATCTCAATTGCTACAATCCAACTTGCTGCAGCAAG
+ATCACACTGTCGTGTTGCCTGCATGGTACA
+>C_II_XVIIIb-53/1
+CACTCGGATACCCTCATTTGATATGAGTACTACTCACTACTGTTATACTCACAATGTTAT
+ACTATCTGGTTGTAGAGATCACTCACACTCACACCAATACTTAGCACTTGGTGTGCTTCG
+GGTATCTGCGACAGGGAGGGTATTCTTTTC
+>C_II_XVIIIb-37/1
+GGCTGCATTGTTGGCACATGACATAACAGTGTTTTCACTAAGGTCACCTGATATTAGCAC
+TAATTTAGATGAGTTTTTGAGGACCTGACTGAGTATTGCTCCGTCTTTGAATATTCGTTT
+GCTATATATGAAGAATGTGTCTGACCTGAT
+>C_II_XVIIIb-25/1
+ACAAATTTCGGAAGATTGAAAGGAAGATTCAGAGTCACAATACAAGGTATGGAGAATTGT
+TCACAAAGCTGTGCACACATATTGAAAAGAAATTGCTAGGATCATCTCGATCTAATAATG
+TCCCACGGTCAGAAGAATTCAGTAGCATCC
+>C_II_XVIIIb-23/1
+TCACACTGCTGGAGTTTCGATTTTGGGTGTCATCTAAATTGATAGAACGCAGAGTAGAAA
+AGAATACCCTCCCTGTCGCAGATACCCGAAGCACACCAAGTGCTAAGTATTGGTGTGAGT
+GTGAGTGATCTCTACAACCAGATAGTATAA
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k1.e10386.children.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k1.e10386.children.fa Wed Mar 15 18:18:07 2023 +0000
b
b'@@ -0,0 +1,252 @@\n+>C_II_XVIIIb-101/1\n+AGTGACATAGTAATGGGGGAGAGGAGGCGCATTACAGTCTTAGGAACGCAACTTACTTCC\n+TGAAGGGATTGTACTTAGCAATGGTATGTCTCTTCTCTATCTTAGTAGAGGTGACTTCAT\n+GATCAGCGGTCACTGCGACAGCGCGCTGAG\n+>C_II_XVIIIb-99/1\n+ATCCAATCACCCATTAGTCTGCTCTATTAAGTGTTCGCTGACGCTAGCAGACTATGCCCG\n+GAACAGAAGCTGGTCACCTTTGACAGGGGGTAGGAAGATACTAGGCGTATCTAATCCTGA\n+CACCATAGAACTTGTGGAAGGAGAGATTCT\n+>C_II_XVIIIb-97/1\n+TGAGTCATGCATAGAGTACGACCCCATTACCAATCTAAGCATGTTCTTAAAGGACAAGGC\n+AATCGCACATCCGAAAGATAACTGGCTCGCCTCATTCAGGCGGAACCTCCTCTCTGAGGA\n+CCAAAAGAAACAGATAAAAGAGGCGACCTC\n+>C_II_XVIIIb-95/1\n+GTCATATTCCACACTGGAGCTCATGAAAATTCTTTCAATCTCTAGCGGAAAGTTGATCGG\n+CCAATCAGTGGTTTCTTATGATGAAGATACCTCTATAAAGAATGATGCTATAATAGTGTA\n+TGACAACACACGAAATTGGATCAGTGAGGC\n+>C_II_XVIIIb-93/1\n+GTCGCAAGTTGATCTGATAGCCTTTGGTTCGACCCCGGAGGAATCTCAATGTCGCAATGC\n+AGTAATGCTACAGACCTGTGGGGCACTACAGATGTGATATATCCCACTGCTCTATCTGAA\n+GTCAGGTCGCTTTCTTCTGCATTCTCTCTC\n+>C_II_XVIIIb-91/1\n+CTTCGCCTGCCTTCCTGCATTACTATCTAACTTGGTAATAGTCTATCTTGATTGCTTATA\n+GCTAGTTCACCTGTCTATCTAATTAGAAAAAACACGGGTAGAAGAGTCTGGAACCTAACT\n+AGCACATTCTGAACGCAATATGGGCTCCAA\n+>C_II_XVIIIb-89/1\n+GGCTATTGCATGTCGGCCTGATATGGCAGTGGAGAAGGACACTGTCCGTGCATTAATCAC\n+CTCTCGCCCAATGCATCCAAGCTCCTCAGCTAAGCTCCTGAGCAAGCTGGATGCAGCCAG\n+GTCGATTGAAGAGATCAGGAAGATCAAACG\n+>C_II_XVIIIb-87/1\n+TTTACCAAGAGTCCCAAATCATGATATATATATTGTCAGTGCTAGGGGGGGCATTGAGGG\n+ATTATGCCAGAAGCTATGGACAATGATCTCAATTGCTACAATCCAACTTGCTGCAGCAAG\n+ATCACACTGTCGTGTTGCCTGCATGGTACA\n+>C_II_XVIIIb-85/1\n+TAGGTTGAACAAGCTGCAAGGGCACCCAGACGTCGGGGATGGAGAGACCCAATTCCTAGA\n+TTTTATGAGAGCAGTGGCGAACAGCGTGCGAGAAGCGTCAAATCCCGCACAAAACGTCAC\n+CCCTCTAGAGCCTCCCCCAACTCCTGGGCC\n+>C_II_XVIIIb-83/1\n+AGGTAGTTTCCCTTGGCAAGAATCAGCAAGATGCCCTAACTCATGCGTCACTGGAGTCTA\n+TACTGATCCGTATCCCTTAGTCTTCAATCGGAATCATACCCTACGAGGGGTATTCGGGAC\n+AATGCTTGATGATGGACAGGCAAGACTCAA\n+>C_II_XVIIIb-81/1\n+GGCCGCCGGCCCACTGAGTGGGGCCGCCTTCTCGGGTGCGGACGGCTGTTTGTCTGGTCT\n+GGCTTGTTGGTCGGGGGCGTCCTGGCCGGCATTTGATTGGACACTCCCGTGCTTCTCCCA\n+TGCAGTGCTCAGAGCCTTGGTCTTGCCTTG\n+>C_II_XVIIIb-79/1\n+GCTTGTAAAGCATGGACAGAAGATAGTTGCTTGCTCCGGTATTGAGCTGTATGTCGCCGG\n+CCTCGCCTGCAGCCTGGGTGGGAGAGGACTCGGTGGATGTGGCCGGCGGGTCACTGTGTT\n+GGGCCGCCTGCTCGGGTGCGGACGGCTGTT\n+>C_II_XVIIIb-77/1\n+TGCTAGGAAGATGAAAGCACGTAACCTTCTCATCTCCACTGTCACACTTTGTGCACCCTC\n+CACTGACACTGAGAATCTCTCCTTCCACAAGTTCTATGGTGTCAGGATTAGATACGCCTA\n+GTATCTTCCTACCCCCTGTCAAAGGTGACC\n+>C_II_XVIIIb-75/1\n+CAATCTGGATGATAACGGACTAGCGCCCGGAGATGACTTAAGGATGAAATGTTGGCGCAA\n+CCGGGGTCCAGAATTTTCATCATGCCTAAGTTGGCTTCCATGATCGCAACAGATGTTTTA\n+AGCTGTTGGATTTCAGACCGCATCATAGGG\n+>C_II_XVIIIb-73/1\n+CAGACGCGGAATGCTCAGCAGTCACTCAATACCTATTGTCGGATGCTGTGAGACCATTGC\n+TTAGGTCGGAACAAGTGAACTCTATCATATCTCCCAACATAATCACATTTCCCGCCAATC\n+TATATTACATGTCCAGGAAGAGTCTTAATC\n+>C_II_XVIIIb-71/1\n+GGGTTAACTCTGAGGCTCAGGGGAGAATGTGATGCAACTTATCAAAAGGATATCTCAATA\n+TTAGATTCTCAAGTCATCGTGACCGGCAATCTCGTTATATCAACTGAGCTTGGAAATGTC\n+AACAATTCAATAAGCAATGCCTTGGATAAG\n+>C_II_XVIIIb-69/1\n+CTTAAAGCTAGTTCACCTGTCTATCTAATTAGAAAAAACACGGGTAGAAGAGTCTGGAAC\n+CTAACTAGCACATTCTGAACGCAATATGGGCTCCAAACCTTCTACCAGGATTCCAGTACC\n+TCTGATGCTGATCATTCGGACTATGCTGAT\n+>C_II_XVIIIb-67/1\n+CTATCTCGGCATCTGTAGCGGTGGCCATCTTCACACTAATCTACTAGGTAGAAGAGAGAA\n+CAGAGACTTTGGTGACTCGTCCTGGTCTCCAGGTGTCTCCCTTCTACCCGTATTTTTTTC\n+TAATACCTCGGGCTCTGAGTGGATTGTTGG\n+>C_II_XVIIIb-65/1\n+ACCCCCACCTGTTGAGTAGGAATATCCATGCAGCCAATTTCTTCAGACACTCGTTGGGCA\n+GCAGCCGCCAGACCTCTCCTTGCTGCCGGAGTCAGTTTTAGTTCAGCAGCCATGTCTTCA\n+TTGATGCTACTTCCCTGAGCCTGAGCATAC\n+>C_II_XVIIIb-63/1\n+GTATGATCCTAGGCCTACGTCAGAAAAGGATCGTGCAAGATTCGACTTAGCTATCTTCAA\n+GAGTTATGAGCTTAATTTGGAGTCATATTCCACACTGGAGCTGATGAACATTCTTTCAAT\n+CTCTAGCGGAAAGTTGATCGGCCAATCAGT\n+>C_II_XVIIIb-61/1\n+GCTGATGAACATTTTGACAATCCCTATCGGAACGTTGATCGGCCAATCAGTGGTTTCTTA\n+TGATGAAGATCCCTCTATAAAGAATGATGCTATAATAGTGTATGACAACACACGAAATTG\n+GATCAGTGAGGCCCAGAATTCAGATGTGGT\n+>C_II_XVIIIb-59/1\n+TACCACAAGAGTTCATCTGCCCTCGAATTTGAGACATGCATAGAGTATGACCCCGTTACC\n+AATCTAAGCATGTTCTTAAAGGACAAGGCAATCGCACATCCGAAAGATAACTGGCTCGCC\n+TCATTCAGGCGGAACCTCCTCTCTGAGGAC\n+>C_II_XVIIIb-57/1\n+CGGGAGCCCAGTTAATTTCCAGTAATAGAGCAACTTGTGCTTGACCAATGGAGAGGATAG\n+ATGTGACTCTGGTAGGATAATCTGGTGCTCTGCCCTTTCGGGACCAGA'..b'TATAA\n+>C_II_XVIIIb-21/1\n+TTCTCTGAACGTTCGACCCGTAGATTGCTGCTTAGTTTCTTATACCCGTCAGTTTTTTCT\n+TGAGTTTTCGATCTGACAGGCTGTCAGAATTGACATACAGGCTGCTGCTGTATGTGCATC\n+CAATTCTCATTTGCGGTCATGCTCTTGTAT\n+>C_II_XVIIIb-19/1\n+GGCACAGACGCGAACAGAGCTTGTCATGGACTATGGAAGGAGTCACAACCATCAGCTGGT\n+GCAACCCCTCATGTGCCCCAGTCAGGGCAGAACCAAGACAATACTCCTGCACCTGTGGAT\n+CATGTCCAGCTACCTGTCGACTTTGTGCAA\n+>C_II_XVIIIb-17/1\n+GAATTTGAATCTGTGAGGCCCAGAATTCAGATGTGGTCCGCTTATTCGAATATGCGGCAC\n+TCGAGGTGCTCCTTGACTGTTCTTACCAACTCTACTATCTGAGAGTGAGAGGCCTAAACA\n+ACATCGTCTTATACATGAGTGACCTATACA\n+>C_II_XVIIIb-15/1\n+TGAGCCGCGATCTTAAAAGCCGAGGATAAGTTCACGACTTTGCTCAGTTCTAGCCAAGAT\n+GTGACAGATAGGATATATAAACAGGTGGCTCTTGAATCTCCGCTAGCGCTACTAAACACC\n+GAATCTATAATTATGAATGCAATAACATCC\n+>C_II_XVIIIb-13/1\n+CCTCACAGTAGCGTCTTGTGAACGACATTCTATCCAGTCCCCTACCGAGGCTGATGAAGT\n+TCTTGGGAGAAAATATTGATGCTGCACTGATCGCAGCAGGGGGACAGCCCGTCCGTCCCT\n+TTTGTGCAGAAAGTTTGGTGAGCACACTAA\n+>C_II_XVIIIb-11/1\n+TGTGCACGGTTTAATTGACGTCATCTTTTTCCCGTCTGTAGAGAGGTTATAAGGAGTAAA\n+TAAGAATACTGTGTCAGCAAGGTCCCCCTCGGCATCCATGTAGATCACAGATCTAATGAC\n+TGTGTCAATGTGGCTGGCAATTATCTGAGT\n+>C_II_XVIIIb-9/1\n+GACTAATCATGTCAGTGTCGTTAGACTGTGCCTACCCGGCAGCACGGGATCTGGGCTTGT\n+GTTGATATTCAGTTTGACTCTTGGCATCGCGAAAGATATTCAGTGCCACATCATTCTTCC\n+AGCCGTCCTATGGCTGACTCAATTGGCTGG\n+>C_II_XVIIIb-7/1\n+GCCATGGTGATTAAGCAGTTCCTAATCCTCTGAACAAGTATCGTACTAAGTAATCTTCCT\n+TTGGGTTTGGTAGTGTGTGCTCGGGGCAAATCTCGCTAAGCGTAAATGCATCATACCTTG\n+CTGGAGCAGTCAAATCTAGCTCTTGTATAA\n+>C_II_XVIIIb-5/1\n+TAGAAACTATTATAGCATCATTGTTTTTATAGGTATCTTCATCATAAGAAACCACTGATT\n+GGCCGATCAACTTTCCGCTAGAGATTGAAAGAATGTTCATCAGCTCCAGTGTGGAATATG\n+ACTCCAAATTAAGCTCATAACTCTTGAAGA\n+>C_II_XVIIIb-3/1\n+CAGATAAAAGAGACGACCCTAACTTATCGCCTCTTGATTGAGTTTTTAGAATCAAATGAT\n+TTTGATCCATACAAGGAGATGGAATACTTGACAACCCTTGAGTACCTAAGAGATGACAGT\n+GTGGCAGTATCATATTCACTTAAAGAGAAG\n+>C_II_XVIIIb-1/1\n+CTTGCTCATTGGCTCATGATTGCGGTTTGTGGAAACTCTCTCTTTGCAGTCGGTGATACG\n+TTTCTTATTGCTGTTGAAGGACAGTTGACTCATTGCCAACATACTCTTAGTCAAGGATAT\n+GCTATCTTGAATGACCCCATTTCCTTGGAA\n+>HQ630064-957/1\n+CAAAGCAAATTTTGTAGCTTAAGATAAAGGGGGCCCTAAAGAAGAAGCAACATCGTCTCG\n+GGCTGTGTTTCCGACGTTGGGGTGCATTATGCGATTTTACTACCAGTGTTTGGTGTGGTG\n+TGGTAAGGGGAAATGAGACGAGAATTTTTC\n+>HQ630064-665/1\n+CCCCGGGAGCTTTTTAGCGCGGCGTCACGTATAAATGCACAGAAGCTTCGTCAGTAACGA\n+AGCGCCTCTGCCAGAGTTATTAACATTAGGACTCGAGAGAGTCTCGAGAAGAGACTTCCC\n+AAAGAGAGTGCCGCTTTGATCTCCAGATCA\n+>HQ630064-403/1\n+TGTCTACCTGTCTGTCAAATTGTGACTACCACGTACCAACTCCTCCCCAAAGTTCTACCC\n+TTTTCCTTAAAAACACATTTAAAATGTCAAAATACCGATCCTTATATTTGTGTGTTATGT\n+TTTCTTTTGCCTTCTAGGATTGATTGAATT\n+>HQ630064-343/1\n+ACGAAAGAAAGAAAGGAAAAGGCCTCAAAGTTAATTTTGTAGCTTAAAAAAGAGGGGGCC\n+CTAAAGAAGAAGCAACATCGTCTCGGGCTGTGTTTCCGACGTTGGGGTGCATTATGCGAT\n+TTTACTACCAGTGTTTGGTGTGGTGTGGTA\n+>HQ630064-307/1\n+TCGTGTGTTTGCTCGCGTGCTTTGGATGTTGGGATGAGGCCTTGGTGTTGCCAGGCTCCG\n+GTTTCTTAAAGGTTATCCGAGGGGGAAGTTCCACGTGTTCTTTTAGCTTTTCTCCATATT\n+CCGGTCTACGGGTACAAATGTCCTCTGGAT\n+>HQ630064-241/1\n+CGTGACCTCATTTACAAGATGTCGGTGGCACCCAGCTACGTCACGTTCGCACCGAGGGTC\n+CAACATTTTCCTAAAATGGCAAAGCGGGACTTGGGTCCGTGGGTAAATTTTCTAATTTTT\n+TCCCACGGTCTGTCGCAAGAGCGTATTCTC\n+>HQ630064-217/1\n+AACTTTGGGGAGGAATTGGTACGTGGTAGACACAATTTGACAGACAGGTAGACAGACAAC\n+GAATTACAAGCGGCCTAGGGTAATAGCGCTGAGCCCATGATGTGTCCGTGCATATGAGAC\n+GACTAATAAGGTAGTCACGTAAACCCCCGT\n+>HQ630064-199/1\n+AGAAGCAGTATTGTGATCTTTCTCTCTTTCGCATGAATCGCTAGGAAGAGGAATTTTGGA\n+AGCATCTGGATACCTTACAGGAGATTTAGAAATCTCGGTAGTTATGACATGTTCATCAGA\n+TGAGGAAGATAGCAAGCGTGACGGTGTTTT\n+>HQ630064-111/1\n+GTGTCTCCGAGTTGTGCTCTCTCGAAACTCACTTTGCCCTCACTTGTTGGCTTTATGGAG\n+GTCGTAATGATCATCGATGGAGTACTTATATAATATTGTGCTCTGCCGGCCATCTCTCCG\n+TGGTCGCAGAACTCCAACACCTTGTTGGCG\n+>HQ630064-105/1\n+AGTTTAAACTTTGTCCGCGCGCCAATCTTGAAGTGATATTAGGACATAGTCAACCAACCG\n+ACGTACGTTGAACCATTTGCGGGCAGCGCTGAGGGAAGACGAGATGTATGGACAACACAA\n+CGACCGGAGAAGAGGGCGTAAAGAACATCA\n+>HQ630064-11/1\n+GAAAGAAAGAAAGGAAAAGGCCTCAAAGCAAATTTTGTAGCTTAAAAAAGAGGGGGCCCT\n+AAAGAAGAAGCAACATCGTCTCGGGCTGTGTTTCCGACGTTGGGGTGCATTATGCGATTT\n+TACTACCAGTGTTTGGTGTGGTGTGGTAAG\n+>HQ630064-3/1\n+AGGACTGTCGCTCATATCTTGCGTACTGTGTCTTATTTACTCGTGTTGTCTTTCCGCAGG\n+GAATGCGCGGTTCGAGTAGGCCGTGGTTTGCGGCCGTTCTCGGATAAGCACAAAATATCG\n+CCCTACTCGTGTGTCAAAAAGGAAAAGTCT\n'
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k2.11176.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k2.11176.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,28 @@
+>C_II_XVIIIb-65/1
+ACCCCCACCTGTTGAGTAGGAATATCCATGCAGCCAATTTCTTCAGACACTCGTTGGGCA
+GCAGCCGCCAGACCTCTCCTTGCTGCCGGAGTCAGTTTTAGTTCAGCAGCCATGTCTTCA
+TTGATGCTACTTCCCTGAGCCTGAGCATAC
+>C_II_XVIIIb-57/1
+CGGGAGCCCAGTTAATTTCCAGTAATAGAGCAACTTGTGCTTGACCAATGGAGAGGATAG
+ATGTGACTCTGGTAGGATAATCTGGTGCTCTGCCCTTTCGGGACCAGAGCCTGCCATGTC
+CTACCCGTGCTTTCTTTTGTGTGTTGGTTT
+>C_II_XVIIIb-49/1
+ATTTAGATCCATACAATGAGATGGAATCTTTGACAACCCTTGAGTACCTAAGAGATGACA
+GTGTGGCAGTATCATATTCACTTAAAGAGAAGGAGGTGAAAGTGAATGGGCGAATTTTTG
+CTAAGTTGACAAAGAAGCTAAGGAACTGCC
+>C_II_XVIIIb-45/1
+GAGAAGATAGTTGATTGGTGTGTGATTGCCTGTAAGAGTCATTGTTACTGTAATATCCCT
+TTGCTGCATTACCTATGGTTTTCATGTAGAATTTTTGTTGAGCACGAGTTAAATAGAGTG
+AAGAGGTATCTGTAAATATTTCTTTGAGTT
+>C_II_XVIIIb-17/1
+GAATTTGAATCTGTGAGGCCCAGAATTCAGATGTGGTCCGCTTATTCGAATATGCGGCAC
+TCGAGGTGCTCCTTGACTGTTCTTACCAACTCTACTATCTGAGAGTGAGAGGCCTAAACA
+ACATCGTCTTATACATGAGTGACCTATACA
+>C_II_XVIIIb-9/1
+GACTAATCATGTCAGTGTCGTTAGACTGTGCCTACCCGGCAGCACGGGATCTGGGCTTGT
+GTTGATATTCAGTTTGACTCTTGGCATCGCGAAAGATATTCAGTGCCACATCATTCTTCC
+AGCCGTCCTATGGCTGACTCAATTGGCTGG
+>C_II_XVIIIb-1/1
+CTTGCTCATTGGCTCATGATTGCGGTTTGTGGAAACTCTCTCTTTGCAGTCGGTGATACG
+TTTCTTATTGCTGTTGAAGGACAGTTGACTCATTGCCAACATACTCTTAGTCAAGGATAT
+GCTATCTTGAATGACCCCATTTCCTTGGAA
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k2.11176.fq
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k2.11176.fq Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,28 @@
+@C_II_XVIIIb-65/1
+ACCCCCACCTGTTGAGTAGGAATATCCATGCAGCCAATTTCTTCAGACACTCGTTGGGCAGCAGCCGCCAGACCTCTCCTTGCTGCCGGAGTCAGTTTTAGTTCAGCAGCCATGTCTTCATTGATGCTACTTCCCTGAGCCTGAGCATAC
++
+G9;G*GGG3GGG*GGF+G#GFG*<G#GD,EB#GGGGFGGGGC8GCG:FGDGFC5G7AGGFGGGGFEGGGGCGDE@G>DGGGGFGGGGGG77FGGG,GGGGFGCGGFG>GGGGCFGGAGGG7GGGGGGGGGGGEGGEGGF<EFGGGF,GFG
+@C_II_XVIIIb-57/1
+CGGGAGCCCAGTTAATTTCCAGTAATAGAGCAACTTGTGCTTGACCAATGGAGAGGATAGATGTGACTCTGGTAGGATAATCTGGTGCTCTGCCCTTTCGGGACCAGAGCCTGCCATGTCCTACCCGTGCTTTCTTTTGTGTGTTGGTTT
++
+FGGGF#GFEGGGG,FE,GGGFGCGGGEG+FCGE#?G:EEGA9FGEGDGGCDFGGGGFF<GGGC:GGGGGD7GFGGGFG*A7EGGGFG@GGGFGGGGEFGGGGGGGGFGGFGGGCGFGGGGGGGGGGGGGGGGGG:GGGGGGGGGGGFGGG
+@C_II_XVIIIb-49/1
+ATTTAGATCCATACAATGAGATGGAATCTTTGACAACCCTTGAGTACCTAAGAGATGACAGTGTGGCAGTATCATATTCACTTAAAGAGAAGGAGGTGAAAGTGAATGGGCGAATTTTTGCTAAGTTGACAAAGAAGCTAAGGAACTGCC
++
+CGG9#G<:AG:GG?G@#0E#6:EGGGF##G=GGGGGG=GG<CFGGGDGCGGGG9GGGGGG7CFGGCEGFCFGFGFGGGGGD+C9=CGGGGGGGGFGGGFFGGGGCGGGG9GG<GGCGGGFGGGEEC,EGGG@GCGGCGGGGGGGGCGFGG
+@C_II_XVIIIb-45/1
+GAGAAGATAGTTGATTGGTGTGTGATTGCCTGTAAGAGTCATTGTTACTGTAATATCCCTTTGCTGCATTACCTATGGTTTTCATGTAGAATTTTTGTTGAGCACGAGTTAAATAGAGTGAAGAGGTATCTGTAAATATTTCTTTGAGTT
++
+GGGGFGFDG>:FGFF*;GCGG#?GGEGFEGC#G<G<GEGGGGFGFGGGGGFGGGFG1C+G9FGGGFGGGGFECGGGGCGC+GGGG9G+DGGFGGGGGGGGCCDFGGGCGGFG:GG4*F8EGGGGGG7GDGGG7GGGGFGGGCGG=E,=G<
+@C_II_XVIIIb-17/1
+GAATTTGAATCTGTGAGGCCCAGAATTCAGATGTGGTCCGCTTATTCGAATATGCGGCACTCGAGGTGCTCCTTGACTGTTCTTACCAACTCTACTATCTGAGAGTGAGAGGCCTAAACAACATCGTCTTATACATGAGTGACCTATACA
++
+GEG#*#F#GGG*G=EG9GG*FGE#GGGGCFGCGFE9GGCGG+GFF5EGG*GDG@EGGGEGG6FFGGGGGGGEG?GGCGGGGF<GGFGGFGGGGFACGG@GCGGGGFGG4FGGGGGGGGGGGG9C?CFGGGGGGEGGGGFGGGF:GGGGGG
+@C_II_XVIIIb-9/1
+GACTAATCATGTCAGTGTCGTTAGACTGTGCCTACCCGGCAGCACGGGATCTGGGCTTGTGTTGATATTCAGTTTGACTCTTGGCATCGCGAAAGATATTCAGTGCCACATCATTCTTCCAGCCGTCCTATGGCTGACTCAATTGGCTGG
++
+ECGGGF1CG:5GCG?)G>#GG6EFE#GG58#GGGFGGGGCGCGFBCGGGGGFGDG@FGGGGGCFGGGCFGFD8CGDFGDGGGACDGFGFGGGG:FGG>GGGG>GGGCFGG4GGCBGGGGFFGEGGGGGGGGGGGGGGGGGGGG:GGGGGG
+@C_II_XVIIIb-1/1
+CTTGCTCATTGGCTCATGATTGCGGTTTGTGGAAACTCTCTCTTTGCAGTCGGTGATACGTTTCTTATTGCTGTTGAAGGACAGTTGACTCATTGCCAACATACTCTTAGTCAAGGATATGCTATCTTGAATGACCCCATTTCCTTGGAA
++
+GFFGG@G#G#GG#GGGC;EGGGGGGE#FG#+GGGGG*GGGGDGGGFCFGG8CCGG1FGGEGGGFGGGGGF*GGGCGGFGGGGG/GF:EGGGGFF@GGG*GE+DGGGFGGCGGGGGGGG@GDGGGGG,GF,CG4FFGGGGFGGGGCFGGGG
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k2.11176.max2.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k2.11176.max2.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,8 @@
+>C_II_XVIIIb-65/1
+ACCCCCACCTGTTGAGTAGGAATATCCATGCAGCCAATTTCTTCAGACACTCGTTGGGCA
+GCAGCCGCCAGACCTCTCCTTGCTGCCGGAGTCAGTTTTAGTTCAGCAGCCATGTCTTCA
+TTGATGCTACTTCCCTGAGCCTGAGCATAC
+>C_II_XVIIIb-57/1
+CGGGAGCCCAGTTAATTTCCAGTAATAGAGCAACTTGTGCTTGACCAATGGAGAGGATAG
+ATGTGACTCTGGTAGGATAATCTGGTGCTCTGCCCTTTCGGGACCAGAGCCTGCCATGTC
+CTACCCGTGCTTTCTTTTGTGTGTTGGTTT
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k2.11176.parents.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k2.11176.parents.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,76 @@
+>C_II_XVIIIb-99/1
+ATCCAATCACCCATTAGTCTGCTCTATTAAGTGTTCGCTGACGCTAGCAGACTATGCCCG
+GAACAGAAGCTGGTCACCTTTGACAGGGGGTAGGAAGATACTAGGCGTATCTAATCCTGA
+CACCATAGAACTTGTGGAAGGAGAGATTCT
+>C_II_XVIIIb-97/1
+TGAGTCATGCATAGAGTACGACCCCATTACCAATCTAAGCATGTTCTTAAAGGACAAGGC
+AATCGCACATCCGAAAGATAACTGGCTCGCCTCATTCAGGCGGAACCTCCTCTCTGAGGA
+CCAAAAGAAACAGATAAAAGAGGCGACCTC
+>C_II_XVIIIb-93/1
+GTCGCAAGTTGATCTGATAGCCTTTGGTTCGACCCCGGAGGAATCTCAATGTCGCAATGC
+AGTAATGCTACAGACCTGTGGGGCACTACAGATGTGATATATCCCACTGCTCTATCTGAA
+GTCAGGTCGCTTTCTTCTGCATTCTCTCTC
+>C_II_XVIIIb-87/1
+TTTACCAAGAGTCCCAAATCATGATATATATATTGTCAGTGCTAGGGGGGGCATTGAGGG
+ATTATGCCAGAAGCTATGGACAATGATCTCAATTGCTACAATCCAACTTGCTGCAGCAAG
+ATCACACTGTCGTGTTGCCTGCATGGTACA
+>C_II_XVIIIb-85/1
+TAGGTTGAACAAGCTGCAAGGGCACCCAGACGTCGGGGATGGAGAGACCCAATTCCTAGA
+TTTTATGAGAGCAGTGGCGAACAGCGTGCGAGAAGCGTCAAATCCCGCACAAAACGTCAC
+CCCTCTAGAGCCTCCCCCAACTCCTGGGCC
+>C_II_XVIIIb-75/1
+CAATCTGGATGATAACGGACTAGCGCCCGGAGATGACTTAAGGATGAAATGTTGGCGCAA
+CCGGGGTCCAGAATTTTCATCATGCCTAAGTTGGCTTCCATGATCGCAACAGATGTTTTA
+AGCTGTTGGATTTCAGACCGCATCATAGGG
+>C_II_XVIIIb-65/1
+ACCCCCACCTGTTGAGTAGGAATATCCATGCAGCCAATTTCTTCAGACACTCGTTGGGCA
+GCAGCCGCCAGACCTCTCCTTGCTGCCGGAGTCAGTTTTAGTTCAGCAGCCATGTCTTCA
+TTGATGCTACTTCCCTGAGCCTGAGCATAC
+>C_II_XVIIIb-57/1
+CGGGAGCCCAGTTAATTTCCAGTAATAGAGCAACTTGTGCTTGACCAATGGAGAGGATAG
+ATGTGACTCTGGTAGGATAATCTGGTGCTCTGCCCTTTCGGGACCAGAGCCTGCCATGTC
+CTACCCGTGCTTTCTTTTGTGTGTTGGTTT
+>C_II_XVIIIb-49/1
+ATTTAGATCCATACAATGAGATGGAATCTTTGACAACCCTTGAGTACCTAAGAGATGACA
+GTGTGGCAGTATCATATTCACTTAAAGAGAAGGAGGTGAAAGTGAATGGGCGAATTTTTG
+CTAAGTTGACAAAGAAGCTAAGGAACTGCC
+>C_II_XVIIIb-45/1
+GAGAAGATAGTTGATTGGTGTGTGATTGCCTGTAAGAGTCATTGTTACTGTAATATCCCT
+TTGCTGCATTACCTATGGTTTTCATGTAGAATTTTTGTTGAGCACGAGTTAAATAGAGTG
+AAGAGGTATCTGTAAATATTTCTTTGAGTT
+>C_II_XVIIIb-43/1
+CCAGATAGTATAACATTGTGAGGATAACAGTAGTGAGTAGTACTCATATCAAATGAGGGT
+ATCCGAGTGCAACCAGACCCTGTAGTAGGGGCCGGGATAAAATTCAAGTGTTCTTGATAT
+GCGGAAGGATAAAATGATGTGACATCACTG
+>C_II_XVIIIb-33/1
+TCCGAACCCAGCATATTCTGGAATATGGAAACAGTCACAACCATCAGCTGGTGCAACCCC
+TCATGTGCCCCAGTCAGGGCAGAACCAAGACAATACTCCTGCACCTGTGGATCATGTCCA
+GCTACCTGTCGACTTTGTGCAAATGATGAT
+>C_II_XVIIIb-29/1
+TTTCCATGCTTGACAAGCTGAGCAATAATTCGTCCAATGCTAAAAAGGGCCCAAGGTCAG
+GCCCCCAAGAAGGGCACCATCAACCTCCGGTCCCACAACAAGGGAACCAGCCGAGCCATG
+GAAGCAACCAGGAAAGACCGCAGCACCAGG
+>C_II_XVIIIb-19/1
+GGCACAGACGCGAACAGAGCTTGTCATGGACTATGGAAGGAGTCACAACCATCAGCTGGT
+GCAACCCCTCATGTGCCCCAGTCAGGGCAGAACCAAGACAATACTCCTGCACCTGTGGAT
+CATGTCCAGCTACCTGTCGACTTTGTGCAA
+>C_II_XVIIIb-17/1
+GAATTTGAATCTGTGAGGCCCAGAATTCAGATGTGGTCCGCTTATTCGAATATGCGGCAC
+TCGAGGTGCTCCTTGACTGTTCTTACCAACTCTACTATCTGAGAGTGAGAGGCCTAAACA
+ACATCGTCTTATACATGAGTGACCTATACA
+>C_II_XVIIIb-9/1
+GACTAATCATGTCAGTGTCGTTAGACTGTGCCTACCCGGCAGCACGGGATCTGGGCTTGT
+GTTGATATTCAGTTTGACTCTTGGCATCGCGAAAGATATTCAGTGCCACATCATTCTTCC
+AGCCGTCCTATGGCTGACTCAATTGGCTGG
+>C_II_XVIIIb-5/1
+TAGAAACTATTATAGCATCATTGTTTTTATAGGTATCTTCATCATAAGAAACCACTGATT
+GGCCGATCAACTTTCCGCTAGAGATTGAAAGAATGTTCATCAGCTCCAGTGTGGAATATG
+ACTCCAAATTAAGCTCATAACTCTTGAAGA
+>C_II_XVIIIb-3/1
+CAGATAAAAGAGACGACCCTAACTTATCGCCTCTTGATTGAGTTTTTAGAATCAAATGAT
+TTTGATCCATACAAGGAGATGGAATACTTGACAACCCTTGAGTACCTAAGAGATGACAGT
+GTGGCAGTATCATATTCACTTAAAGAGAAG
+>C_II_XVIIIb-1/1
+CTTGCTCATTGGCTCATGATTGCGGTTTGTGGAAACTCTCTCTTTGCAGTCGGTGATACG
+TTTCTTATTGCTGTTGAAGGACAGTTGACTCATTGCCAACATACTCTTAGTCAAGGATAT
+GCTATCTTGAATGACCCCATTTCCTTGGAA
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out1.k2.exclude_both.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out1.k2.exclude_both.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,176 @@
+>C_II_XVIIIb-101/1
+AGTGACATAGTAATGGGGGAGAGGAGGCGCATTACAGTCTTAGGAACGCAACTTACTTCC
+TGAAGGGATTGTACTTAGCAATGGTATGTCTCTTCTCTATCTTAGTAGAGGTGACTTCAT
+GATCAGCGGTCACTGCGACAGCGCGCTGAG
+>C_II_XVIIIb-95/1
+GTCATATTCCACACTGGAGCTCATGAAAATTCTTTCAATCTCTAGCGGAAAGTTGATCGG
+CCAATCAGTGGTTTCTTATGATGAAGATACCTCTATAAAGAATGATGCTATAATAGTGTA
+TGACAACACACGAAATTGGATCAGTGAGGC
+>C_II_XVIIIb-91/1
+CTTCGCCTGCCTTCCTGCATTACTATCTAACTTGGTAATAGTCTATCTTGATTGCTTATA
+GCTAGTTCACCTGTCTATCTAATTAGAAAAAACACGGGTAGAAGAGTCTGGAACCTAACT
+AGCACATTCTGAACGCAATATGGGCTCCAA
+>C_II_XVIIIb-89/1
+GGCTATTGCATGTCGGCCTGATATGGCAGTGGAGAAGGACACTGTCCGTGCATTAATCAC
+CTCTCGCCCAATGCATCCAAGCTCCTCAGCTAAGCTCCTGAGCAAGCTGGATGCAGCCAG
+GTCGATTGAAGAGATCAGGAAGATCAAACG
+>C_II_XVIIIb-83/1
+AGGTAGTTTCCCTTGGCAAGAATCAGCAAGATGCCCTAACTCATGCGTCACTGGAGTCTA
+TACTGATCCGTATCCCTTAGTCTTCAATCGGAATCATACCCTACGAGGGGTATTCGGGAC
+AATGCTTGATGATGGACAGGCAAGACTCAA
+>C_II_XVIIIb-81/1
+GGCCGCCGGCCCACTGAGTGGGGCCGCCTTCTCGGGTGCGGACGGCTGTTTGTCTGGTCT
+GGCTTGTTGGTCGGGGGCGTCCTGGCCGGCATTTGATTGGACACTCCCGTGCTTCTCCCA
+TGCAGTGCTCAGAGCCTTGGTCTTGCCTTG
+>C_II_XVIIIb-79/1
+GCTTGTAAAGCATGGACAGAAGATAGTTGCTTGCTCCGGTATTGAGCTGTATGTCGCCGG
+CCTCGCCTGCAGCCTGGGTGGGAGAGGACTCGGTGGATGTGGCCGGCGGGTCACTGTGTT
+GGGCCGCCTGCTCGGGTGCGGACGGCTGTT
+>C_II_XVIIIb-77/1
+TGCTAGGAAGATGAAAGCACGTAACCTTCTCATCTCCACTGTCACACTTTGTGCACCCTC
+CACTGACACTGAGAATCTCTCCTTCCACAAGTTCTATGGTGTCAGGATTAGATACGCCTA
+GTATCTTCCTACCCCCTGTCAAAGGTGACC
+>C_II_XVIIIb-73/1
+CAGACGCGGAATGCTCAGCAGTCACTCAATACCTATTGTCGGATGCTGTGAGACCATTGC
+TTAGGTCGGAACAAGTGAACTCTATCATATCTCCCAACATAATCACATTTCCCGCCAATC
+TATATTACATGTCCAGGAAGAGTCTTAATC
+>C_II_XVIIIb-71/1
+GGGTTAACTCTGAGGCTCAGGGGAGAATGTGATGCAACTTATCAAAAGGATATCTCAATA
+TTAGATTCTCAAGTCATCGTGACCGGCAATCTCGTTATATCAACTGAGCTTGGAAATGTC
+AACAATTCAATAAGCAATGCCTTGGATAAG
+>C_II_XVIIIb-69/1
+CTTAAAGCTAGTTCACCTGTCTATCTAATTAGAAAAAACACGGGTAGAAGAGTCTGGAAC
+CTAACTAGCACATTCTGAACGCAATATGGGCTCCAAACCTTCTACCAGGATTCCAGTACC
+TCTGATGCTGATCATTCGGACTATGCTGAT
+>C_II_XVIIIb-67/1
+CTATCTCGGCATCTGTAGCGGTGGCCATCTTCACACTAATCTACTAGGTAGAAGAGAGAA
+CAGAGACTTTGGTGACTCGTCCTGGTCTCCAGGTGTCTCCCTTCTACCCGTATTTTTTTC
+TAATACCTCGGGCTCTGAGTGGATTGTTGG
+>C_II_XVIIIb-63/1
+GTATGATCCTAGGCCTACGTCAGAAAAGGATCGTGCAAGATTCGACTTAGCTATCTTCAA
+GAGTTATGAGCTTAATTTGGAGTCATATTCCACACTGGAGCTGATGAACATTCTTTCAAT
+CTCTAGCGGAAAGTTGATCGGCCAATCAGT
+>C_II_XVIIIb-61/1
+GCTGATGAACATTTTGACAATCCCTATCGGAACGTTGATCGGCCAATCAGTGGTTTCTTA
+TGATGAAGATCCCTCTATAAAGAATGATGCTATAATAGTGTATGACAACACACGAAATTG
+GATCAGTGAGGCCCAGAATTCAGATGTGGT
+>C_II_XVIIIb-59/1
+TACCACAAGAGTTCATCTGCCCTCGAATTTGAGACATGCATAGAGTATGACCCCGTTACC
+AATCTAAGCATGTTCTTAAAGGACAAGGCAATCGCACATCCGAAAGATAACTGGCTCGCC
+TCATTCAGGCGGAACCTCCTCTCTGAGGAC
+>C_II_XVIIIb-55/1
+AAGTGATCTCCGGGCGATAGTCCGATCTCATCAAGTTTTAATTTCGGGCCCTGGAGACCC
+ATCTCCTTATGTGACACAAGGGGGTGAAATGACACTCAATAAACTCTCACAACCGGTACA
+GCACCCTTCAGAGCTAGTTAAATCTGCTAC
+>C_II_XVIIIb-53/1
+CACTCGGATACCCTCATTTGATATGAGTACTACTCACTACTGTTATACTCACAATGTTAT
+ACTATCTGGTTGTAGAGATCACTCACACTCACACCAATACTTAGCACTTGGTGTGCTTCG
+GGTATCTGCGACAGGGAGGGTATTCTTTTC
+>C_II_XVIIIb-51/1
+GATCCGATCTCGTGTGTTGTCATACACTATTATAGCATCATTCTTTATAGAGGTATCTTC
+ATCATAAGAGACCACTGATTGGCCGATCAACTTTCCGCTAGAGATTGAAAGAATGTTCAT
+CAGCTCCAGTGTGGAATATGACTCCAAATT
+>C_II_XVIIIb-47/1
+ATAACCTGTGGCTCTGCCCTGTCGGGACTAGAGCCTGCCATGTCCTACCCGTGCTTTCTT
+TTGTGTGTTGGTTTCCTTTATATCTCAATGCTTTTTTTATTTTTTCTTAATAAAGTGGCC
+AAAAACAGGACTCAGACTAATCATGTCAGT
+>C_II_XVIIIb-41/1
+ACTTAGGCATGATGGAAAATCGCGACGCCGGCGGCGCCAACATTTCATCCTTAAGTGATC
+TCCGGGCGGTAGCCCGATCTCATCCAGTTTTAATTTCGGGCCCTGGAGACCCATCTCCTT
+ATGTGACACAAGGGGGTGAAATGACACTCA
+>C_II_XVIIIb-39/1
+TCCTCAGGGAGGAGGTTTCGCCTGAATGAGGCGAGCCAGTTATCTTTCGGATGTGCGATT
+GCCTTGTCCTTTAAGAACATGCTTAGATTGGTAACGGGGTCATACTCTATGCATGGCTCA
+AATTCGAGGGCAGATAAACTCTTGTACTCT
+>C_II_XVIIIb-37/1
+GGCTGCATTGTTGGCACATGACATAACAGTGTTTTCACTAAGGTCACCTGATATTAGCAC
+TAATTTAGATGAGTTTTTGAGGACCTGACTGAGTATTGCTCCGTCTTTGAATATTCGTTT
+GCTATATATGAAGAATGTGTCTGACCTGAT
+>C_II_XVIIIb-35/1
+GCGCTAGGCCGAAGCTCGAATCTTAGGGCACCTTCTACCGATATGTTGTCTGTCTTCGAC
+GAATACGAGCAGCTCCTCGCTGCCCAGACCCGTCCTAACGGAACTCACGGAGGAGGAGAG
+AAAGGGAGCACTTTAAAAGTTGAAGTCCCA
+>C_II_XVIIIb-31/1
+GCAGAAGTCAGGCGAATAGAAGAAGTGGGGGTGCTAAGTCCTAGTATTATGACCAACATC
+TTAACTAGGCCACCAGGCAATGGAGACTGGGCCAGCCTGTGCAACGACCCTTATTCCTTT
+AACTTTGAGACTGTTGCAAGCCCAAATATT
+>C_II_XVIIIb-27/1
+CACAGAGACAGATTAGCAGCCAGCTGAGCTAATAGACTTTGGTTAGACCCCGGAGGAATC
+TCAACGTCGCAATGCAGTAATGCTACAGACCTGTGGGGCACTACAGATGTGATATATCCC
+ACTGCTCTATCTGAAGTCAGGTCGCTTTCT
+>C_II_XVIIIb-25/1
+ACAAATTTCGGAAGATTGAAAGGAAGATTCAGAGTCACAATACAAGGTATGGAGAATTGT
+TCACAAAGCTGTGCACACATATTGAAAAGAAATTGCTAGGATCATCTCGATCTAATAATG
+TCCCACGGTCAGAAGAATTCAGTAGCATCC
+>C_II_XVIIIb-23/1
+TCACACTGCTGGAGTTTCGATTTTGGGTGTCATCTAAATTGATAGAACGCAGAGTAGAAA
+AGAATACCCTCCCTGTCGCAGATACCCGAAGCACACCAAGTGCTAAGTATTGGTGTGAGT
+GTGAGTGATCTCTACAACCAGATAGTATAA
+>C_II_XVIIIb-21/1
+TTCTCTGAACGTTCGACCCGTAGATTGCTGCTTAGTTTCTTATACCCGTCAGTTTTTTCT
+TGAGTTTTCGATCTGACAGGCTGTCAGAATTGACATACAGGCTGCTGCTGTATGTGCATC
+CAATTCTCATTTGCGGTCATGCTCTTGTAT
+>C_II_XVIIIb-15/1
+TGAGCCGCGATCTTAAAAGCCGAGGATAAGTTCACGACTTTGCTCAGTTCTAGCCAAGAT
+GTGACAGATAGGATATATAAACAGGTGGCTCTTGAATCTCCGCTAGCGCTACTAAACACC
+GAATCTATAATTATGAATGCAATAACATCC
+>C_II_XVIIIb-13/1
+CCTCACAGTAGCGTCTTGTGAACGACATTCTATCCAGTCCCCTACCGAGGCTGATGAAGT
+TCTTGGGAGAAAATATTGATGCTGCACTGATCGCAGCAGGGGGACAGCCCGTCCGTCCCT
+TTTGTGCAGAAAGTTTGGTGAGCACACTAA
+>C_II_XVIIIb-11/1
+TGTGCACGGTTTAATTGACGTCATCTTTTTCCCGTCTGTAGAGAGGTTATAAGGAGTAAA
+TAAGAATACTGTGTCAGCAAGGTCCCCCTCGGCATCCATGTAGATCACAGATCTAATGAC
+TGTGTCAATGTGGCTGGCAATTATCTGAGT
+>C_II_XVIIIb-7/1
+GCCATGGTGATTAAGCAGTTCCTAATCCTCTGAACAAGTATCGTACTAAGTAATCTTCCT
+TTGGGTTTGGTAGTGTGTGCTCGGGGCAAATCTCGCTAAGCGTAAATGCATCATACCTTG
+CTGGAGCAGTCAAATCTAGCTCTTGTATAA
+>HQ630064-957/1
+CAAAGCAAATTTTGTAGCTTAAGATAAAGGGGGCCCTAAAGAAGAAGCAACATCGTCTCG
+GGCTGTGTTTCCGACGTTGGGGTGCATTATGCGATTTTACTACCAGTGTTTGGTGTGGTG
+TGGTAAGGGGAAATGAGACGAGAATTTTTC
+>HQ630064-665/1
+CCCCGGGAGCTTTTTAGCGCGGCGTCACGTATAAATGCACAGAAGCTTCGTCAGTAACGA
+AGCGCCTCTGCCAGAGTTATTAACATTAGGACTCGAGAGAGTCTCGAGAAGAGACTTCCC
+AAAGAGAGTGCCGCTTTGATCTCCAGATCA
+>HQ630064-403/1
+TGTCTACCTGTCTGTCAAATTGTGACTACCACGTACCAACTCCTCCCCAAAGTTCTACCC
+TTTTCCTTAAAAACACATTTAAAATGTCAAAATACCGATCCTTATATTTGTGTGTTATGT
+TTTCTTTTGCCTTCTAGGATTGATTGAATT
+>HQ630064-343/1
+ACGAAAGAAAGAAAGGAAAAGGCCTCAAAGTTAATTTTGTAGCTTAAAAAAGAGGGGGCC
+CTAAAGAAGAAGCAACATCGTCTCGGGCTGTGTTTCCGACGTTGGGGTGCATTATGCGAT
+TTTACTACCAGTGTTTGGTGTGGTGTGGTA
+>HQ630064-307/1
+TCGTGTGTTTGCTCGCGTGCTTTGGATGTTGGGATGAGGCCTTGGTGTTGCCAGGCTCCG
+GTTTCTTAAAGGTTATCCGAGGGGGAAGTTCCACGTGTTCTTTTAGCTTTTCTCCATATT
+CCGGTCTACGGGTACAAATGTCCTCTGGAT
+>HQ630064-241/1
+CGTGACCTCATTTACAAGATGTCGGTGGCACCCAGCTACGTCACGTTCGCACCGAGGGTC
+CAACATTTTCCTAAAATGGCAAAGCGGGACTTGGGTCCGTGGGTAAATTTTCTAATTTTT
+TCCCACGGTCTGTCGCAAGAGCGTATTCTC
+>HQ630064-217/1
+AACTTTGGGGAGGAATTGGTACGTGGTAGACACAATTTGACAGACAGGTAGACAGACAAC
+GAATTACAAGCGGCCTAGGGTAATAGCGCTGAGCCCATGATGTGTCCGTGCATATGAGAC
+GACTAATAAGGTAGTCACGTAAACCCCCGT
+>HQ630064-199/1
+AGAAGCAGTATTGTGATCTTTCTCTCTTTCGCATGAATCGCTAGGAAGAGGAATTTTGGA
+AGCATCTGGATACCTTACAGGAGATTTAGAAATCTCGGTAGTTATGACATGTTCATCAGA
+TGAGGAAGATAGCAAGCGTGACGGTGTTTT
+>HQ630064-111/1
+GTGTCTCCGAGTTGTGCTCTCTCGAAACTCACTTTGCCCTCACTTGTTGGCTTTATGGAG
+GTCGTAATGATCATCGATGGAGTACTTATATAATATTGTGCTCTGCCGGCCATCTCTCCG
+TGGTCGCAGAACTCCAACACCTTGTTGGCG
+>HQ630064-105/1
+AGTTTAAACTTTGTCCGCGCGCCAATCTTGAAGTGATATTAGGACATAGTCAACCAACCG
+ACGTACGTTGAACCATTTGCGGGCAGCGCTGAGGGAAGACGAGATGTATGGACAACACAA
+CGACCGGAGAAGAGGGCGTAAAGAACATCA
+>HQ630064-11/1
+GAAAGAAAGAAAGGAAAAGGCCTCAAAGCAAATTTTGTAGCTTAAAAAAGAGGGGGCCCT
+AAAGAAGAAGCAACATCGTCTCGGGCTGTGTTTCCGACGTTGGGGTGCATTATGCGATTT
+TACTACCAGTGTTTGGTGTGGTGTGGTAAG
+>HQ630064-3/1
+AGGACTGTCGCTCATATCTTGCGTACTGTGTCTTATTTACTCGTGTTGTCTTTCCGCAGG
+GAATGCGCGGTTCGAGTAGGCCGTGGTTTGCGGCCGTTCTCGGATAAGCACAAAATATCG
+CCCTACTCGTGTGTCAAAAAGGAAAAGTCT
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out2.k1.11176.children.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out2.k1.11176.children.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,20 @@
+>C_II_XVIIIb-87/2
+AGAAGTTGTCACAAACTTTGTGAAGCTAACTTATCACCATATCCGGAGAATCATCTGACC
+TTACCTCTCTCGATACGGCTATTACTTGAATGTCACCTTGTACCATGCAGGCAAAACGAC
+AGTGTGATCTTGCTGCAGCAAGTTGGATTG
+>C_II_XVIIIb-53/2
+CTCTGAGACTTTAGTGCCCAGAACGTCACAACCTAAAGGAGTTGCACTCACACTGCAGGA
+CTTCCGACTTTGGGTGTCATCTAAGTTGATAGAACGCAGAGTAGAAAAGAATACCCTCCC
+TGTCGCAGATACCCGAAGCACACCAAGTGC
+>C_II_XVIIIb-37/2
+GATAGGTAAGGTCTGATGATTCTGCGGATCTGGTGTTAACTCAGCTGCACCAAGCTAGGG
+ATAACTTCTTTAAGGAATTAATTCATGTCAATCATCTGATTGGTCATAACTTGAAAGATC
+GTGAAACCATCAGGTCAGACACATTCTTCA
+>C_II_XVIIIb-25/2
+CAAAGACGTGGGGCATCTTATGAACTAACGTTACTAGCTTGGCGACTGCGGATCTTGTTC
+TTGCTGCTACAATCAGATGCCTTTGGACCTGTTTTATATGGAGCCATGCAAACTTGGCTC
+TGGACCACTTTGAGGGAAACCAGAATGCCG
+>C_II_XVIIIb-23/2
+ATGACCAGCATTATATCCGAGGGCTAGGCCGAGAACTCATAGTAAACGACACCAGTGATG
+TCACATCATTTTATCCTTCCGCATATCAAGAACACTTGAATTTTATCCCGGCGCCTTCTA
+CAGGGTCTGGTTGCACTCGGATACCCTCAT
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out2.k1.e10386.children.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out2.k1.e10386.children.fa Wed Mar 15 18:18:07 2023 +0000
b
b'@@ -0,0 +1,252 @@\n+>C_II_XVIIIb-101/2\n+GGGACCGCGTGCTAGCCCATAGCAAATGCCTCTCCACAGGTTGCTAAGATACTCTGGAGC\n+CAAACCGCGCACCTGCGGAGCGTGAAAGTCATCATTCAAGCTGGCACCCAGCGCGCTGTC\n+GCAGTGACCGCTGATCATGAAGTCACCTCT\n+>C_II_XVIIIb-99/2\n+TAGTGTTATCAGTCAGCTCCATGTTCCTTGGAAGATGAAAGCAAGTAAACTGCTCATCTC\n+CACTGTCACACTTTGTGCACCCTCCACTGACACTGAGAATCTCTCCTTCCACAAGTTCTA\n+TGGTGTCAGGATTAGAGACGCCTAGTATCT\n+>C_II_XVIIIb-97/2\n+GTCCCTAGCGTCTTCGTCTACTGAGGAAATATTCGCCCATTCACTCTCACCTCCTTCTCT\n+TTAAGTGGATATGATACTGCCACACTGTCATCTCTTAGGTACTCAAGGGTTGTCAAGTAT\n+TCCATCTCCTTGTATGGATCAAAATCATTT\n+>C_II_XVIIIb-95/2\n+GTTCTTGTACAGTTCAGTCATGTATAGGACGGTGTTGTTTAGGCCTCTCACTCTCAGATA\n+GTAGAGTTGGTAAGAACAGTCAAGGAGCACCTCGAGTGCCGCATATTCGAATAAGCGGAC\n+CACATCTGAATTCTGGGCCTCACTGATCCA\n+>C_II_XVIIIb-93/2\n+TAAAATTCGCTTCTTTCGGATGAGATGAACCCGGCACGGCGACATTTCGGACCAACACCA\n+ACACAGTTTCTAAATTCGGTCGTCTATAGGAATCTACAAGCGGAAGTTCCATGTAAAGAA\n+GGATTTGTCCAGGAGTTCTACCCATTATGG\n+>C_II_XVIIIb-91/2\n+ATCAACAACATAGGTACGCGAAACCTGGTAGAAGGTTTGGAGCCCATATTGCGTTCATAA\n+TGTGCTAGTTAGGTTCCAGACTCTTCTACACGTGTTTTTTCTAATTAGATAGACAGGTGA\n+ACTAGCTATAAGCAATCAAGATAGACTATT\n+>C_II_XVIIIb-89/2\n+TGTGTGATTTTACAAGGGTGCGAGCCTGTAGAGGGTTGTGATGGCAATCAACCATTCAGT\n+GCAAGGCGTTTGATCTTCCTGATCTCTTCAATCTACCTGGCTGCATCCAGCTTGCTCAGG\n+AGCTTAGCTGAGGAGCTTGGATGCATTGGG\n+>C_II_XVIIIb-87/2\n+AGAAGTTGTCACAAACTTTGTGAAGCTAACTTATCACCATATCCGGAGAATCATCTGACC\n+TTACCTCTCTCGATACGGCTATTACTTGAATGTCACCTTGTACCATGCAGGCAAAACGAC\n+AGTGTGATCTTGCTGCAGCAAGTTGGATTG\n+>C_II_XVIIIb-85/2\n+GGAAGTGCACTTGTCGGCTAGGTTGTGCAGTGGGGTGGGTCGGGGAGCGTAGGAGAGTTC\n+TGGGAGGATGGTTTGGGGTGACCTTATGGGGGCAGGCTGGGTGTTGCCGTTCAGTACCCC\n+CAGTCAGTGTCATTATCTTGGGATGGCCCA\n+>C_II_XVIIIb-83/2\n+TTAAGGCATGTTGATGCCGTGTATGCTGCCTTGATGCTGCTTGAACTCACTGGGGTTACT\n+CGACTGCGGGATATGTTAACGAATACTGCAGACACGGGATTGAGTCTTGCCTGTCCATCA\n+TCAAGCATTGTCCCGAATACCCCTCGTAGG\n+>C_II_XVIIIb-81/2\n+ATTCAGCAGAGACCGCGTGAAGTATCGCAGTTCCACAAGGCAAGACCAAGGCACTGAGCA\n+CTGCATGGGAGAAGCACGGAAGTGTCCAATCAAATGTCGGCCAGGACGCCCCCGACCAAC\n+AAGCCAGACCAGACAAACAGCCGTCCGCAC\n+>C_II_XVIIIb-79/2\n+AGAAGCACGGGCGTGTCCAATCAAATGCAGGCCAGGACGCCCCCGACCAACAAGCCAGAC\n+CAGACAAACAGCCGTCCGCACCCGAGCAGGCGGCCCAACACAGTGACCCGCCGGCCACAT\n+CCACCGAGTCCTCTCCCGCCCAGGCTGCAG\n+>C_II_XVIIIb-77/2\n+GCAATTTTGTTCAGAGATGATGTTTTCTTGTCTAGTAGATCCAATCACCCATTAGTCTCC\n+TCTAATATGTGTTCGCTGACGCTAGCAGACTATGCCCGGAACAGAAGCTGGTTACCTTTG\n+ACAGGGGGTAGGAAGATACTAGGCGTATCT\n+>C_II_XVIIIb-75/2\n+CATTGCTTCTACACCTGTGGCTCATGTCCAGCTACCTGTCGACTTTGTGCAAATGATGAT\n+GTCTATGATGGAGGCATTATCACAGAAGGTAAGTAAAGTTGACTATCAACTAGATCTAGT\n+CTTAAAACAGACATCTTCTATCCCTATGAT\n+>C_II_XVIIIb-73/2\n+TACAGTCAAGGAACGGTACTTGAGGGAACAACGGTGCCAGGATAGTGTCTCTGTCCTCTC\n+TTTCCCTGATCAGATTAAGACTCTTCCTGGACATGTAATATAGATTGGCGGGAAATGTGA\n+TTATGTTGGGAGATATGATAGAGTTCACTT\n+>C_II_XVIIIb-71/2\n+GTAATAAGAGCAGATATGCCGGTCAGTCTTACATTGACTTTGTCTAGCTTGCTGTTACTT\n+TTTGCCAACTTATCCAAGGCATTGCTTATTGAATTGTTGACATTTCCAAGCTCAATTGAT\n+ATATCGAGATTGCCGGTCACGATGACTTGA\n+>C_II_XVIIIb-69/2\n+CGAGGTCTATATGTTGACTGCCTTATCTCCTGTTATTACAATCCCCGCAGCTGCAAGGGG\n+CCTGCCATCAAGGGAGTTTGTCAGACAGACGCAGCTTAATATCAGCATAGTCCGAATGAT\n+CAGCATCAGAGGTACTGGAATCCTGGTAGA\n+>C_II_XVIIIb-67/2\n+CCCTTTCACCTGCCTTCAGCCATGTAACCCCATCCGACCGAGCCAACACAGGCGTAACTC\n+GGTTCGCCAACAATCCACTCAGAGCCCGAGGTATTAGAAAAAAATACGGGTAGAAGGGAG\n+ACACCTGGAGACCAGGACGAGTCACCAAAG\n+>C_II_XVIIIb-65/2\n+GCAAATACCACTTTGCCAGGGCCTTTATGTGCACATCAGTCTGGAGGCTTGGAGTAGAGT\n+ATGCTCAGGCTCAGGGAAGTAGCATCAATGAAGACATGGCTGCTGAACTAAAACTGACTC\n+CGGCAGCAAGGAGAGGTCTGGCGGCTGCTG\n+>C_II_XVIIIb-63/2\n+CACCTCGATTTCCGCATATTCGAATAAGCAGACAACATCTGAATTCTGGGCCTCACTGAT\n+CCAATTTCGTGTGTTGTCATACACTATTATAGCATCATTCTTGATAGAGGTATCTTCATC\n+ATAAGAAACCACTGATTGGCCGATCAACTT\n+>C_II_XVIIIb-61/2\n+GTTAGTCCTCTGACTGTCCGATAGTAGAGTTGGTAAGAACAGTCAAGGAGCACCTCGAGT\n+GCCGCATATTCGAAAAAGCGGTCCACATCTGAATTCTGGGCCTCACTGATCCAATTTCGT\n+GTGTTGTCATACACTATTATAGCATCATTC\n+>C_II_XVIIIb-59/2\n+CCTTGTATGTATAAAAATCATTTGTTTCTAATATCTCAATCAAGAGGCGATTAGTTGAGG\n+TCGCCTCTTTTATCTGTTTCTTTTGGTCCTCAGAGAGGAGGTTCCGCCTGAATGAGGCGA\n+GCCAGTTATCTTTCGGATGTGCGATTGCCT\n+>C_II_XVIIIb-57/2\n+CTGGTATGGGCCACTTTATTAAGAAAAACTAAAAAAAGCATTGTAATATAAAGGAAACCA\n+ACACACAAAAGAAAGCACGGGTAGGACATGGCAGGCTCTGGTCCCGAA'..b'CTCAT\n+>C_II_XVIIIb-21/2\n+CACAAGCAAAACAAAACACATCACTAGGGCTCGGAAACAATACCCTTGATCAGATGAGGG\n+CCAATACAAGAGCATGACCTCAAATGAGAATTGGATGCACGTACAGCAGCAGCCTGTATG\n+TCAATTCTGACAGCCTGTCAGATCGAAAAC\n+>C_II_XVIIIb-19/2\n+ACTAGATCTAGTTGATAGTCAAGTTTACTTACCTTCTGTGATAATGCCTCCATCATAGAC\n+ATCATCATTTGCACAAAGCCGACAGGTAGCTGGACATGATCCACAGGTGCAGGAGTATTG\n+TCTTGGTTCTGCCCTGACTGGGGCACATGA\n+>C_II_XVIIIb-17/2\n+GTACATATTCCCTGGATATGAGCCTGAGACCACGGGTCTGGTGCATGAGACTAATAGTTT\n+CGCAGACATTTCGATGAAATCTATATCTGCCAGCTGGTGTGACCCGTCATGGTTAATTAG\n+GCCTACTGCATTCAACCTTGGGTGTATAAT\n+>C_II_XVIIIb-15/2\n+AACTGGCGCCCCACATCCACTAGTATTCACAGACCCGTTGATTTGATAAGAGAGGGATGT\n+TATTGCATTCATAATTATAGATTCGGTGTTTAGTAGCGCTAGCGGAGATTCAAGAGCCAC\n+CTGTTTATATATCCTATCTGTCACATCTTG\n+>C_II_XVIIIb-13/2\n+TGACCTCCAGGATCTGTAGTGTGCACTGTTTAAGTGATGTTCTCTTTTTCCCGTCTGTAG\n+AGAGGTTATAAGGAGTAAATAAGAATACTGTGTCAGCAAGGTCCCCCTCGGCATCCATGT\n+AGATCACAGATCTAATGACTGTGTCAATGT\n+>C_II_XVIIIb-11/2\n+CACCTCATGCCATCGTCGTGTCACGGAAATCCTATCTAGTCCCCTACCGAGGCTGATGAA\n+GTTCTTGGGAGAAAATATTGATGCTGCACTGATCGCAGCAGGGGGACAGCCCGTCCGTCC\n+CTTTTGTGCAGAAAGTTTGGTGAGCACACT\n+>C_II_XVIIIb-9/2\n+ACAATGCCCGGTTCGGGGGATTTAGAATCGTTCCATTACTAGTTGAGATCCTCAAGGATG\n+ATAGAGTTTAAGGGATTAGACCCAGGCAATTGAATCAGCCATAGGACGGCTGGAAGAATG\n+ATATGGCACTGAATATCTTTCGCGATGCCA\n+>C_II_XVIIIb-7/2\n+CTCATGAACCGTTGCTTGAGCTTCATCCAGTACAAGACATCGGTGCTCGGGTGAAAGACC\n+CGTTAACCCGACAACCTGCATCATTTATACAAGAGCTAGATTTGAGTGCTCCAGCAAGGT\n+ATGATGCATTTACGCTTAGCGAGATTTGCC\n+>C_II_XVIIIb-5/2\n+ATCTCCACAGCAATTTTTGCTTGTGTATCCGGGAAGCGCCTGTTGCGGTTCCTTTCGAGC\n+TACTTGGATTGGCACCAGAATTAAGGATGGTAACCTCAAATAAGTTCATGTATGATCCTA\n+GCCCTATATCAGAGAAGGATTTTGCAAGAC\n+>C_II_XVIIIb-3/2\n+TCAGCTTGTATCCGTTCTGCCATTACCTGGCGGGTCCTTAGCTTCTTTGTCAACTTAGCA\n+AAAATTCGCCCATTCACTTTCACCTCCTTCTCTTTAAGTGAATATGATACTGCCACACTG\n+TCATCTCTTAGGTACTCAAGGGTTGTCAAG\n+>C_II_XVIIIb-1/2\n+AAAAGTGAGGGGGCGAACTTTTGCTAAGTTTACACAGAAGCTAAGGAACTGCCAGGTAAT\n+GGCAGAAGGGATTCTAGCTGACCAGATTGCTCCTTTCTTCCAAGGAAATGGGGTCATTCA\n+AGATAGCATATCCTTGACTAAGAGTATGTT\n+>HQ630064-957/2\n+CTGATCCAGTGGCCCTTCGTTACTCACGAAGTATCTGTGCATTTATACGTGACGCCACGC\n+TAACAAGGAAGCGCCGCGCAATTTGCTCACGGTTGGGCGTAGCTCTGTGCTCGAGCGCAC\n+GCGAAAAATTCTCGTCTCATTTCCCCTTAC\n+>HQ630064-665/2\n+TCGAATAGGCCGTGGTTTGCGGCCGTACTCGTAGAAGCACAAAATATCGCCCTACTCGTG\n+TGTCAAAAAGGAAAAGTCTTAATGGGCTGTGAAGTGATCTGGAGATCAAAGCGGCACTCT\n+CTTTGAGAAGTCTCTTCTCGAGACTCTCTC\n+>HQ630064-403/2\n+TACGGTGGGGGTCACAACAATTTAACCAATGAAATGTTTCTTTATTATAAGAGCGGTATA\n+GCAGATTTAATAAAATACCTTTTTTGAAAAAATATTGACCTGTCTCGCGTTGTAGATTCT\n+TCGCATCCGCGTTCCGTCATAATTCCAAGA\n+>HQ630064-343/2\n+CTCCCGTGCATGTATACGTGACGCCACGCTAACAAGGAAGCGCCGCGCAATTTGCTCACG\n+GTTGGGCGTAGCTCTGTCCTCGAGCGCACGTGAAAAATTCTCGTCTCATTTCCCCTTACC\n+ACACCACACCAAACACTGGTAGTAAAACCG\n+>HQ630064-307/2\n+TTTTTCACAGACGGGGACGTTCTCATCTAACTTCGGGGATCCTGGGGATATATTTTATAG\n+CGATGTAGTAGGGGGAGATGGCCCAATATTCGAAGATCCAGAGGACATTTGTACCCGTAG\n+ACCGGAATATGGAGAAAAGCTAAAAGAACA\n+>HQ630064-241/2\n+CATAATCATTTACAAAACGCGGGACCTCGACAATATGCGGGAAAATCTAACAGCAATCTA\n+TCTTGGAGAATACGCTCTTCCGAAAGACCGTGGGAAAAAATTAGAAAATTTACCCACGGA\n+CCCAAGTCCCGCTTTGCCATTTTAGGAAAA\n+>HQ630064-217/2\n+AAATCGCATGTCACACCCCAGCGTCGGAAAAACCGCCCGAGACGAAGTTGCTTCTTCTTT\n+GGGGCCCCCTCTTTTTTAAGCTACAAAATTTGCTTTGAGGCCTTTTCCTTTCTTTCTTTC\n+TTTCGCGTGTGATTGACGGGGGTTTACGTG\n+>HQ630064-199/2\n+CCGTTAGCACTTAAGTAACCGGAGCTTGGCAACACCAAGGCCTCATCCCAACATCCAAAG\n+CACTCGAGCAAACGCAGGAGTAGCAAATCTCGCTATTTCAAAGCTTCAGAATCGTGTAAT\n+ACACCGTCACGCTTGCTATCTTTCTCATCT\n+>HQ630064-111/2\n+GAAACTTGAATAAGAGCTCCCCAACTACGCTGCTAGAACTAGCAAATTTTATGCTCGCCA\n+ACAAGGTGTTGGAGTTCTGCGACCACGGAGAGATGGCCGGCAGAGCACAATATTATATAA\n+GTACTCCATCGATGATCATTACGACCTCCA\n+>HQ630064-105/2\n+TTATTAAATTGTTATTTCCACCACCGCAGAGCAATACCGTGATGTTCTTTACGCCCTCTT\n+CTCCGGTCGTTGTGTTGTCCATACATCTCGTCTTCCCTCAGCGCTGCCCGCAAATGGTTC\n+AACGTATGTCGGTTGGTTGACTATGTCCTA\n+>HQ630064-11/2\n+TGGGCGTAGTTCTGTCCTCGAGCGCACGCGAAAAATTCTCGTCTCATTTCCCCTTACCAC\n+ACCACACCAAACACTGGTAGTAAAATCGCATAATGCACCCCAACGTCGGAAACACAGCCC\n+GAGACGATGTTGCTTCTTCTTTAGGGCCCC\n+>HQ630064-3/2\n+CCACAGTTATTTGCAATAGGACTCGAGTGAGTTTTAAGAAGAGACTTCTCAAAGAGAGAG\n+CCGCTTTGATCTCCAGATCACTTCACAGCCCATTAAGACTTTTCCTTTTTGACACACGAG\n+TAGGGCGATATTTTGTGCTTATCCGAGAAC\n'
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out2.k2.11176.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out2.k2.11176.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,28 @@
+>C_II_XVIIIb-65/2
+GCAAATACCACTTTGCCAGGGCCTTTATGTGCACATCAGTCTGGAGGCTTGGAGTAGAGT
+ATGCTCAGGCTCAGGGAAGTAGCATCAATGAAGACATGGCTGCTGAACTAAAACTGACTC
+CGGCAGCAAGGAGAGGTCTGGCGGCTGCTG
+>C_II_XVIIIb-57/2
+CTGGTATGGGCCACTTTATTAAGAAAAACTAAAAAAAGCATTGTAATATAAAGGAAACCA
+ACACACAAAAGAAAGCACGGGTAGGACATGGCAGGCTCTGGTCCCGAAAGGGCAGAGCAC
+CAGCTTATCCTACCAGAGTCACATCTATCC
+>C_II_XVIIIb-49/2
+GACCCGATTTCCTTGGAAGAACGGTGCAATCTGTTCAGCTAGAATCCCTTCTGCCATTAC
+CTGGCAGTTCCTTAGCTTCTTTGTCAACTTAGCAAAAATTCTCCCATTCACTTTCACCTC
+CTTCTCGTTAAGTGAATAGGATACTGCCAC
+>C_II_XVIIIb-45/2
+GATCTTGGATGTTACAATACTAGGGCCAAGAGTCAAAGATATCAACAAAATAGGTGATGT
+AATTAGCCTAGTACTCAGAGGTAAGGTTTCCCTGGAGGACCTCATCCCACTAAGGACATA
+CCTGAAGCGCAGTACCTGCCCTAAGTATCT
+>C_II_XVIIIb-17/2
+GTACATATTCCCTGGATATGAGCCTGAGACCACGGGTCTGGTGCATGAGACTAATAGTTT
+CGCAGACATTTCGATGAAATCTATATCTGCCAGCTGGTGTGACCCGTCATGGTTAATTAG
+GCCTACTGCATTCAACCTTGGGTGTATAAT
+>C_II_XVIIIb-9/2
+ACAATGCCCGGTTCGGGGGATTTAGAATCGTTCCATTACTAGTTGAGATCCTCAAGGATG
+ATAGAGTTTAAGGGATTAGACCCAGGCAATTGAATCAGCCATAGGACGGCTGGAAGAATG
+ATATGGCACTGAATATCTTTCGCGATGCCA
+>C_II_XVIIIb-1/2
+AAAAGTGAGGGGGCGAACTTTTGCTAAGTTTACACAGAAGCTAAGGAACTGCCAGGTAAT
+GGCAGAAGGGATTCTAGCTGACCAGATTGCTCCTTTCTTCCAAGGAAATGGGGTCATTCA
+AGATAGCATATCCTTGACTAAGAGTATGTT
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out2.k2.11176.max2.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out2.k2.11176.max2.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,8 @@
+>C_II_XVIIIb-65/2
+GCAAATACCACTTTGCCAGGGCCTTTATGTGCACATCAGTCTGGAGGCTTGGAGTAGAGT
+ATGCTCAGGCTCAGGGAAGTAGCATCAATGAAGACATGGCTGCTGAACTAAAACTGACTC
+CGGCAGCAAGGAGAGGTCTGGCGGCTGCTG
+>C_II_XVIIIb-57/2
+CTGGTATGGGCCACTTTATTAAGAAAAACTAAAAAAAGCATTGTAATATAAAGGAAACCA
+ACACACAAAAGAAAGCACGGGTAGGACATGGCAGGCTCTGGTCCCGAAAGGGCAGAGCAC
+CAGCTTATCCTACCAGAGTCACATCTATCC
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out2.k2.11176.parents.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out2.k2.11176.parents.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,76 @@
+>C_II_XVIIIb-99/2
+TAGTGTTATCAGTCAGCTCCATGTTCCTTGGAAGATGAAAGCAAGTAAACTGCTCATCTC
+CACTGTCACACTTTGTGCACCCTCCACTGACACTGAGAATCTCTCCTTCCACAAGTTCTA
+TGGTGTCAGGATTAGAGACGCCTAGTATCT
+>C_II_XVIIIb-97/2
+GTCCCTAGCGTCTTCGTCTACTGAGGAAATATTCGCCCATTCACTCTCACCTCCTTCTCT
+TTAAGTGGATATGATACTGCCACACTGTCATCTCTTAGGTACTCAAGGGTTGTCAAGTAT
+TCCATCTCCTTGTATGGATCAAAATCATTT
+>C_II_XVIIIb-93/2
+TAAAATTCGCTTCTTTCGGATGAGATGAACCCGGCACGGCGACATTTCGGACCAACACCA
+ACACAGTTTCTAAATTCGGTCGTCTATAGGAATCTACAAGCGGAAGTTCCATGTAAAGAA
+GGATTTGTCCAGGAGTTCTACCCATTATGG
+>C_II_XVIIIb-87/2
+AGAAGTTGTCACAAACTTTGTGAAGCTAACTTATCACCATATCCGGAGAATCATCTGACC
+TTACCTCTCTCGATACGGCTATTACTTGAATGTCACCTTGTACCATGCAGGCAAAACGAC
+AGTGTGATCTTGCTGCAGCAAGTTGGATTG
+>C_II_XVIIIb-85/2
+GGAAGTGCACTTGTCGGCTAGGTTGTGCAGTGGGGTGGGTCGGGGAGCGTAGGAGAGTTC
+TGGGAGGATGGTTTGGGGTGACCTTATGGGGGCAGGCTGGGTGTTGCCGTTCAGTACCCC
+CAGTCAGTGTCATTATCTTGGGATGGCCCA
+>C_II_XVIIIb-75/2
+CATTGCTTCTACACCTGTGGCTCATGTCCAGCTACCTGTCGACTTTGTGCAAATGATGAT
+GTCTATGATGGAGGCATTATCACAGAAGGTAAGTAAAGTTGACTATCAACTAGATCTAGT
+CTTAAAACAGACATCTTCTATCCCTATGAT
+>C_II_XVIIIb-65/2
+GCAAATACCACTTTGCCAGGGCCTTTATGTGCACATCAGTCTGGAGGCTTGGAGTAGAGT
+ATGCTCAGGCTCAGGGAAGTAGCATCAATGAAGACATGGCTGCTGAACTAAAACTGACTC
+CGGCAGCAAGGAGAGGTCTGGCGGCTGCTG
+>C_II_XVIIIb-57/2
+CTGGTATGGGCCACTTTATTAAGAAAAACTAAAAAAAGCATTGTAATATAAAGGAAACCA
+ACACACAAAAGAAAGCACGGGTAGGACATGGCAGGCTCTGGTCCCGAAAGGGCAGAGCAC
+CAGCTTATCCTACCAGAGTCACATCTATCC
+>C_II_XVIIIb-49/2
+GACCCGATTTCCTTGGAAGAACGGTGCAATCTGTTCAGCTAGAATCCCTTCTGCCATTAC
+CTGGCAGTTCCTTAGCTTCTTTGTCAACTTAGCAAAAATTCTCCCATTCACTTTCACCTC
+CTTCTCGTTAAGTGAATAGGATACTGCCAC
+>C_II_XVIIIb-45/2
+GATCTTGGATGTTACAATACTAGGGCCAAGAGTCAAAGATATCAACAAAATAGGTGATGT
+AATTAGCCTAGTACTCAGAGGTAAGGTTTCCCTGGAGGACCTCATCCCACTAAGGACATA
+CCTGAAGCGCAGTACCTGCCCTAAGTATCT
+>C_II_XVIIIb-43/2
+GGGGCAGCGAATGAGAGAGGATATGGGGCGCCTGTTCATGACCCAGATTATATCGGAGGG
+ATAGGCAGAGAACTCATAGTAGACGACACCAGTGATGTCACATCATTTTATCCTTCCGCA
+TATCAATAACACTTGAATTTTATCCCGGCG
+>C_II_XVIIIb-33/2
+CATAGGGATAGAACATGTATGTTTTAAGACGAGATCTAGTTGATAGTCAACTTTACTTAC
+CTTCTGTGATAATGCCTCCATCATAGACATCATCATTTGCACAAAGTCGACAGGTAGCTG
+GACATGATCCACAGGTGCAGGAGCATTGTC
+>C_II_XVIIIb-29/2
+CTCAGTTGCACAAACTCGACAGGTAGCTGCATATGATCCACAGGTGCAGGAGTATTGTCT
+TGGTTCTGCCCTGACTGGGGCACATGAGGGGTTGCACCAGCTGATGGTTGTGACTCCTTC
+CATAGTCCATGATATGCTGTGTTCGCGTCT
+>C_II_XVIIIb-19/2
+ACTAGATCTAGTTGATAGTCAAGTTTACTTACCTTCTGTGATAATGCCTCCATCATAGAC
+ATCATCATTTGCACAAAGCCGACAGGTAGCTGGACATGATCCACAGGTGCAGGAGTATTG
+TCTTGGTTCTGCCCTGACTGGGGCACATGA
+>C_II_XVIIIb-17/2
+GTACATATTCCCTGGATATGAGCCTGAGACCACGGGTCTGGTGCATGAGACTAATAGTTT
+CGCAGACATTTCGATGAAATCTATATCTGCCAGCTGGTGTGACCCGTCATGGTTAATTAG
+GCCTACTGCATTCAACCTTGGGTGTATAAT
+>C_II_XVIIIb-9/2
+ACAATGCCCGGTTCGGGGGATTTAGAATCGTTCCATTACTAGTTGAGATCCTCAAGGATG
+ATAGAGTTTAAGGGATTAGACCCAGGCAATTGAATCAGCCATAGGACGGCTGGAAGAATG
+ATATGGCACTGAATATCTTTCGCGATGCCA
+>C_II_XVIIIb-5/2
+ATCTCCACAGCAATTTTTGCTTGTGTATCCGGGAAGCGCCTGTTGCGGTTCCTTTCGAGC
+TACTTGGATTGGCACCAGAATTAAGGATGGTAACCTCAAATAAGTTCATGTATGATCCTA
+GCCCTATATCAGAGAAGGATTTTGCAAGAC
+>C_II_XVIIIb-3/2
+TCAGCTTGTATCCGTTCTGCCATTACCTGGCGGGTCCTTAGCTTCTTTGTCAACTTAGCA
+AAAATTCGCCCATTCACTTTCACCTCCTTCTCTTTAAGTGAATATGATACTGCCACACTG
+TCATCTCTTAGGTACTCAAGGGTTGTCAAG
+>C_II_XVIIIb-1/2
+AAAAGTGAGGGGGCGAACTTTTGCTAAGTTTACACAGAAGCTAAGGAACTGCCAGGTAAT
+GGCAGAAGGGATTCTAGCTGACCAGATTGCTCCTTTCTTCCAAGGAAATGGGGTCATTCA
+AGATAGCATATCCTTGACTAAGAGTATGTT
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/out2.k2.exclude_both.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/out2.k2.exclude_both.fa Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,176 @@
+>C_II_XVIIIb-101/2
+GGGACCGCGTGCTAGCCCATAGCAAATGCCTCTCCACAGGTTGCTAAGATACTCTGGAGC
+CAAACCGCGCACCTGCGGAGCGTGAAAGTCATCATTCAAGCTGGCACCCAGCGCGCTGTC
+GCAGTGACCGCTGATCATGAAGTCACCTCT
+>C_II_XVIIIb-95/2
+GTTCTTGTACAGTTCAGTCATGTATAGGACGGTGTTGTTTAGGCCTCTCACTCTCAGATA
+GTAGAGTTGGTAAGAACAGTCAAGGAGCACCTCGAGTGCCGCATATTCGAATAAGCGGAC
+CACATCTGAATTCTGGGCCTCACTGATCCA
+>C_II_XVIIIb-91/2
+ATCAACAACATAGGTACGCGAAACCTGGTAGAAGGTTTGGAGCCCATATTGCGTTCATAA
+TGTGCTAGTTAGGTTCCAGACTCTTCTACACGTGTTTTTTCTAATTAGATAGACAGGTGA
+ACTAGCTATAAGCAATCAAGATAGACTATT
+>C_II_XVIIIb-89/2
+TGTGTGATTTTACAAGGGTGCGAGCCTGTAGAGGGTTGTGATGGCAATCAACCATTCAGT
+GCAAGGCGTTTGATCTTCCTGATCTCTTCAATCTACCTGGCTGCATCCAGCTTGCTCAGG
+AGCTTAGCTGAGGAGCTTGGATGCATTGGG
+>C_II_XVIIIb-83/2
+TTAAGGCATGTTGATGCCGTGTATGCTGCCTTGATGCTGCTTGAACTCACTGGGGTTACT
+CGACTGCGGGATATGTTAACGAATACTGCAGACACGGGATTGAGTCTTGCCTGTCCATCA
+TCAAGCATTGTCCCGAATACCCCTCGTAGG
+>C_II_XVIIIb-81/2
+ATTCAGCAGAGACCGCGTGAAGTATCGCAGTTCCACAAGGCAAGACCAAGGCACTGAGCA
+CTGCATGGGAGAAGCACGGAAGTGTCCAATCAAATGTCGGCCAGGACGCCCCCGACCAAC
+AAGCCAGACCAGACAAACAGCCGTCCGCAC
+>C_II_XVIIIb-79/2
+AGAAGCACGGGCGTGTCCAATCAAATGCAGGCCAGGACGCCCCCGACCAACAAGCCAGAC
+CAGACAAACAGCCGTCCGCACCCGAGCAGGCGGCCCAACACAGTGACCCGCCGGCCACAT
+CCACCGAGTCCTCTCCCGCCCAGGCTGCAG
+>C_II_XVIIIb-77/2
+GCAATTTTGTTCAGAGATGATGTTTTCTTGTCTAGTAGATCCAATCACCCATTAGTCTCC
+TCTAATATGTGTTCGCTGACGCTAGCAGACTATGCCCGGAACAGAAGCTGGTTACCTTTG
+ACAGGGGGTAGGAAGATACTAGGCGTATCT
+>C_II_XVIIIb-73/2
+TACAGTCAAGGAACGGTACTTGAGGGAACAACGGTGCCAGGATAGTGTCTCTGTCCTCTC
+TTTCCCTGATCAGATTAAGACTCTTCCTGGACATGTAATATAGATTGGCGGGAAATGTGA
+TTATGTTGGGAGATATGATAGAGTTCACTT
+>C_II_XVIIIb-71/2
+GTAATAAGAGCAGATATGCCGGTCAGTCTTACATTGACTTTGTCTAGCTTGCTGTTACTT
+TTTGCCAACTTATCCAAGGCATTGCTTATTGAATTGTTGACATTTCCAAGCTCAATTGAT
+ATATCGAGATTGCCGGTCACGATGACTTGA
+>C_II_XVIIIb-69/2
+CGAGGTCTATATGTTGACTGCCTTATCTCCTGTTATTACAATCCCCGCAGCTGCAAGGGG
+CCTGCCATCAAGGGAGTTTGTCAGACAGACGCAGCTTAATATCAGCATAGTCCGAATGAT
+CAGCATCAGAGGTACTGGAATCCTGGTAGA
+>C_II_XVIIIb-67/2
+CCCTTTCACCTGCCTTCAGCCATGTAACCCCATCCGACCGAGCCAACACAGGCGTAACTC
+GGTTCGCCAACAATCCACTCAGAGCCCGAGGTATTAGAAAAAAATACGGGTAGAAGGGAG
+ACACCTGGAGACCAGGACGAGTCACCAAAG
+>C_II_XVIIIb-63/2
+CACCTCGATTTCCGCATATTCGAATAAGCAGACAACATCTGAATTCTGGGCCTCACTGAT
+CCAATTTCGTGTGTTGTCATACACTATTATAGCATCATTCTTGATAGAGGTATCTTCATC
+ATAAGAAACCACTGATTGGCCGATCAACTT
+>C_II_XVIIIb-61/2
+GTTAGTCCTCTGACTGTCCGATAGTAGAGTTGGTAAGAACAGTCAAGGAGCACCTCGAGT
+GCCGCATATTCGAAAAAGCGGTCCACATCTGAATTCTGGGCCTCACTGATCCAATTTCGT
+GTGTTGTCATACACTATTATAGCATCATTC
+>C_II_XVIIIb-59/2
+CCTTGTATGTATAAAAATCATTTGTTTCTAATATCTCAATCAAGAGGCGATTAGTTGAGG
+TCGCCTCTTTTATCTGTTTCTTTTGGTCCTCAGAGAGGAGGTTCCGCCTGAATGAGGCGA
+GCCAGTTATCTTTCGGATGTGCGATTGCCT
+>C_II_XVIIIb-55/2
+ATGCAGTGGGCGGGAGGTGATTAATGCACGGGCCGTTTCCTTCTCCACTGCCATATCAGG
+CCCACCTGCAGTAGCAGATTTAACTAGGTCTGAAGGCTGCTGTACCGGTTGTGAGAGTTT
+ATTGAGTGTCATTTCACCCCCTTGTGTCAC
+>C_II_XVIIIb-53/2
+CTCTGAGACTTTAGTGCCCAGAACGTCACAACCTAAAGGAGTTGCACTCACACTGCAGGA
+CTTCCGACTTTGGGTGTCATCTAAGTTGATAGAACGCAGAGTAGAAAAGAATACCCTCCC
+TGTCGCAGATACCCGAAGCACACCAAGTGC
+>C_II_XVIIIb-51/2
+TTAAGGATGGTAACCTCAATTAAGTTCATGTATAATCCTAGCCCTATATCAGAGAAGGAT
+TTTGCAAGACTCGACTTAGCTATCTTCAAGAGTTATGAGCTTAATTTGGAGTCATATTCC
+ACACTGGAGCTGATGAACATTCTTTCAATC
+>C_II_XVIIIb-47/2
+TCGCACTCAATATCTTTCGCGATGCCCGGTGTTAAACTGAATAACAACACAAGCCCAGAT
+CCCGTGCTGCCGGGTAGCCACAATCTAACAACACTGACATGATTAGTCTGAGTCCTGTTT
+TTGGCCACTTTAGTAAGAAAAAATAAAAAA
+>C_II_XVIIIb-41/2
+GTCCTTCTCCACTGCCATATCAGGCCCACCTGCAGTAGCAGATTTAACTAGCTCTGAAGG
+GTGCTGTACCGGTTGTGAGAGTTTATTGAGTGTCATTTCACCCCCTTGTGTCACATAAGG
+AGATGGGTCTCCAGGGCCCGAAATTAAAAC
+>C_II_XVIIIb-39/2
+AAGGCCGCGTGTCAAGGCCGCTACAATATGCGGGAAGGTCATTGGGCAGCTTCATGCCGA
+TTCAGCGGAGATCTCACACGATGTCATGGTGAGAGAGTGCAAGAGTTTATCTGCCCTCGA
+ATTTGAGCCATGCATAGAGTATGACCCCGT
+>C_II_XVIIIb-37/2
+GATAGGTAAGGTCTGATGATTCTGCGGATCTGGTGTTAACTCAGCTGCACCAAGCTAGGG
+ATAACTTCTTTAAGGAATTAATTCATGTCAATCATCTGATTGGTCATAACTTGAAAGATC
+GTGAAACCATCAGGTCAGACACATTCTTCA
+>C_II_XVIIIb-35/2
+CTGAGTGGTTCGTTAGCATCTTCGCTAACAGCAATCCGAAGACACAATACCGCAAAATTC
+CATCTATCCTCTGGATCATCACTGTTAAGAGTAAATACTGGGACTTCAACTTTTAAAGTG
+CTCCCTTTCTCTCCTCCTCCGTGAGTTCCG
+>C_II_XVIIIb-31/2
+GTGCATCAATACTCGAATGAGTAAGGGTACAGCCAATGCTTTCTCCTCTGTCTCATTGTC
+CTCCGTATGCACTCCGGATAATAAAGGGTTTGAACAAGTCTCAAATAAGACCTTTTGTGT
+ATGCTTCTTGAGGACAATATTTGGGCTTGC
+>C_II_XVIIIb-27/2
+TATAGGAATCTACAAGCGGAAGTGTAATGTAAATAAGGATTTGTCCAGGAGTTCTACCCA
+TTATGGAGAGAGAATGCAGAAGAAAGCGACCTGACTTCAGATAGAGCAGTGGGATATATC
+ACATCTGTAGTGCCCCACAGGTCTGTAGCA
+>C_II_XVIIIb-25/2
+CAAAGACGTGGGGCATCTTATGAACTAACGTTACTAGCTTGGCGACTGCGGATCTTGTTC
+TTGCTGCTACAATCAGATGCCTTTGGACCTGTTTTATATGGAGCCATGCAAACTTGGCTC
+TGGACCACTTTGAGGGAAACCAGAATGCCG
+>C_II_XVIIIb-23/2
+ATGACCAGCATTATATCCGAGGGCTAGGCCGAGAACTCATAGTAAACGACACCAGTGATG
+TCACATCATTTTATCCTTCCGCATATCAAGAACACTTGAATTTTATCCCGGCGCCTTCTA
+CAGGGTCTGGTTGCACTCGGATACCCTCAT
+>C_II_XVIIIb-21/2
+CACAAGCAAAACAAAACACATCACTAGGGCTCGGAAACAATACCCTTGATCAGATGAGGG
+CCAATACAAGAGCATGACCTCAAATGAGAATTGGATGCACGTACAGCAGCAGCCTGTATG
+TCAATTCTGACAGCCTGTCAGATCGAAAAC
+>C_II_XVIIIb-15/2
+AACTGGCGCCCCACATCCACTAGTATTCACAGACCCGTTGATTTGATAAGAGAGGGATGT
+TATTGCATTCATAATTATAGATTCGGTGTTTAGTAGCGCTAGCGGAGATTCAAGAGCCAC
+CTGTTTATATATCCTATCTGTCACATCTTG
+>C_II_XVIIIb-13/2
+TGACCTCCAGGATCTGTAGTGTGCACTGTTTAAGTGATGTTCTCTTTTTCCCGTCTGTAG
+AGAGGTTATAAGGAGTAAATAAGAATACTGTGTCAGCAAGGTCCCCCTCGGCATCCATGT
+AGATCACAGATCTAATGACTGTGTCAATGT
+>C_II_XVIIIb-11/2
+CACCTCATGCCATCGTCGTGTCACGGAAATCCTATCTAGTCCCCTACCGAGGCTGATGAA
+GTTCTTGGGAGAAAATATTGATGCTGCACTGATCGCAGCAGGGGGACAGCCCGTCCGTCC
+CTTTTGTGCAGAAAGTTTGGTGAGCACACT
+>C_II_XVIIIb-7/2
+CTCATGAACCGTTGCTTGAGCTTCATCCAGTACAAGACATCGGTGCTCGGGTGAAAGACC
+CGTTAACCCGACAACCTGCATCATTTATACAAGAGCTAGATTTGAGTGCTCCAGCAAGGT
+ATGATGCATTTACGCTTAGCGAGATTTGCC
+>HQ630064-957/2
+CTGATCCAGTGGCCCTTCGTTACTCACGAAGTATCTGTGCATTTATACGTGACGCCACGC
+TAACAAGGAAGCGCCGCGCAATTTGCTCACGGTTGGGCGTAGCTCTGTGCTCGAGCGCAC
+GCGAAAAATTCTCGTCTCATTTCCCCTTAC
+>HQ630064-665/2
+TCGAATAGGCCGTGGTTTGCGGCCGTACTCGTAGAAGCACAAAATATCGCCCTACTCGTG
+TGTCAAAAAGGAAAAGTCTTAATGGGCTGTGAAGTGATCTGGAGATCAAAGCGGCACTCT
+CTTTGAGAAGTCTCTTCTCGAGACTCTCTC
+>HQ630064-403/2
+TACGGTGGGGGTCACAACAATTTAACCAATGAAATGTTTCTTTATTATAAGAGCGGTATA
+GCAGATTTAATAAAATACCTTTTTTGAAAAAATATTGACCTGTCTCGCGTTGTAGATTCT
+TCGCATCCGCGTTCCGTCATAATTCCAAGA
+>HQ630064-343/2
+CTCCCGTGCATGTATACGTGACGCCACGCTAACAAGGAAGCGCCGCGCAATTTGCTCACG
+GTTGGGCGTAGCTCTGTCCTCGAGCGCACGTGAAAAATTCTCGTCTCATTTCCCCTTACC
+ACACCACACCAAACACTGGTAGTAAAACCG
+>HQ630064-307/2
+TTTTTCACAGACGGGGACGTTCTCATCTAACTTCGGGGATCCTGGGGATATATTTTATAG
+CGATGTAGTAGGGGGAGATGGCCCAATATTCGAAGATCCAGAGGACATTTGTACCCGTAG
+ACCGGAATATGGAGAAAAGCTAAAAGAACA
+>HQ630064-241/2
+CATAATCATTTACAAAACGCGGGACCTCGACAATATGCGGGAAAATCTAACAGCAATCTA
+TCTTGGAGAATACGCTCTTCCGAAAGACCGTGGGAAAAAATTAGAAAATTTACCCACGGA
+CCCAAGTCCCGCTTTGCCATTTTAGGAAAA
+>HQ630064-217/2
+AAATCGCATGTCACACCCCAGCGTCGGAAAAACCGCCCGAGACGAAGTTGCTTCTTCTTT
+GGGGCCCCCTCTTTTTTAAGCTACAAAATTTGCTTTGAGGCCTTTTCCTTTCTTTCTTTC
+TTTCGCGTGTGATTGACGGGGGTTTACGTG
+>HQ630064-199/2
+CCGTTAGCACTTAAGTAACCGGAGCTTGGCAACACCAAGGCCTCATCCCAACATCCAAAG
+CACTCGAGCAAACGCAGGAGTAGCAAATCTCGCTATTTCAAAGCTTCAGAATCGTGTAAT
+ACACCGTCACGCTTGCTATCTTTCTCATCT
+>HQ630064-111/2
+GAAACTTGAATAAGAGCTCCCCAACTACGCTGCTAGAACTAGCAAATTTTATGCTCGCCA
+ACAAGGTGTTGGAGTTCTGCGACCACGGAGAGATGGCCGGCAGAGCACAATATTATATAA
+GTACTCCATCGATGATCATTACGACCTCCA
+>HQ630064-105/2
+TTATTAAATTGTTATTTCCACCACCGCAGAGCAATACCGTGATGTTCTTTACGCCCTCTT
+CTCCGGTCGTTGTGTTGTCCATACATCTCGTCTTCCCTCAGCGCTGCCCGCAAATGGTTC
+AACGTATGTCGGTTGGTTGACTATGTCCTA
+>HQ630064-11/2
+TGGGCGTAGTTCTGTCCTCGAGCGCACGCGAAAAATTCTCGTCTCATTTCCCCTTACCAC
+ACCACACCAAACACTGGTAGTAAAATCGCATAATGCACCCCAACGTCGGAAACACAGCCC
+GAGACGATGTTGCTTCTTCTTTAGGGCCCC
+>HQ630064-3/2
+CCACAGTTATTTGCAATAGGACTCGAGTGAGTTTTAAGAAGAGACTTCTCAAAGAGAGAG
+CCGCTTTGATCTCCAGATCACTTCACAGCCCATTAAGACTTTTCCTTTTTGACACACGAG
+TAGGGCGATATTTTGTGCTTATCCGAGAAC
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/read1.fastq
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/read1.fastq Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,16 @@
+@SRR5680996.16159603/1
+TGGGTCAAGCCTTGGACCCCAAACTCTTCGAAACTGACGGTCGTGGTTGCC
++
+@CCFDFFFHHHHHJJIJJJJJJJJJJJIJJJJIJIIGIIIIIIGIIGHIJJ
+@SRR5680996.2860859/1
+GACTGCTGGCTTTTGTGGACTCTGGATTCCAGGCTTTGCAACCCTAGCGGC
++
+@<@DDDFFHHFFFHEFHIIGIIIIIEIIIIGHIHGGIIIIIIIGEIIIIII
+@SRR5680996.18976829/1
+TAGAGAACAGGGCCCACCAGTCGTACTGGCTCCCCCAAAGGTCCCGACTGT
++
+CCCFFFFFHHHHHJJJJJJJIJJJIJJJJJIJJJJJEHJGIDGIJIJJIJC
+@SRR5680996.16094632/1
+GCAATTGTGATAAAAACTGCATGGTACTGGTATAGAGACAGACCAGTAAAC
++
+CCCFFFFFHHHHGJJJJJJJJJJJHHHIJJHHIJJJJJJJJIIGHIHJJJJ
\ No newline at end of file
b
diff -r c26130696a4e -r e0d3f5969458 test-data/extract_kraken_reads/read2.fastq
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extract_kraken_reads/read2.fastq Wed Mar 15 18:18:07 2023 +0000
b
@@ -0,0 +1,16 @@
+@SRR5680996.16159603/2
+CAATCTCTGGATAGATCTTGGGGGGTGCTGATGGGCTGGACTTTCCGGCTG
++
+CCCFFFFFHHFHFIIIJJJIJIIJD<BDDDDDDDDDDDDDDDDDDDDDDDD
+@SRR5680996.2860859/2
+CTTTCAGACATGTGGGCCATCCACTAGCAACAGGATCTAATTTTTTTGACA
++
+;==D??:B:D??AE<<1;AA<ACF+CEGC@@E?8D3?FDG><F>BF:FFE=
+@SRR5680996.18976829/2
+AGACTGACCAGACAGTCCTCTGTTCCTAGGGGCGGTAGGGAAGCCGAGGTA
++
+CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJ@EHHHHFFFFDCDD;@
+@SRR5680996.16094632/2
+ACTGGATGGTTTTAGGTCCCTTGTCGAAGATCAAGTGACCATAGGTGTGTG
++
+@CCFFDFFHDFHHIJJIJJJJJIIJJJJJIIJJJGFIIIJHIHIJDGFHFD
\ No newline at end of file