This is a test version of Biostars. For the public version, visit https://www.biostars.org.
svtyper tool: ValueError: mate not found

Hi there,

I'm using svtyper in order to get the genotype of some structural variants called using Lumpy. Svtyper starts running and everything seems OK until it arrives to a specific variant of the input vcf file. At this record, svtyper crashes with the following error:

File "pysam/calignmentfile.pyx", line 836, in pysam.calignmentfile.AlignmentFile.mate (pysam/calignmentfile.c:10945)
ValueError: mate not found

svtyper calls pysam at some point to get the mate pairs. I've read the code of pysam, and this is what I've found:

Throws a ValueError if read is unpaired or the mate is unmapped.

So, I think that svtyper crashes due to the presence of a read with unmapped mate supporting the SV. If this is the true reason of the error, does this make sense? I mean, I'm interested also in those reads which their mate is unmapped, this is a good pointer about the presence of a SV in that position, isn't it?

Any ideas?

lumpy svtyper structural-variants
• 3,665 views • link
updated by Ram Ram • written by iraun iraun

2 answers

Well, I did not know where the problem comes from but I found a "temporary" solution just in case that anyone is struggling to get this.
Using samblaster tool (available here:https://github.com/GregoryFaust/samblaster) with --addMateTags, mate tags (MC, MQ) are added to the bam file. Without them, svtyper must seek to the corresponding mate position in the bam file. Using this approach, svtyper and therefore pysam, works nice without any error, and also, faster.

• 0 views • link
written by iraun iraun

Now just add -M flag in svtyper!

James

• 0 views • link
written by wjar6718 • wjar6718

Log in to answer this question.

More posts like this