This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Generate head on vcf file from freebayes

well after medhat told me my vcf doesn't have header I try to figure out how to get it, my vcf file I got it with two ways

bin/freebayes --fasta-reference /media/cri/CRIS_DATA/ULT/OTRO2/otro/hg19.fa /media/cri/CRIS_DATA/ULT/OTRO2/otro/input.header.realigned.bam >/media/cri/CRIS_DATA/ULT/OTRO2/otro/freebayessample3.vcf



bin/freebayes --min-coverage 10 -V -f /media/cri/CRIS_DATA/ULT/OTRO2/otro/hg19.fa /media/cri/CRIS_DATA/ULT/OTRO2/otro/input.header.realigned.bam >/media/cri/CRIS_DATA/ULT/OTRO2/otro/freebayessample.vcf

Is any command to add to get my head on the file?

Thanks

vcf
• 2,042 views • link
updated by GenoMax GenoMax • written by cristina_sabiers cristina_sabiers

FreeBayes outputs VCF headers by default (at least, the current version does). Are you sure the VCF with missing header is direct output from FreeBayes, or after additional processing (e.g., BEDtools 'subtract' removes headers)?

• 1 views • link
written by harold.smith.tarheel harold.smith.tarheel

from the documentation page the command is

freebayes -f ref.fa aln.bam >var.vcf

and it should have the header by default what was the command you used for the alignment and is it was to the same reference

/media/cri/CRIS_DATA/ULT/OTRO2/otro/hg19.fa

• 1 views • link
written by Medhat Medhat

0 answers

No answers yet.

Log in to answer this question.

More posts like this