SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Inquiry: minimum length of reads for referece-based assembly or de novo assembly sunfuhui Bioinformatics 1 10-04-2013 09:28 AM
CASIM: Assembly Metrics for Assembly quality casim UK - Cambridge 2 09-09-2013 01:34 AM
How I find not assembly read in a reference assembly??? matiasfreired Bioinformatics 1 04-05-2012 12:13 PM
How much rRNA sequence is reasonable in an RNA-seq dataset? Nigel Saunders General 6 02-07-2012 04:53 AM

Reply
 
Thread Tools
Old 06-02-2013, 02:53 PM   #1
lkral
Member
 
Location: Carrollton, GA

Join Date: May 2011
Posts: 27
Default Is this assembly reasonable?

I have completed some de novo assemblies of MiSeq generated 250PE reads with Minia and SSPACE. I've posted a description of the outcomes and some quality checks at the following site:

http://www.dartergenomics.org/tallapoosa-darter-genome

For those of you with experience in assemblies, could you please look this over and let me know if the scaffold distributions I ended up with are reasonable given the MiSeq data?

Thanks.
lkral is offline   Reply With Quote
Old 06-03-2013, 12:55 PM   #2
rchikhi
Member
 
Location: France

Join Date: Jan 2013
Posts: 13
Default

As an author of Minia, I would say that it looks good, but I would be a bit biased Thanks for using the software.

The k=73 min_abundance=2 assembly appears to be better than the k=63 min_abundance=3, consistently over all metrics. Thus you could focus on the k=73 one.

I have not tested them, but two tools can help you further evaluate accuracy metrics of an assembly without a reference:
https://github.com/vezzi/FRC_align
http://www.sanger.ac.uk/resources/software/reapr/

Also you mentioned that there are ~53M paired reads, yielding 13x coverage. Does each of the 53M paired reads consist of two mates of length 250bp? if so the coverage should be multiplied by two.
rchikhi is offline   Reply With Quote
Old 06-04-2013, 03:56 AM   #3
lkral
Member
 
Location: Carrollton, GA

Join Date: May 2011
Posts: 27
Default

Thanks for the references.

And to be clear, there are 53 million individual reads.
lkral 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 01:37 AM.


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