Changes

From Genome Analysis Wiki
Jump to: navigation, search

VcfCodingSnps

1,782 bytes added, 14:32, 14 May 2010
no edit summary
The annotating result will be added to the entry "INFO" of the input VCF SNP file and outputted together with other information. If a SNP is annotated differently with respect to different genes (or different isoforms of the same gene), all the annotated results will be added into the entry "INFO". If the SNP is NOT in any gene coding region, then the original "INFO" will be outputted. Here is an example of input and output VCF file headlines:
##format=VCFv3.2
##NA12891=../GLF/NA12891.chrom8.SLX.SRP000032.2009_07.glf
##NA12892=../GLF/NA12892.chrom8.SLX.SRP000032.2009_07.glf
##NA12878=../merged/NA12878.chrom8.merged.glf
##minTotalDepth=0
##maxTotalDepth=1000
##minMapQuality=40
##minPosterior=0.9990
##program=glfTrio
##versionDate=Thu Aug 27 18:23:18 2009
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12891 NA12892 NA12878
8 146284 . c a 54 . depth=29;duples=hets;mac=2;tdt=0/2 GT:GQ:GD 1/0:31:12 1/0:32:3 0/0:28:14
8 146703 . c t 92 . depth=41;mac=1;tdt=0/1 GT:GQ:GD 1/1:42:14 0/1:54:9 1/1:24:18
8 151532 . t c 100 . depth=131 GT:GQ:GD 0/0:8:37 1/0:100:26 1/0:100:68
8 151573 . g t 72 . depth=113;mac=1;tdt=1/1 GT:GQ:GD 0/1:48:35 0/0:39:26 0/1:100:52
8 151638 . a c 100 . depth=124;duples=hets;mac=2;tdt=1/2 GT:GQ:GD 0/1:100:55 0/1:100:58 0/1:87:11
8 151651 . c g 100 . depth=124;duples=hets;mac=2;tdt=1/2 GT:GQ:GD 0/1:87:56 0/1:100:56 0/1:24:12
8 151763 . t a 100 . depth=127;duples=hets;mac=2;tdt=1/2 GT:GQ:GD 1/0:100:49 1/0:100:54 1/0:100:24
8 151936 . a g 32 . depth=105;duples=hets;mac=2;tdt=0/2 GT:GQ:GD 0/1:42:44 0/1:23:47 0/0:39:14
8 152578 . c t 87 . depth=108 GT:GQ:GD 1/1:95:31 1/1:89:30 1/1:100:47
76
edits

Navigation menu