Converting blastn output to VCF (for SNPs): what is the quickest way?

Question

[Matteo Schiavinato]

Hello guys,

Title is pretty self-explanatory: what is the quickest way to convert the output of command line blast (say xml format, but I'm open to any suggestion) to a pseudo-VCF file?

I am interested in getting all the positions where I have a SNP (mismatch) between query and subject. Doesn't have to be a VCF or pseudo-VCF by force, I just need scaffold + position + variant.

I guess this question has been asked a lot, but I can't seem to find any answer. Would you suggest me something / link me something that I might have missed?

[Pierre Lindenbaum]

duplicate/see Find Mismatch And Gap Positions