SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   De novo discovery (http://seqanswers.com/forums/forumdisplay.php?f=27)
-   -   Does velvet generate different assemblies from the same input data and same velvet pa (http://seqanswers.com/forums/showthread.php?t=36006)

rahularjun86 11-20-2013 03:43 AM

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

cliffbeall 11-20-2013 05:43 AM

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

SES 11-21-2013 06:43 AM

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.


All times are GMT -8. The time now is 07:03 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.