This is a test version of Biostars. For the public version, visit https://www.biostars.org.
LiftoverVcf no output

Hey all,

I am trying to run Picard's LiftoverVcf to lift over a hg19 to hg38, but there is no output. It runs with no errors, but it doesn't make the output files. Any help would be appreciated.

Here is the relevant section of the script used. I even tried unzipping all input files, but no output files are produced.

# $1 is the input vcf to the script    
java -jar /home/groups/prime-seq/pipeline_tools/bin/picard.jar LiftoverVcf \
    I=$1 \
    O=$output_vcf \
    CHAIN=/home/path/common_files/hg19ToHg38.over.chain \
    REJECT=$output_reject_vcf \
    R=/home/path/genomes/GCF_000001405.39_GRCh38.p13_genomic.fna

Here is the output.

    INFO    2020-06-01 17:22:33 LiftoverVcf Loading up the target reference genome.
[Mon Jun 01 17:22:33 PDT 2020] picard.vcf.LiftoverVcf done. Elapsed time: 0.01 minutes.
Runtime.totalMemory()=2039414784
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:287)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:282)
    at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:103)
    at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:113)
assembly genome liftover sequencing picard

0 answers

No answers yet.

Log in to answer this question.