SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie output to BED format polsum Bioinformatics 11 02-19-2013 11:00 PM
Bowtie/Tophat SAM format detail question traeki Bioinformatics 0 06-12-2011 03:03 PM
Bowtie output format kevlim83 Bioinformatics 5 01-26-2010 01:28 PM
bowtie output file format Ally Bioinformatics 4 08-01-2009 08:40 PM
Bowtie output file format questioner Bioinformatics 0 07-02-2009 12:06 PM

Reply
 
Thread Tools
Old 12-15-2010, 09:00 AM   #1
fahim
Member
 
Location: louisville

Join Date: Oct 2010
Posts: 10
Question Bowtie output format question

Hi
I recently downloaded Bowtie and tying to align affy sequences. I am still using BLAT for longer sequences, but for short sequences, Bowtie is having clear advantage.

I have a doubt about the default output of bowtie. Field number FOUR(4) gives 0-based offset of the start location of reads on the genome. But it is not showing the location where mapping ends. For example, the BLAT output tells you everything including start, end , length of reads matched, mismatched etc. etc.

What I am assuming that since bowtie does not allow gapped alignment, the end should be equal to ( start+ lenghtOfReads-1). Is it true for the alignment results?
Thanks
Fahim
fahim is offline   Reply With Quote
Old 12-16-2010, 12:22 AM   #2
xinwu
Member
 
Location: Beijing

Join Date: Jul 2010
Posts: 33
Default

Quote:
Originally Posted by fahim View Post
Hi
I recently downloaded Bowtie and tying to align affy sequences. I am still using BLAT for longer sequences, but for short sequences, Bowtie is having clear advantage.

I have a doubt about the default output of bowtie. Field number FOUR(4) gives 0-based offset of the start location of reads on the genome. But it is not showing the location where mapping ends. For example, the BLAT output tells you everything including start, end , length of reads matched, mismatched etc. etc.

What I am assuming that since bowtie does not allow gapped alignment, the end should be equal to ( start+ lenghtOfReads-1). Is it true for the alignment results?
Thanks
Fahim
You have the answer and should consider the strand.
xinwu is offline   Reply With Quote
Old 10-26-2012, 11:21 AM   #3
kjlee
Member
 
Location: Atlanta, GA

Join Date: Jun 2011
Posts: 12
Default

From the bowtie manual:

Quote:
0-based offset into the forward reference strand where leftmost character of the alignment occurs.
Leftmost of which sequence: the reference or the query? I believe that it does matter. I presume it is the leftmost of the reference. If this is the case, then the 5' end of the query would be at (offset+length-1).
kjlee is offline   Reply With Quote
Old 10-27-2012, 02:25 PM   #4
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quote:
Originally Posted by kjlee View Post
From the bowtie manual:



Leftmost of which sequence: the reference or the query? I believe that it does matter. I presume it is the leftmost of the reference. If this is the case, then the 5' end of the query would be at (offset+length-1).
The reference. But check the CIGAR string to see if the read has been hard or soft trimmed.
swbarnes2 is offline   Reply With Quote
Old 10-29-2012, 05:41 AM   #5
kjlee
Member
 
Location: Atlanta, GA

Join Date: Jun 2011
Posts: 12
Default

Say the query is 36nts long. The sequence aligns to the negative strand and the last (most 3') nucleotides (let's say 5nts) are soft-clipped, will the position reported be the first position that has a properly-mapped (non-soft-clipped) nucleotide, or will it be the 3' most nucleotide of the 36?
kjlee is offline   Reply With Quote
Reply

Tags
bowtie

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:31 PM.


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