Keyboard shortcuts

Press or to navigate between chapters

Press S or / to search in the book

Press ? to show this help

Press Esc to hide this help

VCF Fields

Rastair's output follows the VCFv4.5 specification.

Filters

NameDescription
PASSAll filters pass
lowDpLow read depth
dnCpG_lowDpLow read depth for de-novo CpG candidate
dnCpG_bqLow base quality for de-novo CpG candidate
dnCpG_mapqLow mapping quality for de-novo CpG candidate
dnCpG_vafLow variant allele frequency for de-novo CpG candidate
dnCpG_adjIncluded as adjacent position for de-novo CpG candidate, but other position did not pass filters
m_vafLow variant allele frequency for methylation candidate
m_bq_ratioLow quality ratio for methylation candidate
m_posAlt allele evidence from read edges for methylation candidate
m_highDpExcessive coverage for methylation candidate
pre_mlLow amount of usable evidence, skipping ML
low_ml_scoreMachine Learning module prediction below threshold

Info Fields

NameDescriptionVCF TypeRust TypeOccurance
ADTotal read depth for each alleleIntegerusizeR
BQRMS base qualityFloatRootMeanSquare1
DPCombined depth across samplesIntegerusize1
MQRMS mapping qualityFloatRootMeanSquare1
MQ0Number of MAPQ == 0 readsIntegerusize1
NSNumber of samples with dataIntegerusize1
AS_SB_OTOT counts per alleleIntegeru32R
AS_SB_OBOB counts per alleleIntegeru32R
SC55-base sequence context centered on the variant positionStringSmolStr1
AFAllele frequency for each ALT allele in the same order as listed (estimated from primary data, not called genotypes)Floatf64A
ABQRMS Base quality per alleleFloatf64R
AMQRMS Map quality per alleleFloatf64R
AS_SS_BQ_OTStrand-specific RMS of base quality per allele on the original top strandFloatf32R
AS_SS_BQ_OBStrand-specific RMS of base quality per allele on the original bottom strandFloatf32R
AS_SS_MQ_OTStrand-specific RMS of mapping quality per allele on the original top strandFloatf32R
AS_SS_MQ_OBStrand-specific RMS of mapping quality per allele on the original bottom strandFloatf32R
PIRRMS of relative position in readFloatf64R
ENT100Shannon entropy of 100bp sequence context around variant position. Value range (0..2)Floatf641
NABRMS of number of aligned basesFloatf64R
NOIRMS of number of indelsFloatf64R
M5mC_StrandsNumber of reads that are evidence for unmodified, modified, no SNP, snpIntegeru324
CPGIs this a CpG site?Flagbool0
CPGnovoDe-novo CPG candidate: Could the alt alleles create a new CpG site?Flagbool0

Format Fields

NameDescriptionVCF TypeRust TypeOccurance
GTGenotypeStringGenotypeAllele1
GLGenotype likelihoods, Phred-scaledIntegerPhred>G
GCGenotype confidence, Phred-scaledIntegerPhred>G
DPRead depthIntegerusize1
M5mCMethylation level at CpG sitesFloatOption<f64>M
MLPrediction of methylation/variant likelihood by rastair's machine learning modelFloatf64A