Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • ruchi1621
    Junior Member
    • Oct 2013
    • 9

    #1

    Error using vcf-stats

    I am trying to figure out the variant statistics from a vcf file using the command:
    Code:
     ./vcf-stats NA11881roi_variants.vcf
    For which I am getting the following output:

    Code:
    The version "4.1" not supported, assuming VCFv4.0
    Expected GT as the first genotype field at chr4:46036725
    Use of uninitialized value $gtype in pattern match (m//) at /N/dc/scratch/bansalr/bowtie2-2.1.0/vcftools_0.1.4a/perl/Vcf.pm line 1060, <__ANONIO__> line 28.
    Use of uninitialized value $gtype in concatenation (.) or string at /N/dc/scratch/bansalr/bowtie2-2.1.0/vcftools_0.1.4a/perl/Vcf.pm line 1060, <__ANONIO__> line 28.
    Could not parse gtype string []
    
     at /N/dc/scratch/bansalr/bowtie2-2.1.0/vcftools_0.1.4a/perl/Vcf.pm line 167, <__ANONIO__> line 28.
            Vcf::throw('Vcf4_0=HASH(0x260f630)', 'Could not parse gtype string []\x{a}') called at /N/dc/scratch/bansalr/bowtie2-2.1.0/vcftools_0.1.4a/perl/Vcf.pm line 1060
            VcfReader::parse_alleles('Vcf4_0=HASH(0x260f630)', 'HASH(0x25221d8)', 'NA11881roi.sorted.bam') called at ./vcf-stats line 83
            main::vcf_stats('HASH(0x2237de8)') called at ./vcf-stats line 9
    Can somebody please help me regarding this? I am new to this field

    Thank you!
  • yl01
    Member
    • Aug 2012
    • 25

    #2
    It looks like the vcftools you are using is of version 0.1.4a. So the first thing i would try is to install the updated version, which is 0.1.11.

    Comment

    • bliu1
      Junior Member
      • Feb 2011
      • 3

      #3
      VCFtools error

      Originally posted by yl01 View Post
      It looks like the vcftools you are using is of version 0.1.4a. So the first thing i would try is to install the updated version, which is 0.1.11.
      Hi, I have same error when I run vcf-isec even I installed the updated version.

      Use of uninitialized value $attr_value in pattern match (m//) at /home/../vcftools_0.1.11/perl/Vcf.pm line 2971, <__ANONIO__> line 6.
      at /home/bliu1/vcftools_0.1.11/perl/Vcf.pm line 2971.

      Could you help? thanks

      Comment

      • Nnogueira
        Junior Member
        • Oct 2011
        • 9

        #4
        I've just started using bcftools stats, so I still don't know much, but for me, I get an error when I don't give a bgzip compressed file, could that be your case? You can get the bgzip in the htslib package, after that, you need to index the bgzip file with bcftools index. Do you think that would help?

        By the way, just getting a ride on your question, does anyone know where to find the documentation for the bcftools stats output? I've got some weird results that I would like to understand better. I've asked the same question here: http://seqanswers.com/forums/showthread.php?t=57687

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
          by SEQadmin2



          Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

          Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
          ...
          10-14-2026, 09:55 AM
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Yesterday, 11:13 AM
        0 responses
        16 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, Yesterday, 10:48 AM
        0 responses
        16 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        35 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        42 views
        0 reactions
        Last Post SEQadmin2  
        Working...