CollectVariantCallingMetrics Error
1
0
Entering edit mode
7 months ago
Sd • 0

Hello, I am trying to get metrics table for a GVCF file using CollectVariantCallingMetrics and I am getting the following error message at the end of run and does not generate metrics table. I was wondering how can I solve this issue.

[Thu Feb 29 16:50:50 EST 2024] picard.vcf.CollectVariantCallingMetrics done. Elapsed time: 3.42 minutes. Runtime.totalMemory()=2013265920 To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp java.lang.NullPointerException: Cannot invoke "htsjdk.samtools.SAMSequenceRecord.getSequenceLength()" because the return value of "htsjdk.samtools.SAMSequenceDictionary.getSequence(String)" is null at picard.util.DbSnpBitSetUtil.loadVcf(DbSnpBitSetUtil.java:235) at picard.util.DbSnpBitSetUtil.createSnpAndIndelBitSets(DbSnpBitSetUtil.java:203) at picard.vcf.CollectVariantCallingMetrics.doWork(CollectVariantCallingMetrics.java:110) at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:289) at org.broadinstitute.hellbender.cmdline.PicardCommandLineProgramExecutor.instanceMain(PicardCommandLineProgramExecutor.java:37) at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:160) at org.broadinstitute.hellbender.Main.mainEntry(Main.java:203) at org.broadinstitute.hellbender.Main.main(Main.java:289)

I use the following commands:

$ /home/gatk-4.4.0.0/./gatk CollectVariantCallingMetrics \
--DBSNP /home/bundle/Homo_sapiens_assembly38.dbsnp138.vcf \
--INPUT /home/sample.g.vcf.gz \
--OUTPUT /home/sample.table \
--GVCF_INPUT true
CollectVariantCallingMetrics GATK GVCF • 713 views
ADD COMMENT
0
Entering edit mode

you asked: Individual vs. joint call VCFs How to slice a CRAM file into the 50kb regions padded with 1kb? Hwo to identify that BQSR is performed on CRAM file How to Split 3000 WGS CRAM files into 1Mbp length chunks

Don't forget to follow up on your threads, that is bad etiquette. If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they all work. If an answer was not really helpful or did not work, provide detailed feedback so others know not to use that answer.

Upvote|Bookmark|Accept

ADD REPLY
0
Entering edit mode
7 months ago

your dbsnp VCF is missing a dictionary in the header (lines starting with "##contig" )

see also: https://gatk.broadinstitute.org/hc/en-us/articles/360042477832-UpdateVcfSequenceDictionary-Picard

ADD COMMENT
0
Entering edit mode

I did UpdateVcfSequenceDictionary for dbsnp VCF file and updated the dbsnp contigs based on hg38 reference file. In the dbsnp VCF header I see the ##contig lines. However, I still have the error in CollectVariantCallingMetrics. What other solutions can help solve this issue?

ADD REPLY
0
Entering edit mode

do you have the very same chromosome notation between the VCF and the DBSNP VCF ? eg. '1' vs 'chr1' ?

ADD REPLY
0
Entering edit mode

Yes. The notations are the same.

ADD REPLY

Login before adding your answer.

Traffic: 879 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6