detect and combine indels into one line like Mutect2 format
0
0
Entering edit mode
2.3 years ago
whb ▴ 60

Hi,

I would like to get a list of common SNVs and indels from two variant callers (Mutect2 and lofreq). but lofreq output indels as two SNVs whereas Mutect2 output indels as one line (shown below) and bcftools -isec does not recognise they are the same.

How can I combine the indels SNVs in the lofreq vcfs into one line like the format in Mutect2?

lofreq vcf:

chr5    117401004   .   G   A
chr5    117401005   .   G   T

Mutect vcf

chr5    117401004   .   GG  AT

Thank you!

indels snvs vcf • 437 views
ADD COMMENT

Login before adding your answer.

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