How to make new VCF files myself?
1
0
Entering edit mode
7.8 years ago
Amirosein ▴ 70

Hi

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

VCF Variant VariantAnnotation • 5.9k views
ADD COMMENT
1
Entering edit mode

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 REPLY
1
Entering edit mode
7.7 years ago
chen ★ 2.5k

Take a look at VCF spec doc, https://samtools.github.io/hts-specs/VCFv4.2.pdf

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

ADD COMMENT

Login before adding your answer.

Traffic: 2148 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6