SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SNP Calling for 454 reads empyrean Bioinformatics 1 11-17-2011 07:23 AM
Samtools mpileup. Which reads are excluded during SNP calling? wwmm933 Bioinformatics 5 05-27-2011 08:12 AM
PubMed: Control-free calling of copy number alterations in deep-sequencing data using Newsbot! Literature Watch 0 04-08-2011 01:10 AM
SNP calling elondin Bioinformatics 3 01-24-2011 09:59 AM
What are the SNP calling programs? foxyg Bioinformatics 2 10-08-2010 05:19 AM

Reply
 
Thread Tools
Old 11-24-2011, 08:36 PM   #1
wanguan2000
Member
 
Location: shanghai

Join Date: Nov 2010
Posts: 24
Default SNP calling reads number

I followed the GATK pipeline.
when I use the IGV to see the bam file I found the question:
position chr17:7578524 : G(51 reads)-A(20 reads) I think this is a snp, but the position doesn't exist the result VCF files. GATK don't think this is a snp,even in lowQual snp.
why? the position have so high sequencing deepth.Who can explain it
wanguan2000 is offline   Reply With Quote
Old 11-25-2011, 08:16 AM   #2
Robby
Member
 
Location: Germany

Join Date: Mar 2011
Posts: 68
Default

That could have many reasons. How are the mapping qualities of the reads at that position and how are the quality scores of the bp supporting the SNP?
Robby is offline   Reply With Quote
Reply

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 03:41 PM.


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