Read summary

Contamination

These two plots show mapping summaries that can highlight potential contamination issues.

The first plot shows the percentage of reads that either map to any combined reference sequence or are unmapped.

The second plot breaks down the the alignment numbers into the specific references (host, helper plasmid, Rep-Cap plasmid, and transgene plasmid).

Truncations

This plot illustrates the frequency of start and end positions of alignments that map completely within the transgene plasmid ITR-ITR region, helping to identify potential truncation hotspots.

ITR-ITR coverage

For each transgene reference, sequencing depth is calculated for both forward and reverse mapping alignments.

AAV Structures

The numbers of different of the various AAV transgene genome types identified in the sample(s) are summarised here.

A detailed report containing more granular genome type assignments per read can be found at: `output/<sample_id>/<sample_id>_aav_per_read_info.tsv`

Assigned_genome_type count percentage sample_id
Backbone contamination 6840 68.4 sample_1
Partial ssAAV 2620 26.2 sample_1
Full_ssAAV 511 5.11 sample_1
Unknown 29 0.29 sample_1
Assigned_genome_type count percentage sample_id
Backbone contamination 6840 68.4 sample_2
Partial ssAAV 2620 26.2 sample_2
Full_ssAAV 511 5.11 sample_2
Unknown 29 0.29 sample_2

Metadata

Key Value
sample sample_1
type test_sample
barcode sample_1
Key Value
sample sample_2
type test_sample
barcode sample_2

Software versions

Name Version
pysam 0.21.0
minimap2 2.26-r1175
seqkit v2.5.1
samtools 1.20
bcftools 1.16
bedtools v2.31.0

Workflow parameters

Key Value
fastq wf-aav-qc/data/test_data/fastq
bam None
ref_host wf-aav-qc/data/test_data/cell_line.fasta.gz
ref_helper wf-aav-qc/data/test_data/helper.fasta
ref_rep_cap wf-aav-qc/data/test_data/repcap.fasta
ref_transgene_plasmid wf-aav-qc/data/test_data/transgene.fasta
transgene_bed None
itr1_start 11
itr1_end 156
itr2_start 2156
itr2_end 2286
non_transgene_refs None
threads 4
out_dir wf-aav-qc
sample None
sample_sheet None
analyse_unclassified False
itr_fl_threshold 100
itr_backbone_threshold 20
symmetry_threshold 10
output_genometype_bams False
override_basecaller_cfg None
igv False
store_dir wf-aav-qc/store_dir