GATK之VariantAnnotator

简介: VariantAnnotator简要说明用途: 利用上下文信息注释识别的变异位点(variant calls)分类: 变异位点操作工具概要: 根据变异位点的背景信息(与功能注释相对)进行注释。

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
目录
相关文章
|
9月前
|
算法 关系型数据库 数据挖掘
Sentieon | 每周文献-Tumor Sequencing(肿瘤测序)-第三期
本期重点关注Sentieon在甲状腺癌和CRISPR/Cas9技术中的应用
86 0
Sentieon | 每周文献-Tumor Sequencing(肿瘤测序)-第三期
|
3月前
|
机器学习/深度学习 测试技术 决策智能
【论文速递】CVPR2022 - 全局跟踪Transformers
【论文速递】CVPR2022 - 全局跟踪Transformers
|
8天前
|
数据挖掘
R语言预测期货波动率的实现:ARCH与HAR-RV与GARCH,ARFIMA模型比较
R语言预测期货波动率的实现:ARCH与HAR-RV与GARCH,ARFIMA模型比较
26 0
R语言预测期货波动率的实现:ARCH与HAR-RV与GARCH,ARFIMA模型比较
|
1月前
单细胞分析|在 Seurat 中使用 sctransform
单细胞分析|在 Seurat 中使用 sctransform
19 0
|
2月前
|
算法 搜索推荐 数据挖掘
Sentieon | 每周文献-Liquid Biopsy-第三十二期
Sentieon | 每周文献-Liquid Biopsy-第三十二期
15 1
|
3月前
|
算法 安全 数据挖掘
Sentieon | 每周文献-Tumor Sequencing-第十四期
Sentieon | 每周文献-Tumor Sequencing-第十四期
41 0
|
7月前
|
存储 JSON Java
GATK4重测序数据怎么分析?
GATK4重测序数据怎么分析?
|
Serverless 测试技术 TensorFlow
ResNetV2模型复现--CVPR2016
关于残差结构的不同尝试
285 0
ResNetV2模型复现--CVPR2016
|
传感器 编解码 算法
Google Earth Engine ——GCOM-C 进行长期和持续的陆地地表温度数据集(JAXA/GCOM-C/L3/LAND/LAI/V2)
Google Earth Engine ——GCOM-C 进行长期和持续的陆地地表温度数据集(JAXA/GCOM-C/L3/LAND/LAI/V2)
176 0
Google Earth Engine ——GCOM-C 进行长期和持续的陆地地表温度数据集(JAXA/GCOM-C/L3/LAND/LAI/V2)
|
机器学习/深度学习 自然语言处理 测试技术
CVPR2021| TimeSformer-视频理解的时空注意模型
transformer在视频理解方向的应用主要有如下几种实现方式:Joint Space-Time Attention,Sparse Local Global Attention 和Axial Attention。这几种方式的共同点是采用ViT中的方式将图像进行分块,而它们之间的区别在于如何用self attention来处理这些块。
CVPR2021| TimeSformer-视频理解的时空注意模型