Open jflot opened 2 years ago
Hi @jflot,
What happened here is that no reads were mapped, as suspected by you. I have added an informative error message.
What I find suspicious is the line INFO, skch::Sketch::build, minimizers picked from reference = 0
- are your reference contigs very short or very repetitive?
You could also look at the file classification_reads20000
, which may be informative.
Best
Alex
Hello, I am trying to run MetaMaps to classify Nanopore reads that I know originate from a mix of four genomes. I built a custom database with these four (unpublished) genomes:
then
I did the mapping step:
but got an the following error message during the classification step:
It looks like none of my reads mapped on my reference genomes, but when I try mapping them using minimap2 for instance they almost all map on at least one of my four reference genomes... Thanks a lot in advance for your help.