This is a test version of Biostars. For the public version, visit https://www.biostars.org.
A question about SPAdes

Hello everyone,

I have been using SPAdes using this command.

spades.py -k 111,113,127 --careful -o spades_output -1 s1-1.fastq -2 s1-2.fastq

If I do not misunderstand, SPAdes will choose the best k-mer based on the longest N50.

How can I know that which k-mer was chosen for the longest N50?

Thank you in advance.

spades assembly
• 1,684 views • link
updated by Istvan Albert Istvan Albert • written by Tong • Tong

don't the log state that? there are lots of files created there, just look into those

• 0 views • link
written by Istvan Albert Istvan Albert

I already looked for it in the log file, but I could not find it.

By the way, thank you for your answer.

• 0 views • link
written by Tong • Tong

I don't find this in the log file either. It says the kmers that were run but not which one was chosen. However is it always the last kmer tested?

A quick trick you can do to be sure, is just grab a part of the first header from the fastg file placed in the main output directory, e.g. "EDGE_20078450_length_31386" , and search for it within all the kmer fastg files. It is specific enough to find the original file.

• 0 views • link
written by samuel.a.odonnell samuel.a.odonnell

I looked at a SPADES output we had. Another way that might work is that the scaffolds.fasta file is created only for the chosen kmer

find . -name 'scaff*'

will tell you what directories contain this file.

• 0 views • link
written by Istvan Albert Istvan Albert

0 answers

No answers yet.

Log in to answer this question.

More posts like this