align_trim

Stand alone version of ARTIC's fieldbioinfomatics align_trim.py

Installation

From conda

conda install bioconda::align_trim 

from pypi

pip install align_trim

from source

git clone https://github.com/artic-network/align_trim.git
cd align_trim
uv sync
uv run align_trim --help

Command Line Interface

Basic Usage

aligntrim [OPTIONS] BEDFILE

The tool reads alignment data from either a SAM file or stdin and outputs trimmed alignments to stdout in SAM format by default.

Required Arguments

• BEDFILE: BED file containing the amplicon primer scheme in v3 format.

Optional Arguments

Input/Output Options

• --samfile, -s : Sorted SAM/BAM file containing the aligned reads, if this is not provided (or '-') then 'align_trim' will read from stdin.
• --output, -o : Output file path. Format determined by extension (.sam/.bam). If not provided or '-', writes SAM to stdout

Processing Options

• --normalise, -n : Subsample to N coverage per amplicon. Use 0 for no normalisation (default: 0)
• --min-mapq, -m : Minimum mapping quality to keep an aligned read (default: 20)
• --primer-match-threshold, -p : Fuzzy match primer positions within this threshold (default: 35)

Primer and Read Handling

• --no-trim-primers : Do not trim primers from reads (by default, primers are trimmed)
• --allow-incorrect-pairs : Allow reads to be assigned to amplicons even if primers are not correctly paired
• --require-full-length : Require all reads to start and stop in primer sites (do not use with rapid barcoding)

Output and Reporting

• --report, -r : Output detailed report TSV to specified filepath
• --amp-depth-report, -a : Output mean depth for each amplicon as TSV to specified filepath
• --no-read-groups : Do not divide reads into pool-based read groups in SAM/BAM output

General Options

• --verbose, -v : Enable debug mode with detailed logging to stderr
• --version : Show version information
• --help : Show help message

Examples

Basic trimming with primer removal

aligntrim primers.bed --bamfile input.bam --output trimmed.bam

Normalize coverage and generate reports

aligntrim primers.bed --bamfile input.bam --normalise 100 \
  --report alignment_report.tsv --amp-depth-report depth_report.tsv \
  --output normalized.bam

Process from stdin with verbose output

samtools view -h input.bam | aligntrim primers.bed --verbose > trimmed.sam 2> verbose.out.txt

Strict full-length read filtering

aligntrim primers.bed --bamfile input.bam --require-full-length \
  --min-mapq 30 --output filtered.bam

Allow mismatched primer pairs with custom threshold

aligntrim primers.bed --bamfile input.bam --allow-incorrect-pairs \
  --primer-match-threshold 50 --output relaxed.bam

Output Formats

The tool supports multiple output formats based on file extension:

• .sam - SAM format (text)
• .bam - BAM format (binary, compressed)
• No extension or - - SAM format to stdout

Report Files

When using --report, a tab-separated file is generated with the following columns:

• chrom: Reference chromosome/contig
• QueryName: Read name
• ReferenceStart/ReferenceEnd: Alignment coordinates
• PrimerPair: Primer pair assignment
• Primer1/Primer2: Individual primer information
• CorrectlyPaired: Boolean indicating proper primer pairing
• Additional alignment metrics

The --amp-depth-report generates a summary of coverage depth per amplicon.