This is a test version of Biostars. For the public version, visit https://www.biostars.org.
verifyBamID generating error 'VcfFileException'

Hello. I was trying to run this too called verifybamID to detect sample contamination with the command below, but for some reason it generates the error below termed as 'VcfFileException'. Someone please explain me the issue here?

command:

 ./verifyBamID --vcf /media/SNPs/AOGC-Genotyping.output.recalibrated.filtered.vcf --bam /media/BAM/AOGC-02-0010_C0J43ACXX-4-13.ReCal.sort.bam --out test --verbose --ignoreRG

Error:

Opening VCF file /media/SNPs/AOGC-Genotyping.output.recalibrated.filtered.vcf
Total of 1 out of 1 individuals successfully matched IDs
WARNING : Skipping marker chr1:17172 with multiple alternative alleles
terminate called after throwing an instance of 'VcfFileException'
  what():  # values = ./.:.:45 do not match with # fields in FORMAT field = 5 at sampleIndex = 133 See line 228.
Aborted (core dumped)
verifybamid

1 answer

Hi MAPK, I am not sure if you solved this.

But it appears to me that the VCF is not sorted. Is your file with multiple samples?

Log in to answer this question.