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 Mapped48524697.319%
Wrongly Mapped133142.67%
Not Mapped540.011%
Total498614100.0%
Read Failure Statistics
Not mapped540.011%
Missing in mapper output00.0%
Mapped to wrong chromosome53361.07%
Mapped to wrong position79781.6%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.986400
Precision0.973300
Recall0.999900
Timing
Raw Mapping Time36.960s
Effective Mapping Time33.660s
Effective Init Time3.300s
Effective Time MeasureCPU
Mapping Time (Wall)23.226s
Mapping Time (CPU)36.960s
Mapping Time (CPU User)34.076s
Mapping Time (CPU System)2.884s
Init Time (Wall)3.123s
Init Time (CPU)3.300s
Init Time (CPU User)2.364s
Init Time (CPU System)0.936s
Additional Information
Mapper Memory Usage2013 MB
Total Test Runtime (Wall)142.564s
Mapper Command Line: 
/project/teaser/genometeaser/software/ngm/ngm --qry2 reads2.fastq --output out_ngm.sam --qry1 reads1.fastq --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --no-progress
Errors and Warnings
No problems were encountered.
Subprocess Log

/project/teaser/genometeaser/software/ngm/ngm --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas

command/project/teaser/genometeaser/software/ngm/ngm --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas
return1
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-10-16.14:40:36
[CONFIG] Parameter: --affine 0 --argos_min_score 0 --bin_size 2 --block_multiplier 2 --bs_cutoff 6 --bs_mapping 0 --cpu_threads 1 --dualstrand 1 --fast 0 --fast_pairing 0 --force_rlength_check 0 --format 1 --gap_extend_penalty 5 --gap_read_penalty 20 --gap_ref_penalty 20 --hard_clip 0 --keep_tags 0 --kmer 13 --kmer_min 0 --kmer_skip 2 --match_bonus 10 --match_bonus_tc 2 --match_bonus_tt 10 --max_cmrs 2147483647 --max_equal 1 --max_insert_size 1000 --max_polya -1 --max_read_length 0 --min_identity 0.650000 --min_insert_size 0 --min_mq 0 --min_residues 0.500000 --min_score 0.000000 --mismatch_penalty 15 --mode 0 --no_progress 0 --no_unal 0 --ocl_threads 1 --overwrite 1 --pair_score_cutoff 0.900000 --paired 0 --parse_all 1 --pe_delimiter / --qry_count -1 --qry_start 0 --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas --ref_mode -1 --sensitive 0 --silent_clip 0 --skip_mate_check 0 --skip_save 0 --slam_seq 0 --step_count 4 --strata 0 --topn 1 --trim5 0 --update_check 0 --very_fast 0 --very_sensitive 0
[NGM] Wrinting output (SAM) to stdout
[SEQPROV] Reading encoded reference from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-enc.2.ngm
[SEQPROV] Reading 119 Mbp from disk took 5.19s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 42.52s
[PREPROCESS] Max. k-mer frequency set so 100!
[MAIN] Finished building hash table.
[MAIN] No qry file specified. If you want to map single-end data use -q/--qry. If you want to map paired-end data, either use -q/--qry and -p or --qry1 and --qry2.
[UPDATE_CHECK] Your version of NGM is more than 6 months old - a newer version may be available. (For performing an automatic check use --update-check)
systime0
time48.3760159016
usrtime0
working_directory/project/teaser/genometeaser/tests_generated/340a758089e3620a129e06c436f2f4af

/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --qry reads_base.fastq --no-progress

