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 04:35 AM
BAM index binning rpandya Bioinformatics 0 07-01-2013 09:09 AM
De novo amplicon reassembly? Elizabeth Bent Bioinformatics 3 07-19-2012 11:33 AM

Thread Tools
Old 05-04-2016, 08:58 AM   #1
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

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:42 PM.

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