No Sequence Segments for secondary alignments in minimap2 while trying to run alignment on ONT data?
0
0
Entering edit mode
9 months ago
Adarsh • 0

Hello all, I was trying to align a set of ONT reads that I have against my reference and I noticed something strange in the generated SAM file. If a read has any secondary alignment, there are no sequence segments present in the output for it. Is this some default behavior for minimap2? I really need the secondary sequence segment for some further analysis. This is the command that I ran.

minimap2 -t 10 -ax map-ont reference.fasta combines_reads.fastq > alignment_out.sam 2> log.txt

I even tried the above command with --secondary=yes command but I still get an empty seq segment for secondary alignments. Could please help me with this?

minimap2 nanopore alignment ont • 461 views
ADD COMMENT

Login before adding your answer.

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