Bowtie output

statsteam · 10-02-2009, 07:39 AM

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

krobison · 10-02-2009, 09:44 AM

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

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

10-02-2009, 07:39 AM	#1
statsteam Member Location: Californica Join Date: Sep 2009 Posts: 19	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

10-02-2009, 09:44 AM	#2
krobison Senior Member Location: Boston area Join Date: Nov 2007 Posts: 747	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