Mapping Quality Thresholds

The plot below shows the percentages of correctly and wrongly mapped reads for all mapping quality thresholds for this mapper. The values at threshold 0 therefore correspond to the unfiltered results.

Basic Statistics
Correctly Mapped12399751.255%
Wrongly Mapped109904.543%
Not Mapped10693444.202%
Total241921100.0%
Read Failure Statistics
Not mapped10693444.202%
Missing in mapper output00.0%
Mapped to wrong chromosome59422.456%
Mapped to wrong position50482.087%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.677700
Precision0.918600
Recall0.536900
Timing
Raw Mapping Time33.701s
Effective Mapping Time32.771s
Effective Init Time0.930s
Effective Time MeasureWall clock
Mapping Time (Wall)33.701s
Mapping Time (CPU)124.997s
Mapping Time (CPU User)123.470s
Mapping Time (CPU System)1.527s
Init Time (Wall)0.930s
Init Time (CPU)0.894s
Init Time (CPU User)0.042s
Init Time (CPU System)0.852s
Additional Information
Mapper Memory Usage3827 MB
Total Test Runtime (Wall)5676.894s
Mapper Command Line: 
/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads.fastq -S out_bowtie2.sam
Errors and Warnings
No problems were encountered.
Subprocess Log

/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2-build /project2/NextGenMapTest/Teaser/references/m_musculus.fasta /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt

