Question: kmer for sample list
0
gravatar for ARich
22 months ago by
ARich60
United States
ARich60 wrote:

Dear All,

I am working with paired end metagenomic illumina data. 

I need to know that on providing list file with multiple sample entries does kmergenie decide kmer for all the samples provide in the list or how it treats the sample list?

I means for metagenomic assembly is feasible to have one kmer for all samples or do you recommend to run kmergenie with all the samples separately.

Thanks

Best,

AR

kmergenie • 747 views
ADD COMMENTlink modified 19 months ago • written 22 months ago by ARich60
0
gravatar for Rayan Chikhi
22 months ago by
Rayan Chikhi1.2k
France, Lille, CNRS
Rayan Chikhi1.2k wrote:

Hi,

 

Kmergenie most likely won't work for metagenomic data, I'm sorry. It expects a single genome. Hence, if your data contains more than one genome, it is quite likely that it won't work. (the model won't be able to fit the kmer histograms that have been generated.)

If you have separate samples and each sample contains a single genome, then it is recommended to run kmergenie on samples separately. However, if each sample contains multiple genomes, then kmergenie most likely won't work.

 

The good news is that for metagenomics assembly, you could try multi-k assemblers, where you do not need to specify a best k value: SPAdes or Megahit.

ADD COMMENTlink written 22 months ago by Rayan Chikhi1.2k
0
gravatar for ARich
19 months ago by
ARich60
United States
ARich60 wrote:

Dear Rayan,

Thanks Yes i tested it with Spaded it does have metagenomic option optimizing to the best kmer.

ADD COMMENTlink written 19 months ago by ARich60
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: 588 users visited in the last hour