Question: convert vcf file from GRCh38 to GRCh37 using picard LiftoverVcf
1
gravatar for jonessara770
3.4 years ago by
jonessara770200
jonessara770200 wrote:

Hello,

I am trying to convert vcf file from GRCh38 to GRCh37 using picard LiftoverVcf. I added chr to my vcf files and ran the software. However, I get the following errors. would you please advise? btw, I am using java openjdk-jdk1.8.0_72.

INFO    2016-12-26 11:52:02 LiftoverVcf Loading up the target reference genome.
INFO    2016-12-26 11:52:19 LiftoverVcf Lifting variants over and sorting.
[Mon Dec 26 11:52:19 EST 2016] picard.vcf.LiftoverVcf done. Elapsed time: 0.30 minutes.
Runtime.totalMemory()=5999427584
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" java.lang.NullPointerException
    at java.lang.String.checkBounds(String.java:384)
    at java.lang.String.<init>(String.java:324)
    at htsjdk.samtools.util.StringUtil.bytesToString(StringUtil.java:301)
    at picard.vcf.LiftoverVcf.doWork(LiftoverVcf.java:196)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:209)
    at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:95)
    at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:105)
dnaseq • 3.4k views
ADD COMMENTlink modified 3.4 years ago by Pierre Lindenbaum128k • written 3.4 years ago by jonessara770200

? btw, I am using java openjdk-jdk1.8.0_72.

and what is the version of picard ?

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by Pierre Lindenbaum128k
1

from the trace StringUtil is called at line 196:

at htsjdk.samtools.util.StringUtil.bytesToString(StringUtil.java:301)
at picard.vcf.LiftoverVcf.doWork(LiftoverVcf.java:196)

But there is no call to StringUtil in the current source of picard: https://github.com/broadinstitute/picard/blob/master/src/main/java/picard/vcf/LiftoverVcf.java#L196 , I suspect you're using an old version of picard.

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by Pierre Lindenbaum128k

and what was your command line ?

ADD REPLYlink written 3.4 years ago by Pierre Lindenbaum128k
1

Thanks Pierre,

I ran it without adding chr to the original vcf file and it runs but obviously can not liftover most of the variants. I don't know what is wrong with my vcf file when I am adding chr. any thoughts?

ADD REPLYlink written 3.4 years ago by jonessara770200

../openjdk-jdk1.8.0_72/bin/java -Dsamjdk.buffer_size=4194304 -Xmx6g -Xms2g -jar ../picard-tools-2.0.1/picard.jar LiftoverVcf TMP_DIR=../tmp I=myvcf.hg19.vcf O=myvcf.hg19.lift.flt.vcf CHAIN=hg38ToHg19.over.chain REJECT= myvcf.hg19.rejected.vcf R= Homosapiens.GRCh37.fa

ADD REPLYlink written 3.4 years ago by jonessara770200

I'm also getting the same error, but I'm not quite making it to the ' LiftoverVcf Lifting variants over and sorting' step. (java-jdk/1.8.0_92, picard/2.8.1) My vcf doesn't have chr prefixes and I removed the 'chr' in the reference.

java -Xmx1G -jar ${PICARD} LiftoverVcf I=./17.2057_passedvariants_T99.9.vcf O=./17.2057_passedvariants_T99.9_chr38.p10.vcf CHAIN=./hg19ToHg38.over.chain REJECT=./17.2057_passedvariants_T99.9_chr38.p10.rejected.vcf R=./Homo_sapiens_assembly38_WithoutHaplotypes.fasta


 [Mon Jul 17 14:07:46 CDT 2017] picard.vcf.LiftoverVcf INPUT=./17.2057_passedvariants_T99.9.vcf OUTPUT=./17.2057_passedvariants_T99.9_chr38.p10.vcf CHAIN=./hg19ToHg38.over.chain REJECT=./17.2057_passedvariants_T99.9_chr38.p10.rejected.vcf REFERENCE_SEQUENCE=./Homo_sapiens_assembly38_WithoutHaplotypes.fasta    WARN_ON_MISSING_CONTIG=false WRITE_ORIGI  2 [Mon Jul 17 14:07:46 CDT 2017] Executing as jdilts1@cri16cn083 on Linux 2.6.32-573.12.1.el6.x86_64 amd64; Java HotSpot(TM) 64-Bit Server VM 1.8.0_92-b14; Picard version: 2.8.1-SNAPSHOT
INFO    2017-07-17 14:07:47     LiftoverVcf     Loading up the target reference genome.
[Mon Jul 17 14:07:47 CDT 2017] picard.vcf.LiftoverVcf done. Elapsed time: 0.01 minutes.
Runtime.totalMemory()=1029177344
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" java.lang.NullPointerException
     at picard.vcf.LiftoverVcf.doWork(LiftoverVcf.java:158)
     at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:208)
     at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:95)
     at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:105)
ADD REPLYlink modified 2.9 years ago • written 2.9 years ago by diltsjeri440
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1242 users visited in the last hour