Mapping to custom reference minimap2/IGV
1
0
Entering edit mode
3.1 years ago
sendhelp • 0

I'm trying to map nanopore reads to a mtDNA reference using minimap2 for the purpose of determining how many mitochondrial reads are present in the nanopore data.

The minimap command I'm using is;

~/apps/minimap2/minimap2 -ax map-ont mtref.fa reads.fq  > aln.sam

I know there are mtDNA reads present as I can align them to the reference, however when I look at the resulting indexed and sorted sam file in IGV there's no coverage?

Thanks!

mapping minimap2 igv nanopore • 2.2k views
ADD COMMENT
0
Entering edit mode

Have you made a custom genome in IGV using the same reference? By default IGV shows you the entire reference so you may need to zoom in before you start seeing the reads.

ADD REPLY
0
Entering edit mode

I have. I solved the initial problem by actually creating a sorted bam file and using that in IGV rather than just the same file. Thanks!

ADD REPLY
1
Entering edit mode
3.1 years ago
ATpoint 64k

You can simply check the full alignment using samtools idxstats for the number of reads mapping to chrM.No need to align against a specific subset (which is also not recommended).

ADD COMMENT
0
Entering edit mode

Thank you! Can I ask why it isn't recommended to align against a subset (I think you're referring to aligning genomic reads to a single chromosome reference?)

ADD REPLY
1
Entering edit mode

The aligner will always try to find the best match for each read in the library. If your library contains reads from the entire genome but you include only one chromosome (or one region) with the true origin of the read missing the aligner might false-positively put it somewhere else. I never quantified what the extend of false-positives is but this opinion is pretty much what I read here at Biostars from many of the more experienced folks.

ADD REPLY
0
Entering edit mode

Ah, makes sense! There isn't a reference genome available for me currently so I guess this is the only method?

ADD REPLY
1
Entering edit mode

Well I am spoiled by working with mouse and human only (and I often forget that when writing answers :-D ) but yeah if you have no reference I would include as much of the draft genome or the contigs you have and then filter for the regions/contigs/chromosomes you need.

ADD REPLY

Login before adding your answer.

Traffic: 1859 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