A strange problem with annotation of vcf file
1
1
Entering edit mode
3.7 years ago
seta ★ 1.6k

Hi friends,

I did variant calling of a bam file (it's related to whole-genome sequencing of a human sample) with GATK (version 4) and then tried to annotate it with dbsnp, but it was not successful. I download human reference genome and dbsnp files from the following addresses:

For annotation with dbsnp, firstly, I used the below command from GATK

Gatk VariantAnnotator –R genom.fa –I file1.bam –V file1.vcf –A coverage –D All_20180418.vcf

GATK4 said

Warning: VariantAnnotator is a BETA tool and is not yet ready for use in production

and didn’t return me any output. So, I used bcftools for annotation with the below command:

bcftools annotate –c ID –a All_20180418.vcf.gz file1.vcf.gz > file1_annotate.vcf

But, bcftools also didn’t annotate my vcf file without returning any error.

Could you please help me what’s wrong and kindly tell me how I can annotate my vcf file with dbsnp?

Many thanks

genome annotation gatk4 bcftools vcf • 2.1k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
3
Entering edit mode
3.7 years ago

SnpSift Annotate does a good job of annotating the ID field of a VCF. You will have to download the dbSNP VCF, though, which is > 10GB last time that I checked.

ADD COMMENT
0
Entering edit mode

Yes, I downloaded dbsnp vcf file with 15.2 GB in gzip. Could you please tell me what is wrong in my work to annotate the vcf file?

ADD REPLY
1
Entering edit mode

I cannot comment on the GATK command because I have never used it. The warning message issue by that command speaks for itself, too.

Regarding BCFtools, you likely have to first unset your VCF IDs with

bcftools annotate --remove ID file1.vcf.gz | \
bcftools annotate -c ID -a All_20180418.vcf.gz file1_annotate.vcf -Ov \
> file1_annotate.vcf

Of course, SnpSift also works.

Nota bene: also check that your contig names are the same as those in the dbSNP VCF. One may have a 'chr' prefix, while, the other, not

ADD REPLY
0
Entering edit mode

Many thanks for your points. I'll try them

ADD REPLY
0
Entering edit mode

Good luck my friend.

ADD REPLY

Login before adding your answer.

Traffic: 693 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