Variant Sites Not Listed In Vcf
hi all hope someone can help. The .vcf file generated after SNP calling contains variants with genotypes of 0/1 or 1/0 or 1/1.
does this mean that variants not detected in the vcf are homozygous wildtype.
is there a way to extract all such variants from the bam file an generate a file that has the QC score?
thanks in advance. ashwin
1 answer
use samtools mpileup to call all position in your BAM -> http://samtools.sourceforge.net/samtools.shtml
Log in to answer this question.
More posts like this
-
how to extract variants from the vcf.gz files linked below?
written by Mahad •
-
How To Get Count Of 0|0 From A Gvcf
written by win
-
Can Igv View Blast+ Output File?
written by shl198
-
Removing Phix Data From Ngs Data For Analysis
written by crivenster
-
Can I Use The Gatk Snp Calling On Sam Files With Missing Quality Values
written by Varun Gupta
-
How To Extract Consensus Sequence Of Mapped Regions From A Bam File?
written by lin.barnum •
-
How To Interpret The Flag Field In The Sam Specification
written by deepthithomaskannan
-
Samtools And Region List
written by win
-
Create A Tbi File From Vcf
written by win
-
Filtering A Sam File For Quality Scores
written by Varun Gupta
Ashwin -- what method are you using to call variants? Most callers have an option to enable reporting of confident reference calls as well. For instance, in GATK UnifiedGenotyper use
EMIT_ALL_CONFIDENT_SITEShttp://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_genotyper_UnifiedGenotyper.html#--output_mode