This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Minimap2 giving 2 TB sam file as an output

Hello,

I have been trying to run Minimap2 with fasta files of ~200Mb. I have 3 such files which I am trying to map onto a reference genome. The output sam file of first 2 files are about 500MB. But the third file gives a 2TB sam file and it has not completed running. Is there any reason of why this could be happening ? Any help would be appreciated.

The command I used is as follows:

minimap2 -a reference.mmi query.fasta -o output.sam -t 20
sequence minimap2 alignment
• 6,305 views • link
updated by h.mon h.mon • written by Inquisitive8995 Inquisitive8995

You could consider avoiding creating a sam intermediate:

minimap2 -a -t 20 reference.mmi query.fasta | samtools sort -o output.bam
• 0 views • link
written by WouterDeCoster WouterDeCoster

Thanks for your answer. But I am still confused as the other 2 files of the same size gave the output of around 500Mb and this gives 2TB.

• 0 views • link
written by Inquisitive8995 Inquisitive8995

Perhaps you have a lot of multi-mapping or secondary alignments?

• 0 views • link
written by GenoMax GenoMax

Thanks for your reply. Is there anyway with which I can restrict the mapping to only the top hit or something like that ? Are there any parameters in Minimap to do so ?

• 0 views • link
written by Inquisitive8995 Inquisitive8995

You can set -N 0 to not retain secondary alignments.

• 0 views • link
written by GenoMax GenoMax

Why do you tag this as "assembly" while you are aligning? Why do you think this is a "software error"? Wouldn't "minimap2" be the best tag here?

• 0 views • link
written by WouterDeCoster WouterDeCoster

0 answers

No answers yet.

Log in to answer this question.

More posts like this