Heads up:
We’re moving the GATK website, docs and forum to a new platform. Read the full story and breakdown of key changes on this blog.
Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Base Recalibrator Issue

AshuAshu Member
edited November 2012 in Ask the GATK team

Hi
I got the following error with GenomeAnalysisTK-2.2-2-gf44cc4e's Base Recalibrator.

##### ERROR MESSAGE: Key 2006 is too large for dimension 2 (max is 2001)

I also ran the picard's validateSamFile to validate my BAM file and it says NO ERRORs.
What exactly does this error mean? what key is it talking about? And how can I fix it?
Thanks,
Ashu

Post edited by Geraldine_VdAuwera on

Answers

  • droazendroazen Cambridge, MAMember, Broadie, Dev ✭✭

    Hi Ashu,

    Could you please post the full command line you're using to run the BaseRecalibrator, and describe any previous processing steps in your pipeline involving the BAM file you're using as input (eg., did you run it through ReduceReads first?)

    David

  • AshuAshu Member
    edited November 2012

    Hi David,

    Here's the command I used to run BaseRecalibrator.

    java -jar GenomeAnalysisTK.jar -T BaseRecalibrator -I aligned_reads.bam -R reference.fasta -knownSites dbsnp.vcf -o recal.grp
    

    And, I did not know if I have to run it through ReduceReads? I have not used that before. The older versions of GATK were able to process my alignment files without using "ReduceReads".

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    No Ashu, you don't need to run ReduceReads. That is a new tool that can be used after processing, just before calling variants. I think David was asking because if you had used ReduceReads before BaseRecalibrator, that might explain your problem.

    So, did you use any other tools before running BaseRecalibrator? Like maybe IndelRealigner? Telling David all the steps you performed on your dataset will help him find out what the problem might be...

  • AshuAshu Member

    Dear David/Geraldine,

    I got these alignment(BAM) files(for a bacteria) from PacificBiosciences Secondary Analysis Software. They come with a mapping quality of 255. With the older version of GATK, I used to change the mapping quality to 60 using PrintReads and then do a quality recalibration on the alignment files and then call snps on these recalibrated files. The whole process worked fine. Now, I wanted to repeat this process with the new version of GATK as it says you can call indels with the new Unified Genotyper.

    Ashu

  • droazendroazen Cambridge, MAMember, Broadie, Dev ✭✭
    edited November 2012

    Hi Ashu,

    I believe your problem is with the CycleCovariate, which is very platform-dependent in the way it works. To test this hypothesis, could you please try re-running the BaseRecalibrator without the CycleCovariate and report back as to whether or not you get the same error? You can do this by adding the following options to your command line:

    --no_standard_covs -cov ContextCovariate

  • AshuAshu Member

    It worked, thank you so much. I wanted to ask you if running without Cycle Covariate would change things for me from the way they were working before with the old version of GATK.

  • ebanksebanks Broad InstituteMember, Broadie, Dev ✭✭✭✭

    Hi Ashu,

    This will be fixed in the next version of the GATK (2.3). For now you should run BQSR without the cycle covariate.

    That being said, you will have a very (very) tough time calling indels from PacBio data. The most significant error mode for PacBio is their high indel error rates and it makes accurate indel calling nearly impossible.

  • AshuAshu Member

    Ok...that means, no version of GATK can call indels for Pacbio as of now? Could you may be give me a timeline as to when would the new version that could call indels be released?

  • ebanksebanks Broad InstituteMember, Broadie, Dev ✭✭✭✭

    That's not exactly what it means. To be specific:

    • The GATK can call indels in PacBio data.
    • However, no version of any tool in the world can call indels in PacBio data well.
    • We have no intention of supporting development of indel calling in PacBio data.
  • AshuAshu Member

    Yes I know no other tool can call indels well from PacBio data.... I was hoping this new version could...but thanks for the help.
    I really appreciate it.

    Ashu

Sign In or Register to comment.