The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

#### Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

#### ☞ Did you remember to?

1. Search using the upper-right search box, e.g. using the error message.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

#### ☞ Formatting tip!

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks (  ) each to make a code block as demonstrated here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

# No plots generated by the BaseRecalibrator walker

Member Posts: 4

I cannot produce BQSR plots, although I can open the grp file with gsa.read.gatkreport.

Here's the command:

java -Xmx1g -jar $shares/GenomeAnalysisTK-2.3-6-gebbba25/GenomeAnalysisTK.jar \ -T BaseRecalibrator \ -I ./0.reorder.bam \ -R$shares/ftp.broadinstitute.org/bundle/2.3/hg19/ucsc.hg19.fasta \
-knownSites $shares/ftp.broadinstitute.org/bundle/2.3/hg19/dbsnp_137.hg19.vcf \ -BQSR ./0.reorder.bam.recal.grp \ -o ./0.reorder.bam.post_recal.grp \ --plot_pdf_file ./0.reorder.bam.post_recal.grp.pdf \ -L chr1:1-1000 \ -l DEBUG \ --intermediate_csv_file ./0.reorder.bam.post_recal.grp.csv ##### ERROR stack trace java.lang.NullPointerException at org.broadinstitute.sting.utils.Utils.join(Utils.java:286) at org.broadinstitute.sting.utils.recalibration.RecalUtils.writeCSV(RecalUtils.java:450) at org.broadinstitute.sting.utils.recalibration.RecalUtils.generateRecalibrationPlot(RecalUtils.java:394) at org.broadinstitute.sting.gatk.walkers.bqsr.BaseRecalibrator.generatePlots(BaseRecalibrator.java:474) at org.broadinstitute.sting.gatk.walkers.bqsr.BaseRecalibrator.onTraversalDone(BaseRecalibrator.java:464) at org.broadinstitute.sting.gatk.walkers.bqsr.BaseRecalibrator.onTraversalDone(BaseRecalibrator.java:112) at org.broadinstitute.sting.gatk.executive.Accumulator$StandardAccumulator.finishTraversal(Accumulator.java:129)
at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:91)

It looks like the csv file is not being produced.

Thanks!

Tagged:

Alright, what's failing is the method that creates the data lines to be written to the csv file. It's an extremely simple operation so there's got to be something wrong with your data. Considering the interval that you're running on is extremely short, there may not even be any valid data at all in the part of the table that fails to write. You should try running again from the first step of recalibration with a much longer interval, eg 20:10000000-20000000.

Geraldine Van der Auwera, PhD

All the data needed for the plots should be in the gatkreport; what do you see in that file? Have you tried running without specifying the intermediate csv file?

Geraldine Van der Auwera, PhD

• Member Posts: 4

The grp file looks good. I tried to use BQSR.R to make the plot, but it seems to require a csv file.
If I specify the intermediate csv file, then it contains only the header line. If not, then it is not produced.
Thanks

What error do you get when you try using BQSR.R to make the plot?

Geraldine Van der Auwera, PhD

• Member Posts: 4

After making some changes for ggplot2 v0.9.3:
opts( -> theme( and theme_ -> element_

I get:

Error in distributeGraphRows(list(a, b, c), c(1, 1, 1)) :
object 'a' not found

It looks like a, b and c are created based on the lines in the csv file, and I don't have any.

Thanks
Alex

OK, that makes sense. Can you please run your original command again with -l DEBUG but without specifying the intermediate csv file filename, then post the console output? I need to know if with the internal default, the csv file gets created properly for you and if not, at what point it fails exactly.

Geraldine Van der Auwera, PhD

• Member Posts: 4

Here is the output for the command above minus --intermediate_csv_file ./0.reorder.bam.post_recal.grp.csv`
Thanks!

Alright, what's failing is the method that creates the data lines to be written to the csv file. It's an extremely simple operation so there's got to be something wrong with your data. Considering the interval that you're running on is extremely short, there may not even be any valid data at all in the part of the table that fails to write. You should try running again from the first step of recalibration with a much longer interval, eg 20:10000000-20000000.

Geraldine Van der Auwera, PhD

• ShanghaiMember Posts: 14

GATK 3.4 is used, but the option "--intermediate_csv_file" is not difined. How to generate the file *.grp.csv? Thanks a lot.