Option to output insertion/deletion coordinates in query sequence?

[cizydorczyk]

Is there an option to output the coordinates of an insertion or deletion in the query sequence?

For example, if I wanted to investigate the genes present in a 5kb insertion, I would have to go back to my query contig and look at the annotations in that region. Is there a way to get svim-asm to output such coordinates? I see it can output query names, which I believe would at least tell me which contig if the assembly consists of multiple, but I have not found anything re: coordinates.

Thank you, Conrad

[clemgoub]

Hi there, I am asking myself the same question. I am calling SVs with svim-asm and would like to be able to project the VCF coordinate to each query (alternative) genome. The goal is to make a gff for each genome.

I see the question was also asked here https://github.com/eldariont/svim-asm/issues/11 and here https://github.com/eldariont/svim-asm/issues/14. But I don't think this was resolved. I understand this may not be a priority for the devs, but if someone in the community as come with a nice solution I'll be happy to try!

Before I reinvent the wheel (I'm thinking cross-mapping the flanking regions), I was wondering if we could get this info from SVIM itself.

Cheers,

Clément