BAM_QC_example
A modular tool to aggregate results from bioinformatics analyses across many samples into a single report.
example-out/03_results
General Statistics
| Sample Name | ≥ 1X | ≥ 5X | ≥ 10X | ≥ 30X | ≥ 50X | Median | Mean Cov. | Min Cov. | Max Cov. | Mb Total Coverage Bases | Genome length | Error rate | Non-primary | Reads mapped | % Mapped | % Proper pairs | % MapQ 0 reads | Total seqs | Mean insert | Reads | Reads mapped | % Reads mapped | Reads | Bases | Coverage | Mean depth | Mean BQ | Mean MQ |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| tiny | 0.00% | 0.0M | 0.0M | 100.0% | 0.0% | 0.0% | 0.0M | 0.0bp | 0.0M | 0.0M | 100.0% | |||||||||||||||||
| tiny.coverage | 0.0M | 0.0Mb | 20.0% | 0.2x | 40.0 | 60.0 | ||||||||||||||||||||||
| tiny.genome | 20.0% | 0.0% | 0.0% | 0.0% | 0.0% | 0X | 0.2X | 0.0X | 1.0X | 0.0Mb | 200 | |||||||||||||||||
| tiny.regions | 51.0% | 0.0% | 0.0% | 0.0% | 0.0% | 1X | 0.5X | 0.0X | 1.0X | 0.0Mb | 79 |
Mosdepth
Fast BAM/CRAM depth calculation for WGS, exome, or targeted sequencing.https://github.com/brentp/mosdepthDOI: 10.1093/bioinformatics/btx699
Cumulative coverage distribution
Proportion of bases in the reference genome with, at least, a given depth of coverage. Note that for 2 samples, a BED file was provided, so the data was calculated across those regions. For 1 samples, it's calculated across the entire genome length. 1 samples have both global and region reports, and we are showing the data for regions
For a set of DNA or RNA reads mapped to a reference sequence, such as a genome or transcriptome, the depth of coverage at a given base position is the number of high-quality reads that map to the reference at that position, while the breadth of coverage is the fraction of the reference sequence to which reads have been mapped with at least a given depth of coverage (Sims et al. 2014).
Defining coverage breadth in terms of coverage depth is useful, because sequencing experiments typically require a specific minimum depth of coverage over the region of interest (Sims et al. 2014), so the extent of the reference sequence that is amenable to analysis is constrained to lie within regions that have sufficient depth. With inadequate sequencing breadth, it can be difficult to distinguish the absence of a biological feature (such as a gene) from a lack of data (Green 2007).
For increasing coverage depths (1×, 2×, …, N×), coverage breadth is calculated as the percentage of the reference sequence that is covered by at least that number of reads, then plots coverage breadth (y-axis) against coverage depth (x-axis). This plot shows the relationship between sequencing depth and breadth for each read dataset, which can be used to gauge, for example, the likely effect of a minimum depth filter on the fraction of a genome available for analysis.
Average coverage per contig
Average coverage per contig or chromosome
Samtools
1.23.1
Toolkit for interacting with BAM/CRAM files.http://www.htslib.orgDOI: 10.1093/bioinformatics/btp352
Percent mapped
Alignment metrics from samtools stats; mapped vs. unmapped reads vs. reads mapped with MQ0.
For a set of samples that have come from the same multiplexed library, similar numbers of reads for each sample are expected. Large differences in numbers might indicate issues during the library preparation process. Whilst large differences in read numbers may be controlled for in downstream processings (e.g. read count normalisation), you may wish to consider whether the read depths achieved have fallen below recommended levels depending on the applications.
Low alignment rates could indicate contamination of samples (e.g. adapter sequences), low sequencing quality or other artefacts. These can be further investigated in the sequence level QC (e.g. from FastQC).
Reads mapped with MQ0 often indicate that the reads are ambiguously mapped to multiple locations in the reference sequence. This can be due to repetitive regions in the genome, the presence of alternative contigs in the reference, or due to reads that are too short to be uniquely mapped. These reads are often filtered out in downstream analyses.
Alignment stats
This module parses the output from samtools stats. All numbers in millions.
Flagstat
This module parses the output from samtools flagstat
Flagstat: Percentage of total
This module parses the output from samtools flagstat
Mapped reads per contig
The samtools idxstats tool counts the number of mapped reads per chromosome / contig. Chromosomes with < 0.1% of the total aligned reads are omitted from this plot.
Coverage: global stats
Stats parsed from samtools coverage output, and summarized (added up or weighted-averaged) across all regions.
| Sample Name | Reads | Bases | Coverage | Mean depth | Mean BQ | Mean MQ |
|---|---|---|---|---|---|---|
| tiny.coverage | 0.0M | 0.0Mb | 20.0% | 0.2x | 40.0 | 60.0 |
Coverage: stats per region
Per-region stats parsed from samtools coverage output.
Software Versions
Software Versions lists versions of software tools extracted from file contents.
| Software | Version |
|---|---|
| Samtools | 1.23.1 |