To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at

How do we generate BQSR plotsduring 1st phase i.e while running BaseRecalibrator?

There is an option to generate pdf plots of BQSR results, but I am not getting how to use it.
I gave command line:
java -jar GenomeAnalysisTK-2.3-9-ge5ebf34/GenomeAnalysisTK.jar -T BaseRecalibrator -ql 95 -I /storage/BQSR_testing/test.bam -R /storage/BQSR_testing/MISEQ/basic_files/hg19.fa -knownSites /storage/BQSR_testing/dbsnp137.vcf -o /storage/BQSR_testing/NA12878_with_dbsnp_137_recal_data_exome_for_plots.grp -plots -nct 3

In this case nothing was generated. Then I gave another command line mentioned below:

java -jar GenomeAnalysisTK-2.3-9-ge5ebf34/GenomeAnalysisTK.jar -T BaseRecalibrator -ql 95 -I /storage/BQSR_testing/test.bam -R /storage/BQSR_testing/MISEQ/basic_files/hg19.fa -knownSites /storage/BQSR_testing/dbsnp137.vcf -o /storage/BQSR_testing/NA12878_with_dbsnp_137_recal_data_exome_for_plots.grp -plots file.pdf -nct 3

Then also nothing was generated in the plots form.In logs it gave some thing that I have posted below:

INFO 12:36:02,166 BaseRecalibrator - Calculating quantized quality scores...
INFO 12:36:02,259 BaseRecalibrator - Writing recalibration report...
INFO 12:36:04,310 BaseRecalibrator - ...done!
**INFO 12:36:04,311 BaseRecalibrator - Generating recalibration plots...
WARN 12:36:04,374 RScriptExecutor - Skipping: Rscript (resource)org/broadinstitute/sting/utils/recalibration/BQSR.R /tmp/BQSR829189491498949821.csv /storage/BQSR_testing/NA12878_with_dbsnp_137_recal_data_exome_for_plots.grp /storage/BQSR_testing/file.pdf **
INFO 12:36:04,374 BaseRecalibrator - Processed: 9974 reads
INFO 12:36:04,376 ProgressMeter - done 9.97e+03 8.7 s 14.5 m 0.0% 3.1 w 3.1 w


Sign In or Register to comment.