SEQanswers

Go Back   SEQanswers > Core Facilities



Similar Threads
Thread Thread Starter Forum Replies Last Post
searching 454 (single read) + Illumina paired-end or matepair data boetsie General 0 10-04-2010 01:45 AM
Paired end Short read data SS1234 Bioinformatics 6 06-09-2010 01:16 PM
help! what is a paired-end read? hitdavid Bioinformatics 1 01-14-2010 11:42 AM
Difference in paired-end and single-end read ? darshan Bioinformatics 1 09-30-2009 11:44 PM

Reply
 
Thread Tools
Old 10-14-2011, 05:37 AM   #1
lisa1102
Junior Member
 
Location: UK

Join Date: Oct 2011
Posts: 4
Default Average Read Coverage for 454 paired end read data

Hi
I have just received some paired end read data from 454 flx sequencer. The coverage is 6X for three of my samples and 3X for one of my samples.
These samples are microbes with no reference strain.
Can anyone tell me if this coverage is the norm?
Thanks,
Lisa
lisa1102 is offline   Reply With Quote
Old 10-14-2011, 06:29 AM   #2
MissDNA
Senior Member
 
Location: Brazil

Join Date: Nov 2010
Posts: 146
Default

Coverage depends on the size of your genomes. Do you have any idea of sizes?
Also, what kind of PE was it?

In our last 3kb PE run, we sequenced 2 bacteria in a 2 lane PTP. The 1st one had a genome around 4.5M and the coverage was 56X. The other one is larger than 7M and we got 16X coverage. We knew the 2nd library was not very good but decided to sequence anyway.

So if your organisms were small, I think 6x and 3x coverage is not very good.
MissDNA is offline   Reply With Quote
Old 10-14-2011, 07:16 AM   #3
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,184
Default

Hi Lisa,
You will have to give us more to go on. Like how many and what type of regions were used on the 454 sequencer? Also what genome size are your microbes?

--
Phillip
pmiguel is offline   Reply With Quote
Old 10-18-2011, 08:10 AM   #4
vamosia
Member
 
Location: New York

Join Date: Mar 2009
Posts: 15
Default

On my previous microbes projects, 20X seems to be the depth that works well with Newbler. More depth will typically break the contig more since the errors accumulate and less depth does not generate enough signal.

But this is not absolute and it depends on your genome.
vamosia is offline   Reply With Quote
Old 10-18-2011, 08:25 AM   #5
lisa1102
Junior Member
 
Location: UK

Join Date: Oct 2011
Posts: 4
Default

Hi
Thanks for your comments. This is the only information I have at present. The size of the genomes is 2.7MB for the 6X coverage and 1.5MB for the 3X coverage.
If people think that 20X is average then is anything I can do to improve the analyses or is it best to repeat the runs?
lisa1102 is offline   Reply With Quote
Old 10-18-2011, 08:28 AM   #6
MissDNA
Senior Member
 
Location: Brazil

Join Date: Nov 2010
Posts: 146
Default

Don't you know how the PTP was divided: 2, 4, 8 or 16 regions? If you had a sample per region in a 2-lane PTP the coverage you got is way too low. How was the library prep?
MissDNA is offline   Reply With Quote
Old 10-18-2011, 08:36 AM   #7
lisa1102
Junior Member
 
Location: UK

Join Date: Oct 2011
Posts: 4
Default

I don't have that information of the sample per region at the moment. I am waiting for the information. The only information I have on the library prep is from the protocol "Paired end lib prep method manual - 20kb and 8kb span"
lisa1102 is offline   Reply With Quote
Old 10-18-2011, 08:39 AM   #8
MissDNA
Senior Member
 
Location: Brazil

Join Date: Nov 2010
Posts: 146
Default

I never worked with 8 and 20kb libraired. My experience is with 3kb.
My recomendation is you need to gather all the info: about quality of the library, number of reads generated, and how the PTP was divided before deciding what to do next.
MissDNA is offline   Reply With Quote
Old 10-18-2011, 08:40 AM   #9
lisa1102
Junior Member
 
Location: UK

Join Date: Oct 2011
Posts: 4
Default

Yes Thanks for your comments. This will help me pose more specific questions.
lisa1102 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 10:40 PM.


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