Changes

From Genome Analysis Wiki
Jump to navigationJump to search
no edit summary
Line 88: Line 88:     
; If Nuclear Families are Available, Assess Mendelian Consistency
 
; If Nuclear Families are Available, Assess Mendelian Consistency
: When reporting rates of Mendelian inconsistencies, report these not on a per site basis but, instead, based on the number of sites with at least one non-reference call in the trio.
+
: When reporting rates of Mendelian inconsistencies, report these not as a fraction of all sites, but as a fraction of sites with at least one non-reference call in the trio.
    
== Variant Filters ==
 
== Variant Filters ==
   −
Initial sets of SNP calls invariable will include many false positives. As a proportion of total variants called, the number of false positives is likely to increase as more and more samples are sequenced. To keep the number of false positives under control, it is important to both apply an increasingly strict set of quality control filters but also to experimentally validate some newly discovered variants. Many of these filters are currently implemented in [[GATK]].
+
Initial sets of SNP calls invariably will include many false positives. The fraction of false positives among all variants called is likely to increase as more and more samples are sequenced. To keep the fraction of false positives under control, it is important to both apply an increasingly strict set of quality control filters but also to experimentally validate some newly discovered variants. Many of these filters are currently implemented in [[GATK]].
    
; Mapping Quality Filter
 
; Mapping Quality Filter
: Consider removing variants at sides with low mapping quality scores. Even if the average mapping quality score for a site is high, consider removing variants at sites where a large fraction of reads have low mapping quality scores.
+
: Consider removing variants at sites with low mapping quality scores. Even if the average mapping quality score for a site is high, consider removing variants at sites where a noticeable fraction of reads have low mapping quality scores.
    
; Allele Balance Filter
 
; Allele Balance Filter
 
: Among individuals who are assigned an heterozygous genotype, check the proportion of reads supporting each allele. Consider filtering out variants where one of the alleles accounts for <30% of reads.
 
: Among individuals who are assigned an heterozygous genotype, check the proportion of reads supporting each allele. Consider filtering out variants where one of the alleles accounts for <30% of reads.
 +
 +
; Local Realignment Filter
 +
: Consider removing single nucleotide variants at sites where local realignment of all covering sequence reads suggests that a short indel polymorphism is present in the population.
    
; Read Depth
 
; Read Depth
: Consider filtering out variants at sites where total read depth is unusually low. Unfortunately, capture protocols introduce very large amounts of variation in sequencing depth and it is usually not possible to also accurately filter out sites sequenced at very high depth.  
+
: Consider filtering out variants at sites where total read depth is unusually low. Unfortunately, capture protocols introduce very large amounts of variation in sequencing depth and it is usually not possible to accurately filter out sites sequenced at very high depth.  
    
= Special Considerations for Admixed Samples =
 
= Special Considerations for Admixed Samples =
Line 112: Line 115:  
= Initial Association Analyses =
 
= Initial Association Analyses =
   −
We anticipate that, at least early on, the initial association analysis of whole genome datasets in the context of complex trait association studies will focus on identifying and resolving quality control issues that might result in unexpected artifacts.
+
We anticipate that, at least early on, the initial association analysis of whole exome datasets in the context of complex trait association studies will focus on identifying and resolving quality control issues that might result in unexpected artifacts.
    
== Initial Single SNP Tests ==
 
== Initial Single SNP Tests ==
29

edits

Navigation menu