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

Answers

Sign In or Register to comment.