SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie 2 explanation of the pair-end mode report wariobrega Bioinformatics 12 Yesterday 10:57 PM
bowtie: extra parameter error feralBiologist Bioinformatics 0 11-21-2013 05:05 AM
Bowtie -e (--maqerr) parameter nemesis Bioinformatics 0 09-29-2011 08:47 AM
bowtie parameter setting nunu_ping RNA Sequencing 0 09-28-2010 01:59 PM
question about bowtie -e parameter ranel Bioinformatics 3 08-17-2010 07:55 AM

Reply
 
Thread Tools
Old 01-03-2018, 01:01 AM   #1
Gorgon_VZ
Member
 
Location: Germany

Join Date: Oct 2016
Posts: 10
Default bowtie 2 dpad parameter explanation

Can someone provide a detailed explanation, what the bowtie2 parameter "--dpad" does? The explanation given in the manual is as follows: "--dpad <int> include <int> extra ref chars on sides of DP table (15)".

I would interprete it like this:
1) A sequence string of length <int> is attached to every 5' and 3' end of a read after seeding and global alignment.
2) Attached sequences are used to solve mapping issues at readends (e.g. deletions present at readends that have been falsley softclipped in previous steps).

!!!But I am really not confident with this!!!

I am asking because I realised that increasing dpad to an integer of 100, solves issues I have to face with large deletions in a merged 350bp Illumina amplicon dataset. I mean solving this issue is fine, but I do not have a clue what the possible drawbacks of this setting might be. Furthermore, as I do not exactly understand the parameter and its behaviour, I do not know how to test for possible side-effects.

Here are my option settings:

bowtie2 --rdg 5,1 --rfg 5,1 --dpad 100 --score-min L,0,-0.75 --very-sensitive
Gorgon_VZ is offline   Reply With Quote
Old 01-03-2018, 03:11 AM   #2
Gorgon_VZ
Member
 
Location: Germany

Join Date: Oct 2016
Posts: 10
Default

Sorry the first explanation I posted was from bowties implemented command line help. The following description was actually taken from the manual:
--dpad <int>
“Pads” dynamic programming problems by <int> columns on either side to allow gaps. Default: 15.
Gorgon_VZ is offline   Reply With Quote
Old 01-10-2018, 11:25 AM   #3
sdriscoll
I like code
 
Location: San Diego, CA, USA

Join Date: Sep 2009
Posts: 438
Default

I sometimes use bowtie2 for the alignment stage of running RSEM (RNA-Seq gene expression). RSEM didn't allow INDELS in the alignments so as part of tweaking bowtie2 to not report those I was advised to set '--dpad 0'.
__________________
/* Shawn Driscoll, Gene Expression Laboratory, Pfaff
Salk Institute for Biological Studies, La Jolla, CA, USA */
sdriscoll is offline   Reply With Quote
Reply

Tags
bowtie, indel, large gap, mapping, parameters

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 09:18 AM.


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