Unconfigured Ad

Collapse
X
Collapse
 
  • Page of 1
  • Time
  • Show
Clear All
new posts
Previous template Next
  • josk
    Junior Member
    • Oct 2011
    • 1
    #1

    Insertions follow reference in vcf file

    Hi,
    I'm working with gvcf files, containing also non-variant positions. I have noticed some insertions that I'm not able to correctly interpret or handle. An example (simplified):
    Code:
    NC_013991.2 1403 . G GAA 176.32 . GT:AD:DP:GQ:PGT:PID:PL:PS ./. 0/0 0/0 ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. 1|1 0/0 ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. 0/0 ./. ./. 1/1 ./. ./. 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 ./. 0/0 0/0 0/0 0/0 ./. 0/0 ./. 0/0 ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./.
NC_013991.2 1404 . A . . . GT:AD:DP:RGQ ./. 0/0 0/0 0/0 ./. ./. ./. ./. ./. ./. ./. 0/0 ./. ./. ./. ./. 0/0 ./. 0/0 0/0 0/0 ./. ./. ./. ./. ./. ./. 0/0 ./. ./. ./. ./. ./. ./. 0/0 0/0 0/0 0/0 ./. ./. 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 ./. 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 ./. 0/0 0/0 0/0 0/0 ./. ./. ./. 0/0 ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. 0/0 ./. ./. ./.
NC_013991.2 1405 . A . . . GT:AD:DP:RGQ ./. 0/0 0/0 0/0 ./. ./. ./. ./. ./. ./. ./. 0/0 ./. ./. ./. ./. 0/0 ./. 0/0 0/0 0/0 ./. ./. ./. ./. ./. ./. 0/0 ./. ./. ./. ./. ./. ./. 0/0 0/0 0/0 0/0 ./. ./. 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 ./. 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 ./. 0/0 0/0 0/0 0/0 ./. ./. ./. 0/0 ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. ./. 0/0 ./. ./. ./.
    Here, the reference sequence is "GAA". Why, at site 1403, do I get an alternative allele "GAA"? There seem to be two samples at this position that are "1/1", thus homozygous for GAA. But they are also homozygous for A at the next positions. What thus is their genotype ? GAAAA, thus containing an insertion of AA?
    I've seen multiple cases like this, also with longer alternative alleles, and each time the alternative allele sequence exactly recapitulates the reference sequence. This seems like an error to me, but how to correct?
    The mapping was done either with bwa aln or bwa mem, and we've observed this behaviour both when using GATK4 and bcftools for genotyping.

    I would be very grateful if someone could clarify this.
    Thanks,
    Jos
    Tags: bcftools, bwa, gatk, mapping, vcf
    Previous template Next

    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
    View All

    ad_right_rmr

    Collapse

    News

    Collapse
    Topics Statistics Last Post
    Exeter-Led Study Maps What an Average Human Genome Looks Like by SEQadmin2
    Started by SEQadmin2, 10-22-2026, 11:13 AM
    0 responses
    35 views
    0 reactions
    		 Last Post SEQadmin2
    by SEQadmin2
    10-22-2026, 11:13 AM  
    Long-Read Sequencing Uncovers Cancer Mutations Conventional Testing Missed by SEQadmin2
    Started by SEQadmin2, 10-22-2026, 10:48 AM
    0 responses
    25 views
    0 reactions
    		 Last Post SEQadmin2
    by SEQadmin2
    10-22-2026, 10:48 AM  
    Mayo Clinic-Led Study Detects DNA Mutations Across Noncancerous Tissues by SEQadmin2
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    41 views
    0 reactions
    		 Last Post SEQadmin2
    by SEQadmin2
    10-15-2026, 10:52 AM  
    Study Finds Chromatin Motion Holds Steady Until Cells Face Stress by SEQadmin2
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    45 views
    0 reactions
    		 Last Post SEQadmin2
    by SEQadmin2
    10-14-2026, 05:39 AM  
    View All
    Working...