Question: How to make new VCF files myself?
gravatar for Amirosein
3.4 years ago by
UBC, Vancouver, CA.
Amirosein70 wrote:


I'm doing some computations on sanger seq. data.

for example i have a read, i aligned it to the reference genome and found an SNP, now i want to make a new VCF file for using in other tools like Annovar or SIFT, how to do this? (specially in Rstudio/Bioconductor)

"VcfR" is not published yet and i could not find any related function in "VariantAnnotation" package

Edit: finally i decide to write a function myself, it seems not hard, but first is there anybody wrote a function before?

thanks all

variantannotation variant vcf • 1.4k views
ADD COMMENTlink modified 3.3 years ago by chen1.9k • written 3.4 years ago by Amirosein70

It should not be difficult to make a VCF from scratch unless you want to make a multisample VCF. Can be done with simple Awk. What information do you have and how the file looks like ?

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by geek_y10k
gravatar for chen
3.3 years ago by
chen1.9k wrote:

Take a look at VCF spec doc,

Follow the spec, write the file directly by yourself. It is not difficult.

ADD COMMENTlink written 3.3 years ago by chen1.9k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 869 users visited in the last hour