SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
plasmid annotation biocc Bioinformatics 3 07-10-2011 05:06 PM
how to detect ori of plasmid biocc General 2 02-24-2011 08:15 PM
Plasmid DNA Extraction Kit Orochem Vendor Forum 1 07-29-2010 09:54 AM
Plasmid contamination in Long Tag Paired End library kmcarr 454 Pyrosequencing 3 03-11-2009 04:29 PM

Reply
 
Thread Tools
Old 04-03-2012, 10:20 AM   #1
*#1*
Junior Member
 
Location: East coast

Join Date: Jul 2010
Posts: 8
Default How to circularize plasmid assembly

Hi,
My assembly resulted in a single contig, which is linear. How do I check where is the complimentary sequences on the ends so that I could circularize the plasmid? Is there a tool that can do that easily?

The thing is I have ran several assemblies and got several one-contigs, but they differ slightly in size, so I expect the difference to be in tails. I have aligned the contigs in Mauve and see that 95% of contigs overlap, but I still cannot figure out where to trim the tails so that I could circularize the plasmid.

I would appreciate any help and pointers.
Thanks!
*#1* is offline   Reply With Quote
Old 04-04-2012, 02:58 AM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,541
Default

You can use tools like mummerplot to compare the contig to itself, and thus deduce the circle.

Since (AFAIK) no read mappers support a circular reference, you can check the circle like so: Once you've made the circle, try mapping the reads to it. Then shift the origin (say half way along) and repeat this mapping. Both times you should get nice uniform coverage. Or, take two copies of the circle as a single long reference, and map the reads to that - again you're looking for nice uniform coverage.
maubp is offline   Reply With Quote
Old 04-04-2012, 07:14 AM   #3
ians
Member
 
Location: St. Louis, MO

Join Date: Aug 2011
Posts: 53
Default

What you could do is chop your sequence in half (at what ever break point is appropriate) and then run it through minimus2 or another OLC. If you get one contig back, you now have a perfect linear representation of a circular plasmid. Now you know that the ends can be connected bluntly.
ians 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:34 PM.


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