SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Is bowtie2 really this unstable? kgulukota Bioinformatics 5 02-22-2012 08:27 AM
Bowtie2 --al and --un: no output arnav Bioinformatics 3 02-14-2012 01:08 PM
Bowtie2 and 2 GB limit mscholz Bioinformatics 1 01-03-2012 01:08 PM
bowtie2 fun mscholz Bioinformatics 3 11-18-2011 07:30 PM
Bowtie2 with Tophat plassaaw Bioinformatics 4 11-11-2011 06:12 AM

Reply
 
Thread Tools
Old 03-01-2012, 06:53 AM   #1
LP_SEP23
Junior Member
 
Location: Philadelphia

Join Date: Feb 2012
Posts: 9
Default Anyone has experience with Bowtie2?

Hi All,
my name is Lucia, and I've been lurking around for a while but this is my first post.
I have a pretty big ChIP-seq experiment to analyze, and initially did the alignments with Bowtie. The stats were not that great, and since I have 100bp reads, I decided to redo the alignments with Bowtie2
I do get a lot more mapped reads, but the quality scores seem a lot lower
The programs compute the scores differently, so I am not sure if the scores are really "lower".
Does Bowtie2 increase mapping at the expense of quality?
Or is the scoring scheme so completely different that the quality scores are not comparable at all?
thanks
LP_SEP23 is offline   Reply With Quote
Old 03-01-2012, 07:05 AM   #2
epi
Member
 
Location: USA

Join Date: Jan 2012
Posts: 38
Default

Sorry can't help you with your specific problem.

Though in my experience, bowtie works much better than bowtie2. Apart from some bugs, many options are missing (strata/unique alignments). I shared my opinion with authors as well.
Must mention that I have 50bp reads and my final goal (ChIP-Seq) may be different from you.
epi is offline   Reply With Quote
Old 03-01-2012, 07:08 AM   #3
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

Maybe try bwa
NicoBxl is offline   Reply With Quote
Old 12-05-2012, 03:42 PM   #4
sdriscoll
I like code
 
Location: San Diego, CA, USA

Join Date: Sep 2009
Posts: 438
Default

Depends on what score you are talking about. The MAPQ score from bowtie2 has to do with the uniquness of the alignment and not necessarily the quality of the alignment. You want to look at the AS attribute in the extended attributes of the SAM files. Also read up on this section of the bowtie2 manual:

http://bowtie-bio.sourceforge.net/bo...r-more-similar

By default bowtie2 is allowing no mismatches in the seed (22 bases for default end-to-end alignment) and a minimum alignment score of -60.6 (with a max of 0) for 100bp reads. Default penalties are something like -6 for a mismatch at a high quality base an -11 for a length-2 gap. Gaps in reads and reference are scored the same. Specifically you get a -5 for a gap open and -3 for each base. So a length-1 gap is -8, length-2 is -11, etc. So with a min score of about -60 you could have as many as 6 mismatches past the seed or all kinds of combinations of gaps and mismatches. You can control all of those parameters. The min score setting is defined as -0.6 + -0.6*L where L is you read length. If you want to set that value you do so with --min-score L,-0.6,-0.6. So you can figure out what you want your min score to be, like say -40, and work backward from that to figure out what values to put in the setting.
sdriscoll is offline   Reply With Quote
Old 12-05-2012, 05:03 PM   #5
adaptivegenome
Super Moderator
 
Location: US

Join Date: Nov 2009
Posts: 437
Default

Mapping quality is valuable for variant calling. If the quality score does not correspond to accuracy of alignment, this is a big problem I think.
adaptivegenome is offline   Reply With Quote
Old 12-05-2012, 10:34 PM   #6
sdriscoll
I like code
 
Location: San Diego, CA, USA

Join Date: Sep 2009
Posts: 438
Default

Quote:
Originally Posted by adaptivegenome View Post
Mapping quality is valuable for variant calling. If the quality score does not correspond to accuracy of alignment, this is a big problem I think.
You're right about that. That's probably why it's better to use BWA for those alignments. I don't know of any advantage other than speed for using bowtie2 over BWA. I haven't done any hardcore comparisons between the two. I do like BWA's alignments more though...not sure why.
sdriscoll 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:40 AM.


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