SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Genome Res De novo bacterial genome sequencing: millions of very short reads assembly b_seite Literature Watch 1 10-05-2017 12:26 AM
Question about next step in a de novo bacterial genome assembly JodyFranke General 3 07-28-2017 03:55 PM
De novo genome assembly (Hybrid) - bacterial genome MuG General 2 03-28-2017 05:13 PM
Automated pipeline for de novo Bacterial genome assembly Morgane_AUS Bioinformatics 0 10-22-2014 06:25 PM
DNASTAR’s De Novo Bacterial Genome Assembly App Now Available on BaseSpace DNASTAR Vendor Forum 0 11-13-2012 09:22 AM

Reply
 
Thread Tools
Old 11-14-2018, 12:58 AM   #1
SDH
Junior Member
 
Location: Western Australia

Join Date: Sep 2016
Posts: 5
Default Issues with bacterial genome assembly (mapping vs. de novo)

I am attempting to assemble a genome from an organism I am working on. Initial BLAST analysis of some assembled contigs has led me to believe I am dealing with previously described bacterium “X”. Read mapping in geneious supports this as the majority of my reads are used and I get good coverage (bar a few regions of no coverage).

My de novo assembly however is giving me mixed messages…Though the assembled contigs have a high ANI, >95% when compare to the reference genome, the contigs will not map to the reference. I have a handful of contigs between 300 – 600 kb however, none will map back to the reference. I initially thought this was due to a poor assembly however my read quality is good and I can more-or-less (give or take a few base pairs) replicate the large contigs across two different assemblers (Spades and Abyss) and two different data sets from illumina and ion torrent reads generated from the same organism. DNA was extracted from a known pure culture so it is unlikely to be due to contaminating reads.

Has anyone come across this before or have any idea of where I might have gone wrong?
SDH is offline   Reply With Quote
Old 11-16-2018, 04:39 AM   #2
Bukowski
Senior Member
 
Location: Aberdeen, Scotland

Join Date: Jan 2010
Posts: 362
Default

What are you mapping the contigs back to the reference with?

You could perhaps try reordering the contigs with progressive Mauve against the 'reference' genome, which might not be that good a reference though..
Bukowski is offline   Reply With Quote
Old 11-19-2018, 12:13 AM   #3
SDH
Junior Member
 
Location: Western Australia

Join Date: Sep 2016
Posts: 5
Default

Good point - I was lazy and just ran the mapping within geneious though I did adjust the thresholds.

I will try bbmap - any other suggestions for bacterial genomes?

The data was trimmed/filtered prior to assembly and the fact I get very similar large contigs from two different assemblers I would like to think they are reasonably accurate.
SDH is offline   Reply With Quote
Reply

Tags
bacterial genome assembly, de novo, read mapping

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 11:10 AM.


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