Entering edit mode
6.3 years ago
NaskoMak
•
0
I have a library of 50,000 short sequences (20 bases each). I also have a fastq file from Illumina with short reads (23 bases each). I need to map the sequenced reads (first 20 bases) to the catalog, allowing for at most 1 mismatch and no indels. Could anyone point me to a suitable aligner? If you could also give me some advice on what parameters of that aligner would be suitable given this problem, that would be awesome!
Bowtie has an explicit option to control for mismatches (-n). You may check out the manual.