SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
question about denovo assembly kenietz Bioinformatics 26 05-13-2013 06:12 PM
Denovo assembly problem huma Asif Illumina/Solexa 1 03-27-2013 10:20 PM
denovo assembly nagaraj Bioinformatics 5 07-11-2012 07:13 AM
SOAP denovo assembly results bioinfosm Bioinformatics 9 06-13-2012 08:22 AM
DeNovo Transcriptome Assembly (human) Uwe Appelt Bioinformatics 0 10-28-2009 02:42 AM

Reply
 
Thread Tools
Old 12-17-2012, 08:46 AM   #1
asunen
Junior Member
 
Location: Michigan

Join Date: Dec 2012
Posts: 2
Default single gene denovo assembly?

Hello,

I'm looking at population level variation in a 100kb region in Drosophila. When I use standard alignment programs (bwa and bowtie) on reads from some sequenced populations, I keep losing insertions that I know are present (and appear in the raw reads), but are missing in the reference genome.

Someone suggested I do a denovo assembly on my 100kb region (since I'm not really interested in the rest of the genome), to get around this problem.

Is there a way to filter reads for a specific region using denovo assembly? Is this a good approach to the problem, or am I better off tweaking my methods with alignment to the reference genome?

Thanks!
asunen is offline   Reply With Quote
Old 12-17-2012, 03:38 PM   #2
jgibbons1
Senior Member
 
Location: Worcester, MA

Join Date: Oct 2009
Posts: 133
Default

I've actually done something similar. I had whole genome data but for a particular piece of analysis I was only interested in a highly variable 35kb region. I also knew that I was missing indels in mapping results.

I ended up compiling a database with this 35kb region (plus 5kb flanking on each side) from all of the sequenced genomes/genes in my region I could find. I then blasted my reads against this database. Reads which significantly hit my region were extracted and then independently de novo assembled. It worked very well.

Hope this helps.

John
jgibbons1 is offline   Reply With Quote
Old 12-18-2012, 08:21 AM   #3
asunen
Junior Member
 
Location: Michigan

Join Date: Dec 2012
Posts: 2
Default

Thank you, I'll try that
asunen is offline   Reply With Quote
Old 12-18-2012, 11:20 AM   #4
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

The only issue is, how do you know which reads are in your region? If your sample has an indel as compared to the reference, the reads might not map there, so you won't know to include in the de novo. You'd have to denovo everything, or do another experiment where you long PCR around your region of interest.
swbarnes2 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 01:35 AM.


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