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!

Did you remember to?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
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.

Did we ask for a bug report?


Then follow instructions in Article#1894.

Formatting tip!


Surround blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block.
Powered by Vanilla. Made with Bootstrap.
Picard 2.9.0 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

Odd variant call quality distribution

BlueBlue Member Posts: 48

Basically I have an odd-looking distribution of my variant quality scores (see attached png), and was wondering how concerned should I be and how can I rectify it.

The input data from the graph is from UnifiedGenotyper vcf output file, QUAL values.
The four samples in the vcf file are one Drosophila reference line, and three more which are outcrosses of the reference line and thus are heterozygous for the reference allele.

My fastq read-mapping pipeline includes adapter and low-quality base removal, and local re-alignment. I've also attached a pdf showing read quality distribution from one of the samples which also looks a bit odd.

qualscoredistr.png
1261 x 804 - 67K
pdf
pdf
TH1_CRSNR.multiplemetrics.quality_distribution.pdf
5K

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Administrator, Dev Posts: 11,163 admin

    Unfortunately we don't have the resources to provide detailed troubleshooting of your results, sorry. If you think the distributions look odd, perhaps you can look in more detail at subsets of variants from different tranches in IGV, to get an idea of whether they look real or not. You may need to play around with parameters of variant recalibration to find the settings that suit your data best. Good luck!

    Geraldine Van der Auwera, PhD

  • BlueBlue Member Posts: 48

    This QUAL distribution anomaly exists at a similar level in all chromosomes.
    I'm not performing variant recalibration because most of the variants are likely to be novel.

    In the UnifiedGenotyper, what metrics go into the calculation which generates the QUAL value for each variant ?

    I cannot see this information in the documentation.

    I assume that total depth, 'DP' is one of them.

    My DP distribution is normal so there must be something else affecting my QUAL.

Sign In or Register to comment.