Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Haplotype calls from Amplicon Sequencing empyrean Bioinformatics 3 06-25-2012 02:27 PM
dubious SNP calls from GATK caswater Bioinformatics 0 04-24-2012 05:21 PM
Complete Genomics Variant Calls quicksand21 Bioinformatics 2 10-11-2011 07:21 AM
None reference heterozygous calls Chiel Bioinformatics 0 01-03-2011 04:34 AM
Hi all! [SNP calls from Illumina w/o reference] gdilla Introductions 0 08-13-2010 12:29 PM

Thread Tools
Old 07-19-2012, 11:05 AM   #1
Location: Cupertino, CA

Join Date: Nov 2011
Posts: 59
Default Homozygous Ref calls and no calls

Hi all, I was wondering if someone could help.

I am generating a vcf file based on the following command chain.

samtools sort eg2.bam eg2.sorted.bam samtools mpileup -uf referencegenomefile eg2.sorted.bam | bcftools view -bvcg - > eg2.raw.bcf bcftools view eg2.raw.bcf | perl /usr/share/samtools/ varFilter -D100 > results.vcf

What is being observed is that the VCF file only contains "1/1" or "0/1" or "1/0" calls. I don't see a single "0/0" call.

In bcftools i have tried to eliminate a the "v" parameter but then i dont get any genotype calls.

What I am really after is: is there a way that i can generate another VCF file (or include in the same file) that also has "0/0" genotypes. Also, is is possible to include no call, i.e. a way to indicate if a site was not sequenced at all?

Thanks in advance. A
ashkot is offline   Reply With Quote
Old 08-09-2012, 10:01 AM   #2
Location: Tennessee

Join Date: Jun 2012
Posts: 10

You definitely don't want -v in your bcftools view call, it hides those invariant sites.

samtools mpileup -uf referencegenomefile eg2.sorted.bam | bcftools view -cg - | perl /usr/share/samtools/ varFilter -D100 > results.vcf

When you do this, though, they homozygous ref aren't written as 0/0, they appear as follows:
21 9411856 . G . 60 . DP=61;VDB=0.0124;AF1=0;AC1=0;DP4=19,30,0,12;MQ=39;FQ=-57;PV4=0.012,0.041,2.2e-06,1 PL 0
21 9411868 . C . 71 . DP=65;VDB=0.0124;AF1=0;AC1=0;DP4=23,30,0,12;MQ=39;FQ=-68;PV4=0.0053,1.6e-18,1.2e-06,1 PL 0
21 9411871 . G . 65 . DP=62;VDB=0.0124;AF1=0;AC1=0;DP4=22,28,0,12;MQ=39;FQ=-62;PV4=0.0051,5.2e-26,1.1e-06,1 PL 0
21 9411931 . G . 79 . DP=42;VDB=0.0003;AF1=0;AC1=0;DP4=18,20,0,4;MQ=48;FQ=-76;PV4=0.12,1,0.00024,0.24 PL 0
etorstenson is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:12 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO