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 Mapped6310.421%
Wrongly Mapped480.032%
Not Mapped14931499.547%
Total149993100.0%
Read Failure Statistics
Not mapped14931499.547%
Missing in mapper output00.0%
Mapped to wrong chromosome00.0%
Mapped to wrong position480.032%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment7
Secondary Alignments0
F-Measure0.008400
Precision0.929300
Recall0.004200
Timing
Raw Mapping Time113.644s
Effective Mapping Time106.768s
Effective Init Time6.876s
Effective Time MeasureCPU
Mapping Time (Wall)60.311s
Mapping Time (CPU)113.644s
Mapping Time (CPU User)108.748s
Mapping Time (CPU System)4.896s
Init Time (Wall)7.122s
Init Time (CPU)6.876s
Init Time (CPU User)4.932s
Init Time (CPU System)1.944s
Additional Information
Mapper Memory Usage1819 MB
Total Test Runtime (Wall)141.444s
Mapper Command Line: 
/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/E_coli.fasta -t 4 --qry reads.fastq --no-progress
Errors and Warnings
No problems were encountered.
Subprocess Log

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

command/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/E_coli.fasta -t 4 --qry reads_base.fastq --no-progress
memory721752000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-08-30.13:40:47
[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/E_coli.fasta --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/E_coli.fasta-enc.2.ngm
[SEQPROV] Reading 5 Mbp from disk took 0.27s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/E_coli.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 34.70s
[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.033s
[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 (0 reads mapped (0.00%), 1 reads not mapped, 1 lines written)(elapsed: 44.709793s)
[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.272
time44.8535540104
usrtime4.616
working_directory/project/teaser/genometeaser/tests_generated/3f0952a92d640dea6a7f5d449db65923

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

command/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/E_coli.fasta -t 4 --qry reads_base.fastq --no-progress
memory721780000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-08-30.13:41:32
[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/E_coli.fasta --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/E_coli.fasta-enc.2.ngm
[SEQPROV] Reading 5 Mbp from disk took 0.02s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/E_coli.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 1.08s
[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 (0 reads mapped (0.00%), 1 reads not mapped, 1 lines written)(elapsed: 6.936268s)
[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)
systime1.944
time7.121560812
usrtime4.932
working_directory/project/teaser/genometeaser/tests_generated/3f0952a92d640dea6a7f5d449db65923

/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/E_coli.fasta -t 4 --qry reads.fastq --no-progress

command/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/E_coli.fasta -t 4 --qry reads.fastq --no-progress
memory1819524000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2018-08-30.13:41:40
[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.fastq --qry_count -1 --qry_start 0 --ref /project/teaser/genometeaser/references/E_coli.fasta --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/E_coli.fasta-enc.2.ngm
[SEQPROV] Reading 5 Mbp from disk took 0.01s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/E_coli.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 0.58s
[PREPROCESS] Max. k-mer frequency set so 100!
[INPUT] Input is single end data.
[INPUT] Opening file reads.fastq for reading
[INPUT] Input is Fastq
[INPUT] Estimating parameter from data
[INPUT] Reads found in files: 150000
[INPUT] Average read length: 172 (min: 100, max: 324)
[INPUT] Corridor width: 30
[INPUT] Average kmer hits pro read: 1.245136
[INPUT] Max possible kmer hit: 53
[INPUT] Estimated sensitivity: 0.300000
[INPUT] Estimating parameter took 0.979s
[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: 149898
[MAIN] Done (679 reads mapped (0.45%), 149321 reads not mapped, 150000 lines written)(elapsed: 48.766285s)
[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)
systime4.896
time60.3112120628
usrtime108.748
working_directory/project/teaser/genometeaser/tests_generated/3f0952a92d640dea6a7f5d449db65923