SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Emergent self-organizing maps from 4mer frequencies (binning metagenomic contigs) rhinoceros Bioinformatics 0 05-16-2014 03:35 AM
BAM index binning rpandya Bioinformatics 0 07-01-2013 08:09 AM
De novo amplicon reassembly? Elizabeth Bent Bioinformatics 3 07-19-2012 10:33 AM

Reply
 
Thread Tools
Old 05-04-2016, 07:58 AM   #1
jtamames
Junior Member
 
Location: Madrid

Join Date: May 2014
Posts: 7
Default shorter contigs after binning and reassembly

Hello friends, maybe you can help me with this.

I have a coassembly of two metagenomes, from which I can retrieve the contigs that belong to a particular taxon of interest. My criteria for assigning a taxon to a contig is rather strict, asking for 80% of the genes in the contig to be undoubtedly assigned to the taxon. In this way I get >500 contigs summing up to 5 Mb, with some big contigs >100 Kb. Then I map the original reads to these contigs using Bowtie2, and extract the mapping reads. When I reassemble this set of reads using Spades, to my surprise I get considerably shorter contigs. Inspection with Mauve or nucmer shows extensive splitting of the original contigs, exactly the opposite of my expectations.

Could you please hint me on this? Could it be that the mapping step is producing an unbalanced coverage and as a consequence the contigs split?

Thanks a lot!
jtamames is offline   Reply With Quote
Reply

Tags
assembly, binning, metagenome

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 06:20 PM.


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