Help regg. tools for comparing two sets of contigs from de-novo assembly

jsreddy82 · 02-11-2013, 02:41 PM

Hello

I have two sets of fastq files (control & treatment) from DNA-sequencing (illumina) of a bacterial strain. The strain is new and has no (very) closely related strains. I have used the CLC Genomics Workbench to build contigs de-novo and then used Sequencher to improve the quality and length of contigs generated from the CLC Workbench. I have to compare these two sets of contigs (fasta sequences) to search for similarities and differences. The treatment was an insertion in the genome and I wish to study how this affects the genome and/or genes. Since I don't have a way to perform a genome wide study, is there a way to perform Control_contigs vs Treatment_contigs pairwise search to study the difference and how easy or hard would it be to parse the result files from the tools suggested?
Statistical info: Size of contigs range from 128Kb to 100bases in both sets, majority of them between 5K-15Kbases. Each set has around 140 contigs. Majority of contigs have a 40% G/C content. Some of the contigs have N's. Average length of contigs:12Kbases.
PS. I have tried to use CD-HIT-EST-2D and standalone blast+ for which the results have been inconclusive or jobs failed.

Thanks,

Joel_91 · 03-04-2014, 06:13 AM

Hello,

I am interested in doing the same comparison between two contigs files as you are, did you solve it? CD-HIT-EST-2D is not working for me because I get an error of segmentation fault ('core' dumped), any help?

Thanks

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Comparing sets of contigs	LizBent	Bioinformatics	1	05-17-2012 02:39 AM
Comparative studies of de novo assembly tools for next-generation sequencing technolo	strob	Literature Watch	2	07-25-2011 03:58 PM
PubMed: Local De Novo Assembly of RAD Paired-End Contigs Using Short Sequencing Reads	Newsbot!	Literature Watch	0	05-05-2011 11:40 PM
How to allow mapped contigs to be used in de novo assembly projects	edleloth	454 Pyrosequencing	1	03-08-2011 07:13 AM
De-novo assembly of bacteria genomes - which tools?	maasha	Bioinformatics	8	08-12-2009 01:41 PM

02-11-2013, 02:41 PM	#1
jsreddy82 Junior Member Location: Starkville, MS Join Date: May 2011 Posts: 1	Help regg. tools for comparing two sets of contigs from de-novo assembly Hello I have two sets of fastq files (control & treatment) from DNA-sequencing (illumina) of a bacterial strain. The strain is new and has no (very) closely related strains. I have used the CLC Genomics Workbench to build contigs de-novo and then used Sequencher to improve the quality and length of contigs generated from the CLC Workbench. I have to compare these two sets of contigs (fasta sequences) to search for similarities and differences. The treatment was an insertion in the genome and I wish to study how this affects the genome and/or genes. Since I don't have a way to perform a genome wide study, is there a way to perform Control_contigs vs Treatment_contigs pairwise search to study the difference and how easy or hard would it be to parse the result files from the tools suggested? Statistical info: Size of contigs range from 128Kb to 100bases in both sets, majority of them between 5K-15Kbases. Each set has around 140 contigs. Majority of contigs have a 40% G/C content. Some of the contigs have N's. Average length of contigs:12Kbases. PS. I have tried to use CD-HIT-EST-2D and standalone blast+ for which the results have been inconclusive or jobs failed. Thanks, __________________ Joseph Last edited by jsreddy82; 02-12-2013 at 06:32 AM. Reason: Re-post

03-04-2014, 06:13 AM	#2
Joel_91 Junior Member Location: Barcelona Join Date: Oct 2013 Posts: 2	Hello, I am interested in doing the same comparison between two contigs files as you are, did you solve it? CD-HIT-EST-2D is not working for me because I get an error of segmentation fault ('core' dumped), any help? Thanks