Why does `bowtie` only find perfect matches even when specifying mismatches? (ie. `-n 3`)
I tried a few methods to align a set of 32-mers against a reference genome, to find where they may potentially align with at most 3 mismatches.
bowtie --sam -v 3 -p 8 --all -f --norc hg38 kmers.fa kmers_aln.sam
When using bowtie, even though I specify the number of mismatches -v 3 and to capture all possible alignments -a, the resultant matches are all perfect matches (CIGAR = 32M). Is there a different parameter I am missing that increases the aligner's sensitivity? Additionally, for personal learning, why isn't bowtie able to capture the mismatches - is it because it only ranks its matches and saves only the exact matches?
1 answer
Log in to answer this question.
More posts like this
-
Which locus should be considered the true ortholog after liftOver disagreement between chrUn and ch…
written by jhonatan.henrique •
-
phage sequence no alignment to genome
written by Luwell •
-
Tutorial: How to align RNA and DNA 3D structures using ARTEMIS
written by Eugene F. Baulin •
-
Mapping all 32-mers within Hamming Distance = 3
written by giova34 •
-
Mitochondrial light/heavy strand allocation of variants
written by caearthworms •
-
Aligning SNP array probes to human genome with BLAT -- how to prioritize matches?
written by croberts377 •
-
De novo RNA-seq
written by blawney •
-
Exact Matching With Bowtie, Blat And Blast+
written by Alex Reynolds
-
Bowtie Inaccuracy Limit
written by Arpssss •
-
Tutorial: Mapping Reads With Bwa And Bowtie
written by MSU course 2011 •
