7.4 VCFDiagnostics

Diagnosing a VCF file

VCFDiagnostics estimates the fraction of reference to the alternate allele. This task requires the input VCF file to have allelic depth field. None of the VCF files generated by any of the ATLAS tasks have allelic depth field.

7.4.1 Input

Required inputs :

--vcf example1.vcf.gz VCF file with allelic depth field.

Optional inputs :

  • None

Specific Parameters :

  • None

  • See Filter parameters to apply specific filters for bases, reads and parsing window setting.

Engine parameters that are common to all tasks can be found here.

7.4.2 Output

*.txt Text file with ratio of reference to the alternate allele.

7.4.3 Usage Example

#! /bin/bash

# `--fixedSeed = N` is needed to have reproducable results in regression test

. $(dirname $0)/find_atlas
. $(dirname $0)/simulate_vcf --fixedSeed 3 --sampleSize 2 --chrLength 2321 --ploidy 2

out="VCFDiagnostics"
$atlas --task VCFDiagnostics --vcf simulate.vcf.gz \
       --fixedSeed 2 --out $out --logFile $out.out 2> $out.eout

out="VCFFixInt"
$atlas --task VCFDiagnostics --fixInt --vcf simulate.vcf.gz \
       --fixedSeed 1 --out $out --logFile $out.out 2> $out.eout

out="VCFToInvariantBed"
$atlas --task VCFDiagnostics --writeBED --vcf simulate.vcf.gz \
       --fixedSeed 0 --out $out --logFile $out.out 2> $out.eout