Question: How to decide Gap/mismatch penalties ?
gravatar for KVC_bioinfo
20 months ago by
KVC_bioinfo380 wrote:


I am a new student of Bioinformatics. I am trying to learn how to set optimal the gap or mismatch penalties for aligning the reads. So, If I am trying to align the reads(RNA-seq) How to decide what is the best gap penalty or what is the best mismatch penalty? particularly while using bowtie2.

Thank you in advance.

#gappenalty #bowtie2 • 943 views
ADD COMMENTlink modified 4 weeks ago by Biostar ♦♦ 20 • written 20 months ago by KVC_bioinfo380
gravatar for Devon Ryan
20 months ago by
Devon Ryan89k
Freiburg, Germany
Devon Ryan89k wrote:

Firstly, don't use bowtie2 for spliced data.

Secondly, use the defaults. As a rule of thumb, always use the defaults with tools unless you have a good reason to do otherwise (the exception is the number of threads to use). It's generally the case that the person who wrote a tool put more thought into the various implications of the settings than you'll want to. Have a look at the results and judge then if some of the settings need to be tweaked.

ADD COMMENTlink written 20 months ago by Devon Ryan89k

I am aligning it to reference transcriptome.

ADD REPLYlink written 20 months ago by KVC_bioinfo380
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 927 users visited in the last hour