Correct Annovar input
Entering edit mode
15 months ago
storm1907 ▴ 30

Hello, I tried to convert my raw unfiltered Deepvariant vcf samples to Annovar format. First, I runned script with hg19 build, but got this error

WARNING: 1884 exonic SNPs have WRONG reference alleles specified in your input file!
WARNING: An example input line is <chr11        11964514        11964514        T       C       hom     54.5    171     chr11   11964514        .       T       C       54.5    PASS    .       GT:GQ:DP:AD:VAF:PL      1/1:54:171:0,171:1:54,62,0>
WARNING: ANNOVAR can still annotate exonic_variant_function for the mutation correctly!
WARNING: you may have used wrong -buildver, or specified incorrect reference allele, or used outdated mRNA FASTA file!

I switched to hg38 build, so I downloaded additional packages. Then I got this error message

A total of 27174 sequences cannot be found in humandb/hg38_refGeneMrna.fa

here is my command line

./ -format vcf4 $file -outfile $outfile -includeinfo -withzyg -withfreq

So which reference genome should I use then? Or there are certain file types I can only use for Annovar?

Thank you!

software error • 455 views

Login before adding your answer.

Traffic: 1616 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6