SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
hg19 genome reference for short read mapping yh253 Bioinformatics 4 12-29-2013 10:11 PM
Problem with BWA mapping second read pbluescript Bioinformatics 3 10-12-2011 08:52 AM
BWA behaviour with Mate Pair data + Multi read mapping apratap Bioinformatics 5 06-23-2011 05:38 PM
velvet assembler contigs into gbrowse Zimbobo Bioinformatics 0 04-15-2010 02:36 PM
number of contigs in velvet bioenvisage Bioinformatics 6 03-24-2010 09:10 PM

Reply
 
Thread Tools
Old 01-19-2010, 09:34 AM   #1
apratap
Member
 
Location: Bay Area

Join Date: Jan 2009
Posts: 58
Unhappy BWA and Velvet : Anamolies in contigs and read mapping to reference

Hi All

I believe I have an interesting problem to look at or may be its just me who is experiencing this for the first time.

I have reads from GA for a bacterial genome. I used BWA to map them to known reference. Only 4% reads map back. Seeing the low % I decided to de novo assemble the reads using Velvet .


1. Total of 2150 contigs
2. 601 align to reference using 52 % reads from the original set.
3. k-mer size I used was 29. (read length : 50 bp PE)

I am not sure what to infer from these very contrasting numbers according to me. With Velvet contigs I am able to use 52% of the reads that aligns well to reference whereas with BWA only 4% reads align.

Looking forward for your thoughts.

Thanks!
-Abhi

Last edited by apratap; 01-19-2010 at 09:35 AM. Reason: typo in subject line
apratap is offline   Reply With Quote
Old 01-29-2010, 01:44 PM   #2
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

What is the percent identity between the contig and the known reference?
lh3 is offline   Reply With Quote
Reply

Tags
bwa, velvet

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 12:48 PM.


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