Question: Best tools to search for and report all alignments
0
gravatar for hersh
3 days ago by
hersh0
hersh0 wrote:

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.

sequence alignment • 104 views
ADD COMMENTlink modified 3 days ago by Biogeek310 • written 3 days ago by hersh0
2
gravatar for h.mon
3 days ago by
h.mon21k
Brazil
h.mon21k wrote:

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? .

ADD COMMENTlink written 3 days ago by h.mon21k

Thanks for the response. I will take a look at segemehl for my analysis, and I appreciate the additional links and info.

ADD REPLYlink written 2 days ago by hersh0
0
gravatar for Biogeek
3 days ago by
Biogeek310
Biogeek310 wrote:

Subread align or STAR with custom settings to your desire? . These two are pretty fast and accurate.subread slightly pipping STAR according to literature

ADD COMMENTlink written 3 days ago by Biogeek310
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1051 users visited in the last hour