Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie Paired end 100% failed to align rahul.m.dhodapkar Bioinformatics 15 09-12-2013 03:05 PM
Bowtie and reads that failed to align: (100.00%) michy Bioinformatics 7 02-08-2011 07:42 PM
soap2 PE results do not match bowtie or novoalign axiom7 Bioinformatics 0 01-26-2011 01:43 PM
454IsotigsLayout.txt - contigs don't match to any isotig? dschika Bioinformatics 9 10-11-2010 05:58 AM
BOWTIE - match with empty chr and 0 position dukevn Bioinformatics 6 12-01-2009 12:22 PM

Thread Tools
Old 02-05-2012, 08:37 PM   #1
Junior Member
Location: Caltech

Join Date: May 2011
Posts: 9
Default Can bowtie align mates individually if they don't match?


I am trying to use Bowtie for paired end read alignment, using SAM as the output format. If mates align within a specified insert size on the same chromosome, then alignments of two mates are reported with flags 99 and 147, or 163 and 83.

I was wondering if there is way to make bowtie to report alignments of mates even if they don't match (i.e. they don't fall into the insert size, potentially on different chromosomes)? I see that there are flag combinations that could allow such report in SAM, but I wasn't able to find settings that could make Bowtie do that.

"Letís start with the three fundamental Rules of Robotics...."
history_of_robots is offline   Reply With Quote
Old 02-06-2012, 06:01 AM   #2
Senior Member
Location: USA, Midwest

Join Date: May 2008
Posts: 1,172

Better handling of discordant read pairs is a feature of Bowtie 2. From the FAQ:

Bowtie 2's paired-end alignment mode is more flexible than Bowtie 1's. For example, for pairs that do not align in a paired fashion, it will attempt to find unpaired alignments for each mate.
kmcarr 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 05:45 AM.

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