From Genome Analysis Wiki
Jump to navigationJump to search
2,981 bytes removed
, 21:52, 15 June 2014
| Line 292: |
Line 292: |
| | | | |
| | Also remember to index this file and extract the sites. | | Also remember to index this file and extract the sites. |
| − |
| |
| − | ==Insertion/Deletion ratios, Coding Regions and Overlap analysis==
| |
| − |
| |
| − | You can obtain measure of insertion deletion ratios, coding region indels and sensitivity analysis by using the profile_indels analysis.
| |
| − |
| |
| − | vt profile_indels -g indel.reference.txt -r ~/ref/vt/grch37/hs37d5.fa mills.normalized.sites.bcf
| |
| − |
| |
| − | The indel.reference.txt file contains the required reference to perform the overlap analysis.
| |
| − |
| |
| − | data set
| |
| − | No Indels : 8904 [0.93] //#variants in your data set [ins/del ratio]
| |
| − | FS/NFS : 0.66 (67/35) //Proportion of frameshift Indels. (#Frameshift Indels/#Nonframeshift Indels)<br>
| |
| − | dbsnp //A represents the data set you input, B represents dbsnp
| |
| − | A-B 2975 [1.06] //#variants in A only [ins/del ratio]
| |
| − | A&B 5929 [0.86] //#variants in A and B
| |
| − | B-A 2059845 [1.51]
| |
| − | Precision 66.6% //A&B/A this represents how novel your data set is in the variants represented.
| |
| − | Sensitivity 0.3% //A&B/B this represents sensitivity somewhat if dbsnp is considered a high quality Indel
| |
| − | //set and the sample are the same in both data sets. (which they usually are not, this is still
| |
| − | //nonetheless a useful indicator)<br>
| |
| − |
| |
| − |
| |
| − | Ins/Del ratios: Reference alignment based methods tend to be biased towards the detection of deletions. This provides a useful measure for discovery Indel sets to show the varying degree of biasness. It also appears that as coverage increases, the ins/del ratio tends to 1.
| |
| − |
| |
| − | Coding region analysis: Coding region Indels may be categorised as Frame shift Indels and Non frameshift Indels. A lower proportion of Frameshift Indels may indicate a better quality data set but this depends also on the individuals sequenced.
| |
| − |
| |
| − | Overlap analysis: overlap analysis with other data sets is an indicator of sensitivity.
| |
| − |
| |
| − | * dbsnp: contains Indels submitted from everywhere, I am not sure what does this represent exactly. But assuming most are real, then precision is a useful estimated quantity from this reference data set.
| |
| − | * Mills: contains doublehit common indels from the Mills. et al paper and is a relatively good measure of sensitivity for common variants. Because not all Indels in this set is expected to be present in your sample, this actually gives you an underestimate of sensitivity.
| |
| − | * Mills chip: This is a subset of the Mills data set. There are genotypes here that are useful for subsetting polymophic subsets of variants that are present in samples common with your data set, this can potentially provide a better estimate of sensitivity. In general not very useful unless you happen to be working on 1000 Genomes data or any data set who's individuals are commonly studied.
| |
| − | * Affy Exome Chip: This contains somewhat rare variants in exonic regions and is useful for exome chip analysis. You should subset your exome data to exome region Indels before comparing against this data set.
| |
| − |
| |
| − | This analysis supports filters too.
| |
| | | | |
| | ==to document== | | ==to document== |
