Question: Kmergenie fails to fit models with --diploid settings
0
gravatar for pepijn.kooij
20 months ago by
Royal Botanic Gardens, Kew
pepijn.kooij0 wrote:

I’ve been trying to run kmergenie on a trimmed dataset of MiSeq 2x300 paired end reads (trimmed with trimmomatic), using a list with the forward and reverse, paired and unpaired output from trimmomatic. When I run kmergenie with the normal haploid settings it works fine, but running it with --diploid many of the many of the kmer models are not able to fit the on the histograms (see screenshot of terminal+html files). Any idea why this is?

screenshot of terminal

results without --diploid

results with --diploid

I really would like to be able to fit the diploid models, which I think would improve the kmer estimations, based on the fact that you can clearly see two humps in the histograms.

Hope you can help me out with this.

miseq kmergenie paired-end kmer • 691 views
ADD COMMENTlink modified 8 months ago by Rayan Chikhi1.2k • written 20 months ago by pepijn.kooij0
0
gravatar for Rayan Chikhi
8 months ago by
Rayan Chikhi1.2k
France, Lille, CNRS
Rayan Chikhi1.2k wrote:

for the record: Pepijn and I communicated by email, and the diploid model has been improved in recent kmergenie versions, hopefully addressing this issue.

ADD COMMENTlink written 8 months ago by Rayan Chikhi1.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 576 users visited in the last hour