Question: Bowtie Seed Length (-L) And Non-Unique Reads Mode (M=-1)
gravatar for Eric Guo
6.7 years ago by
Eric Guo60
Eric Guo60 wrote:
Hi there, I have two questions regarding alignment using Bowtie: 1. Is there a way to set the Seed Length (-l) to the full length of each read instead of using a single Seed Length for all reads? 2. When using m = -1 mode (Suppress all alignments for a read if more than n reportable alignments exist (-m): OFF), will a read be randomly assigned to one of the alignment positions in the genome? Or, will all the alignments for one single read be reported in the final output .bam file.  Thanks in advance for your help. -Eric
alignment bowtie • 2.0k views
ADD COMMENTlink modified 6.7 years ago by Jennifer Hillman Jackson25k • written 6.7 years ago by Eric Guo60
gravatar for Jennifer Hillman Jackson
6.7 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hi Eric, The seed length is a single input value for all sequences in any particular job. If alignment are not suppressed, then all reportable hits are in the output. More about Bowtie can be found here: Best, Jen Galaxy team
ADD COMMENTlink written 6.7 years ago by Jennifer Hillman Jackson25k
Hello Eric, I apologize, but I made a mistake, please see below This is incorrect. If alignments are not suppressed, then reporting "-a" (all) is not the default. Rather, there are several 'reporting options', and they are used in combination and with 'alignment options' to achieve different output results. It is probably best to show you the documentation, so that you can compare different parameters and output they produce. Please see: Hopefully this clears up any confusion my original reply may have caused! Jen Galaxy team
ADD REPLYlink written 6.7 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.


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