Input data and parameters 

QualiMap command line

qualimap bamqc -bam ../mtuberculosis2/2-alignment/pathogen/bam/mycobacteriumTuberculosis/SRR25787976_hisat2.bam -nw 400 -hm 3

Alignment

Command line: "/home/dguevara/pipeline/src/../tools/hisat2-2.1.0/hisat2-align-s --wrapper basic-0 -x ../mtuberculosis2/2-alignment/pathogen/indices/hisat2/mycobacteriumTuberculosis -S ../mtuberculosis2/2-alignment/pathogen/sam/mycobacteriumTuberculosis/SRR25787976_hisat2.sam -p 16 -1 ../mtuberculosis2/2-alignment/host/fastq/SRR25787976_1.fastq -2 ../mtuberculosis2/2-alignment/host/fastq/SRR25787976_2.fastq"
Draw chromosome limits: no
Analyze overlapping paired-end reads: no
Program: hisat2 (2.1.0)
Analysis date: Thu Feb 01 04:51:15 CST 2024
Size of a homopolymer: 3
Skip duplicate alignments: no
Number of windows: 400
BAM file: ../mtuberculosis2/2-alignment/pathogen/bam/mycobacteriumTuberculosis/SRR25787976_hisat2.bam

Summary 

Globals

Reference size 4,411,532
Number of reads 33,341,480
Mapped reads 32,302,443 / 96.88%
Unmapped reads 1,039,037 / 3.12%
Mapped paired reads 32,302,443 / 96.88%
Mapped reads, first in pair 16,162,819 / 48.48%
Mapped reads, second in pair 16,139,624 / 48.41%
Mapped reads, both in pair 32,027,228 / 96.06%
Mapped reads, singletons 275,215 / 0.83%
Read min/max/mean length 100 / 100 / 100
Duplicated reads (estimated) 30,043,140 / 90.11%
Duplication rate 63.86%
Clipped reads 2,667,928 / 8%

ACGT Content

Number/percentage of A's 630,135,338 / 19.57%
Number/percentage of C's 946,950,944 / 29.41%
Number/percentage of T's 583,281,943 / 18.11%
Number/percentage of G's 1,059,747,664 / 32.91%
Number/percentage of N's 180,006,970,651 / 5590.08%
GC Percentage 62.32%

Coverage

Mean 41,533.6853
Standard Deviation 64,951.3363

Mapping Quality

Mean Mapping Quality 23.44

Insert size

Mean 10,824.18
Standard Deviation 146,498.52
P25/Median/P75 139 / 169 / 199

Mismatches and indels

General error rate 0%
Mismatches 3,546,495
Insertions 8,380
Mapped reads with at least one insertion 0.03%
Deletions 69,331
Mapped reads with at least one deletion 0.21%
Homopolymer indels 55.34%

Chromosome stats

Name Length Mapped bases Mean coverage Standard deviation
AL123456.3 4411532 183227181772 41,533.6853 64,951.3363

Coverage across reference 

Coverage Histogram 

Coverage Histogram (0-50X) 

Genome Fraction Coverage 

Duplication Rate Histogram 

Mapped Reads Nucleotide Content 

Mapped Reads GC-content Distribution 

Mapped Reads Clipping Profile 

Homopolymer Indels 

Mapping Quality Across Reference 

Mapping Quality Histogram 

Insert Size Across Reference 

Insert Size Histogram