command/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --qry reads_base.fastq --no-progress
memory928620000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-10-16.14:41:24
[CONFIG] Parameter: --affine 0 --argos_min_score 0 --bin_size 2 --block_multiplier 2 --bs_cutoff 6 --bs_mapping 0 --cpu_threads 4 --dualstrand 1 --fast 0 --fast_pairing 0 --force_rlength_check 0 --format 1 --gap_extend_penalty 5 --gap_read_penalty 20 --gap_ref_penalty 20 --hard_clip 0 --keep_tags 0 --kmer 13 --kmer_min 0 --kmer_skip 2 --match_bonus 10 --match_bonus_tc 2 --match_bonus_tt 10 --max_cmrs 2147483647 --max_equal 1 --max_insert_size 1000 --max_polya -1 --max_read_length 0 --min_identity 0.650000 --min_insert_size 0 --min_mq 0 --min_residues 0.500000 --min_score 0.000000 --mismatch_penalty 15 --mode 0 --no_progress 1 --no_unal 0 --ocl_threads 1 --output out_ngm.sam --overwrite 1 --pair_score_cutoff 0.900000 --paired 0 --parse_all 1 --pe_delimiter / --qry reads_base.fastq --qry_count -1 --qry_start 0 --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas --ref_mode -1 --sensitive 0 --silent_clip 0 --skip_mate_check 0 --skip_save 0 --slam_seq 0 --step_count 4 --strata 0 --topn 1 --trim5 0 --update_check 0 --very_fast 0 --very_sensitive 0
[NGM] Opening for output (SAM): out_ngm.sam
[SEQPROV] Reading encoded reference from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-enc.2.ngm
[SEQPROV] Reading 119 Mbp from disk took 0.06s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 0.64s
[PREPROCESS] Max. k-mer frequency set so 100!
[INPUT] Input is single end data.
[INPUT] Opening file reads_base.fastq for reading
[INPUT] Input is Fastq
[INPUT] Estimating parameter from data
[INPUT] Reads found in files: 1
[INPUT] Average read length: 60 (min: 60, max: 62)
[INPUT] Corridor width: 14
[INPUT] Not enough reads to estimate parameter
[INPUT] Sensitivity parameter neither set nor estimated. Falling back to default.
[INPUT] Estimating parameter took 0.004s
[INPUT] Input is Fastq
No protocol specified
[OPENCL] Available platforms: 1
[OPENCL] AMD Accelerated Parallel Processing
[OPENCL] Selecting OpenCl platform: AMD Accelerated Parallel Processing
[OPENCL] Platform: OpenCL 1.2 AMD-APP (1214.3)
[OPENCL] 1 CPU device found.
[OPENCL] Device 0: AMD Opteron(tm) Processor 6348 (Driver: 1214.3 (sse2,avx,fma4))
[OPENCL] 48 CPU cores available.
[MAIN] Alignments computed: 1
[MAIN] Done (1 reads mapped (100.00%), 0 reads not mapped, 1 lines written)(elapsed: 6.702757s)
[UPDATE_CHECK] Your version of NGM is more than 6 months old - a newer version may be available. (For performing an automatic check use --update-check)
systime0.956
time6.82838606834
usrtime2.38
working_directory/project/teaser/genometeaser/tests_generated/340a758089e3620a129e06c436f2f4af

/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --qry reads_base.fastq --no-progress

command/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --qry reads_base.fastq --no-progress
memory930708000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-10-16.14:41:31
[CONFIG] Parameter: --affine 0 --argos_min_score 0 --bin_size 2 --block_multiplier 2 --bs_cutoff 6 --bs_mapping 0 --cpu_threads 4 --dualstrand 1 --fast 0 --fast_pairing 0 --force_rlength_check 0 --format 1 --gap_extend_penalty 5 --gap_read_penalty 20 --gap_ref_penalty 20 --hard_clip 0 --keep_tags 0 --kmer 13 --kmer_min 0 --kmer_skip 2 --match_bonus 10 --match_bonus_tc 2 --match_bonus_tt 10 --max_cmrs 2147483647 --max_equal 1 --max_insert_size 1000 --max_polya -1 --max_read_length 0 --min_identity 0.650000 --min_insert_size 0 --min_mq 0 --min_residues 0.500000 --min_score 0.000000 --mismatch_penalty 15 --mode 0 --no_progress 1 --no_unal 0 --ocl_threads 1 --output out_ngm.sam --overwrite 1 --pair_score_cutoff 0.900000 --paired 0 --parse_all 1 --pe_delimiter / --qry reads_base.fastq --qry_count -1 --qry_start 0 --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas --ref_mode -1 --sensitive 0 --silent_clip 0 --skip_mate_check 0 --skip_save 0 --slam_seq 0 --step_count 4 --strata 0 --topn 1 --trim5 0 --update_check 0 --very_fast 0 --very_sensitive 0
[NGM] Opening for output (SAM): out_ngm.sam
[SEQPROV] Reading encoded reference from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-enc.2.ngm
[SEQPROV] Reading 119 Mbp from disk took 0.05s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 0.59s
[PREPROCESS] Max. k-mer frequency set so 100!
[INPUT] Input is single end data.
[INPUT] Opening file reads_base.fastq for reading
[INPUT] Input is Fastq
[INPUT] Estimating parameter from data
[INPUT] Reads found in files: 1
[INPUT] Average read length: 60 (min: 60, max: 62)
[INPUT] Corridor width: 14
[INPUT] Not enough reads to estimate parameter
[INPUT] Sensitivity parameter neither set nor estimated. Falling back to default.
[INPUT] Estimating parameter took 0.004s
[INPUT] Input is Fastq
No protocol specified
[OPENCL] Available platforms: 1
[OPENCL] AMD Accelerated Parallel Processing
[OPENCL] Selecting OpenCl platform: AMD Accelerated Parallel Processing
[OPENCL] Platform: OpenCL 1.2 AMD-APP (1214.3)
[OPENCL] 1 CPU device found.
[OPENCL] Device 0: AMD Opteron(tm) Processor 6348 (Driver: 1214.3 (sse2,avx,fma4))
[OPENCL] 48 CPU cores available.
[MAIN] Alignments computed: 1
[MAIN] Done (1 reads mapped (100.00%), 0 reads not mapped, 1 lines written)(elapsed: 3.052332s)
[UPDATE_CHECK] Your version of NGM is more than 6 months old - a newer version may be available. (For performing an automatic check use --update-check)
systime0.936
time3.12333989143
usrtime2.364
working_directory/project/teaser/genometeaser/tests_generated/340a758089e3620a129e06c436f2f4af

