Line 54: |
Line 54: |
| Some possible annotating results for a single SNP with the meanings of their output format are listed below: | | Some possible annotating results for a single SNP with the meanings of their output format are listed below: |
| | | |
− | 5'UTR=A26C2[-] means the SNP is in the 5'UTR region of gene A26C2 with a minus strand. | + | 5'UTR=A26C2[-] means the SNP is in the 5'UTR region of gene A26C2 with a minus strand. |
− | INTRONIC=POTEG[-] means the SNP is in the intronic region of gene POTEG with a minus strand. | + | INTRONIC=POTEG[-] means the SNP is in the intronic region of gene POTEG with a minus strand. |
| SYNONYMOUS_CODING=BARD1(uc002veu.2):His506His[-] means the SNP is synonymous coding at the 506th codon in gene BARD1 with a minus strand and it keeps amino-acid His unchanged. | | SYNONYMOUS_CODING=BARD1(uc002veu.2):His506His[-] means the SNP is synonymous coding at the 506th codon in gene BARD1 with a minus strand and it keeps amino-acid His unchanged. |
| NON_SYNONYMOUS_CODING=BARD1(uc002veu.2):Arg658Cys[-] means the SNP is non_synonymous coding at the 658th codon in gene BARD1 (ucsc gene name uc002veu.2)with a minus strand and it changes amino-acid from Arg to Cys. | | NON_SYNONYMOUS_CODING=BARD1(uc002veu.2):Arg658Cys[-] means the SNP is non_synonymous coding at the 658th codon in gene BARD1 (ucsc gene name uc002veu.2)with a minus strand and it changes amino-acid from Arg to Cys. |
| SPLICE_SITE=FARP2(uc002wbi.1)[+] means the SNP is in the SPLICE_SITE (5 bp within exon start or end positions in the coding region) of gene FARP2 (ucsc gene name uc002wbi.1) with a plus strand. | | SPLICE_SITE=FARP2(uc002wbi.1)[+] means the SNP is in the SPLICE_SITE (5 bp within exon start or end positions in the coding region) of gene FARP2 (ucsc gene name uc002wbi.1) with a plus strand. |
| STOP_GAINED=C2orf83(uc002vph.1):Trp141stop[-] means the SNP is the 141th codon in gene MAPK12 (ucsc gene name uc002vph.1) with a minus strand and it changes amino-acid Trp to a stop codon. | | STOP_GAINED=C2orf83(uc002vph.1):Trp141stop[-] means the SNP is the 141th codon in gene MAPK12 (ucsc gene name uc002vph.1) with a minus strand and it changes amino-acid Trp to a stop codon. |
| + | STOP_LOST=OR2M3(uc001ieb.1):stop313Arg[+] means the SNP is the 313th codon in gene OR2M3 (ucsc gene name uc001ieb.1) with a plus strand and it changes a stop codon to amino-acid Arg. |
| | | |
− | The annotating result will be added to the entry "INFO" of the input VCF SNP file and outputted together with other information. Here is an example of output VCF file headlines: | + | 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 output VCF file headlines: |
| | | |
− | ##format=VCFv3.2
| + | ##format=VCFv3.2 |
− | ##NA12891=../GLF/NA12891.chrom22.SLX.SRP000032.2009_07.glf
| + | ##NA12891=../GLF/NA12891.chrom8.SLX.SRP000032.2009_07.glf |
− | ##NA12892=../GLF/NA12892.chrom22.SLX.SRP000032.2009_07.glf
| + | ##NA12892=../GLF/NA12892.chrom8.SLX.SRP000032.2009_07.glf |
− | ##NA12878=../merged/NA12878.chrom22.merged.glf
| + | ##NA12878=../merged/NA12878.chrom8.merged.glf |
− | ##minTotalDepth=0
| + | ##minTotalDepth=0 |
− | ##maxTotalDepth=1000
| + | ##maxTotalDepth=1000 |
− | ##minMapQuality=40
| + | ##minMapQuality=40 |
− | ##minPosterior=0.9990
| + | ##minPosterior=0.9990 |
− | ##program=glfTrio
| + | ##program=glfTrio |
− | ##versionDate=Thu Aug 27 18:23:18 2009
| + | ##versionDate=Thu Aug 27 18:23:18 2009 |
− | #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12891 NA12892 NA12878
| + | #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12891 NA12892 NA12878 |
− | 22 15464609 . a g 100 . depth=109;mac=1;tdt=0/1;3'UTR=psiTPTE22[+] GT:GQ:GD 0/1:100:44 1/1:81:28 1/1:100:37
| + | 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 |
− | 22 15464609 . a g 100 . depth=109;mac=1;tdt=0/1;3'UTR=psiTPTE22[+] GT:GQ:GD 0/1:100:44 1/1:81:28 1/1:100:37
| + | 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 |
− | 22 15464609 . a g 100 . depth=109;mac=1;tdt=0/1;3'UTR=psiTPTE22[+] GT:GQ:GD 0/1:100:44 1/1:81:28 1/1:100:37
| + | 8 151532 . t c 100 . depth=131;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] GT:GQ:GD 0/0:8:37 1/0:100:26 1/0:100:68 |
− | 22 15464609 . a g 100 . depth=109;mac=1;tdt=0/1;3'UTR=psiTPTE22[+] GT:GQ:GD 0/1:100:44 1/1:81:28 1/1:100:37
| + | 8 151573 . g t 72 . depth=113;mac=1;tdt=1/1;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] GT:GQ:GD 0/1:48:35 0/0:39:26 0/1:100:52 |
− | 22 15482433 . a g 38 . depth=21;3'UTR=psiTPTE22[+] GT:GQ:GD 1/1:34:11 1/1:14:3 1/1:35:7
| + | 8 151638 . a c 100 . depth=124;duples=hets;mac=2;tdt=1/2;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] GT:GQ:GD 0/1:100:55 0/1:100:58 0/1:87:11 |
− | 22 15644565 . g t 77 . depth=140;NON_SYNONYMOUS_CODING=XKR3:His15644565Asn[-] GT:GQ:GD 1/1:100:49
| + | 8 151651 . c g 100 . depth=124;duples=hets;mac=2;tdt=1/2;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] 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;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] 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;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] GT:GQ:GD 0/1:42:44 0/1:23:47 0/0:39:14 |
| + | 8 152578 . c t 87 . depth=108;5'UTR=RPL23A_20_869(uc010lra.1)[-];5'UTR=RPL23A_20_869(uc003woq.2)[-];5'UTR=RPL23A_20_869(uc010lrb.1)[-] GT:GQ:GD 1/1:95:31 1/1:89:30 1/1:100:47 |