Entering edit mode
10 months ago
Art1ess
•
0
I am building the database of kraken2. The previous add-to-library used my custom taxID and wrote names.dmp and nodes.dmp myself.
When running kraken2-build --build --db
., the following prompt appears:
Creating sequence ID to taxonomy ID map (step 1)...
Sequence ID to taxonomy ID map already present, skipping map creation.
Estimating required capacity (step 2)...
Estimated hash table requirement: 21298865004 bytes
Capacity estimation complete. [11m0.264s]
Building database files (step 3)...
Taxonomy parsed and converted.
CHT created with 13 bits reserved for taxid.
Completed processing of 0 sequences, 0 bp
Writing data to disk... complete.
Database files completed. [2m59.229s]
Database construction complete. [Total: 14m0.492s]
Obviously, running kraken2 after the build is complete does not produce the correct results:
U CL100096888L1C013R048_557162/1 unclassified (taxid 0) 100 A:29 0:12 A:29
U CL100096888L1C012R064_96906/1 unclassified (taxid 0) 100 0:20 A:31 0:19
U CL100096888L1C004R035_314813/1 unclassified (taxid 0) 100 A:70
U CL100096888L1C017R037_497775/1 unclassified (taxid 0) 100 A:29 0:10 A:31
U CL100096888L1C012R006_314272/1 unclassified (taxid 0) 100 0:70
U CL100096888L1C017R067_125180/1 unclassified (taxid 0) 100 A:59 0:11
U CL100096888L1C017R067_147919/1 unclassified (taxid 0) 100 A:59 0:11
U CL100096888L1C003R044_122470/1 unclassified (taxid 0) 100 0:7 A:52 0:4 A:7
U CL100096888L1C009R091_47094/1 unclassified (taxid 0) 100 0:70
U CL100096888L1C011R005_253537/1 unclassified (taxid 0) 100 0:70
What do you get instead? Error/no output? Have you looked at the log files?
I have 2 files output, no error logs...
the .kreport file:
head .result.kraken file:
this is names.dmp:
and nodes.dmp
Were you able to fix this issue? I'm also trying to add external species to krakendb