This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Comparing variants in vcf files that have been mapped to GRCh37 and GRCh38

I have a VCF format file, which contains variants found for file mapped to hg38 (reference from GATK hg38 bundle) and I would like to compare this variants with VCF file from 1000 genome project, which is mapped to GRCh37. By comparing this files I mean finding variants shared by both files. Are the coordinates in this files the same or should I somehow convert them?

hg38 hg19 vcf
• 5,000 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by cl10101 • cl10101

1 answer

convert the file with picard LiftoverVcf : https://broadinstitute.github.io/picard/command-line-overview.html#LiftoverVcf then compare How to compare 2 VCF files

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you for response. I removed char prefix in hg19ToHg38.over.chain file as you suggested and now I have error:

"Exception in thread "main" htsjdk.tribble.TribbleException: Badly formed variant context at location chr1:789016; getEnd() was 789016 but this VariantContext contains an END key with value 724396"

I would be grateful for any suggestion how to solve this problem.

• 0 views • link
written by cl10101 • cl10101

there is no output for 789016

• 0 views • link
written by cl10101 • cl10101

try "grep -Fw 724396 your.vcf"

I got the same problem:

Exception in thread "main" htsjdk.tribble.TribbleException: Badly formed variant context at location 13:32707645; getEnd() was 32707645 but this VariantContext contains an END key with value 33281782

and the corresponding record is:

13  33281782    esv3631727  C   <INS:MT>    100 PASS    AC=154;AF=0.0307508;AN=5008;CIEND=-1,0;CIPOS=0,1;CS=NUMT_umich;END=33281782;NS=2504;SVTYPE=INS;IMPRECISE;DP=15932;EAS_AF=0.0129;AMR_AF=0.0058;AFR_AF=0.0908;EUR_AF=0.001;SAS_AF=0.0164;VT=SV    GT  0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0
• 0 views • link
written by Min Dai Min Dai

Log in to answer this question.

More posts like this