Package: amplimap
amplimap: amplicon mapping and analysis pipelineamplimap is a command-line tool to automate the processing and analysis of data from targeted next-generation sequencing (NGS) experiments with PCR-based amplicons or capture-based enrichment systems.From raw sequencing reads, amplimap generates a variety of output files including read alignments, per-basepair nucleotide counts, target coverage data and annotated variant calls.In addition to its focus on user-friendliness and reproducibility, amplimap supports advanced features such as the generation of consensus base calls for read families based on molecular identifiers/barcodes (UMIs) and the detection of chimeric reads caused by amplification of off-target loci.InstallationWe recommend that you install amplimap through Conda:wget https://raw.githubusercontent.com/koelling/amplimap/master/environment.yml
conda env create --file environment.ymlWe also have aDocker imageavailable.
Please see ourfull installation instructionsfor additional details.OverviewTo run amplimap you create a directory containing a small set of input files:A subdirectory with FASTQ.GZ or BAM files representing your different samples (tested with Illumina MiSeq, HiSeq and NextSeq)Optionally: Files describing the targeted genomic regions, the primers you used or other custom configuration parametersThen you can runamplimapto generate a variety of different output files, depending on your experiment.
These include, for example:A target coverage table, showing you how well-covered each target region was in each sample.A table of germline variants in your samples, annotated with gene, impact, population frequencies, deleteriousness scores, etc.A per-basepair “pileup” table telling you how often each nucleotide was seen in each sample at each position.Built on top ofSnakemakeand Python 3, amplimap is entirely
automated and can be run on a single machine as well as on an HPC cluster
(e.g. LSF, SGE).Supported experimental protocolsamplimap is compatible with most targeted sequencing protocols that generate paired-end short read data.For protocols utilising PCR or smMIPs each read should start with a known primer (or targeting arm) sequence, followed by the amplified target DNA.
Reads can optionally contain a unique molecular identifier (UMI) sequence in front of the primer, which can be used to group reads into families.
Data should be available as demultiplexed FASTQ.GZ files, with each pair of files representing a different sample.For capture-based protocols data can be provided in FASTQ.GZ or unmapped/mapped BAM format, which may contain UMIs as BAM tags.Some of the protocols we have analyzed with amplimap include:PCR-based targeted resequencing (single/multiplex)smMIPs with and without UMIsProbe based target enrichment, for example:IDT xGen Lockdown probesTwist Bioscience Custom PanelsTutorialsCalling germline variants in amplicon-based resequencing dataIdentifying low-frequency somatic mutations in FGFR2 with UMI-tagged smMIPsQuantifying allele-specific expressionLinksPackage:https://pypi.org/project/amplimap/Code:https://github.com/koelling/amplimap/Documentation:https://amplimap.readthedocs.io/Citation and LicenseLicensed under the Apache License, version 2.0.
Copyright 2020 Nils Koelling.
When you use amplimap,
please cite theamplimap paperin your work:Nils Koelling, Marie Bernkopf, Eduardo Calpena, Geoffrey J Maher, Kerry A Miller, Hannah K Ralph, Anne Goriely, Andrew O M Wilkie, amplimap: a versatile tool to process and analyze targeted NGS data, Bioinformatics, Volume 35, Issue 24, 15 December 2019, Pages 5349–5350,https://doi.org/10.1093/bioinformatics/btz582