GATK Showcase in Terra

Check out these fully configured workspaces to test drive the Best Practices pipelines and workshop tutorials with zero installation required!

Get email notifications!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements on the blog, by following the instructions given here.

Got a problem?

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

StrandArtifect fields missed in normal sample in Mutect2 output

I used the latest version of Mutect2 to call somatic variants by comparing a tumor sample versus a normal sample. There is inconsistency of the fields between columns "FORMAT", "NORMAL" and "TUMOR". Taking one variant as an example:

The last two fields "SA_MAP_AF" and "SA_POST_PROB" were calculated for only tumor sample (last column), but missed in normal sample (10th column), which resulted in discrepancy between FORMAT column and the normal sample column. Is it true that Mutect2 only reports them for tumor sample? If so, it'd be best to let Mutect2 print missing values (eg. a period) to fill the last two fields for normal sample.