GATK no "SB" annotation after GenotypeGVCFs
0
1
Entering edit mode
13 months ago
wang.yiguan ▴ 10

Hi,

I run the following command to generate a gvcf file.

gatk HaplotypeCaller -R dmel-all-chromosome-r6.36.fasta \
    -I sample_43.sort.noDup.realign.bam \
    -ERC GVCF \
    -L 4 \
    -O sample_43_chr4_test.g.vcf.gz

This gvcf file contains the StrandBiasBySample or SB annotation. And then I run GenotypeGVCFs as following:

gatk GenotypeGVCFs -R dmel-all-chromosome-r6.36.fasta \
    -V sample_43__chr4_test.g.vcf.gz \
    -A StrandBiasBySample \
    -O sample_43_chr4_test.vcf.gz

The SB annotation disappeared in the final output vcf file sample_43_chr4_test.vcf.gz ! I want to keep the SB annotation in the final vcf.

I was trying to google a solution, but failed.

Thanks in advance!

gatk strandbias • 934 views
ADD COMMENT
0
Entering edit mode

Hi,

There is a post on GATK forum that may help: Missing annotations in the output callset VCF.

Have you tried to remove the "-A StrandBiasBySample" option from your GenotypeGVCFs command?

ADD REPLY
0
Entering edit mode

Yes, I've tried to remove -A StrandBiasBySample, but still no SB annotation in the output vcf.

ADD REPLY

Login before adding your answer.

Traffic: 1096 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6