Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

FilterByOrientationBias in GATK4 hasn't filtered out any artefacts in FFPE sample

Hi,

I am currently using GATK4 to identify somatic mutations from FFPE WGS data. Everything else works fine apart from the second filtering step trough FilterByOrientationBias. Since the given sample is FFPE, I assume there should be some artifacts (C>T) being filtered out in this step. However, it doesn't filter any.

The command lines I used:

java -jar $picard CollectSequencingArtifactMetrics I=tumour.bam O=tumour_artifact FILE_EXTENSION=.txt R=$genome

gatk FilterByOrientationBias -AM 'G/T' -AM 'C/T' -V "$sample"_somatic_oncefiltered.vcf.gz -P tumour_artifact.pre_adapter_detail_metrics.txt -O "$sample"_twicefiltered.vcf.gz

Part of the summary of the filtering:

Sample Artifact_Mode Num_Artifact_Mode Num_Artifact_Mode_Filtered
Tumour C>T 11175 0
Normal C>T 0 0
Tumour G>T 3494 0
Normal G>T 0 0

As you can see, there are over 10K artifacts from FFPE and 3K from OxoG. But none of them have been filtered. Is it because the tool is still under development or is it because I did something wrong?

Looking forward to hearing from you soon!

Best wishes,
Qingli

Answers

Sign In or Register to comment.