HISAT2: does the parameter -k 1 means unique mapping?
1
0
Entering edit mode
4.1 years ago
biolab ★ 1.3k

Hi, everyone,

The -k <int> parameter in HISAT2 is explained by the following in the manual: It searches for at most <int> distinct, primary alignments for each read.

If I use -k 1, does it mean unique mapping? Thank you very much!

hisat2 • 2.3k views
ADD COMMENT
1
Entering edit mode

See this for more: Does Bowtie2 Preform Unique Mapping
If you are interested in removing multi-mapped reads then see @Devon's answer here: How to extract unique mapped results from Bowtie2 bam results?

ADD REPLY
0
Entering edit mode

Thank you very much, genomax, these posts are very helpful.

ADD REPLY
1
Entering edit mode
4.1 years ago

As far as I understood it will then only produce 1 alignment per read - even if there are multiple locations in the genome to which that read would fit. So you would indeed get unique mappings, but that's cheating and incorrect.

ADD COMMENT
0
Entering edit mode

Hi, WouterDeCoster, thank you very much for your explanations.

ADD REPLY
0
Entering edit mode

Is it incorrect to include all reads even those that map >1 times? I have always kept uniquely mapped. Also is there an option to keep only uniquely mapped using HISAT2?

ADD REPLY

Login before adding your answer.

Traffic: 1407 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6