Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • alternative to vcftools for vcf v4.1

    Hi everyone,

    I need to retrieve some basic information from a vcf file of whole genome reads mapped to a reference, such as coverage, nucleotide diversity, etc.. I found vcftools to do all I needed, but my vcf files are in version 4.1 which is not supported by the vcftools binary executable, and I can't seem to find a linux alternative for this. Anyone has any suggestions?

    Thanks!

  • #2
    The latest versions of VCFtools do support v4.1 - what version are you using? I suggest you update.

    Comment


    • #3
      Hi, thanks for your help. I'm using version 0.1.9. To be fair though, the problem is that I tried to retrieve average depth from my vcf file with the "vcftools input.vcf --depth" command and although everything looked fine with the analysis, the output had "nan" in the column when I expected the average coverage do be. I the assumed this might be due to the lack of support for vcf v4.1 (as I understand it is only supported in the perl scripts and not in the binary executable), but I might be wrong of course.

      Comment


      • #4
        Does your VCF file contain depth information (i.e. is there a a DP flag in the FORMAT fields)? If not, it's unlikely that VCFtools will be able to calculate depth.

        Comment


        • #5
          Yes, my vcf files have the DP field in the INFO column. That's why I got confused and assumed t was a problem of version incompatibility.

          Comment


          • #6
            GATK has a whole suite of walkers that can get info from VCF files.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Genetic Variation in Immunogenetics and Antibody Diversity
              by seqadmin



              The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
              11-06-2024, 07:24 PM
            • seqadmin
              Choosing Between NGS and qPCR
              by seqadmin



              Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
              10-18-2024, 07:11 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 11-08-2024, 11:09 AM
            0 responses
            59 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 11-08-2024, 06:13 AM
            0 responses
            38 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 11-01-2024, 06:09 AM
            0 responses
            35 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-30-2024, 05:31 AM
            0 responses
            23 views
            0 likes
            Last Post seqadmin  
            Working...
            X