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 Mapped50022495.671%
Wrongly Mapped205833.937%
Not Mapped20530.393%
Total522860100.0%
Read Failure Statistics
Not mapped20530.393%
Missing in mapper output00.0%
Mapped to wrong chromosome134072.564%
Mapped to wrong position71761.372%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.977900
Precision0.960500
Recall0.995900
Timing
Raw Mapping Time263.324s
Effective Mapping Time256.328s
Effective Init Time6.996s
Effective Time MeasureCPU
Mapping Time (Wall)148.605s
Mapping Time (CPU)263.324s
Mapping Time (CPU User)250.492s
Mapping Time (CPU System)12.832s
Init Time (Wall)7.140s
Init Time (CPU)6.996s
Init Time (CPU User)0.180s
Init Time (CPU System)6.816s
Additional Information
Mapper Memory Usage4892 MB
Total Test Runtime (Wall)695.581s
Mapper Command Line: 
/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads1.fastq -t 4 > out_bwa.sam_1.sai; /project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads2.fastq -t 4 > out_bwa.sam_2.sai; /project/teaser/genometeaser/software/bwa sampe /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam_1.sai out_bwa.sam_2.sai reads1.fastq reads2.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
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.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: 266.951 sec; CPU: 3.552 sec
[bwa_aln_core] convert to sequence coordinate... 4.50 sec
[bwa_aln_core] refine gapped alignments... 0.86 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: 202.973 sec; CPU: 5.360 sec
systime8.76
time469.99259901
usrtime0.164
working_directory/project/teaser/genometeaser/tests_generated/24f21f076cb23e5eaf20f9023e376370

/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
memory4598880000
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.535 sec; CPU: 2.504 sec
[bwa_aln_core] convert to sequence coordinate... 3.81 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.535 sec; CPU: 4.480 sec
systime6.816
time7.13957095146
usrtime0.18
working_directory/project/teaser/genometeaser/tests_generated/24f21f076cb23e5eaf20f9023e376370

/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads1.fastq -t 4 > out_bwa.sam_1.sai; /project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads2.fastq -t 4 > out_bwa.sam_2.sai; /project/teaser/genometeaser/software/bwa sampe /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam_1.sai out_bwa.sam_2.sai reads1.fastq reads2.fastq > out_bwa.sam

command/project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads1.fastq -t 4 > out_bwa.sam_1.sai; /project/teaser/genometeaser/software/bwa aln /project/teaser/genometeaser/references/H_sapiens.fa reads2.fastq -t 4 > out_bwa.sam_2.sai; /project/teaser/genometeaser/software/bwa sampe /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam_1.sai out_bwa.sam_2.sai reads1.fastq reads2.fastq > out_bwa.sam
memory4892336000
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... 89.94 sec
[bwa_aln_core] write to the disk... 0.15 sec
[bwa_aln_core] 261430 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 reads1.fastq
[main] Real time: 33.416 sec; CPU: 93.308 sec
[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... 96.13 sec
[bwa_aln_core] write to the disk... 0.15 sec
[bwa_aln_core] 261430 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 reads2.fastq
[main] Real time: 33.352 sec; CPU: 99.536 sec
[bwa_sai2sam_pe_core] convert to sequence coordinate...
[infer_isize] (25, 50, 75) percentile: (80, 100, 121)
[infer_isize] low and high boundaries: 60 and 203 for estimating avg and std
[infer_isize] inferred external isize from 225153 pairs: 102.609 +/- 27.658
[infer_isize] skewness: 0.558; kurtosis: -0.232; ap_prior: 1.00e-05
[infer_isize] inferred maximum insert size: 298 (7.05 sigma)
[bwa_sai2sam_pe_core] time elapses: 59.23 sec
[bwa_sai2sam_pe_core] changing coordinates of 7722 alignments.
[bwa_sai2sam_pe_core] align unmapped mate...
[bwa_paired_sw] 3253 out of 4614 Q17 singletons are mated.
[bwa_paired_sw] 523 out of 28619 Q17 discordant pairs are fixed.
[bwa_sai2sam_pe_core] time elapses: 2.17 sec
[bwa_sai2sam_pe_core] refine gapped alignments... 4.11 sec
[bwa_sai2sam_pe_core] print alignments... 3.27 sec
[bwa_sai2sam_pe_core] 261430 sequences have been processed.
[main] Version: 0.7.16a-r1181
[main] CMD: /project/teaser/genometeaser/software/bwa sampe /project/teaser/genometeaser/references/H_sapiens.fa out_bwa.sam_1.sai out_bwa.sam_2.sai reads1.fastq reads2.fastq
[main] Real time: 81.703 sec; CPU: 70.420 sec
systime12.832
time148.605185986
usrtime250.492
working_directory/project/teaser/genometeaser/tests_generated/24f21f076cb23e5eaf20f9023e376370