command/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2-build /project2/NextGenMapTest/Teaser/references/m_musculus.fasta /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt
return1
status1
stderr
None
stdout
Settings:
Output files: "/project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt.*.bt2l"
Line rate: 7 (line is 128 bytes)
Lines per side: 1 (side is 128 bytes)
Offset rate: 4 (one in 16)
FTable chars: 10
Strings: unpacked
Max bucket size: default
Max bucket size, sqrt multiplier: default
Max bucket size, len divisor: 4
Difference-cover sample period: 1024
Endianness: little
Actual local endianness: little
Sanity checking: disabled
Assertions: disabled
Random seed: 0
Sizeofs: void*:8, int:4, long:8, size_t:8
Input files DNA, FASTA:
/project2/NextGenMapTest/Teaser/references/m_musculus.fasta
Building a LARGE index
Reading reference sizes
Time reading reference sizes: 00:00:33
Calculating joined length
Writing header
Reserving space for joined string
Joining reference sequences
Time to join reference sequences: 00:00:28
bmax according to bmaxDivN setting: 664765538
Using parameters --bmax 498574154 --dcv 1024
Doing ahead-of-time memory usage test
Passed! Constructing with these parameters: --bmax 498574154 --dcv 1024
Constructing suffix-array element generator
Building DifferenceCoverSample
Building sPrime
Building sPrimeOrder
V-Sorting samples
V-Sorting samples time: 00:01:00
Allocating rank array
Ranking v-sort output
Ranking v-sort output time: 00:00:11
Invoking Larsson-Sadakane on ranks
Invoking Larsson-Sadakane on ranks time: 00:00:28
Sanity-checking and returning
Building samples
Reserving space for 12 sample suffixes
Generating random suffixes
QSorting 12 sample offsets, eliminating duplicates
QSorting sample offsets, eliminating duplicates time: 00:00:00
Multikey QSorting 12 samples
(Using difference cover)
Multikey QSorting samples time: 00:00:00
Calculating bucket sizes
Binary sorting into buckets
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Binary sorting into buckets time: 00:01:26
Splitting and merging
Splitting and merging time: 00:00:00
Split 2, merged 7; iterating...
Binary sorting into buckets
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Binary sorting into buckets time: 00:01:19
Splitting and merging
Splitting and merging time: 00:00:00
Split 1, merged 0; iterating...
Binary sorting into buckets
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Binary sorting into buckets time: 00:01:22
Splitting and merging
Splitting and merging time: 00:00:00
Avg bucket size: 3.79866e+08 (target: 498574153)
Converting suffix-array elements to index image
Allocating ftab, absorbFtab
Entering Ebwt loop
Getting block 1 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:23
Sorting block of length 200501613
(Using difference cover)
Sorting block time: 00:02:08
Returning block of 200501614
Getting block 2 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:28
Sorting block of length 400805978
(Using difference cover)
Sorting block time: 00:04:32
Returning block of 400805979
Getting block 3 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:34
Sorting block of length 461882454
(Using difference cover)
Sorting block time: 00:05:24
Returning block of 461882455
Getting block 4 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:40
Sorting block of length 443285039
(Using difference cover)
Sorting block time: 00:05:12
Returning block of 443285040
Getting block 5 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:38
Sorting block of length 333625971
(Using difference cover)
Sorting block time: 00:03:45
Returning block of 333625972
Getting block 6 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:39
Sorting block of length 363005442
(Using difference cover)
Sorting block time: 00:04:09
Returning block of 363005443
Getting block 7 of 7
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:28
Sorting block of length 455955649
(Using difference cover)
Sorting block time: 00:05:23
Returning block of 455955650
Exited Ebwt loop
fchr[A]: 0
fchr[C]: 775105804
fchr[G]: 1329074091
fchr[T]: 1883091452
fchr[$]: 2659062152
Exiting Ebwt::buildToDisk()
Returning from initFromVector
Wrote 894765894 bytes to primary EBWT file: /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt.1.bt2l
Wrote 1329531084 bytes to secondary EBWT file: /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt.2.bt2l
Re-opening _in1 and _in2 as input streams
Returning from Ebwt constructor
Headers:
len: 2659062152
bwtLen: 2659062153
sz: 664765538
bwtSz: 664765539
lineRate: 7
offRate: 4
offMask: 0xfffffffffffffff0
ftabChars: 10
eftabLen: 20
eftabSz: 160
ftabLen: 1048577
ftabSz: 8388616
offsLen: 166191385
offsSz: 1329531080
lineSz: 128
sideSz: 128
sideBwtSz: 96
sideBwtLen: 384
numSides: 6924642
numLines: 6924642
ebwtTotLen: 886354176
ebwtTotSz: 886354176
color: 0
reverse: 0
Total time for call to driver() for forward index: 00:46:16
Reading reference sizes
Time reading reference sizes: 00:00:25
Calculating joined length
Writing header
Reserving space for joined string
Joining reference sequences
Time to join reference sequences: 00:00:28
Time to reverse reference sequence: 00:00:02
bmax according to bmaxDivN setting: 664765538
Using parameters --bmax 498574154 --dcv 1024
Doing ahead-of-time memory usage test
Passed! Constructing with these parameters: --bmax 498574154 --dcv 1024
Constructing suffix-array element generator
Building DifferenceCoverSample
Building sPrime
Building sPrimeOrder
V-Sorting samples
V-Sorting samples time: 00:01:00
Allocating rank array
Ranking v-sort output
Ranking v-sort output time: 00:00:11
Invoking Larsson-Sadakane on ranks
Invoking Larsson-Sadakane on ranks time: 00:00:28
Sanity-checking and returning
Building samples
Reserving space for 12 sample suffixes
Generating random suffixes
QSorting 12 sample offsets, eliminating duplicates
QSorting sample offsets, eliminating duplicates time: 00:00:00
Multikey QSorting 12 samples
(Using difference cover)
Multikey QSorting samples time: 00:00:00
Calculating bucket sizes
Binary sorting into buckets
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Binary sorting into buckets time: 00:01:33
Splitting and merging
Splitting and merging time: 00:00:00
Split 3, merged 6; iterating...
Binary sorting into buckets
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Binary sorting into buckets time: 00:01:25
Splitting and merging
Splitting and merging time: 00:00:00
Split 1, merged 3; iterating...
Binary sorting into buckets
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Binary sorting into buckets time: 00:01:18
Splitting and merging
Splitting and merging time: 00:00:00
Avg bucket size: 3.32383e+08 (target: 498574153)
Converting suffix-array elements to index image
Allocating ftab, absorbFtab
Entering Ebwt loop
Getting block 1 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:26
Sorting block of length 486788261
(Using difference cover)
Sorting block time: 00:05:38
Returning block of 486788262
Getting block 2 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:33
Sorting block of length 371944694
(Using difference cover)
Sorting block time: 00:04:15
Returning block of 371944695
Getting block 3 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:34
Sorting block of length 311926459
(Using difference cover)
Sorting block time: 00:03:31
Returning block of 311926460
Getting block 4 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:39
Sorting block of length 283955416
(Using difference cover)
Sorting block time: 00:03:13
Returning block of 283955417
Getting block 5 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:37
Sorting block of length 335651746
(Using difference cover)
Sorting block time: 00:03:49
Returning block of 335651747
Getting block 6 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:40
Sorting block of length 182814825
(Using difference cover)
Sorting block time: 00:02:00
Returning block of 182814826
Getting block 7 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:38
Sorting block of length 376753567
(Using difference cover)
Sorting block time: 00:04:21
Returning block of 376753568
Getting block 8 of 8
Reserving size (498574154) for bucket
Calculating Z arrays
Calculating Z arrays time: 00:00:00
Entering block accumulator loop:
10%
20%
30%
40%
50%
60%
70%
80%
90%
100%
Block accumulator loop time: 00:00:27
Sorting block of length 309227177
(Using difference cover)
Sorting block time: 00:03:34
Returning block of 309227178
Exited Ebwt loop
fchr[A]: 0
fchr[C]: 775105804
fchr[G]: 1329074091
fchr[T]: 1883091452
fchr[$]: 2659062152
Exiting Ebwt::buildToDisk()
Returning from initFromVector
Wrote 894765894 bytes to primary EBWT file: /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt.rev.1.bt2l
Wrote 1329531084 bytes to secondary EBWT file: /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt.rev.2.bt2l
Re-opening _in1 and _in2 as input streams
Returning from Ebwt constructor
Headers:
len: 2659062152
bwtLen: 2659062153
sz: 664765538
bwtSz: 664765539
lineRate: 7
offRate: 4
offMask: 0xfffffffffffffff0
ftabChars: 10
eftabLen: 20
eftabSz: 160
ftabLen: 1048577
ftabSz: 8388616
offsLen: 166191385
offsSz: 1329531080
lineSz: 128
sideSz: 128
sideBwtSz: 96
sideBwtLen: 384
numSides: 6924642
numLines: 6924642
ebwtTotLen: 886354176
ebwtTotSz: 886354176
color: 0
reverse: 1
Total time for backward call to driver() for mirror index: 00:46:49
systime0
time5586.61908579
usrtime0
working_directory/project2/NextGenMapTest/Teaser/tests_generated/M3_n

