SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 05-09-2013, 04:26 AM   #1
Inma
Member
 
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,

Inma

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

Join Date: Jun 2010
Posts: 34
Default

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, 06:34 PM   #3
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

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
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 06:10 PM.


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