SEQanswers

Go Back   SEQanswers > Applications Forums > De novo discovery



Similar Threads
Thread Thread Starter Forum Replies Last Post
Merging Velvet Assemblies millermr Bioinformatics 20 02-26-2015 12:02 PM
Using bam file as the input for velvet elisadouzi Bioinformatics 1 04-13-2012 04:27 PM
velvet assemblies rahularjun86 Bioinformatics 3 01-31-2012 02:59 AM
Merging Velvet Assemblies millermr General 0 11-23-2010 08:59 PM
Fraction of reads incorporated in Velvet assemblies foram Bioinformatics 3 02-04-2009 05:36 PM

Reply
 
Thread Tools
Old 11-20-2013, 03:43 AM   #1
rahularjun86
Member
 
Location: Frankfurt(M), Germany

Join Date: Jan 2011
Posts: 58
Default Does velvet generate different assemblies from the same input data and same velvet pa

Hi all,

I am using four different libraries of read size 76bps, insert sizes are 300bps, 1kb, 8kb and 12kb. Expected genome size is 80MB.

I am running velvet using these four libraries. Actually I tried velvet first for different k-mers and then for the best k-mer I tried different cov_cutoffs. In all these assemblies from same k-mer (69) with different cov_cutoffs, I used the same Roadmap and Sequences files from the initial velveth run (K-mer 69 and cov_cutoff default).

Surprisingly I got 10MB of N50 and 23MB of largest scaffold size using cov_cutoff of 12 (median coverage is 30.76 in Log file) on the previously generated Roadmap and Sequences files. Then later I tried all new assembly from the same reads, K-mer 69 and cov_cutoff 12, now my N50 is 2MB and largest scaffold size is 6.78 MB.

Later I tried the same input files and same velvet parameters and I figured out that velveth is generating different Roadmap file for the same k-mer 69, for all three runs. What could be the reason behind this? In this case it is not possible to regenerate the results.

My question is, does velvet generates different assemblies from the same data and with same parameters?

I would really appriciate your response on this.

Best regards,
Rahul
__________________
Rahul Sharma,
Ph.D
Frankfurt am Main, Germany
rahularjun86 is offline   Reply With Quote
Old 11-20-2013, 05:43 AM   #2
cliffbeall
Senior Member
 
Location: Ohio

Join Date: Jan 2010
Posts: 144
Default

I saw something about this issue on the velvet-users email list a while back - check out the thread titled "Deterministic?":

http://listserver.ebi.ac.uk/pipermai...er/001734.html
cliffbeall is offline   Reply With Quote
Old 11-21-2013, 06:43 AM   #3
SES
Senior Member
 
Location: Vancouver, BC

Join Date: Mar 2010
Posts: 275
Default

This was answered on biostars, for reference. It is related to OpenMP, and you can be certain of the results if you set the thread level to one.
SES is offline   Reply With Quote
Reply

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 05:18 PM.


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