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 09:11 PM
Problem with BWA mapping second read pbluescript Bioinformatics 3 10-12-2011 07:52 AM
BWA behaviour with Mate Pair data + Multi read mapping apratap Bioinformatics 5 06-23-2011 04:38 PM
velvet assembler contigs into gbrowse Zimbobo Bioinformatics 0 04-15-2010 01:36 PM
number of contigs in velvet bioenvisage Bioinformatics 6 03-24-2010 08:10 PM

Thread Tools
Old 01-19-2010, 08:34 AM   #1
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.


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

Join Date: Feb 2008
Posts: 693

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

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 01:51 AM.

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