Hi, Please check the problem and the solutoin I apply and tell me if it is correct or not.
The problem:
I have a few variances (chr,position,ref_base,alt_base,gene) and also 50 vcf files. I need to check that if my variances exist in any of VCF files.
The solution:
First I use Annovar program to convert vcf files to Annovar Input file (convert2annovar command) and second I user Annovar again (table_annovar command with refGene,esp6500siv2_all,1000g2014oct_all and snp138 protocols) to generate a tab delimited file which contain VCF variances and Lastly I user linux "cat" command and pipe in into "grep" command to check if any of my variaces exist in that annovar output file.
please tell that if this is a common way for searching a varince in a vcf file ?
•
written
by
amir.jebelli •
1 answer
You might take a look at tools like bedtools, bedops, vcftools, or VariantAnnotation/GenomicRanges (Bioconductor).
•
written
by
Sean Davis
Thanks Sean, I use vcfintersect from vcflib to check if my variances exists in vcf file.
vcfintersect -b MY_INPUT_VARIANCES.bed MY_VCF_FILE.vcf.gz
Log in to answer this question.
-
Convert NCBI Downloaded files to ANNOVAR format
written by robert •
-
ANNOVAR Annotation of multiple VCF files error because of the absence of hg38_gnomad211_genome.txt …
written by Can Abdullah •
-
Annotated txt and csv files show empty columns in SIFT,Polyphen2, Mutation tester and other columns
written by Dhara Awasthi •
-
table_annovar.pl script shows error
written by Dhara Awasthi •
-
Modifying Annovar Command Line: table_annovar.pl
written by ghalib.bardai •
-
Is it possible to add all possible VCF annotations in one command using ANNOVAR ?
written by kay
-
Annovar using VCF file
written by Samar •
-
Error in bcf/vcf lifes
written by shakeelbiochemist •
-
count number of variants per gene from .vcf annovar&snpsift
written by User6891
-
Error When Convert Human Snp138.Txt To .Vcf Format
written by asphodel6 •
i think with VCFtools you can define and filter indels...az koja bioinformatics yad gerefti kalak?
I use vcfinterset from vcflib to solve the problem. I also find that I can filter vcf files with vcftools program.
4-5 mahi hast khodam shoro kardam be yad gereftam.
eyval dadash, carry on