Troubleshoot common issues
If you have a question that is not addressed here, please don’t hesitate to ask for help
My country / division does not show up on the map
This is most often a result of the country / division not being present in the file defining the latitude & longitude of each deme. Adding it to that file (and rerunning the Snakemake rules downstream of this) should fix this.
My trait (e.g. division) is grey instead of colored
We generate the colors from the colors
rule in the Snakefile, which uses the ordering TSV to generate these. See Workflow config file guide for more info.
A note about locations and colors: Unless you want to specifically override the colors generated, it’s usually easier to add information to the default ncov
files, so that you can benefit from all the information already in those files.
My genomes aren’t included in the analysis
There are a few steps where sequences can be removed:
During the
filter
step:Samples that are included in the exclude file are removed
Samples that fail the filtering criteria, as defined in your filter config, are removed.
If you do not have any custom filtering criteria, the default filters in the parameters.yaml are applied.
Check the
results/{build_name}/filtered_log.tsv
file to see the filtered reason for each sequence.
Samples may be randomly removed during subsampling; see Workflow config file guide for more info.
During the
refine
step, Augur can drop samples that deviate from the expected clock rate. Inspect the log file named likelogs/refine_{build_name}.txt
to look for samples filtered by this step. See the refine configuration guide, for details on the clock rate filter.
Sequencing and alignment errors
Genome sequencing, bioinformatic processing of the raw data, and alignment of the sequences are all steps were errors can slip in. Such errors can distort the phylogenetic analysis. To avoid sequences with known problems to mess up the analysis, we keep a list of problematic sequences in config/exclude.txt
and filter them out. To facilitate spotting such problematic sequences, we added an additional quality control step that produces the file results/excluded_by_diagnostics.txt
.
This file is the output of scripts/diagnostic.py
and is produced by rule diagnostic
. This file contains only those sequences with diagnostics exceeding thresholds and mirrors the format of config/exclude.txt
. These names could be added to config/exclude.txt
for permanent exclusion. Note, however, that some sequences might look problematic due to alignment issues rather than intrinsic problems with the sequence. The flagged sequences will be excluded from the current run.
To only run the sequence diagnostic, you can specify any of the three above files as target, or use the diagnostic
target:
nextstrain build ... diagnostic