SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Meta assembly Autotroph Metagenomics 1 04-05-2012 12:32 PM
Meta-Velvet mihir.karnik Metagenomics 4 03-27-2012 03:41 AM
[Velvet,assembly] core dumped occured by runnning velvet matador3000 De novo discovery 0 12-17-2011 07:31 AM
velvet assembler contigs into gbrowse Zimbobo Bioinformatics 0 04-15-2010 01:36 PM
number of contigs in velvet bioenvisage Bioinformatics 6 03-24-2010 08:10 PM

Reply
 
Thread Tools
Old 10-24-2012, 05:50 PM   #1
deprekate
Junior Member
 
Location: San Diego

Join Date: Oct 2012
Posts: 4
Default meta-velvet returns nodes instead of contigs in assembly?

Hi I just started using meta-velvet to assemble a small metagenome 100k. The average length of reads are 102, with a min of 29 and a max of 227. I tried running meta-velvet with a k-mer of 21.

After the assembly is finished and I look at the output file "meta-velvetg.contigs.fa" I noticed that there are contigs of length 21 (less then the smallest read). This behavior seems strange, that it is returning just the k-mer overlap (or node) and the rest of the reads did not match, and shouldn't be reported as a contig
deprekate is offline   Reply With Quote
Old 10-25-2012, 11:41 AM   #2
cliffbeall
Senior Member
 
Location: Ohio

Join Date: Jan 2010
Posts: 144
Default

I am not sure if this applies to meta-velvet, I only have used the original. However, with standard velvet, it reports the contig lengths in kmers rather than nt in the fasta headers. So is it possible that you have contigs with 21 kmers, that are actually 41 nt long?
cliffbeall is offline   Reply With Quote
Old 10-25-2012, 01:53 PM   #3
deprekate
Junior Member
 
Location: San Diego

Join Date: Oct 2012
Posts: 4
Default

Ah, yes you are correct. The length of the actual sequence is >41 even though the header says length_21.

Last edited by deprekate; 10-25-2012 at 03:59 PM.
deprekate is offline   Reply With Quote
Reply

Tags
contigs, k-mer, velvet

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:48 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO