I am looking to align a fastq file of 50mer sequences against a reference genome. I am looking to return ALL possible alignments, not just the best alignments. I looked at using Bowtie2 with the -a flag, however in the documenation it says "Some tools are designed with this reporting mode in mind. Bowtie 2 is not! For very large genomes, this mode is very slow.".
Which tools are designed with this in mind? Ideally the output would return alignments in descending order by alignment score, equivalent to Bowtie2. If there is not a better option for this workflow, I will stick to using Bowtie2. Thanks in advance.
2 answers
There are several "all" mappers, segemehl is one I see recommended frequently. Also, see the HTS Mappers (long but even so incomplete) list, search for "A" under "alignment reported" for all-mappers.
Another suggestion at Best aligner for reporting all exact matches of multi-mapping short reads? .
Log in to answer this question.
-
bowtie2 mapping qualities of 0
written by Egelbets •
-
bowtie mismatch 0 alignment output
written by kb_93 •
-
Achieving equivalent of -m mode for Bowtie2
written by hersh •
-
Piping bowtie2 output directly into BAM
written by hersh •
-
Bowtie2 reports no secondary alignments
written by goodez
-
Efficiently performing intersections for multiple subsets of a BED file
written by hersh •
-
which alignment tool is better to use to align genome on short reference sequence ~5000b.p
written by bbb •
-
Bowtie -m and -v parameter equivalent in bowtie2?
written by datascientist28
-
Best aligner for reporting all exact matches of multi-mapping short reads?
written by Ryan Thompson
-
Best way to map paired-end, uniquely mapped reads with Bowtie2?
written by Ian
