Question: Increase the kmer size in IDBA_UD?
0
gravatar for mail2steff
11 months ago by
mail2steff50
Potsdam, Germay
mail2steff50 wrote:

Dear all

Ive pair-end reads from papaya samples with average length of 151 bp. I used IDBA_UD for assembling. I increased the read length and kmer size by changing the default value as follows

static const uint32_t kMaxShortSequence = 200;

static const uint32_t kNumUint64 = 16;

./idba_ud -r reads.fa -o output mink 20 maxk 140

I have two questions.

  1. When I executed IDBA_UD, It gave the contigs for the kmer size 20-100 only. I did not get any contigs for the kmer more than 100. Where I have gone wrong?
  2. My second question is, The program internally decides the best kmer and produces the contig.fa file. How do we know which kmer is used to create this contig files?
kmer idba idba_ud ngs illumina • 463 views
ADD COMMENTlink modified 11 months ago by Sej Modha4.1k • written 11 months ago by mail2steff50
0
gravatar for Sej Modha
11 months ago by
Sej Modha4.1k
Glasgow, UK
Sej Modha4.1k wrote:
  1. It is possible to control the min and max kmer sizes as well as the increment of the kmer steps using following parameters:

    --mink arg (=20) minimum k value (<=504)

    --maxk arg (=100) maximum k value (<=504)

    --step arg (=20) increment of k-mer of each iteration

  2. If program used multiple kmer values than each output would be saved and contigs files labelled contgs.*.fa save the corresponding contigs in the fasta format. Final contigs and scaffolds are saved in the contigs.fa and scaffolds.fa files that contain the consolidated contigs and scaffolds from each kmer assembly.

ADD COMMENTlink written 11 months ago by Sej Modha4.1k

Dear Sej Modha, Thank you for the commands. I increased the kmer value as above mentioned and the set the increment of k-mer of each iteration as 20. Eventhough I got the contigs for the kmwe 20, 40, 60, 80 and 100. Not more than that. Will IDBA_UD produce contigs upto 100 or AM I missing something?

With respect to the 2nd answer, How do I know which kmer is the best?

ADD REPLYlink written 11 months ago by mail2steff50
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: 1709 users visited in the last hour