I'm using bowtie in repetitive regions and need all possible alignments of these reads. My question is about the backtracking mentioned in the original bowtie paper, http://genomebiology.com/2009/10/3/R25. If I understand correctly the algorithm limits the number of backtracking steps to 125. That means that you get all valid alignments for 125 strings? Is this still the case? Is this irrelevant if you are only allowing 0-1 mismatches?
-
-
Assuming you're using the --all option:
In -n mode (quality-aware) you can specify --tryhard to get all (?) the alignments. Not sure if it guarantees all of them, but the manual claims maximum sensitivity with that.
In -v mode (ignore quals) it'll always report all the alignments (implied by the manual).
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...09-01-2026, 05:41 AM
Comment