Question: Reporting of multiple alignments using hisat2
gravatar for
21 months ago by
markus.glass20 wrote:


I would like try HISAT2 for aligning RNA-Seq data but I am uncertain about using the -k parameter for reporting N multiple alignments. The manual says a primary alignment is an alignment whose alignment score is equal or higher than any other alignments. But it also says HISAT2 does not guarantee that the N alignments reported are the best possible in terms of alignment score. Does this mean i) the N reported alignments may not be all alignments with the best alignment score or does it mean ii) there may be alignments with a better score than of those N reported?

i) would be clear, but in case of ii) I don't understand how this is possible if a primary alignment has an alignment score equal or higher than any other alignments.

Thanks in advance!

rna-seq • 1.1k views
ADD COMMENTlink modified 9 months ago by predeus890 • written 21 months ago by markus.glass20
gravatar for predeus
9 months ago by
predeus890 wrote:

This is actually one of the most confusing things about alignment in general.

What this description means (I think) is that if you have -k set to 10, and aligner will find 10 locations where the read aligns with a certain mapping score, it will report them and stop looking further.

What I do not get is that, according to the description, there should only be primary alignments in the resulting SAM. That's certainly not the case: quite a few of them have "secondary" bit-wise flag (256) set.

ADD COMMENTlink written 9 months ago by predeus890
Please log in to add an answer.


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