/project/teaser/genometeaser/software/ngm/ngm --qry2 reads2.fastq --output out_ngm.sam --qry1 reads1.fastq --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --no-progress

command/project/teaser/genometeaser/software/ngm/ngm --qry2 reads2.fastq --output out_ngm.sam --qry1 reads1.fastq --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas -t 4 --no-progress
memory2013200000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-10-16.14:41:35
[CONFIG] Parameter: --affine 0 --argos_min_score 0 --bin_size 2 --block_multiplier 2 --bs_cutoff 6 --bs_mapping 0 --cpu_threads 4 --dualstrand 1 --fast 0 --fast_pairing 0 --force_rlength_check 0 --format 1 --gap_extend_penalty 5 --gap_read_penalty 20 --gap_ref_penalty 20 --hard_clip 0 --keep_tags 0 --kmer 13 --kmer_min 0 --kmer_skip 2 --match_bonus 10 --match_bonus_tc 2 --match_bonus_tt 10 --max_cmrs 2147483647 --max_equal 1 --max_insert_size 1000 --max_polya -1 --max_read_length 0 --min_identity 0.650000 --min_insert_size 0 --min_mq 0 --min_residues 0.500000 --min_score 0.000000 --mismatch_penalty 15 --mode 0 --no_progress 1 --no_unal 0 --ocl_threads 1 --output out_ngm.sam --overwrite 1 --pair_score_cutoff 0.900000 --paired 1 --parse_all 1 --pe_delimiter / --qry1 reads1.fastq --qry2 reads2.fastq --qry_count -1 --qry_start 0 --ref /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas --ref_mode -1 --sensitive 0 --silent_clip 0 --skip_mate_check 0 --skip_save 0 --slam_seq 0 --step_count 4 --strata 0 --topn 1 --trim5 0 --update_check 0 --very_fast 0 --very_sensitive 0
[NGM] Opening for output (SAM): out_ngm.sam
[SEQPROV] Reading encoded reference from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-enc.2.ngm
[SEQPROV] Reading 119 Mbp from disk took 0.06s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/Arabidopsis_thaliana.fas-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 0.68s
[PREPROCESS] Max. k-mer frequency set so 100!
[INPUT] Input is paired end data.
[INPUT] Opening file reads1.fastq for reading
[INPUT] Opening file reads2.fastq for reading
[INPUT] Input is Fastq
[INPUT] Estimating parameter from data
[INPUT] Reads found in files: 249307
[INPUT] Average read length: 60 (min: 60, max: 62)
[INPUT] Corridor width: 14
[INPUT] Average kmer hits pro read: 14.313253
[INPUT] Max possible kmer hit: 16
[INPUT] Estimated sensitivity: 0.894578
[INPUT] Estimating parameter took 1.462s
[INPUT] Input is Fastq
[INPUT] Input is Fastq
No protocol specified
[OPENCL] Available platforms: 1
[OPENCL] AMD Accelerated Parallel Processing
[OPENCL] Selecting OpenCl platform: AMD Accelerated Parallel Processing
[OPENCL] Platform: OpenCL 1.2 AMD-APP (1214.3)
[OPENCL] 1 CPU device found.
[OPENCL] Device 0: AMD Opteron(tm) Processor 6348 (Driver: 1214.3 (sse2,avx,fma4))
[OPENCL] 48 CPU cores available.
[MAIN] Valid pairs found: 99.74%
[MAIN] Estimated insert size: 103 bp
[MAIN] Alignments computed: 498570
[MAIN] Done (498560 reads mapped (99.99%), 54 reads not mapped, 498614 lines written)(elapsed: 12.909911s)
[UPDATE_CHECK] Your version of NGM is more than 6 months old - a newer version may be available. (For performing an automatic check use --update-check)
systime2.884
time23.2261199951
usrtime34.076
working_directory/project/teaser/genometeaser/tests_generated/340a758089e3620a129e06c436f2f4af