SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why edgeR ouput some p values larger than 1 idyll_ty RNA Sequencing 1 02-01-2013 03:14 PM
SOAP denovo assembly results bioinfosm Bioinformatics 9 06-13-2012 07:22 AM
samtools (or other tools) to find larger indels waalkes Bioinformatics 0 06-15-2011 02:53 PM
Ion Torrent de novo assembly results nickloman Ion Torrent 2 05-11-2011 09:27 PM
illumina de novo txome assm: overlapping paired ends or long inserts? Zigster De novo discovery 1 10-15-2009 04:46 AM

Reply
 
Thread Tools
Old 01-22-2012, 02:06 PM   #1
martin_313
Junior Member
 
Location: Germany

Join Date: Dec 2011
Posts: 6
Default assembly with unsatisfying results - use new reads with larger inserts?

We recently finished an assembly of Guignardia bidwellii paired end DNA reads using Abyss and refined the results with SSpace. However, we were not quite satisified with the results (still to many contigs).
We are now considering generating new libraries with larger insert sizes between paired ends to resolve repeats. Is there a genereal guideline for choosing those insert sizes? From what I understood so far, an optimal insert size would be larger than the length of the largest repeat appearing in the genome (is this correct?).
Sadly, we do not have too much information about the genome for estimating the size of repeats in the genome. Is there some way to make an eductaed guess based on the Abyss-output files (like the contig linkage graph) of our last assembly? Or should we simply go for the largest insert size possible?

Thank your for answers
Martin
martin_313 is offline   Reply With Quote
Old 01-22-2012, 04:30 PM   #2
peromhc
Senior Member
 
Location: Durham, NH

Join Date: Sep 2009
Posts: 108
Default

before talking about making new libs, tell us a little about the current libs and sequencing

How many libs do you currently have, what insert sizes?
How deep is your sequencing currently--10X, 20X, more?
How big is your genome.
peromhc is offline   Reply With Quote
Old 01-22-2012, 05:21 PM   #3
MadsAlbertsen
Member
 
Location: Denmark

Join Date: Aug 2010
Posts: 26
Default

The first thing I do after an assembly is always a contig length vs. coverage plot to estimate the repeat content. It's usually quite easy to see the repeats.

rgds
Mads
MadsAlbertsen is offline   Reply With Quote
Old 01-23-2012, 12:37 AM   #4
martin_313
Junior Member
 
Location: Germany

Join Date: Dec 2011
Posts: 6
Default

Thank you for your answers.

peromhc, since we received the libraries in a single file without too much information, i don't know the exact number of libraries. The insert sizes range from 100 to 500 with the majority between 200 and 400.
All reads were 101 bp long.
The average sequencing depth seems to be around 50. The size of the genome is estimated to be around 40 mbp. Our assembly contains around 30 mbp.

MadsAlbertsen, how do you get such a plot? What program/command do you use?
martin_313 is offline   Reply With Quote
Old 01-23-2012, 12:52 AM   #5
MadsAlbertsen
Member
 
Location: Denmark

Join Date: Aug 2010
Posts: 26
Default

Most assemblers have a function to obtain the average coverage for each contig - it's a long time since I've used ABySS so I can't remember if it has an output contig coverage option.

Otherwise map the reads back to your contigs using any short read alligner and calculate the coverage.

You can easy do the coverage vs. contig length in excel. But if you have loads of contigs the smoothScatter function in R is great.

rgds
Mads
MadsAlbertsen 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 12:29 PM.


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