Question: Regarding K mer analysis
0
gravatar for DL
19 months ago by
DL10
India
DL10 wrote:

Hello everyone,

I am going to do genome assembly for that i did kmer analysis to see that which kmer value is good for the assembly. genome is heterozygous so i found two peaks in kmer analysis. i performed kmer analysis using kmer value 21, 32 and 64 now i am confused which kmer value is best for assembly. At kmer value=64, i got large genome size. Please explain which parameter i will see for good assembly. K_mer_analysis

Please Reply as soon as possible

Thanks in advance

sequence next-gen assembly genome • 1.0k views
ADD COMMENTlink modified 19 months ago by WouterDeCoster38k • written 19 months ago by DL10

From my perspective, everything in your figure says to me that the best assembled genome out of the three is the one with kmer=64. I looked at each statistic and kmer=64 has the best result for each, including the narrowest gap between min and max. A larger k-mer is ideal as it essentially provides for more 'uniqueness' in the base sequences that will form the assembled genome. The downside to using a large k-mer is the computational expense associated with it, as the resulting de Bruijn graph that is generated using such a larger k-mer would be complex.

Some questions for you:

  • Why did you choose these k-mers in particular?
  • How are the alignment rates if you re-align your sample(s) to each assembled genome? - I assume that the highest alignment is at k-mer=64
ADD REPLYlink written 19 months ago by Kevin Blighe41k

Thank You Kevin for your reply.

i have one doubt. what should i do if the genome is heterozygous ?? for that 63 kmer is ohk ???

ADD REPLYlink written 19 months ago by DL10

For that, I cannot comment - apologies. I just know that by the standard measures, the k-mer=63 genome looks better.

You may find this thread of use: Recommendations For Heterozygous Genome Assembly Software

I also found this recent (2014) publication: https://www.ncbi.nlm.nih.gov/pubmed/24755901

Apologies that I cannot assist further!

ADD REPLYlink written 19 months ago by Kevin Blighe41k
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: 832 users visited in the last hour