SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about next step in a de novo bacterial genome assembly JodyFranke General 3 07-28-2017 03:55 PM
De novo Assembly k-mer question ronaldrcutler Bioinformatics 1 08-19-2016 05:11 PM
Compare de-novo transcriptome assembly to genome reference guided assembly IdoBar Bioinformatics 1 04-04-2014 01:28 AM
Abyss de novo assembly question hakattack Bioinformatics 4 05-09-2013 08:54 AM

Reply
 
Thread Tools
Old 12-06-2018, 02:56 PM   #1
tarias
Member
 
Location: Missouri

Join Date: Sep 2011
Posts: 10
Default De novo genome assembly question: which assembly should I use?

I have sequenced two individuals and done three libraries per individual each with a different insert size (170, 500, 800).Using flow citometry I know this genome is 332MB. Also I know is diploid and the genetic variation is not to high. I have done different assemblies for a this plant genome and found this two were the best:

One Individual
3 libraries
Insert size 170, 500, 800
Number of clean reads 70,287,710
Number of contigs 134,557
Mean N50 length (bp) 5,050
BUSCO 78.54


Two Individuals
2 individuals
6 libraries
Insert size 170, 500, 800
Number of clean reads 209,022,970
Number of contigs 446,608
Mean N50 length (bp) 4,010.65
BUSCO 76.38

I then ran PEP_SCAFFOLDER in this two alignments and got this

One Individual
3 libraries


Main genome scaffold total: 136339
Main genome scaffold sequence total: 189.946 MB
Main genome contig N/L50: 17091/2.413 KB

Two Individuals
2 individuals
6 libraries


Main genome scaffold total: 444336
Main genome scaffold sequence total: 332.652 MB
Main genome contig N/L50: 83374/706

I have two questions:
1. Would you annotate any of these two genomes?
2. Which one would you annotate and why?

Thanks,
Tatiana
tarias is offline   Reply With Quote
Old 12-07-2018, 05:23 AM   #2
ctseto
Member
 
Location: SE MN

Join Date: Oct 2013
Posts: 44
Default

Annotation may reveal what's not there, esp if you know from exp that a given gene is present (eg someone PCR/subcloned ). Guessing eukaryote and that you may have hit some contigs that are repeat heavy (eg contig is pure repeat), and perhaps that your library isn't traversing those repeats?
ctseto 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 09:07 AM.


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