Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
PCA interpretation gene_x Bioinformatics 0 03-13-2013 12:18 PM
interpretation of output from HTSeq syintel87 Bioinformatics 1 01-06-2013 02:12 PM
Velvet Interpretation nxtgenkid10 De novo discovery 2 02-24-2012 02:43 AM
Help interpretation indel labmtl Illumina/Solexa 0 03-29-2011 08:01 AM
Interpretation of BWA mapping quality christophpale Bioinformatics 0 07-21-2010 04:15 AM

Thread Tools
Old 05-09-2013, 05:26 AM   #1
Location: Spain

Join Date: Feb 2013
Posts: 25
Default BWA interpretation

Hi all,

I have seen in BWA that several reads stop together in the same base. Do you know why it happens?

You can see like a "wall" where the sequencing can't go on and only a few reads cross this "wall".

My genome is a bacterial genome.

Thank you very much,


Last edited by Inma; 05-09-2013 at 07:26 AM.
Inma is offline   Reply With Quote
Old 05-09-2013, 06:41 PM   #2
Location: Sydney, Australia

Join Date: Jun 2010
Posts: 34

Hi Inma,

I too have seen this. BWA can trim reads with low quality, and if there is a difficult-to-sequence region, such as a homopolymer run, then many reads will have low quality at that position and get trimmed.

I think that is what may cause the 'wall' effect.

Either that, or you have not removed duplicate reads?
rbagnall is offline   Reply With Quote
Old 05-09-2013, 07:34 PM   #3
Senior Member
Location: San Diego

Join Date: May 2008
Posts: 912

I doubt the problem is bwa. I bet if you aligned with bowtie, you'd see the same thing.

There are a lot of possibilities. Maybe your reference isn't right. Or maybe your library really looks like that. There are probably certain areas of the genome that fragment poorly.

Are you doing any quality trimming of the reads?
swbarnes2 is offline   Reply With Quote

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 09:26 PM.

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