SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Annotate contigs with BLAST hit names; remove contigs with no hit Bueller_007 Bioinformatics 10 02-27-2013 10:22 AM
Help:scaffold 2 contigs! rogerholmes.novogene Bioinformatics 11 07-28-2011 10:05 PM
Strange contigs Ramet RNA Sequencing 0 06-08-2011 04:53 AM
gsMapper contigs haonmada 454 Pyrosequencing 1 01-22-2010 11:25 AM
blast of the contigs dina Bioinformatics 5 11-02-2009 09:44 AM

Reply
 
Thread Tools
Old 11-13-2013, 02:23 AM   #1
vpi
Member
 
Location: Sweden

Join Date: Oct 2013
Posts: 18
Default SPAdes contigs

Hi,

I am new to SPAdes and I have some questions. I wonder about reads, like how many reads make a contig. I donít see any information in the SPAdes.log file. Is there a way to see the number of reads for each contig?
vpi is offline   Reply With Quote
Old 11-13-2013, 09:03 AM   #2
ctseto
Member
 
Location: SE MN

Join Date: Oct 2013
Posts: 44
Default

Velvet used to store read information differently; which then allowed extraction of that information via a perl script.

It may be worth mapping your original reads back to your contigs anyways-to check for changes in coverage, or to see if there are spurious gaps. Programs like Opera and REAPr rely on this mapping of the starting reads to the finished product as part of their quality control.
ctseto is offline   Reply With Quote
Old 10-19-2015, 12:41 PM   #3
Mario Moreno
Junior Member
 
Location: Chile

Join Date: Nov 2013
Posts: 7
Default

Hello all!
I am new to SPAdes too,
I have the contigs number but I dont have info about the number of read used for make the contigs. there any way to know this information without using external programs?
Thanks
Mario Moreno is offline   Reply With Quote
Old 12-02-2015, 03:49 PM   #4
Dmitry.Antipov
Junior Member
 
Location: St. Petersburg, Russia

Join Date: Feb 2013
Posts: 2
Default

Hi.
SPAdes do not output the exact number of reads mapped to the contig. However you can estimate it with help of contig length and average contig coverage (NB, it is kmer coverage for the last k used for assembly) that can be found in contig's fasta header.
Dmitry.Antipov 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 07:51 PM.


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