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.

# UnifiedGenotyper missed indels calling on 100 base pair reads

Member Posts: 7

Hi, I'm using The Genome Analysis Toolkit (GATK) v2.3-9-ge5ebf34

I was running the UnifiedGenotyper tools on a 100bp paired end read with the the command below:

Program Args: -T UnifiedGenotyper -R ref.fasta -I realignedBwa.bam -glm BOTH -log variantsCalling.log -o bwa_variants.vcf

I noticed that only SNPs were called and not a single indel was called. however, when i used a 150bp paired end read with the similar command above it worked fine, the SNPs and Indels were called.

To make sure that the realignedBWA.bam file that i used was not corrupt/malformed, i used two different other pileup program on this bam file and the SNPs and Indels were called nicely. Hope you can take a look at this. I attached together the log file just in case

Hi there,

Can you tell me what happens if you set the -glm argument to INDEL?

Geraldine Van der Auwera, PhD

• Member Posts: 7

If the -glm argument is set up to INDEL, nothing is called up except for the 150bp which called up the indel

while if i set up the argument to BOTH as well as SNP, only the SNPs were called.

Hmm. Try using the --output_mode EMIT_ALL_SITES` argument and see if the indels get called and with what quality scores.

Geraldine Van der Auwera, PhD

Oh, scratch that, that option doesn't work for indels.

Have you tried calling it with the HaplotypeCaller? It is more advanced than the UG.

Geraldine Van der Auwera, PhD

• Member Posts: 7

I've been trying the HaplotypeCaller, and i'm not sure if it's my procedure/methods are wrong, but not a single SNPs or indel was called. the only output are the progress meter and no calling (for both 100bp and 150bp)

That's very odd. What kind of dataset are you running on?

Geraldine Van der Auwera, PhD

• Member Posts: 7

I simulated mutation read from a reference genome using DWGSIM, which produced reads (used the parameter to generate the read as if coming from illumina) which I align using several aligner including BWA..

• Member Posts: 7

Sorry for the long delay, it's holiday here.. I'll take a quick look and reply back, thanks a lot for helping out

Geraldine Van der Auwera, PhD

• Member Posts: 7

The base qualities are good enough, higher than the default minimum of 17 needed by the UG. A check on BQSR.log showed only 253 reads were filtered out from 20257 (1.25%) and samtools flagstat showed a 100% QC-passed reads.

Should I try to change any parameters involving the minimum base qualities because as for now i'm using the default.

• Member Posts: 7

ok, my bad.. it turns out the DP for indels is lower than the minimum thus making the indels being ommited.

Thanks for helping out pdexheimer and Geraldine.