alcov 1.1.11

Identify frequencies of concerning mutations from aligned reads

Alcov

Abundance learning for SARS-CoV-2 variants. The primary purpose of the tool is:

• Estimating abundace of variants of concern from wastewater sequencing data

You can read more about how Alcov works in the preprint, Alcov: Estimating Variant of Concern Abundance from SARS-CoV-2 Wastewater Sequencing Data

The tool can also be used for:

• Converting nucleotide and amino acid mutations for SARS-CoV-2 such as those found on https://covariants.org/variants/S.N501
• Determining the frequency of mutations of interest in BAM files
• Plotting the depth for each ARTIC amplicon (https://github.com/artic-network/artic-ncov2019/tree/master/primer\_schemes/nCoV-2019/V3)
• Comparing amplicon GC content with its read depth (as a measure of degredation)

The tool is under active development. If you have questions or issues, please open an issue on GitHub or email me (email in setup.py).

Installing

The latest release can be downloaded from PyPI

pip install alcov

This will install the Python library and the CLI.

To install the development version, clone the repository and run

pip install .

Usage example

Preprocessing

Alcov expects a BAM file of reads aligned to the SARS-CoV-2 reference genome. For an example of how to process Illumina reads, check the prep directory for a script named "prep.py" which outlines our current preprocessing pipeline, including the generation of a "samples.txt" file used by alcov "find_lineages" command.

Estimating relative abundance of variants of concern:

alcov find_lineages reads.bam

Finding lineages in BAM files for multiple samples:

alcov find_lineages samples.txt

Where samples.txt looks like:

path/to/reads1.bam	Sample 1 name
path/to/reads2.bam	Sample 2 name
...

Example usage: To estimate the relative abundance of lineages in a list of samples (samples.txt), while considering only positions with a minimum depth of 10 reads, the following command can be used. This will also save the heatmap as a .png image and the corresponding frequencies as a csv file.

alcov find_lineages --min_depth=10 --save_img=True --csv=True samples.txt

Optionally specify which VOCs to look for (Note: This will restrict alcov to only consider the lineages specified in this text file. Do not provide this file if you wish alcov to consider all lineages for which it has constellation files.)

alcov find_lineages reads.bam lineages.txt

Where lineages.txt looks like: Note: These lineages must be chosen from the list of lineages that alcov has constellation files for (updated weekly) found in "alcov/alcov/data/constellations/"

BA.5-like
BQ.1.1-like
XBB-like
XBB.1.5-like
...

Optionally change minimum read depth (default 40)

alcov find_lineages --min_depth=5 reads.bam

Optionally show how predicted mutation rates agree with observed mutation rates

alcov find_lineages --show_stacked=True reads.bam

Use mutations which are found in multiple VOCs (can help for low coverage samples) - This is now the default behaviour.

alcov find_lineages --unique=False reads.bam

Plotting change in lineage distributions over time for multiple sites

alcov find_lineages --ts samples.txt

Where samples.txt looks like:

reads1.bam	SITE1_2021-09-10
reads2.bam	SITE1_2021-09-12
...
reads3.bam	SITE2_2021-09-10
reads4.bam	SITE2_2021-09-12
...

Converting mutation names:

$ alcov nt A23063T
A23063T causes S:N501Y
$ alcov aa S:E484K
G23012A causes S:E484K

Finding mutations in BAM file:

alcov find_mutants reads.bam

Finding mutations in BAM files for multiple samples:

alcov find_mutants samples.txt

Where samples.txt looks like:

reads1.bam	Sample 1 name
reads2.bam	Sample 2 name
...

Running find_mutants will print the number of reads with and without each mutation in each sample and then generate a heatmap showing the frequencies for all samples.

You can also specify a custom mutations file:

alcov find_mutants samples.txt mutations.txt

Where mutations.txt looks like:

S:N501Y
G23012A
...

Getting the read depth for each amplicon

alcov amplicon_coverage reads.bam

or

alcov amplicon_coverage samples.txt

Plotting amplicon GC content against amplicon depth

alcov gc_depth reads.bam

or

alcov gc_depth samples.txt