/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads_base.fastq -S out_bowtie2.sam

command/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads_base.fastq -S out_bowtie2.sam
memory3745320000
return0
status1
stderr
None
stdout
1 reads; of these:
1 (100.00%) were unpaired; of these:
1 (100.00%) aligned 0 times
0 (0.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
0.00% overall alignment rate
systime1.738735
time41.6170399189
usrtime0.043993
working_directory/project2/NextGenMapTest/Teaser/tests_generated/M3_n

/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads_base.fastq -S out_bowtie2.sam

command/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads_base.fastq -S out_bowtie2.sam
memory3743276000
return0
status1
stderr
None
stdout
1 reads; of these:
1 (100.00%) were unpaired; of these:
1 (100.00%) aligned 0 times
0 (0.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
0.00% overall alignment rate
systime0.85187
time0.930344104767
usrtime0.041993
working_directory/project2/NextGenMapTest/Teaser/tests_generated/M3_n

/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads.fastq -S out_bowtie2.sam

command/project2/NextGenMapTest/Teaser/software/bowtie2/bowtie2 -p 4 -x /project2/NextGenMapTest/Teaser/references/m_musculus.fasta_bt -U reads.fastq -S out_bowtie2.sam
memory3827744000
return0
status1
stderr
None
stdout
241921 reads; of these:
241921 (100.00%) were unpaired; of these:
106934 (44.20%) aligned 0 times
118054 (48.80%) aligned exactly 1 time
16933 (7.00%) aligned >1 times
55.80% overall alignment rate
systime1.526767
time33.7012329102
usrtime123.470229
working_directory/project2/NextGenMapTest/Teaser/tests_generated/M3_n