Go Back   SEQanswers > Bioinformatics > Bioinformatics
Similar Threads
Thread Thread Starter Forum Replies Last Post
what is DP4 in VCF (Variant Call Format) and AF1? gcrdb Bioinformatics 8 01-20-2015 09:37 AM
strand bias variant call yinshe Bioinformatics 2 10-10-2011 06:01 AM
Variant call using BAM from Tophat emilyjia2000 Bioinformatics 5 08-01-2011 04:37 PM
ERANGE to call peaks from ChIP-Seq data sonja Bioinformatics 3 08-31-2010 08:47 AM
Using Galaxy to Call Peaks? Giles Bioinformatics 4 02-05-2010 10:43 AM

Thread Tools
Old 07-31-2012, 11:42 AM   #1
Location: Toronto, ON

Join Date: Dec 2009
Posts: 16
Default SNP quality distribution peaks at 222 from variant call pile

I have reference mapped paired end illumina reads and called variants using BWA and Samtools respectively. The resulting vcf was treated to remove high coverage SNPs with
Code: varFilter -D30
and then filtered for low quality SNPs using awk
I graphed the distribution of SNP quality and observed a huge peak at 222., I repeated it with other samples and observed the same peak. Any clues as to why I may be seeing this?
Attached Images
File Type: png Picture1.png (22.4 KB, 10 views)

Last edited by elfuser; 07-31-2012 at 12:02 PM.
elfuser is offline   Reply With Quote

samtools vcf bwa

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 10:02 PM.

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