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 Mapped48641393.029%
Wrongly Mapped311715.962%
Not Mapped52761.009%
Total522860100.0%
Read Failure Statistics
Not mapped52761.009%
Missing in mapper output00.0%
Mapped to wrong chromosome232194.441%
Mapped to wrong position79521.521%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.963900
Precision0.939800
Recall0.989300
Timing
Raw Mapping Time206.752s
Effective Mapping Time199.640s
Effective Init Time7.112s
Effective Time MeasureCPU
Mapping Time (Wall)96.504s
Mapping Time (CPU)206.752s
Mapping Time (CPU User)190.704s
Mapping Time (CPU System)16.048s
Init Time (Wall)7.243s
Init Time (CPU)7.112s
Init Time (CPU User)0.156s
Init Time (CPU System)6.956s
Additional Information
Mapper Memory Usage4730 MB
Total Test Runtime (Wall)639.764s
Mapper Command Line: 
/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads.fastq > out_bwa.sam
Errors and Warnings
No problems were encountered.
Subprocess Log

/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads_base.fastq > out_bwa.sam

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads_base.fastq > out_bwa.sam
memory4598912000
return0
status1
stderr
None
stdout
[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_aln_core] calculate SA coordinate... 0.08 sec
[bwa_aln_core] write to the disk... 0.00 sec
[bwa_aln_core] 1 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa aln -t 4 /project/teaser/genometeaser/references/H_sapiens.fa reads_base.fastq
[main] Real time: 267.029 sec; CPU: 3.664 sec
[bwa_aln_core] convert to sequence coordinate... 4.60 sec
[bwa_aln_core] refine gapped alignments... 0.88 sec
[bwa_aln_core] print alignments... 0.00 sec
[bwa_aln_core] 1 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads_base.fastq
[main] Real time: 203.062 sec; CPU: 5.484 sec
systime9.012
time470.153521061
usrtime0.152
working_directory/project/teaser/genometeaser/tests_generated/73b7a31da9ec873be5a7818478a53cf8

/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads_base.fastq > out_bwa.sam

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads_base.fastq > out_bwa.sam
memory4598884000
return0
status1
stderr
None
stdout
[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_aln_core] calculate SA coordinate... 0.06 sec
[bwa_aln_core] write to the disk... 0.00 sec
[bwa_aln_core] 1 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa aln -t 4 /project/teaser/genometeaser/references/H_sapiens.fa reads_base.fastq
[main] Real time: 2.597 sec; CPU: 2.552 sec
[bwa_aln_core] convert to sequence coordinate... 3.87 sec
[bwa_aln_core] refine gapped alignments... 0.67 sec
[bwa_aln_core] print alignments... 0.00 sec
[bwa_aln_core] 1 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads_base.fastq
[main] Real time: 4.584 sec; CPU: 4.540 sec
systime6.956
time7.2432410717
usrtime0.156
working_directory/project/teaser/genometeaser/tests_generated/73b7a31da9ec873be5a7818478a53cf8

/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads.fastq > out_bwa.sam

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads.fastq > out_bwa.sam
memory4730420000
return0
status1
stderr
None
stdout
[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_aln_core] calculate SA coordinate... 87.58 sec
[bwa_aln_core] write to the disk... 0.15 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 92.61 sec
[bwa_aln_core] write to the disk... 0.14 sec
[bwa_aln_core] 522860 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa aln -t 4 /project/teaser/genometeaser/references/H_sapiens.fa reads.fastq
[main] Real time: 64.136 sec; CPU: 184.600 sec
[bwa_aln_core] convert to sequence coordinate... 6.38 sec
[bwa_aln_core] refine gapped alignments... 2.63 sec
[bwa_aln_core] print alignments... 1.32 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] convert to sequence coordinate... 6.35 sec
[bwa_aln_core] refine gapped alignments... 2.63 sec
[bwa_aln_core] print alignments... 1.28 sec
[bwa_aln_core] 522860 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam.bwa reads.fastq
[main] Real time: 32.236 sec; CPU: 22.116 sec
systime16.048
time96.5041370392
usrtime190.704
working_directory/project/teaser/genometeaser/tests_generated/73b7a31da9ec873be5a7818478a53cf8