This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Saving file as .vcf

Hello,

I really need help to save my file as .vcf.

I have modified my vcf file as I want using VCFR package in R. Now I want to save my new in .vcf format using code write.vcf.

But the saved file is some kind of weird data file. Can someone help me where I am making mistake?

Thank youenter image description here

ngs variant vcf
• 2,192 views • link
updated by Ram Ram • written by Maryam • Maryam

from manual:

The function write.vcf takes an object of either class vcfR or chromR and writes the vcf data to a
vcf.gz file (gzipped text)
• 0 views • link
written by cpad0112 cpad0112

1 answer

This gibberish means that it is probably gzip-compressed and should be vcf.gz. Please don't get used to open OMICS files in plain text editors. Use Unix tools such as less, head or something like that.

• 0 views • link
written by ATpoint ATpoint

Rename your file to .vcf.gz and then try and see whether something like zless can properly open it. Frankly, and I really mean no offense, if you are stuck with these kinds of basics you must spend time with Unix fundamentals before diving into analysis.

• 0 views • link
written by ATpoint ATpoint

Thank you so much. I just uncompressed my files.

• 0 views • link
written by Maryam • Maryam

Why did you delete the post?

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this
  • Understanding IGV
    written by zizigolu zizigolu
    Hi I have targeted sequencing (WGS) I have uploaded bam and vcf file for sample; From these screenshot which one shows/support this mutation ![enter image …
  • singlecell annotation
    written by kayah • kayah
    I want to annotate my file to reference file(downloaded from single cell portal). I want to annotate cell_type! But don't know what packages to use …
  • help with weird PCA? (vcfR)
    written by MaeBH • MaeBH
    Hi everyone, very new to bioinformatics I have a SNP datasheet (.vcf) that I tried to make a PCA graph with using Rstudio (vcfR package) …
  • How can I create more accurate admixture plot like this?
    written by kgwkk2 • kgwkk2
    ![enter image description here][1] ![enter image description here][2] I made simple vcf file processing as ``` plink --vcf Asp_multicall2.vcf.gz --double-id --allow-extra-chr \ --set-missing-var-ids @:# \ …
  • How to modify VCF file?
    written by CHENG • CHENG
    Hi community, I have a question: the SNP position in vcf file is from GRCh37/hg19, I need to change the position to GRCh38. So, I …
  • Weird WGCNA plot
    written by synat.keam • synat.keam
    Dear Seniors and All members, I am wondering whether anyone has done weighted gene co-expression network analysis (WGCNA) from Bulk-RNAseq before. I followed the WGNCA …
  • Help needed in retrieving data from Candida genome database
    written by Trinh • Trinh
    Hi, I want to retrieve data from http://www.candidagenome.org/. For example with Candida albicans, I want to retrieve the data from all genes of Candida albicans …
  • SNP ID
    written by shamsurrahaman132 • shamsurrahaman132
    Hello Community, I'm Shamsur. I was dealing with variant calling. Is it possible to get snp id's while generating vcf file, like this image ![vcf …
  • intersectBed a gtf file with a vcf file
    written by Susmita Mandal Susmita Mandal
    Hey all, I have been trying to intersect a gtf file with a vcf file using intersectBed without using any options. I have edited my …
  • issue while converting vcf file to genind with vcfR
    written by Ayité K • Ayité K
    I am using vcfR for several analysis. It works for a previous vcf file. But actually I got the following message while converting a vcf …