<p>There is some description of homocompressed polymers (-p kmer) in this paper: <a href="https://arxiv.org/pdf/1708.01492.pdf" rel="nofollow noreferrer">https://arxiv.org/pdf/1708.01492.pdf</a></p>

CC BY-SA 4.0
2020-09-01T18:53:49.183Z
2020-09-01T18:53:49.183Z
Andor Kiss
<p>I am using <a href="https://github.com/ruanjue/wtdbg2" rel="noreferrer">wtdbg2</a> 2.3 to assemble a human genome (sequenced on PromethION from a cell line). I filtered out reads with low average quality, and now I am trying to determine the parameters that will optimize overall assembly quality.</p>

<p>The wtdbg2 assembler has two k-mer related settings, documented as follows.</p>

<pre><code> -k &lt;int&gt;    Kmer fsize, 0 &lt;= k &lt;= 25, [0]
 -p &lt;int&gt;    Kmer psize, 0 &lt;= p &lt;= 25, [21]
             k + p &lt;= 25, seed is &lt;k-mer&gt;+&lt;p-homopolymer-compressed&gt;
</code></pre>

<p>The requirement that k + p is &lt;= 25<sup>**</sup> suggests that these parameters can be set simultaneously. However, in all of the presets for the various long read sequencing technologies, one of the two settings is always set to 0.</p>

<p>My questions are as follows.</p>

<ol>
<li>Are <code>k</code> and <code>p</code> ever intended to be used simultaneously?</li>
<li>If so, how does the behavior change with <code>-k 2 -p 19</code> versus <code>-k 19 -p 2</code>?</li>
<li>I'm familiar with some common implications of choosing k size for k-mer based analyses of Illumina reads. (Smaller k-sizes are more sensitive but also more subject to spurious matches. Larger k-sizes are more specific but also more likely to include a sequencing error.) Do similar implications apply for <code>-k</code> and <code>-p</code> with wtdbg2?</li>
<li>BONUS: Is there an explanation anywhere of how <code>-k</code> and <code>-p</code> are used? It looks to me like <code>-k</code> is related to k-mers in the traditional sense, and <code>-p</code> is a homopolymer compressed k-mer. Are seeds made by concatenating the regular k-mer and the compressed k-mer? It's possible I missed some important details in <a href="https://doi.org/10.1101/530972" rel="noreferrer">the preprint</a> but I haven't been able to make the connection yet.</li>
</ol>

<hr>

<p><sup>**</sup>Incidentally, the documentation here is wrong. During runtime wtdbg2 enforces a constraint of &lt;= 23, not 25. This matches the description in the bioRxiv preprint stating that 46 bits (2 bits per base) are used to store each k-mer.</p>

<genome><nanopore><assembly><long-reads>
wtdbg2: practical implications of k-mer fsize and psize choice
# wtdbg2: practical implications of k-mer fsize and psize choice

(Asked by: username_0 on 2019-02-14)

I am using wtdbg2 2.3 to assemble a human genome (sequenced on PromethION from a cell line). I filtered out reads with low average quality, and now I am trying to determine the parameters that will optimize overall assembly quality.
The wtdbg2 assembler has two k-mer related settings, documented as follows.
 -k     Kmer fsize, 0 <= k <= 25, [0]
 -p     Kmer psize, 0 <= p <= 25, [21]
             k + p <= 25, seed is +

The requirement that k + p is <= 25** suggests that these parameters can be set simultaneously. However, in all of the presets for the various long read sequencing technologies, one of the two settings is always set to 0.
My questions are as follows.

Are k and p ever intended to be used simultaneously?
If so, how does the behavior change with -k 2 -p 19 versus -k 19 -p 2?
I'm familiar with some common implications of choosing k size for k-mer based analyses of Illumina reads. (Smaller k-sizes are more sensitive but also more subject to spurious matches. Larger k-sizes are more specific but also more likely to include a sequencing error.) Do similar implications apply for -k and -p with wtdbg2?
BONUS: Is there an explanation anywhere of how -k and -p are used? It looks to me like -k is related to k-mers in the traditional sense, and -p is a homopolymer compressed k-mer. Are seeds made by concatenating the regular k-mer and the compressed k-mer? It's possible I missed some important details in the preprint but I haven't been able to make the connection yet.


**Incidentally, the documentation here is wrong. During runtime wtdbg2 enforces a constraint of <= 23, not 25. This matches the description in the bioRxiv preprint stating that 46 bits (2 bits per base) are used to store each k-mer.

**Tags:** genome, nanopore, assembly, long-reads

---

## 1 answers

### Answer by username_1 on 2020-09-01. Score: 1

There is some description of homocompressed polymers (-p kmer) in this paper: <https://arxiv.org/pdf/1708.01492.pdf>


---

