Question: how to assign kmer value for paired end reads
gravatar for girijakaushal
4.5 years ago by
girijakaushal10 wrote:


I am trying to find optimal kmer value using kmergenie so as to use that kmer value in Metavelevet tool for metagenome assembly.

Q. 1) My first query is related to size of my left and right read's file after trimming sequences with Phred score<20. So, When I did trimming of such reads it reduced the orignal size of my read files as earlier my both left and right reads were 10,014 MiB but after using NGSqc toolkit trimming tool, My left read file was 9412 MiB and right reads file was 9288 MiB. So do it is any problem? why their file size is different as earlier they were same.

Q. 2) As I am having paired end reads of 102bp so I firstly calculated for left reads that gave me best k=25 and then I calculated it for right read and it gave me beat k=21. So what kmer value should I take to assemble these left and right reads?

Please guide me, I would be heartily thankful.

Best regards

assembly • 1.8k views
ADD COMMENTlink modified 4.5 years ago by dbrowne.up70 • written 4.5 years ago by girijakaushal10

I'd use a program like to make sure that reads are properly paired. Following command can do that, it will also remove read shorter than the length 10. -fastq file_R1.fastq -fastq2 file_R2.fastq -min_len 10 -out_bad null -out_good file_clean

You can re-run kmergenie on the prinseq output files. Also you should check out other de novo assembly tools like SPAdes and IDBA-UD that can use multiple k-mer values for the de novo assembly.

ADD REPLYlink written 4.5 years ago by Sej Modha4.8k
gravatar for dbrowne.up
4.5 years ago by
United States
dbrowne.up70 wrote:

Personally, I don't think that a single k-mer value will give you an optimal assembly. Different k-mer sizes will optimally assemble different regions of the genome. Since you're doing metagenomic assembly, check out this program called MeGAMerge:

It utilizes assemblies from multiple k-mer values, as well as long reads, if you have them, to build an improved draft genome assembly.

ADD COMMENTlink written 4.5 years ago by dbrowne.up70
Please log in to add an answer.


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