detect and combine indels into one line like Mutect2 format
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!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Splitting VCF using phase information
written by Kyra •
-
bcftools isec output for multiple vcf files interpretations
written by minoo •
-
drop duplicate insertion deletions in VCF at same position while keeping one
written by curious
-
How to handle VCFs from the same sample but using different aligners and variant callers?
written by aldhairmedico •
-
Tool: all2vcf: a tool to convert non-standard variant outputs (mummer. bcftools isec) to VCF
written by Matteo Schiavinato
-
bcftools isec not producing sites.txt
written by steve
-
bcftools isec: tandem mutations
written by whb •
-
Combine several VCF files
written by Max Ivon
-
VCF output when merging with bcftools merge
written by steven.davis •
-
Venn diagram for vcf files
written by sukesh1411 •
