SEQanswers

Go Back   SEQanswers > Applications Forums > De novo discovery



Similar Threads
Thread Thread Starter Forum Replies Last Post
memory requirments of velvet tool (de novo assembly) bioinfosm Bioinformatics 12 04-19-2012 04:26 AM
de novo assembly of repeat elements HESmith Bioinformatics 12 11-09-2011 04:11 PM
Velvet de novo assembly to amosvalidate canuck Bioinformatics 5 07-17-2011 12:24 PM
de novo assembly (velvet or others) strob Bioinformatics 1 01-20-2010 05:53 AM
Velvet de novo assembly of Solid reads HOWTO KevinLam De novo discovery 1 01-10-2010 01:11 AM

Reply
 
Thread Tools
Old 07-12-2011, 07:50 AM   #1
salmonella
Junior Member
 
Location: texas

Join Date: Feb 2011
Posts: 5
Default how to resolve repeat areas with Velvet when doing de novo assembly

I am sequencing a bacterial genome and have assembled my Illumina reads (40 bp single) using Velvet. I am particularly interested in a particular gene (my gene of interest, MGIS) along with its neighboring genes that is probably located on a low copy number plasmid. One complicating factor is that MGIS is flanked by two genes that are repeated in the genome. I have tried to optimize parameters (k mer length, no gaps allowed, no mismatches allowed) in Velvet, and thus far, the largest contig that contains MGIS is only about 4 kb.
1. Are there any other parameters that I can change that may help to increase my contig size?
2. I have about 6 kb of sequence from MGIS and it flanking regions obtained from manual sequencing. Can I use this as a ‘reference’ to make Velvet or any other assembly program begin the assembly from this point and work outwards from this point to reassemble as much of the region as possible? I believe that I have enough read coverage, but I am guessing that because of the way that the assembly is working, the reads are being sequestered?

If anyone has any suggestions on some alternative approaches to reconstructing the plasmid based on the existing Illumina data, I would very much like to hear them.
salmonella is offline   Reply With Quote
Old 10-24-2011, 09:42 PM   #2
av_d
Member
 
Location: Pune, India

Join Date: Sep 2009
Posts: 12
Default

May be you can the optimize "-exp-cov" and "-cov_cutoff" , also a higher value of K may help you to resolve the repeat structure.
av_d 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 11:21 AM.


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