problem in annotate my vcf file with another vcf file ??

0

Entering edit mode

8.3 years ago

rastogi.suraj93 ▴ 30

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 • 2.3k views

ADD COMMENT • link updated 8.0 years ago by Biostar 20 • written 8.3 years ago by rastogi.suraj93 ▴ 30

0

Entering edit mode

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

ADD REPLY • link 8.3 years ago by GouthamAtla 12k

0

Entering edit mode

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

ADD REPLY • link 8.3 years ago by rastogi.suraj93 ▴ 30

0

Entering edit mode

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

ADD REPLY • link 8.3 years ago by GenoMax 147k

0

Entering edit mode

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

ADD REPLY • link 8.3 years ago by rastogi.suraj93 ▴ 30

Login before adding your answer.