SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bowtie output fancy Bioinformatics 0 09-23-2011 07:26 AM
Bowtie --un output cascoamarillo Bioinformatics 0 06-07-2011 12:46 PM
bowtie output frymor Bioinformatics 2 01-12-2011 11:45 AM
bowtie output fargo Bioinformatics 4 02-19-2010 11:57 PM
Bowtie output pythonlovesbowtie Bioinformatics 2 11-20-2009 05:07 AM

Reply
 
Thread Tools
Old 10-02-2009, 07:39 AM   #1
statsteam
Member
 
Location: Californica

Join Date: Sep 2009
Posts: 19
Default Bowtie output

Hi all,

I got an RNA seq file that contains about 15million reads.
As it has more than 32bp reads, I decided to use Bowtie instead of my usual pick which is Eland.

So I mapped the sequence on the hg18. For the comparison purpose, I also ran Eland using first 32bp of the original read.

As usual, Eland output showed about 92% mapping rate but the Bowtie output contains only 5.1 million rows which is only about 34% mapping rate.

The options I used when running Bowtie was -f and -t.
The exact command line was:

bowtie -f -t bowtieIndexPath inputRNASequence outputFileName

The Bowtie output contiains all the sequences (5.1 million sequences) mapped with at most two mismatches.
Meanwhile, the Eland mapped 11.2 million sequences with at most two mismatches using first 32bp reads.
Since there were 15 million reads initially,
the mapping rate of Eland with at most two mismatches is 11.2/15 = 0.747
the mapping rate of Bowtie with at most two mismatches is 5.1/15 = 0.340

Is there any way to make Bowtie display the result of each read even though it does not have match?

Thank you,
Statsteam
statsteam is offline   Reply With Quote
Old 10-02-2009, 09:44 AM   #2
krobison
Senior Member
 
Location: Boston area

Join Date: Nov 2007
Posts: 747
Default

There is an option in bowtie to trim the sequences - -to get a really fair comparison you could try bowtie with that option set so it is looking only at the first 32 nt.

It would seem unlikely that splice junctions explain things, but you could use TopHat which is the bowtie-based spliced aligner
krobison is offline   Reply With Quote
Reply

Tags
bowtie, rna

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 04:29 AM.


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