From Genome Analysis Wiki
Jump to navigationJump to search
308 bytes added
, 00:01, 15 April 2014
Line 128: |
Line 128: |
| | | |
| == Additional Analysis Options == | | == Additional Analysis Options == |
| + | |
| + | === Generate a VCF File to Annotate Outside RAREMETAL === |
| + | * --writeVCF allows user to write a VCF file including pooled single variants from all studies. Then users can use their favorite annotation tool to annotate the VCF file. After annotating the VCF file, users can use that file as input for '''RAREMETAL''' for further gene-based or region-based meta analysis. |
| + | * The output vcf file will be name as: yourPrefix.pooled.variants.vcf. An example output vcf file is in the following: |
| + | #CHROM POS ID REF ALT QUAL FILTER INFO |
| + | 1 115658497 115658497 G A . . ALT_AF=0.380906; |
| + | 2 74688884 74688884 G A . . ALT_AF=8.33611e-05; |
| + | 3 121414217 121414217 C A . . ALT_AF=0.0747833; |
| + | ===Annotation=== |
| + | * RAREMETAL automatically recognizes the annotation format generated by [[ANNO | '''ANNO''']] or [[EPACTS#Annotating_VCF_file_using_EPACTS | '''EPACTS''']]. |
| + | * To annotate a the VCF generated in previous step, you can use the following command: |
| + | ./anno --in your.in.vcf.gz --out your.out.vcf.gz |
| | | |
| === Group Rare Variants from Annotated VCF === | | === Group Rare Variants from Annotated VCF === |
Line 148: |
Line 160: |
| * 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. | | * 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 ===
| |
− | * --writeVCF allows user to write a VCF file including pooled single variants from all studies. Then users can use their favorite annotation tool to annotate the VCF file. After annotating the VCF file, users can use that file as input for '''RAREMETAL''' for further gene-based or region-based meta analysis.
| |
− | * The output vcf file will be name as: yourPrefix.pooled.variants.vcf. An example output vcf file is in the following:
| |
− | #CHROM POS ID REF ALT QUAL FILTER INFO
| |
− | 1 115658497 115658497 G A . . ALT_AF=0.380906;
| |
− | 2 74688884 74688884 G A . . ALT_AF=8.33611e-05;
| |
− | 3 121414217 121414217 C A . . ALT_AF=0.0747833;
| |
| | | |
| ===Options for Report Generation=== | | ===Options for Report Generation=== |