Skip to content

Latest commit

 

History

History
 
 

1.5.1

Folders and files

NameName
Last commit message
Last commit date

parent directory

..
 
 
 
 

VADR container

Main tool : VADR

VADR is a suite of tools for classifying and analyzing sequences homologous to a set of reference models of viral genomes or gene families. It has been mainly tested for analysis of Norovirus, Dengue, and SARS-CoV-2 virus sequences in preparation for submission to the GenBank database.

You can find additional information on the SARS-CoV-2 models used by VADR here. The models are downloaded from the NCBI FTP server

Additional tools:

  • perl v5.22.1
  • infernal v1.1.4
  • hmmer v3.3.2
  • ncbi-blast+ v2.12.0
  • fasta v36.3.8h (the tool, not the file format)
  • minimap2 v2.24

Available VADR models:

  • sarscov2 v1.3-2
  • Mpox (AKA MPXV, formerly known as "Monkeypox") v1.4.2-1
  • Norovirus and other Caliciviridae
  • Dengue virus and other Flaviviridae
  • RSV v1.5.1-2

FYI

  • Mpox FYIs
    • Note: Support for MonkeyPox genome annotation was added to the VADR software (July 2022) and is under active development. Things may change quickly. See the above documentation ^ to see the latest information on the state of MPXV annotation with VADR.
    • Also be aware that some Mpox sequences may take up to 30 minutes to annotate, depending on how divergent it is from the RefSeq NC_063383 sequence. Some sequences may only take a minute or so.
  • Most of the VADR model files are located at /opt/vadr/vadr-models in the container filesystem and this path is stored in the globally accessible bash variable $VADRMODELDIR. For most applications, there is no need to specify v-annotate.pl --mdir /path/to/model/files since $VADRMODELDIR is set in the environment.
    • The exception is that Dengue and other Flaviviridae model files are located at /opt/vadr/vadr-models-flavi/ within the container filesystem. To use these models, please specify the 2 options: v-annotate.pl --mdir /opt/vadr/vadr-models-flavi/ --mkey flavi. A full example command can be found below.

VADR Documentation

Example Usage

# trim terminal Ns from my input genome (VADR requires this as the first step)
# for MPXV, adjust maxlen to 210000
/opt/vadr/vadr/miniscripts/fasta-trim-terminal-ambigs.pl \
    /data/SRR13957123.consensus.fa \
    --minlen 50 \
    --maxlen 30000 \
    > /data/SRR13957123.consensus.trimmed.fasta

# run v-annotate.pl using the sarscov2 model to annotate a trimmed input genome
v-annotate.pl --noseqnamemax --glsearch -s -r --nomisc \
    --mkey sarscov2 --lowsim5seq 6 \
    --lowsim3seq 6 --alt_fail lowscore,insertnn,deletinn \
    /data/SRR13957123.consensus.trimmed.fasta \
    SRR13957123-vadr-outdir

# run v-annotate.pl using mpxv model to annotate a trimmed input genome
v-annotate.pl --split --cpu 8 --glsearch -s -r --nomisc --mkey mpxv \
    --r_lowsimok --r_lowsimxd 100 --r_lowsimxl 2000 --alt_pass discontn,dupregin \
    --minimap2 --s_overhang 150 \
    mpxv.consensus.trimmed.fasta \
    mpxv-vadr-1.5-test-output

# run v-annotate.pl using Flaviviridae model to annotate a Dengue viral genome
v-annotate.pl --split --cpu 1 --group Dengue --nomisc --noprotid \
    --mdir /opt/vadr/vadr-models-flavi/ --mkey flavi \
    GCF_000862125.1_ViralProj15306_genomic.fna \
    dengue-test-outdir