Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • elfuser
    Member
    • Dec 2009
    • 16

    #1

    Understanding coverage in VCF

    Hello and sorry if this question has been asked before or it is super obvious.
    I'm having a bit of a hard time finding out what the numbers mean in the vcf pileup.
    I created an alignment BAM file of my illumina paired-end reads to a bacterial genome using the following:

    Code:
    bwa aln ref.fa r1.fq > r1.sai
    bwa aln ref.fa r2.fq > r2.sai
    bwa sampe ref.fa r1.sai r2.sai r1.fq r2.fq > aln.sam
    samtools view -bS aln.sam > aln.bam
    Then I create the pileup with
    Code:
    samtools mpileup -uf aln.bam > pile.bcf
    bcftools view pile.bcf > pile.vcf
    I also create a depth file with
    Code:
    samtools depth aln.bam > depth.txt
    Let's take position 214 in the reference genome
    In the vcf
    Code:
    gi|218888746|ref|NC_011770.1|	214	.	G	C	202	.	DP=45;AF1=1;AC1=2;DP4=0,0,25,13;MQ=55;FQ=-141	GT:PL:GQ	1/1:235,114,0:99
    In the depth
    Code:
    gi|218888746|ref|NC_011770.1|	214	53
    I understand that the DP(45) value in the vcf file refers to the raw coverage of the given position with high quality(q13?). An that the coverage number in the depth file (53) refers to the number of bases called at that given position regardless of quality and/or conflict. I also understand that the DP4 values (0,0,25,13) reflect if the reads that cover that position match the reference or an alternative variant. In this case 38 (25+13) high quality reads match the alternative variant and 0 the reference. If this is right what happened to the other 7 (45-38) high quality reads that cover that position.

    Another doubt I have is when I find multiple genotypes at a position. For example:
    In the VCF
    Code:
    gi|218888746|ref|NC_011770.1|	3081	.	C	T,G	213	.	DP=70;AF1=1;AC1=2;DP4=0,0,31,32;MQ=50;FQ=-214	GT:PL:GQ	1/1:246,187,0,247,170,244:99
    My understanding is that in this case that 63 reads support either a T or a G at this position. So what I am trying to understand is how I can tell the support for each genotype. I tried looking at the GT:PL:GQ, but I dont quite understand it.
    Thanks!

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...