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 Mapped52478287.63%
Wrongly Mapped6790211.339%
Not Mapped61741.031%
Total598858100.0%
Read Failure Statistics
Not mapped61741.031%
Missing in mapper output00.0%
Mapped to wrong chromosome401986.712%
Mapped to wrong position277044.626%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.934100
Precision0.885400
Recall0.988400
Timing
Raw Mapping Time109.824s
Effective Mapping Time109.448s
Effective Init Time0.376s
Effective Time MeasureCPU
Mapping Time (Wall)58.975s
Mapping Time (CPU)109.824s
Mapping Time (CPU User)107.020s
Mapping Time (CPU System)2.804s
Init Time (Wall)0.552s
Init Time (CPU)0.376s
Init Time (CPU User)0.024s
Init Time (CPU System)0.352s
Additional Information
Mapper Memory Usage345 MB
Total Test Runtime (Wall)154.330s
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
memory220292000
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: 12.389 sec; CPU: 0.176 sec
[bwa_aln_core] convert to sequence coordinate... 0.22 sec
[bwa_aln_core] refine gapped alignments... 0.05 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: 9.691 sec; CPU: 0.276 sec
systime0.46
time22.1680150032
usrtime0.012
working_directory/project/teaser/genometeaser/tests_generated/9df2b6f30934a6f088552bf04d77f27b

/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
memory220292000
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.150 sec; CPU: 0.120 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.278 sec; CPU: 0.232 sec
systime0.352
time0.552062988281
usrtime0.024
working_directory/project/teaser/genometeaser/tests_generated/9df2b6f30934a6f088552bf04d77f27b

/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
memory345240000
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... 41.41 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.63 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... 10.81 sec
[bwa_aln_core] write to the disk... 0.04 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: 31.696 sec; CPU: 94.124 sec
[bwa_aln_core] convert to sequence coordinate... 2.59 sec
[bwa_aln_core] refine gapped alignments... 2.16 sec
[bwa_aln_core] print alignments... 1.29 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.13 sec
[bwa_aln_core] print alignments... 1.27 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_aln_core] convert to sequence coordinate... 0.86 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: 27.092 sec; CPU: 15.660 sec
systime2.804
time58.9752280712
usrtime107.02
working_directory/project/teaser/genometeaser/tests_generated/9df2b6f30934a6f088552bf04d77f27b