Picard LiftoverVcf warning - input malformated

I'm working on liftover vcf from hg38 to hg19. I got this error message back:

Exception in thread "main" htsjdk.tribble.TribbleException: The provided VCF file is malformed at approximately line number 906: unparsable vcf record with allele -, for input source: MYPATH

Has anyone encountered similar problem? Please share any thoughts, thanks!

vcf picard warning error liftover

• 1,695 views • link

written 5.2 years ago by Yuri •

1 answer

- is not an allele supported by the picard/htsjdk api: https://github.com/samtools/htsjdk/blob/7719274fe370a51a24e6067de21bbe7e18c160a9/src/main/java/htsjdk/variant/variantcontext/Allele.java#L332

you could fix this by converting such alleles by N or <.> (symbolic allele)

• 0 views • link

written 5.2 years ago by Pierre Lindenbaum

Thanks! I will try it :D

• 0 views • link

written 5.2 years ago by Yuri •

More posts like this

How to format "I" and "D" in vcf version 4.2 for liftover analysis in GATK
written 2.2 years ago by Omics data mining

Hello everyone I am facing challenges with liftover of a VCF file from hg19 to hg38 using GATK because of 'I' and 'D' annotations representing …
gatk SelectVariants is giving duplicate allele error while extracting SNPs out of vcf file
written 2.7 years ago by Bushra Bibi

I am trying to extract snps out of merged vcf file using gatk SelectVariants command but it is giving following error: htsjdk.tribble.TribbleException: The provided VCF …
VCF record format error when trying to run Beagle 5.4
written 3.1 years ago by r.tor •

Hi everyone, I am trying to run the Beagle 5.4 software where it needs to read a VCF file. The error message I have is: …
vcf.Reader() in Python doesn't read my VCF
written 3.3 years ago by Najla •

I was working on a building a GUI with wxPython and wanted to import a VCF file with open(file_path, 'r') as file: vcf_reader = vcf.Reader(file) …
htsjdk.tribble.TribbleException: The provided VCF file is malformed
written 4.7 years ago by Egelbets •

I have VCF files that I want to convert to a more readable TSV file using [GATK VariantsToTable][1], and I also want to load in …
No header in VCF file
written 4.7 years ago by Begonia_pavonina •

Hello everyone, I am working with a specific variant calling pipeline, and the output is a VCF file missing headers. It seems there are no …
GATK Haplotype Caller error
written 8.7 years ago by kirti.amity •

While ruuning GATK Haplotype caller for human, java -jar /biobank/seq/kirti/new/NA12878/GATK/GenomeAnalysisTK.jar -T HaplotypeCaller -R /biobank/seq/kirti/new/NA12878/reference/hg19_new.fa -I /biobank/seq/kirti/new/NA12878/nature_paper/"$a"/"$a"sorted2_bwa_aligned.bam --dbsnp /biobank/seq/kirti/new/NA12878/nature_paper/Homo_sapiens.vcf --sample_name "$a" -ploidy 2 -o /biobank/seq/kirti/new/NA12878/nature_paper/"$a"/"$a"bwa_GATK.vcf I …
GATK Indel sorting, vcf error
written 8.9 years ago by hdtms •

Hi I used the GATK pipeline until I got a vcf that had SNPs and Indels, so I used GATK Analysis tools to remove SNPs …
SAM validation error
written 9.5 years ago by fire_water

**Command** picard-tools-1.118/MarkDuplicates.jar I=file.Aligned.sortedByCoord.bam O=file_out.bam METRICS_FILE=file_out.metrics **Error Message** Exception in thread "main" htsjdk.samtools.SAMFormatException: SAM validation error: ERROR: Record 801397832, Read name K00135:24:H3V7TBBXX:1:1101:13078:1332, Mate Alignment start (1651771652) …
VariantRecalibrator Error message
written 12.1 years ago by cvu

Hi, I'm using VariantRecalibrator from GATK. I've generated my vcf files with Mpileup/bcftools. When I am using VariantRecalibrator, with this argument, ``` java -Xmx4g \ …