This is a test version of Biostars. For the public version, visit https://www.biostars.org.
problem in annotate my vcf file with another vcf file ??

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

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

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

0 answers

No answers yet.

Log in to answer this question.