Hi,
I have short reads that I would like to align on a bacterial genome. I need to know the exact identity between the reads and the reference genome, even if there are mismatches or gaps near the 5' or 3' ends.
To be sure to make myself clear: on this 20nt sequence (MM in lowercase) I would expect an identity of 16/20= 80%
GaaGCCGTTCTTATAGTaaT
I tried with BWA but for a 36bp read I can only have up to max of 5 mismatches, it doesnt seem to be able to give the identity for more MM/gaps (is there anyway to do that ?)
I am thinking of using a global aligner (free license) that should also work correctly on relatively short reads (36bp). If this aligner could output an easy-to-parse output file that would be X-mas for me.
Also, I tried Exonerate but I cannot figured out how to make it work the way I want (the options --exhaustive --model affine:global do not seem to work ...).
Any help is appreciated.
Thanks
•
written
by
Pasta •
1 answer
Novoalign can do full Needleman-Wunsch. It is quite sensitive as well.
(Note, Novoalign is free for academic and non-profit use)
Seconding novoalign. I've ran it on deepseq datasets of 10-50 million 36bp reads, it's much slower than bowtie but still finishes within a day or two on a reasonably powerful desktop. The "native" novoalign output is pretty easy to parse, and it can also output sam files.
Log in to answer this question.
-
Exonerate affine:local output
written by Ak •
-
kmer alignment with mismatch
written by Vince
-
Ultrafast alignment of short sequences to a genome (with up to 5 mismatches!)
written by Aurelie MLB
-
Why does Blast report a different alignment when using the -m8 option
written by sebastian.f.ing •
-
Using Agrep With 2 Patterns
written by toshnam
-
Why Do I Get The Exonerate Error: Cannot Perform Heuristic Alignment Using Non-Local Model
written by Woa
-
Convert Bam To Fastq With Pe & Se Reads
written by Pasta •
-
Bwa Specify Number Of Maximum Mismatches And Gaps
written by Nicolas Rosewick
-
Shrimp Short Reads Aligner
written by Joe •
-
Cannot Install Latest Version Of Deseq
written by Pasta •
You will need a more sensitive tool than bwa, but this comes for the price of increased run time. How many reads do you have? You can nowadays run exact algorithms on up to a few thousand reads, but if it is more than that a heuristic approach is required.
I have 13 million reads....