SEQanswers

Go Back   SEQanswers > Applications Forums > De novo discovery



Similar Threads
Thread Thread Starter Forum Replies Last Post
How can be possible? (Velvet) aloliveira Bioinformatics 4 09-28-2012 12:43 AM
[Velvet,assembly] core dumped occured by runnning velvet matador3000 De novo discovery 0 12-17-2011 07:31 AM
Velvet assembler bioinf Bioinformatics 31 08-24-2011 09:19 AM
velvet bioenvisage Bioinformatics 1 02-17-2010 08:16 AM
velvet ANJAN PURKAYASTHA Bioinformatics 0 03-09-2009 07:38 AM

Reply
 
Thread Tools
Old 02-12-2010, 07:12 AM   #1
bioenvisage
Member
 
Location: it

Join Date: Oct 2009
Posts: 40
Default velvet

hi After running velvet assembly I got Estimated Coverage = 25.12000
Estimated Coverage cutoff = 10.250000
Final graph has 11106 nodes and n50 of 498, max 13308, total 243547, using 221137/2500000 reads, so I have 3 questions here.

1) what does this max 13308 ( It is the maximum contig length obtained)
2)what does this total 243547 refers to?
3) should I run the velvet again by using the obtained estimated cov cutoff and estimated coverege (for eg here should I run velvet again using cov_cutoff 10
est_coverege 25)

can any one answer this???

thanks
bioenvisage is offline   Reply With Quote
Old 02-12-2010, 08:42 AM   #2
lcollado
Member
 
Location: Baltimore, MD

Join Date: Jun 2009
Posts: 65
Default

1) Yes, it's the length of the longest contig you got.
2) It's the sum of the length of all your contigs.
3) You should probably try out the velvetOptimizer.pl script that is distributed along Velvet. Or you can join the Velvet mailing list and ask there.
__________________
L. Collado Torres, Ph.D. student in Biostatistics.
lcollado is offline   Reply With Quote
Old 02-12-2010, 11:11 AM   #3
Zigster
(Jeremy Leipzig)
 
Location: Philadelphia, PA

Join Date: May 2009
Posts: 116
Default

Quote:
using 221137/2500000 reads
that is a very low read usage - i would focus on figuring out why so few of your reads are making it into the assembly. Maybe they need to be trimmed further?
__________________
--
Jeremy Leipzig
Bioinformatics Programmer
--
My blog
Twitter
Zigster is offline   Reply With Quote
Old 02-13-2010, 04:05 PM   #4
bioenvisage
Member
 
Location: it

Join Date: Oct 2009
Posts: 40
Default

hi zingster thnx for ur reply ...but my read length is 75 bp, do u think it should be trimmed further.
bioenvisage is offline   Reply With Quote
Old 02-14-2010, 09:25 AM   #5
Zigster
(Jeremy Leipzig)
 
Location: Philadelphia, PA

Join Date: May 2009
Posts: 116
Default

Trim on quality, not length
Joseph Fass's script does soft trimming - try 10, 15, and 20

http://wiki.bioinformatics.ucdavis.e...rimBWAstyle.pl
__________________
--
Jeremy Leipzig
Bioinformatics Programmer
--
My blog
Twitter
Zigster 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 07:01 AM.


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