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 Mapped29251797.506%
Wrongly Mapped50681.689%
Not Mapped24150.805%
Total300000100.0%
Read Failure Statistics
Not mapped24150.805%
Missing in mapper output00.0%
Mapped to wrong chromosome41141.371%
Mapped to wrong position9540.318%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.987400
Precision0.983000
Recall0.991800
Timing
Raw Mapping Time44.500s
Effective Mapping Time44.400s
Effective Init Time0.100s
Effective Time MeasureCPU
Mapping Time (Wall)35.945s
Mapping Time (CPU)44.500s
Mapping Time (CPU User)43.268s
Mapping Time (CPU System)1.232s
Init Time (Wall)0.231s
Init Time (CPU)0.100s
Init Time (CPU User)0.012s
Init Time (CPU System)0.088s
Additional Information
Mapper Memory Usage152 MB
Total Test Runtime (Wall)87.399s
Mapper Command Line: 
/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta 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/Saccharomyces_cerevisiae.fasta reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads_base.fastq > out_bwa.sam

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads_base.fastq > out_bwa.sam
memory48832000
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.00 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/Saccharomyces_cerevisiae.fasta reads_base.fastq
[main] Real time: 1.810 sec; CPU: 0.040 sec
[bwa_aln_core] convert to sequence coordinate... 0.05 sec
[bwa_aln_core] refine gapped alignments... 0.01 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/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads_base.fastq
[main] Real time: 1.035 sec; CPU: 0.072 sec
systime0.104
time2.88079881668
usrtime0.02
working_directory/project/teaser/genometeaser/tests_generated/37932f46d81271bed2cae5ca41c9d972

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

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads_base.fastq > out_bwa.sam
memory48832000
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.00 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/Saccharomyces_cerevisiae.fasta reads_base.fastq
[main] Real time: 0.059 sec; CPU: 0.028 sec
[bwa_aln_core] convert to sequence coordinate... 0.04 sec
[bwa_aln_core] refine gapped alignments... 0.01 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/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads_base.fastq
[main] Real time: 0.131 sec; CPU: 0.060 sec
systime0.088
time0.231396913528
usrtime0.012
working_directory/project/teaser/genometeaser/tests_generated/37932f46d81271bed2cae5ca41c9d972

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

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads.fastq > out_bwa.sam
memory152964000
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... 31.48 sec
[bwa_aln_core] write to the disk... 0.09 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 4.74 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 300000 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa aln -t 4 /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta reads.fastq
[main] Real time: 17.320 sec; CPU: 37.112 sec
[bwa_aln_core] convert to sequence coordinate... 3.20 sec
[bwa_aln_core] refine gapped alignments... 1.16 sec
[bwa_aln_core] print alignments... 0.94 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] convert to sequence coordinate... 0.53 sec
[bwa_aln_core] refine gapped alignments... 0.25 sec
[bwa_aln_core] print alignments... 0.26 sec
[bwa_aln_core] 300000 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/Saccharomyces_cerevisiae.fasta out_bwa.sam.bwa reads.fastq
[main] Real time: 18.514 sec; CPU: 7.344 sec
systime1.232
time35.9451630116
usrtime43.268
working_directory/project/teaser/genometeaser/tests_generated/37932f46d81271bed2cae5ca41c9d972