How to detect real deletion or gaps in sequencing projects

watashi · 08-01-2008, 02:15 PM

Dear all

We have 3 genomes sequenced of 11, 15 and 23 fold coverage by 454 Sequencer. At least more than 10kb of the gap regions (defined by unsuccessful reference mapping) are coding regions containing little repetitive sequences. Interestingly, these gap regions appear quite concordantly among the 3 genomes.

Would anybody know how to utilize the scaffold files or related in order to check whether they are real gaps or not? Automation cannot be done if the process has to be studied by manually comparing the results generated by reference mapper and de novo assembler.

_________ ................ ________
|..............|................|..........|
|.ATTTCC..|---------->| CGCCC |
|_________|................|______|

Say, in reference genome, it is:
ATTTCCTTAGGAACGCCC

can we have a quick way to identify ATTTCCCGCCC (not to do that manually…) so to confirm the genome has deletion of (TTAGGAA) or not?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
large sequencing projects	barak	Service Providers	12	11-11-2013 01:55 PM
large sequencing projects	barak	Bioinformatics	1	11-11-2013 01:47 PM
detecting deletion in amplicon sequencing of mixed population	Josim	Genomic Resequencing	0	11-10-2011 06:24 AM
Designing de-novo sequencing projects	isharon	De novo discovery	0	05-25-2010 01:01 AM
PubMed: Brazilian genome sequencing projects: state of the art.	Newsbot!	Literature Watch	0	12-17-2008 06:03 AM

08-01-2008, 02:15 PM	#1
watashi Junior Member Location: cn Join Date: Aug 2008 Posts: 1	How to detect real deletion or gaps in sequencing projects Dear all We have 3 genomes sequenced of 11, 15 and 23 fold coverage by 454 Sequencer. At least more than 10kb of the gap regions (defined by unsuccessful reference mapping) are coding regions containing little repetitive sequences. Interestingly, these gap regions appear quite concordantly among the 3 genomes. Would anybody know how to utilize the scaffold files or related in order to check whether they are real gaps or not? Automation cannot be done if the process has to be studied by manually comparing the results generated by reference mapper and de novo assembler. _________ ................ ________ \|..............\|................\|..........\| \|.ATTTCC..\|---------->\| CGCCC \| \|_________\|................\|______\| Say, in reference genome, it is: ATTTCCTTAGGAACGCCC can we have a quick way to identify ATTTCCCGCCC (not to do that manually…) so to confirm the genome has deletion of (TTAGGAA) or not? Last edited by watashi; 08-01-2008 at 02:19 PM.