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 Mapped29237997.46%
Wrongly Mapped51521.717%
Not Mapped24690.823%
Total300000100.0%
Read Failure Statistics
Not mapped24690.823%
Missing in mapper output00.0%
Mapped to wrong chromosome41691.39%
Mapped to wrong position9830.328%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.987100
Precision0.982700
Recall0.991600
Timing
Raw Mapping Time63.312s
Effective Mapping Time63.180s
Effective Init Time0.132s
Effective Time MeasureCPU
Mapping Time (Wall)87.469s
Mapping Time (CPU)63.312s
Mapping Time (CPU User)61.956s
Mapping Time (CPU System)1.356s
Init Time (Wall)0.348s
Init Time (CPU)0.132s
Init Time (CPU User)0.012s
Init Time (CPU System)0.120s
Additional Information
Mapper Memory Usage152 MB
Total Test Runtime (Wall)207.791s
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
memory48836000
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.808 sec; CPU: 0.032 sec
[bwa_aln_core] convert to sequence coordinate... 0.06 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: 3.090 sec; CPU: 0.080 sec
systime0.112
time4.95787215233
usrtime0.016
working_directory/project/teaser/genometeaser/tests_generated/316281ab4f84058f3c27332d7303df8f

/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
memory48836000
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.133 sec; CPU: 0.044 sec
[bwa_aln_core] convert to sequence coordinate... 0.06 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.143 sec; CPU: 0.076 sec
systime0.12
time0.347581148148
usrtime0.012
working_directory/project/teaser/genometeaser/tests_generated/316281ab4f84058f3c27332d7303df8f

/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
memory152968000
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... 47.65 sec
[bwa_aln_core] write to the disk... 0.12 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 5.09 sec
[bwa_aln_core] write to the disk... 0.01 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: 46.245 sec; CPU: 53.712 sec
[bwa_aln_core] convert to sequence coordinate... 3.85 sec
[bwa_aln_core] refine gapped alignments... 1.75 sec
[bwa_aln_core] print alignments... 1.82 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.23 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: 41.132 sec; CPU: 9.564 sec
systime1.356
time87.468547821
usrtime61.956
working_directory/project/teaser/genometeaser/tests_generated/316281ab4f84058f3c27332d7303df8f