SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
strange SAM output zhidong Bioinformatics 5 08-23-2011 05:10 AM
Strange/erroneous TopHat output with paired-end data hoho Bioinformatics 7 03-09-2011 12:21 PM
bowtie output frymor Bioinformatics 2 01-12-2011 11:45 AM
strange bowtie index building and mapping problem Gangcai Bioinformatics 0 08-04-2010 05:02 PM
bowtie output fargo Bioinformatics 4 02-19-2010 11:57 PM

Reply
 
Thread Tools
Old 08-17-2011, 07:51 AM   #1
vebaev
Senior Res.
 
Location: Plovdiv, Bulgaria

Join Date: Oct 2008
Posts: 110
Default strange bowtie output

Hi all,
I have used bowtie with the following options as option -a result some reads mapping thousands times:
bowtie -f -m 200 -v 1 -p 3

As I know -m 200 will remove a tag which maps more than 200 times, but if maps 199 times they will be in the output?

As a result from -m 200 my output include only one maping (one alignment for one read? Even reads I know that map 5 times are showed as only 1 alignment?

How is that possible
vebaev is offline   Reply With Quote
Old 08-17-2011, 01:33 PM   #2
dreesbl
Registered Vendor
 
Location: Seattle, WA

Join Date: Nov 2009
Posts: 6
Default

@vebaev

As you found out, option -a is usually trouble! To answer your question, -m 200 means that NO alignments are reported for a read that has 200+ hits. For a read with 1, 2, or 199 alignments, the default behavior of Bowtie is to report a single alignment for the read. You can set the max number of alignments to report for each read with the -k option (default is -k 1).

Another thing to watch out for is that the single alignment reported is not necessarily the best one. Unless you use the --best option Bowtie will report the first valid alignment it encounters. The --best option sets it to find all valid alignments and report them in best to worst order (only the best one at -k 1), which is much more computationally intensive. Many of the "options" in Bowtie are shortcuts made necessary by limited computing resources and don't necessarily improve the quality of your results.

As an aside, why are you using read data in fasta format? What instrument platform is it from? Just curious why you don't have fastq or qseq. Bowtie takes base quality score into account and leaving it out will affect your alignment results. You could try BWA, which doesn't use base quality in its calculations for alignment.

dreesbl is offline   Reply With Quote
Old 08-17-2011, 01:44 PM   #3
vebaev
Senior Res.
 
Location: Plovdiv, Bulgaria

Join Date: Oct 2008
Posts: 110
Default

Thanks for the reply.
I have notised the details that you described.

I succseed outputing all alignments up to 200 for each read and discarding ones with 200+ by combining "-a -m 200"

I'm using fasta since I have already trimmed, filtered, cleaned and unified reads with their counts in fasta.
vebaev 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 08:49 AM.


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