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

Thread Tools
Old 11-14-2018, 12:58 AM   #1
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
Senior Member
Location: UK

Join Date: Jan 2010
Posts: 390

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
Junior Member
Location: Western Australia

Join Date: Sep 2016
Posts: 5

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

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 03:43 AM.

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