Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa perfect match only guil Bioinformatics 4 08-30-2013 09:18 AM
exact match with bowtie2 papori General 0 06-04-2013 02:44 AM
Bowtie not mapping (even) perfect matches ? yvan.wenger Bioinformatics 0 02-05-2013 05:17 AM
perfect alignment match NicoBxl Bioinformatics 3 10-07-2010 04:28 AM
Perfect match disagreement between bowtie and BLAT on human genome danielsbrewer Bioinformatics 3 03-10-2009 09:10 PM

Thread Tools
Old 11-05-2013, 12:15 AM   #1
Location: china

Join Date: Aug 2013
Posts: 14
Default perfect match with bowtie/bowtie2

Hi guys,

I have clean small RNA sequencing data at my hands. I wanna perform a perfect-match of miRNA to hairpin data from miRBase with Bowtie. However, i'm confused with its two mapping modes -n & -v. Can anyone tell me the difference of these two modes? Which option fits my case best? why?

Here is my code, for now.
bowtie -p 2 -v 0 hairpin -f sample_1.fa sample_1.bout
The reason i choose -v mode is that i donnt understand the Maq strategy which -n mode adopt, thus i donnt know best options <-n -l -e> settings in my case.
jason_ARGONAUTE is offline   Reply With Quote
Old 11-05-2013, 01:06 AM   #2
Location: San Francisco, CA

Join Date: Mar 2012
Posts: 21

-n mode only considers mismatches in part of the read sequence (in the seed). Also, it takes the base quality scores into account when deciding what is the best alignment.
(See point 2. here)

-v 0 is right for finding exact matches.
bw. is offline   Reply With Quote
Old 11-05-2013, 02:28 AM   #3
Devon Ryan
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480

With bowtie2, you can either just use "grep AS:i:0" to get only perfect matches or just use the "--score-min 'C,0,-1'" option, which should do effectively the same thing.
dpryan is offline   Reply With Quote

-n & -v option, bowtie perfect match

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 03:04 PM.

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