Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

picard liftovervcf parsing error

$
0
0

I am using picard liftover vcf to algin variants from hg38 to hg19. I am not sure what the error is referring to though it seems to be an issue with the vcf (I have copied a few lines of the format). Thank you :).

  java -jar /home/cmccabe/Desktop/NGS/picard-tools-1.140/picard.jar LiftoverVcf \
  I=/home/cmccabe/Desktop/align/out.vcf \
  O=/home/cmccabe/Desktop/tvc/IDP.lifted_over.vcf \
  CHAIN=/home/cmccabe/Desktop/NGS/picard-tools-1.140/resources/hg38ToHg19.over.chain \
  REJECT=/home/cmccabe/Desktop/tvc/IDP_rejected_variants.vcf \
  R=/home/cmccabe/Desktop/NGS/picard-tools-1.140/resources/ucsc.hg19.fasta

java version on ubuntu 14.04
java version "1.8.0_111"
Java(TM) SE Runtime Environment (build 1.8.0_111-b14)
Java HotSpot(TM) 64-Bit Server VM (build 25.111-b14, mixed mode)

trace error
[Mon Jan 09 08:04:50 CST 2017] picard.vcf.LiftoverVcf INPUT=/home/cmccabe/Desktop/tvc/out.vcf OUTPUT=/home/cmccabe/Desktop/tvc/IDP.lifted_over.vcf CHAIN=/home/cmccabe/Desktop/NGS/picard-tools-1.140/resources/hg38ToHg19.over.chain REJECT=/home/cmccabe/Desktop/tvc/IDP_rejected_variants.vcf REFERENCE_SEQUENCE=/home/cmccabe/Desktop/NGS/picard-tools-1.140/resources/ucsc.hg19.fasta VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false GA4GH_CLIENT_SECRETS=client_secrets.json
[Mon Jan 09 08:04:50 CST 2017] Executing as cmccabe@DTV-A5211QLM on Linux 4.4.0-57-generic amd64; Java HotSpot(TM) 64-Bit Server VM 1.8.0_111-b14; Picard version: 1.140(a81bc82e781dae05c922d1dbcee737334612399f_1444244284) IntelDeflater
INFO 2017-01-09 08:04:51 LiftoverVcf Loading up the target reference genome.
INFO 2017-01-09 08:05:00 LiftoverVcf Lifting variants over and sorting.
[Mon Jan 09 08:05:00 CST 2017] picard.vcf.LiftoverVcf done. Elapsed time: 0.17 minutes.
Runtime.totalMemory()=6771703808
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" java.lang.IllegalStateException: Key CS found in VariantContext field INFO at chr1:11082610 but this key isn't defined in the VCFHeader. We require all VCFs to have complete VCF headers by default.
at htsjdk.variant.vcf.VCFEncoder.fieldIsMissingFromHeaderError(VCFEncoder.java:176)
at htsjdk.variant.vcf.VCFEncoder.encode(VCFEncoder.java:115)
at htsjdk.variant.variantcontext.writer.VCFWriter.add(VCFWriter.java:222)
at picard.vcf.LiftoverVcf.doWork(LiftoverVcf.java:206)
at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:206)
at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:95)
at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:105)

VCF format: (few lines total 370)

fileformat=VCFv4.1

fileDate=20130610

source=ensembl;version=75;url=http://e75.ensembl.org/homo_sapiens

reference=ftp://ftp.ensembl.org/pub/release-75/fasta/homo_sapiens/dna/

INFO=<ID=TSA,Number=0,Type=String,Description="Type of sequence alteration. Child of term sequence_alteration as defined by the sequence ontology project.">

INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency">

INFO=<ID=E_MO,Number=0,Type=Flag,Description="Multiple_observations.http://www.ensembl.org/info/docs/variation/data_description.html#evidence_status">

INFO=<ID=E_1000G,Number=0,Type=Flag,Description="1000Genomes.http://www.ensembl.org/info/docs/variation/data_description.html#evidence_status">

INFO=<ID=E_HM,Number=0,Type=Flag,Description="HapMap.http://www.ensembl.org/info/docs/variation/data_description.html#evidence_status">

INFO=<ID=E_Freq,Number=0,Type=Flag,Description="Frequency.http://www.ensembl.org/info/docs/variation/data_description.html#evidence_status">

INFO=<ID=E_C,Number=0,Type=Flag,Description="Cited.http://www.ensembl.org/info/docs/variation/data_description.html#evidence_status">

INFO=<ID=dbSNP_137,Number=0,Type=Flag,Description="Variants (including SNPs and indels) imported from dbSNP">

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT Individual

chr1 11082610 rs11689432 G A . . TSA=SNV;E_MO;E_Freq;E_HM;E_C;CS=pathogenic;dbSNP_137 GT 1/0
chr1 11107061 rs77977199 T G . . TSA=SNV;E_Freq;dbSNP_137 GT 1/0
chr1 17380507 rs11203289 G C . . TSA=SNV;E_MO;E_Freq;E_1000G;CS=pathogenic;dbSNP_137 GT 1/1


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>