SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Few contigs map to reads a532942 Bioinformatics 1 07-02-2014 08:07 AM
Too little reads mapped back to contigs using Bowtie oliviaxinw Bioinformatics 4 09-28-2012 05:18 AM
Map transcripts from de novo assemblier such as trinity back to the genome ? sptmbr Bioinformatics 5 02-29-2012 09:47 AM
Too few reads mapping back to contigs cerebralrust Bioinformatics 8 02-26-2012 10:25 PM
How to map 454 reads/contigs to a mitochondrial genome? fruktimport Bioinformatics 2 03-28-2011 10:35 AM

Reply
 
Thread Tools
Old 07-24-2014, 07:02 AM   #1
vanillasky
Member
 
Location: Europe

Join Date: Mar 2014
Posts: 42
Default CLC-Reads do not map back to contigs

I have been using CLC genomics to try and assemble paired end reads from an illumina run (this is an environmental metagenome sample). The assembly is set up so that at the end, the reads are mapped back to the contigs to see how many matched and how many did not match the contigs. I see from our output that a large number of reads do not match the contigs. Is this indicative of misassembly or what could this indicate? Any input would be most appreciated.

p.s. the trimming was done with an adapter input file to remove any of these present in the sequences. We see only a small percentage having undergone adapter removal following the trimming process. The majority of sequences were just trimmed based on quality score.
vanillasky is offline   Reply With Quote
Old 07-24-2014, 07:53 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,076
Default

This would be a good candidate to contact CLC Tech support about. They should be able to remote in and see what is going on. I doubt we can give you any good advice here (as far as CLC goes).
GenoMax is offline   Reply With Quote
Old 07-24-2014, 10:36 AM   #3
vanillasky
Member
 
Location: Europe

Join Date: Mar 2014
Posts: 42
Default

Hi GenoMax,

I actually obtain similar results using other assembly programs such as velvet for example. Why would a high number of reads in general not match up with a set of contigs post assembly?
vanillasky is offline   Reply With Quote
Old 07-24-2014, 10:43 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,076
Default

Was that with metavelvet (http://metavelvet.dna.bio.keio.ac.jp/) or regular velvet? Since this is a metagenomic sample you should use an assembler that is designed for that purpose. Ray-meta, MetAMOS look like other options.
GenoMax is offline   Reply With Quote
Old 07-24-2014, 11:54 AM   #5
yzzhang
Member
 
Location: florida

Join Date: Jan 2013
Posts: 67
Default

for metagenome assembly, is it normal that a large percentage of reads can not be assembled to contigs, right?
yzzhang 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 11:31 PM.


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