Hi!
I made a minimap2 index for a custom genome. I was planning on using it for pacbio and nanopore data.
I created the index using minimap2 -x map-pb -d genome.mmi genome.fa. Unfortunately, I forgot to make one specifically for the nanopore data using map-ont... Instead, I mapped both data sets to the same index using
For the nanopore: minimap2 -a -x map-ont genome.mmi sample.fastq.gz
For the pacbio: minimap2 -a -x map-pb genome.mmi sample.fastq.gz
I have limited experience with minimap2 and I am wondering if this is a big deal, since the mapping was using the correct parameters - it's just the index that might not be perfect for the data...
Can someone explain to me how this could affect the performance? If it would? I can't seem to find a thorough explanation about this.
Thanks,
2 answers
Directly from the minimap2 page
Importantly, it should be noted that once you build the index, indexing parameters such as -k, -w, -H and -I can't be changed during mapping. If you are running minimap2 for different data types, you will probably need to keep multiple indexes generated with different parameters. This makes minimap2 different from BWA which always uses the same index regardless of query data types.
So the answer is yes, I don't know how these parameters change between data types but I wouldn't assume they would work interchangeably without any benchmarking.
Why build an index at all? I never build indices for minimap2 mapping, unless forced to by a further tool eg megalodon. Then I can set the mapping type using appropriately using -x map-ont or map-pb as appropriate.
Log in to answer this question.
-
Differentiate CLR and Hifi readings in datasets
written by María José •
-
Ribosomal RNA Depletion in silico
written by CL •
-
Sequencing Depth (Read Depth) Calculations
written by R.L. •
-
How to use minimap2?
written by MobiusT •
-
IGV Sashimi plot from ONT long read sequencing of cDNA amplicons
written by WRUZIC1 •
-
no header in SAM file
written by Maxine •
-
Polishing PACBIO corrected assembly using ONT and ILLUMINA data (.fastq)
written by K •
-
Question regarding minimap2 presets
written by Rogerio Ribeiro •
-
Mapping to custom reference minimap2/IGV
written by sendhelp •
-
Long-read alignment with Minimap2
written by gokberk
compuTE : If you can provide an update in this thread with some stats (
samtools idxstatsmay be good enough) as to how the numbers change with a common index used for ONT and PacBio VS specific indexes made for the two datasets.minimap2is such an excellent aligner that the gross stats may remain very similar (specifics may change some).Rarely do people have both types of data and this insight would provide a very useful datapoint.