Question: Annovar - Mouse Exome Sequencing
0
gravatar for cbatist
2.3 years ago by
cbatist0
cbatist0 wrote:

Hi all,

Currently, I have been trying to annotate a mouse VCF file generated by Strelka using Annovar. I was able to download the gene annotations from UCSC using the following commands, and everything went fine until this step:

perl annotate_variation.pl -buildver mm10 -downdb -webfrom annovar refGene mouse_mm10/ perl annotate_variation.pl --buildver mm10 --downdb seq mouse_mm10/mm10_seq retrieve_seq_from_fasta.pl mouse_mm10/mm10_refGene.txt -seqdir mouse_mm10/mm10_seq -format refGene -outfile mouse_mm10/mm10_refGeneMrna.fa

When I tried to use table_annovar by running the following command, the job runs, but it generates an empty output.

table_annovar.pl/Users/caro mouse_mm10/ --vcfinput --outfile final --buildver mm10 --protocol refGene --operation g

Below is the message that is generated after I run the job. Would someone have any suggestions in what I am doing wrong here?

NOTICE: Running with system command <convert2annovar.pl -includeinfo="" -allsample="" -withfreq="" -format="" vcf4="" users="" carolinabatista="" desktop="" results_853="" strelka="" results="" passed.somaticteste.indels.vcf=""> final.avinput> NOTICE: Finished reading 1 lines from VCF file NOTICE: A total of 0 locus in VCF file passed QC threshold, representing 0 SNPs (0 transitions and 0 transversions) and 0 indels/substitutions NOTICE: Finished writing allele frequencies based on 0 SNP genotypes (0 transitions and 0 transversions) and 0 indels/substitutions for 0 samples NOTICE: Running with system command <table_annovar.pl final.avinput="" mouse_mm10="" -outfile="" final="" --buildver="" mm10="" --protocol="" refgene="" --operation="" g="" -otherinfo="" -nastring="" .="">

NOTICE: Processing operation=g protocol=refGene NOTICE: Running with system command <annotate_variation.pl -geneanno="" -buildver="" mm10="" -dbtype="" refgene="" -outfile="" final.refgene="" -exonsort="" final.avinput="" mouse_mm10=""/> NOTICE: Output files were written to final.refGene.variant_function, final.refGene.exonic_variant_function NOTICE: Reading gene annotation from mouse_mm10/mm10_refGene.txt ... Done with 36864 transcripts (including 5292 without coding sequence annotation) for 24540 unique genes

NOTICE: Multianno output file is written to final.mm10_multianno.txt NOTICE: Reading from final.mm10_multianno.txt

NOTICE: VCF output is written to final.mm10_multianno.vcf Use of uninitialized value $pre_anno_string in concatenation (.) or string at table_annovar.pl line 191, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1. Use of uninitialized value in join or string at table_annovar.pl line 192, <manno> line 1.

ADD COMMENTlink modified 5 months ago by Kevin Blighe48k • written 2.3 years ago by cbatist0
0
gravatar for Kevin Blighe
5 months ago by
Kevin Blighe48k
Kevin Blighe48k wrote:

First convert the VCF to ANNOVAR format with:

perl annovar/convert2annovar.pl -format vcf4 -includeinfo -allsample -withfreq -outfile out.annovar in.vcf
ADD COMMENTlink written 5 months ago by Kevin Blighe48k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1381 users visited in the last hour