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?