Struo2 icon indicating copy to clipboard operation
Struo2 copied to clipboard

0 Bracken added reads with GTDB database for all taxonomic units

Open alimayy opened this issue 2 years ago • 4 comments

Hi there,

Thanks for developing Struo2 and for releasing the GTDB kraken2 dbs prepared with struo2.

I'm using Kraken v2.1.1 and Bracken v2.8 on paired-end animal microbiome data and the struo2 GTDB_release207 kraken2 DB GTDB . The files I have under in my kraken2/bracken DB are as follows with matching md5sums to those in the mentioned DB repository.

-rw-rw-r--  1 ali  ali  287G dec 21 18:53 hash.k2d
-rw-rw-r--  1 ali  ali    64 dec 21 18:53 opts.k2d
-rw-rw-r--  1 ali  ali  9,8M dec 21 18:53 taxo.k2d
-rw-rw-r--  1 ali  ali    59 dec 21 18:53 database150mers.kmer_distrib
-rw-r--r--  1 ali  ali  578M dec 21 18:53 database150mers.kraken

When I run kraken2 and Bracken using this database, I don't get any added reads from Bracken. Any idea why this might be the case, e.g. could I be missing some required files/folders here, or a kraken2 or Bracken version issue?

image

alimayy avatar Dec 21 '22 20:12 alimayy