Question: HISAT2: does the parameter -k 1 means unique mapping?
0
gravatar for biolab
20 months ago by
biolab1.1k
biolab1.1k wrote:

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 • 1.2k views
ADD COMMENTlink modified 20 months ago by WouterDeCoster37k • written 20 months ago by biolab1.1k
1

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 REPLYlink written 20 months ago by genomax64k

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

ADD REPLYlink written 20 months ago by biolab1.1k
1
gravatar for WouterDeCoster
20 months ago by
Belgium
WouterDeCoster37k wrote:

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 COMMENTlink written 20 months ago by WouterDeCoster37k

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

ADD REPLYlink written 20 months ago by biolab1.1k
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: 1576 users visited in the last hour