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 https://software.broadinstitute.org/firecloud/documentation/freecredits

Picard CollectOxoGMetrics empty metrics

Zhenyu_ZhangZhenyu_Zhang University of ChicagoMember
edited February 2017 in Ask the GATK team

I tried to collect OxoG metrics, but less than 2% of my BAMs gave empty metrics as below. I have also updated my picard to 2.8.3. The example file I have is of normal size (14G). Also I saw from another thread asking about read group information, and it does.

@RG ID:0    CN:BCM  DS:Wilms Tumor PI R Gibbs   DT:2012-05-31T12:00:00+0000 LB:ING_WTSRG.PAKGZX01A-1_1pA    PL:Illumina PU:110208_SN601_0075_A81K6MABXX_6   SM:TARGET-50-PAKGZX-01A-01D
## htsjdk.samtools.metrics.StringHeader
# picard.analysis.CollectOxoGMetrics INPUT=TARGET-50-PAKGZX-01A-01D_hg19_Illumina_gdc_realn.bam OUTPUT=temp REFERENCE_SEQUENCE=/home/ubuntu/SCRATCH/GRCh38.d1.vd1.fa    MINIMUM_QUALITY_SCORE=20 MINIMUM_MAPPING_QUALITY=30 MINIMUM_INSERT_SIZE=60 MAXIMUM_INSERT_SIZE=600 INCLUDE_NON_PF_READS=true USE_OQ=tr
ue CONTEXT_SIZE=1 STOP_AFTER=2147483647 VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false GA4GH_CLIENT_SECRETS=client_secrets.json
## htsjdk.samtools.metrics.StringHeader
# Started on: Mon Feb 20 18:44:42 UTC 2017

## METRICS CLASS    picard.analysis.CollectOxoGMetrics$CpcgMetrics
SAMPLE_ALIAS    LIBRARY CONTEXT TOTAL_SITES TOTAL_BASES REF_NONOXO_BASES    REF_OXO_BASES   REF_TOTAL_BASES ALT_NONOXO_BASES    ALT_OXO_BASES   OXIDATION_ERROR_RATE    OXIDATION_Q C_REF_REF_BASES G_REF_REF_BASES C_REF_ALT_BASES G_REF_ALT_BASES C_REF_OXO_ERROR_RATE    C_REF_OXO_Q G_REF_OXO_ERROR_RATE    G_REF_OXO_Q
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   CCA 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   ACC 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   CCC 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   GCA 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   TCT 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   ACG 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   CCG 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   GCC 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   GCG 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   ACT 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   TCA 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   CCT 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   TCC 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   GCT 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   TCG 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
TARGET-50-PAKGZX-01A-01D    ING_WTSRG.PAKGZX01A-1_1pA   ACA 0   0   0   0   0   0   0   ?   -?  0   0   0   0   ?   ?   ?   ?
Post edited by shlee on

Issue · Github
by Sheila

Issue Number
1769
State
open
Last Updated
Assignee
Array
Milestone
Array

Answers

Sign In or Register to comment.