No sample in the VCF file!

Hi guys,

Can you please help me with this. I am working with trying to do:

vcf.fn <- "sequences_covid_SRA_april21_36genomes_sequence_aligned.trimmed.fasta.vcf" 
snpgdsVCF2GDS(vcf.fn, "sarscov2.gds", method="biallelic.only")

and I got this error:

No sample in the VCF file!

Which could be the reasons for this error? I check the file and it seems ok. Thank you in advance.

vcf vcfr

• 1,328 views • link

updated 5.3 years ago by Ram

• written 5.3 years ago by ebhmayra •

what is the output of

grep "#CHROM" -m1 sequences_covid_SRA_april21_36genomes_sequence_aligned.trimmed.fasta.vcf

If it is only #CHROM POS ID REF ALT QUAL FILTER INFO , then there is no sample in the vcf.

• 0 views • link

written 5.3 years ago by Pierre Lindenbaum

What is the output to the shell command:

bcftools query -l sequences_covid_SRA_april21_36genomes_sequence_aligned.trimmed.fasta.vcf

You'll need to have bcftools installed to run the above command. If it's not installed, give us the output to:

grep -m1 -E "^#CHROM" sequences_covid_SRA_april21_36genomes_sequence_aligned.trimmed.fasta.vcf | tr "\t" "\n"

• 0 views • link

written 5.3 years ago by Ram

0 answers

No answers yet.

More posts like this

VCF sample id format for snpgdsVCF2GDS (SNPRelate package)
written 2 months ago by shpak.max •

I need to convert vcfs into gds format as input for another software package, and I've been working with SNPRelate functions to do so. When …
Saving the output of LD pruning from SNPRelate package as a new GDS file
written 2.9 years ago by Patrick •

I successfully ran LD pruning from SNPRelate package and obtained the selected SNPs IDs. I would like to create a new GDS file with only …
Label groups with samples names in snpgdsDrawTree (SNPRelate package)
written 4.7 years ago by Begonia_pavonina •

Hello everyone, I have drawn a dendrogram with the snpgdsDrawTree from the SNPRelate R package. There are 8 samples in my dataset, and I wish …
.gds format for SAIGEgds
written 5.9 years ago by abyousaf •

Hello everyone, I am trying to use SAIGEgds. I originally had chromosome-separated VCF files that I combined into .gds format using [SNPRelate][1]'s `snpgdsVCF2GDS` and successfully …
Problem in LD-based SNP pruning using SNPRelate
written 6.1 years ago by evelyn

Hello, I have an expression data matrix and phenotype data. I started with SNPRelate to get the PCA and IBS analysis using: library(SNPRelate) library(gdsfmt) #Convert …
PCA plot color differentiation
written 6.5 years ago by tothepoint

I followed the post related to [PCA from vcf file][1] followed by library("SNPRelate") vcf.fn<-"~/xxx/tmp.vcf" snpgdsVCF2GDS(vcf.fn, "ccm.gds", method="biallelic.only") genofile <- openfn.gds("ccm.gds") ccm_pca<-snpgdsPCA(genofile) plot(ccm_pca$eigenvect[,1],ccm_pca$eigenvect[,2] ,col=as.numeric(substr(ccm_pca$sample, 1,3) == …
SNPRelate LD Pruning List of Loci
written 7.6 years ago by jessivannatta •

Hi there, I'm having trouble finding which SNPs have been pruned after LD pruning. I want to be able to list those loci that are …
Retrieve list of kept SNP's after LD pruning
written 8.3 years ago by leeandroid

Hi everyone! I followed the [SNPRelate package tutorial][1] in order to filter a VCF by LD and finally pursued PCA. My problem is that by …
file formatting for SNPRelate
written 10.8 years ago by lovestowell •

Hello! I am trying to use the R package SNPRelate to calculate pairwise relatedness between 8 males and 8 females. I have a VCF file …
converting vcf to gds losing some snps
written 11.3 years ago by gulcek •

I have vcf files of many individuals that I converted from 23andme raw data. I need to convert vcf files to gds files. I preprocessed …