This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract Dp From The Vcf File

I want to extract some information from the info column of my VCF file, when I load my VCF file, I do not have access to my depth (DP) and number of alleles (AN), however, the NS(number of samples for this allele works fine), DP and AN shown as "NA". My vcf file clearly show that neither DP nor AN is empty, here is just one line example of my vcf file, thanks!

gi|393925858|gb|AGTA02071966.1|    0000000739    .    G    A    121.20    PASS    NS=74:AN=2:DP=8448    GT:DP:GQ:EC:SG    0/1:262:144:116:R

> vcf <- readVcf("z.vcf", "Genome")
> hdr <- exptData(vcf)[["header"]]
>  info(hdr)

DataFrame with 3 rows and 3 columns
        Number        Type                     Description
   <character> <character>                     <character>
NS           1     Integer     Number of Samples With Data
DP           1     Integer                     Total Depth
AN           1     Integer Number of Alleles in Population

>   info(vcf)
DataFrame with 2648 rows and 3 columns
                                                  NS        DP        AN
                                           <integer> <integer> <integer>
gi|393925858|gb|AGTA02071966.1|:0000000739        74        NA        NA
gi|393925858|gb|AGTA02071966.1|:0000000781        74        NA        NA
gi|393925983|gb|AGTA02071903.1|:0000000957        74        NA        NA
gi|393925983|gb|AGTA02071903.1|:0000000960        74        NA        NA
gi|393925983|gb|AGTA02071903.1|:0000001007        73        NA        NA
...                                              ...       ...       ...
r vcf
• 4,740 views • link
written by mad.cichlids mad.cichlids

Can you be more specific? How are you reading the files? Bioconductor? Package name? Version?

• 0 views • link
written by JC JC

Thanks. I used vcf <- readVcf("z.vcf", "Genome"), z.vcf is my vcf file, "Genome" is the folder of my indexed ref genome. The package is variantannotation in Bioconductor. Here is how i installed it:

source("http://bioconductor.org/biocLite.R")
    biocLite("VariantAnnotation")

According to the archive in the Bioconductor, this should be Version1.8.13. Please let me know if you need additional information that I can provide, I really appreciate your comment. And I am running R in ubuntu

• 0 views • link
written by mad.cichlids mad.cichlids

1 answer

Thanks to Valerie Obenchain's help, it turned out that it is the colon sign in the INFO column caused this, after replacing the colon with semicolon, everything works fine.

• 0 views • link
written by mad.cichlids mad.cichlids

Log in to answer this question.

More posts like this
  • Extract vcf file header and save it in Key-Value csv
    written by tinyskinn • tinyskinn
    Hi everyone, If I have a vcf file with header like this: ``` ##fileformat=VCFv4.0 ##fileDate=20090805 ##source=myImputationProgramV3.1 ##reference=1000GenomesPilot-NCBI36 ##phasing=partial ##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data"> ##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth"> …
  • Get genotype GT from multi-sample VCF (per chromosome)
    written by zx8754 zx8754
    Example input: multi-sample VCF (adapted from [www.internationalgenome.org][1]): **Note:** my actual file is bgzipped and tabixed with ~2Mln variants (rows) and ~1000 samples (columns). ##fileformat=VCFv4.0 ##fileDate=20090805 …
  • Index Vcf File
    written by mad.cichlids mad.cichlids
    <p>Hi, I was trying to index my vcf file, i first sort it and then zip it before indexing as this post: <a href='/p/78444/'>Tabix -p …
  • Vcf File From Mpileup ... What Next
    written by Dataminer Dataminer
    <p>Hi,</p> <p>I have done variant calling on my file using <a href='http://samtools.sourceforge.net/'>samtools</a> mpile up and I have converetd my file from bcf to vcf. </p> …
  • How To Find The Most Frequent Allel Snp In Vcf File
    written by Chen Sun Chen Sun
    <p>When I use <code>samtools mpileup</code>, there are usually several allel SNP in ALT column according to REF column. just as follows:</p> <pre><code>##fileformat=VCFv4.1 ##samtoolsVersion=0.1.19-44428cd ##reference=file:///gpfs/home/cxs1031/backup/ref/chr22.fa ##contig=&lt;ID=chr22,length=51304566&gt; …
  • Variant Effect Predictor On Domain Level
    written by jackuser1979 jackuser1979
    <p>I have non-model organism and I have annotated SNPs using Variant effect Predictor with VCF file using Ensembl assembly and I got the list of …
  • Chrx Genotype Counts In Esp6500 Exomes
    written by Caddymob Caddymob
    <p>Does anyone know the format of how the NHLBI reports the genotypes for chromosome X varriants? If we look at the example <code>rs6525447</code> we get …
  • Error Loading Vcf File In Ucsc Genome Browser
    written by win win
    <p>Hi all, I generated a vcf file using samtools and the wanted to load it into UCSC Genome Browser as suggested at this link</p> <p><a …
  • Importing Varscan Vcf To Igv
    written by jackuser1979 jackuser1979
    Did any one tried importing varscan variant caller ( V.2.3.2) VCF output to IGV (V 2.1.2.4)?. I tried indexing VCF file using igvtools-index & Tabix …
  • Snp Call Within A Given Region And Vcf Result
    written by Zhshqzyc • Zhshqzyc
    <p>Hello, I used the command to call snp. I hope to find something.</p> <pre><code>samtools mpileup -C50 -r chr2:1000000-2000000 -uf ref.fa s1.bam s2.bam s3.bam s4.bam s5.bam …