「GATK 3.8」VariantAnnotator介绍

VariantAnnotator

简要说明

用途: 利用上下文信息注释识别的变异位点(variant calls)
分类: 变异位点操作工具
概要: 根据变异位点的背景信息(与功能注释相对)进行注释。目前有许多的注释模块(见注释模块一节)可供使用。

输入文件

用于注释的VCF文件和可选的BAM文件

输出文件

注释完毕的VCF文件

使用案例

HaplotypeCallerUnifiedGenotyper的结果中增加每个样本的深度和dbSNP ID信息。

java -jar GenomeAnalysisTK.jar \
   -R reference.fasta \
   -T VariantAnnotator \
   -I input.bam \
   -V input.vcf \
   -o output.vcf \
   -A Coverage \
   --dbsnp dbsnp.vcf

参数说明:

-R/--reference_sequence:参考基因组
-T/--analysis_type : 运行的工具
-I/--input_file: 和vcf相应的BAM文件
-o :输出文件
-V/--varaint: 输入的VCF文件
-A/--annotation: 要添加哪些注释项
--dbsnp: 已有的snp信息注释数据库

HaplotypeCaller和MuTect2也有-A选项,并且有些注释模块只能在HaplotypeCaller和MuTect2计算,例如StrandAlleleCountsBySample
如下是 -A可接的内容:

Standard annotations in the list below are marked with a '*'.
Available annotations for the VCF INFO field:
        AS_BaseQualityRankSumTest
        AS_FisherStrand
        AS_InbreedingCoeff
        AS_InsertSizeRankSum
        AS_MQMateRankSumTest
        AS_MappingQualityRankSumTest
        AS_QualByDepth
        AS_RMSMappingQuality
        AS_ReadPosRankSumTest
        AS_StrandOddsRatio
        AlleleBalance
        BaseCounts
        *BaseQualityRankSumTest
        *ChromosomeCounts
        ClippingRankSumTest
        ClusteredReadPosition
        *Coverage
        *ExcessHet
        *FisherStrand
        FractionInformativeReads
        GCContent
        GenotypeSummaries
        *HaplotypeScore
        HardyWeinberg
        HomopolymerRun
        *InbreedingCoeff
        LikelihoodRankSumTest
        LowMQ
        MVLikelihoodRatio
        *MappingQualityRankSumTest
        MappingQualityZero
        NBaseCount
        PossibleDeNovo
        *QualByDepth
        *RMSMappingQuality
        *ReadPosRankSumTest
        SampleList
        SnpEff
        SpanningDeletions
        *StrandOddsRatio
        TandemRepeatAnnotator
        TransmissionDisequilibriumTest
        VariantType
        
Available annotations for the VCF FORMAT field:
        AlleleBalanceBySample
        AlleleCountBySample
        BaseCountsBySample
        BaseQualitySumPerAlleleBySample
        *DepthPerAlleleBySample
        DepthPerSampleHC
        MappingQualityZeroBySample
        OxoGReadCounts
        StrandAlleleCountsBySample
        StrandBiasBySample


Available classes/groups of annotations:
        AS_RMSAnnotation
        AS_RankSumTest
        AS_StandardAnnotation
        AS_StrandBiasTest
        ActiveRegionBasedAnnotation
        BetaTestingAnnotation
        ExperimentalAnnotation
        RMSAnnotation
        RankSumTest
        ReducibleAnnotation
        RodRequiringAnnotation
        StandardAnnotation
        StandardHCAnnotation
        StandardSomaticAnnotation
        StandardUGAnnotation
        StrandBiasTest
        WorkInProgressAnnotation

注释模块

这是官方文档提供的注释模块:

Name Summary
AS_BaseQualityRankSumTest Allele-specific rank Sum Test of REF versus ALT base quality scores
AS_FisherStrand Allele-specific strand bias estimated using Fisher's Exact Test *
AS_InbreedingCoeff Allele-specific likelihood-based test for the inbreeding among samples
AS_InsertSizeRankSum Allele specific Rank Sum Test for insert sizes of REF versus ALT reads
AS_MQMateRankSumTest Allele specific Rank Sum Test for mate's mapping qualities of REF versus ALT reads
AS_MappingQualityRankSumTest Allele specific Rank Sum Test for mapping qualities of REF versus ALT reads
AS_QualByDepth Allele-specific call confidence normalized by depth of sample reads supporting the allele
AS_RMSMappingQuality Allele-specific Root Mean Square of the mapping quality of reads across all samples.
AS_ReadPosRankSumTest Allele-specific Rank Sum Test for relative positioning of REF versus ALT allele within reads
AS_StrandOddsRatio Allele-specific strand bias estimated by the Symmetric Odds Ratio test
AlleleBalance Allele balance across all samples
AlleleBalanceBySample Allele balance per sample
AlleleCountBySample Allele count and frequency expectation per sample
BaseCounts Count of A, C, G, T bases across all samples
BaseCountsBySample Count of A, C, G, T bases for each sample
BaseQualityRankSumTest Rank Sum Test of REF versus ALT base quality scores
BaseQualitySumPerAlleleBySample Sum of evidence in reads supporting each allele for each sample
ChromosomeCounts Counts and frequency of alleles in called genotypes
ClippingRankSumTest Rank Sum Test for hard-clipped bases on REF versus ALT reads
ClusteredReadPosition Detect clustering of variants near the ends of reads
Coverage Total depth of coverage per sample and over all samples.
DepthPerAlleleBySample Depth of coverage of each allele per sample
DepthPerSampleHC Depth of informative coverage for each sample.
ExcessHet Phred-scaled p-value for exact test of excess heterozygosity
FisherStrand Strand bias estimated using Fisher's Exact Test
FractionInformativeReads The fraction of reads deemed informative over the entire cohort
GCContent GC content of the reference around the given site
GenotypeSummaries Summarize genotype statistics from all samples at the site level
HaplotypeScore Consistency of the site with strictly two segregating haplotypes
HardyWeinberg Hardy-Weinberg test for transmission disequilibrium
HomopolymerRun Largest contiguous homopolymer run of the variant allele
InbreedingCoeff Likelihood-based test for the inbreeding among samples
LikelihoodRankSumTest Rank Sum Test of per-read likelihoods of REF versus ALT reads
LowMQ Proportion of low quality reads
MVLikelihoodRatio Likelihood of being a Mendelian Violation
MappingQualityRankSumTest Rank Sum Test for mapping qualities of REF versus ALT reads
MappingQualityZero Count of all reads with MAPQ = 0 across all samples
MappingQualityZeroBySample Count of reads with mapping quality zero for each sample
NBaseCount Percentage of N bases
OxoGReadCounts Count of read pairs in the F1R2 and F2R1 configurations supporting the reference and alternate alleles
PossibleDeNovo Existence of a de novo mutation in at least one of the given families
QualByDepth Variant call confidence normalized by depth of sample reads supporting a variant
RMSMappingQuality Root Mean Square of the mapping quality of reads across all samples.
ReadPosRankSumTest Rank Sum Test for relative positioning of REF versus ALT alleles within reads
SampleList List samples that are non-reference at a given site
SnpEff Top effect from SnpEff functional predictions
SpanningDeletions Fraction of reads containing spanning deletions
StrandAlleleCountsBySample Number of forward and reverse reads that support each allele
StrandBiasBySample Number of forward and reverse reads that support REF and ALT alleles
StrandOddsRatio Strand bias estimated by the Symmetric Odds Ratio test
TandemRepeatAnnotator Tandem repeat unit composition and counts per allele
TransmissionDisequilibriumTest Wittkowski transmission disequilibrium test
VariantType General category of variant
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 221,576评论 6 515
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,515评论 3 399
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 168,017评论 0 360
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,626评论 1 296
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,625评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,255评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,825评论 3 421
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,729评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,271评论 1 320
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,363评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,498评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 36,183评论 5 350
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,867评论 3 333
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,338评论 0 24
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,458评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,906评论 3 376
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,507评论 2 359

推荐阅读更多精彩内容