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 Mapped52349787.416%
Wrongly Mapped6916111.549%
Not Mapped62001.035%
Total598858100.0%
Read Failure Statistics
Not mapped62001.035%
Missing in mapper output00.0%
Mapped to wrong chromosome409466.837%
Mapped to wrong position282154.711%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.932900
Precision0.883300
Recall0.988300
Timing
Raw Mapping Time107.348s
Effective Mapping Time106.964s
Effective Init Time0.384s
Effective Time MeasureCPU
Mapping Time (Wall)62.082s
Mapping Time (CPU)107.348s
Mapping Time (CPU User)104.796s
Mapping Time (CPU System)2.552s
Init Time (Wall)0.511s
Init Time (CPU)0.384s
Init Time (CPU User)0.024s
Init Time (CPU System)0.360s
Additional Information
Mapper Memory Usage345 MB
Total Test Runtime (Wall)148.938s
Mapper Command Line: 
/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/D_melanogaster.fasta reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/D_melanogaster.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/D_melanogaster.fasta reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/D_melanogaster.fasta out_bwa.sam.bwa reads_base.fastq > out_bwa.sam

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/D_melanogaster.fasta reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/D_melanogaster.fasta out_bwa.sam.bwa reads_base.fastq > out_bwa.sam
memory220216000
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/D_melanogaster.fasta reads_base.fastq
[main] Real time: 0.152 sec; CPU: 0.136 sec
[bwa_aln_core] convert to sequence coordinate... 0.19 sec
[bwa_aln_core] refine gapped alignments... 0.03 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/D_melanogaster.fasta out_bwa.sam.bwa reads_base.fastq
[main] Real time: 0.268 sec; CPU: 0.232 sec
systime0.36
time0.455067873001
usrtime0.02
working_directory/project/teaser/genometeaser/tests_generated/58281c12b324089e734fe0553bf2f5c0

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

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/D_melanogaster.fasta reads_base.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/D_melanogaster.fasta out_bwa.sam.bwa reads_base.fastq > out_bwa.sam
memory220216000
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/D_melanogaster.fasta reads_base.fastq
[main] Real time: 0.158 sec; CPU: 0.136 sec
[bwa_aln_core] convert to sequence coordinate... 0.19 sec
[bwa_aln_core] refine gapped alignments... 0.03 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/D_melanogaster.fasta out_bwa.sam.bwa reads_base.fastq
[main] Real time: 0.308 sec; CPU: 0.232 sec
systime0.36
time0.510673999786
usrtime0.024
working_directory/project/teaser/genometeaser/tests_generated/58281c12b324089e734fe0553bf2f5c0

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

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/D_melanogaster.fasta reads.fastq -t 4 > out_bwa.sam.bwa; /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/D_melanogaster.fasta out_bwa.sam.bwa reads.fastq > out_bwa.sam
memory345168000
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... 37.61 sec
[bwa_aln_core] write to the disk... 0.14 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 39.94 sec
[bwa_aln_core] write to the disk... 0.12 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 11.86 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 598858 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa aln -t 4 /project/teaser/genometeaser/references/D_melanogaster.fasta reads.fastq
[main] Real time: 32.043 sec; CPU: 91.628 sec
[bwa_aln_core] convert to sequence coordinate... 2.63 sec
[bwa_aln_core] refine gapped alignments... 2.16 sec
[bwa_aln_core] print alignments... 1.28 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_aln_core] convert to sequence coordinate... 2.62 sec
[bwa_aln_core] refine gapped alignments... 2.14 sec
[bwa_aln_core] print alignments... 1.26 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_aln_core] convert to sequence coordinate... 0.88 sec
[bwa_aln_core] refine gapped alignments... 0.62 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 598858 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa samse /project/teaser/genometeaser/references/D_melanogaster.fasta out_bwa.sam.bwa reads.fastq
[main] Real time: 29.801 sec; CPU: 15.680 sec
systime2.552
time62.0824239254
usrtime104.796
working_directory/project/teaser/genometeaser/tests_generated/58281c12b324089e734fe0553bf2f5c0