-
Description of bugI am working on a pipeline of running a metagenomic assembly with metaSPAdes, followed by metaViralSPAdes. I am attempting to use the assembly graph as input to metaviralspades to save time. The error I get is stating that the kmer length in the assembly graph does not match the cli settings (I guess the default for metaviralSPAdes is 21. How do I retrieve the kmer size from the assembly graph, so I can use it to set k properly for metaViralSpades? I can tell it gets parsed here, but I can't make heads nor tails of it. Thanks in adnvance spades.logparams.txtSPAdes versionSPAdes-3.15.4 Operating SystemLinux-3.10.0-1160.45.1.el7.x86_64-x86_64-with-glibc2.29 Python Version3.8.10 Method of SPAdes installationDocker No errors reported in spades.log
|
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 1 reply
-
Hello You need to use last k-mer length if the input is assembly graph. If you have used default set of k-mers for metaSPAdes (did not alter |
Beta Was this translation helpful? Give feedback.
Hello
You need to use last k-mer length if the input is assembly graph. If you have used default set of k-mers for metaSPAdes (did not alter
-k
option), then it would be 55.