From Genome Analysis Wiki
Jump to navigationJump to search
1,119 bytes removed
, 22:02, 7 August 2013
| Line 103: |
Line 103: |
| | | | |
| | ===== From an Annotated VCF File or Output a VCF File to Annotate ===== | | ===== From an Annotated VCF File or Output a VCF File to Annotate ===== |
| − | * If --groupFile option is '''NOT''' specified, '''rareMETAL''' will look for an annotated vcf file as blue print for variants to group.
| + | If --groupFile option is '''NOT''' specified, '''rareMETAL''' will look for an annotated vcf file as blue print for variants to group. Users are also allowed to generate a vcf file based on the superset of variants from pooled samples, and annotate outside rareMETAL. Then, annotated vcf file can be used as input for rareMETAL for gene-level meta-analysis, or group files can be generated based on the annotated vcf file. Detailed description of these options are available. There are also examples of this usage at the bottom of this page. |
| − | * The annotated VCF file should be specified using --annotatedVcf option.
| |
| − | * --annotation should be used with --annotatedVcf together when specific category of functional variants are of interest to be grouped. For example, if grouping nonsynonymous and splicing variants are of interests, the following should be included in command line:
| |
| − | | |
| − | --annotatedVcf your.annotated.vcf --annotation nonsyn/splicing
| |
| − | Note: this allows you to group variants that are annotated starting with nonsyn or splicing (not case-sensitive).
| |
| − | | |
| − | * Special format for the annotated VCF file is required: all annotation information should be coded in INFO field in VCF file, starting with the key "ANNO=". An example annotated VCF file is in the following:
| |
| − | | |
| − | #CHROM POS ID REF ALT QUAL FILTER INFO
| |
| − | 1 19208194 . G A 100 PASS
| |
| − | AC=3;'''ANNO='''nonsynonymous:ALDH4A1:NM_170726:exon8:c.C866T:p.P289L,ALDH4A1:NM_001161504:exon8:c.C686T:p.P229L,ALDH4A1:NM_003748:exon8:c.C866T:p.P289L,;
| |
| − | '''ANNO='''splicing:ALDH4A1
| |
| − | 1 19208293 . G C 100 PASS AC=7;STUDIES=5;MAC=7;MAF=0.001;DESIGN=TBD_ASSAY;DSCORE=1.00;
| |
| − | '''ANNO='''nonsynonymous:ALDH4A1:NM_170726:exon8:c.C767G:p.P256R,ALDH4A1:NM_001161504:exon8:c.C587G:p.P196R,ALDH4A1:NM_003748:exon8:c.C767G:p.P256R,
| |
| − | | |
| − | * Notice that each variant is allowed to have more than one annotations; but each annotation should start with a new key "ANNO=" followed by annotation:genename:other transcript information.
| |
| | | | |
| | ===== Generate a VCF File to Annotate Outside of Rare Metal ===== | | ===== Generate a VCF File to Annotate Outside of Rare Metal ===== |
