Question: problem in annotate my vcf file with another vcf file ??
0
gravatar for rastogi.suraj93
3.2 years ago by
rastogi.suraj9330 wrote:

hello every one, I am trying to annotate my vcf files which is generate using same reference genome for both vcf files but query reads is different. after alignment and SNP calling i got vcf file and count also variants. now i want to annotate these files to each other. i am using vcf annotate command like this

cat input.vcf.gz | vcf-annotate -a input2.vcf.gz -d column=HASH[0x2],key=INFO,ID=INDEL,Number=0,type=flag,Description=my annotate files -c CHROM,POS,ID,REF,ALT | bgzip -c > output.vcf

this command is working and dont show any error at window but output result i am not getting proper as a vcf format.i appreciate if can some one tell me how to access this command proper and get good result?? Thanking you

snp next-gen genome • 999 views
ADD COMMENTlink modified 2.9 years ago by Biostar ♦♦ 20 • written 3.2 years ago by rastogi.suraj9330

"this command is working and showing error" ? Can you rephrase the question ?

ADD REPLYlink written 3.2 years ago by geek_y9.9k

sorry, I again edit my question...please check it

ADD REPLYlink written 3.2 years ago by rastogi.suraj9330

Should you not be using zcat input.vcf.gz with a gzipped file?

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by genomax73k

sorry, i used but i forgot to mention there ..

ADD REPLYlink written 3.2 years ago by rastogi.suraj9330
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: 3327 users visited in the last hour