SEQanswers

Go Back   SEQanswers > Introductions



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie Paired-end options days369 Bioinformatics 5 12-12-2012 05:38 PM
Multiple alignment why and what for (Bowtie -k options) pmcqt Bioinformatics 0 11-29-2011 06:48 AM
Setting Bowtie options from the Tophat command line GiladZil RNA Sequencing 2 08-02-2011 02:42 PM
Tophat options to report unaligned reads and controlling Bowtie options Siva Bioinformatics 0 10-15-2010 08:38 PM
Using BWA with illumina output, do I need to adjust score foxyg Bioinformatics 0 06-22-2010 08:51 AM

Reply
 
Thread Tools
Old 09-11-2013, 01:42 AM   #1
antoza
Member
 
Location: France

Join Date: Aug 2013
Posts: 18
Default bowtie adjust alignment options

Hi everyone,

I am also new. I have started to work as a PostDoc in a mutant resequencing project and i am using bowtie for mapping paired end reads against a reference assembly. I would be grateful if someone knows how exactly can i refine aligment options in bowtie in aspects of lenght of match (40%) and similarity levels (70%) for each read in order to relax mapping stringency. I am trying to figure out this and i think is probably something with -n/--seedmms and -e/--maqerr but i am not quite sure.
Any help would be so helpful.
Antonio
antoza is offline   Reply With Quote
Old 09-11-2013, 03:01 AM   #2
mastal
Senior Member
 
Location: uk

Join Date: Mar 2009
Posts: 667
Default

I am not very familiar with Bowtie2, but in Bowtie1 the parameters to change would be -n or -v, and -l.
mastal is offline   Reply With Quote
Old 09-11-2013, 03:47 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

In bowtie2, the parameters you might want to look at are:
  • -N number of seed mismatches, default is 0, but this can also be 1, though that will significantly slow things down
  • -L seed length, make it shorter if you want to increase sensitivity (see also the presets)
  • --score-min This effectively sets the minimum alignment score to declare an alignment valid. It's length dependent and the defaults change according to whether you use local or end-to-end alignment. In general, if you relax the --score-min function, you'll get more alignments with greater mismatches. Keep in mind that bowtie2's MAPQ calculation function depends heavily on this (the lower the minimum score is, the higher the MAPQ will be for a given read with a given number of mismatches).
  • -D Number of failing seed extensions, see also the presets
  • -R Number of repeat seeds, see also the presets
dpryan 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:03 PM.


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