Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

CollectSamErrorMetrics error

Answered
0

3 comments

  • Avatar
    Joe Georgeson

    Hi Anthony DiCi

    Thank you for the reply, after adding read groups this function works fine.

    Best,

    Joe

    1
    Comment actions Permalink
  • Avatar
    Anthony DiCi

    Hi Joe Georgeson,

    Thank you for writing to the GATK forum! I hope that we can help you sort this out.

    I brought this issue to our developers and received some feedback and things for you to try. The crash seems to have occurred in a part of Picard where it’s trying to get the read group ID. It could be because either the bam didn’t have it set or it was set to something Picard didn’t like. 

    You can verify whether there is a problem with your bam file by running our GATK tool ValidateSamFile. Assuming the issue is a missing read group, you may be able to repair it by using either --samtools addreplacerg or --picard AddOrReplaceReadGroups.

    I hope this helps! Please let me know if this leads you to success. If you have any other questions in the meantime, please do not hesitate to reach out.

    Best,
    Anthony

    0
    Comment actions Permalink
  • Avatar
    Anthony DiCi

    Hi Joe Georgeson,

    Great! Happy to hear that. Thank you again for contributing to the GATK forum. Please do not hesitate to reach out with any other questions that pop up in the future.

    Best,

    Anthony

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk