This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to use minimap2?

Hi, I am new to bioinformatics and I am tasked to do a read alignment. I have fastq.gz files, and I should run minimap2 over those datasets, which are long reads from Oxford nanopore. In the manual sent me, minimap2 [-x $preset] -d reference_index.mmi reference.fa and minimap2 -ax map-ont [-t $threads] reference_index.mmi ont_reads.fq > output.sam is said to be run over dataset. I understood that first one is creating minimizer index and second for mapping. How does .mmi represent here and how can use fast.gz files for read mapping(lets assume I have file dataset.fastq.gz)? Should I run these two command consequtively?

sequencing indexing minimap2 read-mapping

1 answer

Correct. You create index in step 1 and the align your data in step 2. Two command above represent these two steps.

How does .mmi represent here

That is the index file created from reference.fa.

Log in to answer this question.