Skip to content

rhpvorderman/sequali

Repository files navigation

Sequali

Sequence quality metrics for FASTQ and uBAM files.

Features:

  • MultiQC support since MultiQC version 1.22.
  • Low memory footprint, small install size and fast execution times.
  • Informative graphs that allow for judging the quality of a sequence at a quick glance.
  • Overrepresentation analysis using 21 bp sequence fragments. Overrepresented sequences are checked against the NCBI univec database.
  • Estimate duplication rate using a fingerprint subsampling technique which is also used in filesystem duplication estimation.
  • Checks for 6 illumina adapter sequences and 17 nanopore adapter sequences for single read data.
  • Determines adapters by overlap analysis for paired read data.
  • Insert size metrics for paired read data.
  • Per tile quality plots for illumina reads.
  • Channel and other plots for nanopore reads.
  • FASTQ and unaligned BAM are supported. See "Supported formats".

Example reports:

  • GM24385_1.fastq.gz; HG002 (Genome In A Bottle) on ultra-long Nanopore Sequencing. Sequence file download.
  • GM24385_1_cut.fastq.gz; GM24385_1.fastq.gz processed with cutadapt: cutadapt -o GM24385_1_cut.fastq.gz --cut -64 --cut 64 --minimum-length 500 -Z --max-aer 0.1 GM24385_1.fastq.gz. The resulting file has 64 bp cut off from both its ends and after that filtered for a minimum length of 500 and a maximum average error rate of 0.1.
  • 21C125_R1.fastq.gz; Illumina NovaSeq X paired-end sequencing of Campylobacter jejuni. ENA accession: ERR11204024.

For more information check the documentation.

Supported formats

  • FASTQ. Only the Sanger variation with a phred offset of 33 and the error rate calculation of 10 ^ (-phred/10) is supported. All sequencers use this format today.
    • Paired end sequencing data is supported.
    • For sequences called by illumina base callers an additional plot with the per tile quality will be provided.
    • For sequences called by guppy additional plots for nanopore specific data will be provided.
  • unaligned BAM. Any alignment flags are currently ignored.
    • For uBAM data as delivered by dorado additional nanopore plots will be provided.

Installation

Installation via pip is available with:

pip install sequali

Sequali is also distributed via bioconda. It can be installed with:

conda install -c conda-forge -c bioconda sequali

Quickstart

sequali path/to/my.fastq.gz

This will create a report my.fastq.gz.html and a json my.fastq.gz.json in the current working directory.

For all command line options checkout the usage documentation.

For more extensive information about the module options check the documentation on the module options.

Acknowledgements

  • FastQC for its excellent selection of relevant metrics. For this reason these metrics are also gathered by Sequali.
  • The matplotlib team for their excellent work on colormaps. Their work was an inspiration for how to present the data and their RdBu colormap is used to represent quality score data. Check their writings on colormaps for a good introduction.
  • Wouter de Coster for his excellent post on how to correctly average phred scores.
  • Marcel Martin for providing very extensive feedback.

License

This project is licensed under the GNU Affero General Public License v3. Mainly to avoid commercial parties from using it without notifying the users that they can run it themselves. If you want to include code from Sequali in your open source project, but it is not compatible with the AGPL, please contact me and we can discuss a separate license.