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 Mapped70877787.894%
Wrongly Mapped8236010.213%
Not Mapped152671.893%
Total806404100.0%
Read Failure Statistics
Not mapped152671.893%
Missing in mapper output00.0%
Mapped to wrong chromosome501866.223%
Mapped to wrong position321743.99%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment0
Secondary Alignments0
F-Measure0.935600
Precision0.895900
Recall0.978900
Timing
Raw Mapping Time647.936s
Effective Mapping Time497.172s
Effective Init Time150.764s
Effective Time MeasureCPU
Mapping Time (Wall)287.737s
Mapping Time (CPU)647.936s
Mapping Time (CPU User)420.228s
Mapping Time (CPU System)227.708s
Init Time (Wall)134.559s
Init Time (CPU)150.764s
Init Time (CPU User)4.340s
Init Time (CPU System)146.424s
Additional Information
Mapper Memory Usage7945 MB
Total Test Runtime (Wall)1260.564s
Mapper Command Line: 
/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/M_musculus.fasta -t 4 --qry reads.fastq --no-progress
Errors and Warnings
No problems were encountered.
Subprocess Log

/project/teaser/genometeaser/software/ngm/ngm --ref /project/teaser/genometeaser/references/M_musculus.fasta

command/project/teaser/genometeaser/software/ngm/ngm --ref /project/teaser/genometeaser/references/M_musculus.fasta
return1
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2017-07-12.13:01:29
[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/M_musculus.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] Wrinting output (SAM) to stdout
[SEQPROV] Reading encoded reference from /project/teaser/genometeaser/references/M_musculus.fasta-enc.2.ngm
[SEQPROV] Reading 4838 Mbp from disk took 205.70s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/M_musculus.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 372.45s
[PREPROCESS] Max. k-mer frequency set so 1007!
[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.
systime0
time582.225491047
usrtime0
working_directory/project/teaser/genometeaser/tests_generated/a720be7d510db33409606fe9fee515f7

/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/M_musculus.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/M_musculus.fasta -t 4 --qry reads_base.fastq --no-progress
memory7402520000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2017-07-12.13:11:11
[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/M_musculus.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/M_musculus.fasta-enc.2.ngm
[SEQPROV] Reading 4838 Mbp from disk took 42.34s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/M_musculus.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 86.99s
[PREPROCESS] Max. k-mer frequency set so 1007!
[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: 147.433807s)
systime156.456
time149.153815031
usrtime4.512
working_directory/project/teaser/genometeaser/tests_generated/a720be7d510db33409606fe9fee515f7

/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/M_musculus.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/M_musculus.fasta -t 4 --qry reads_base.fastq --no-progress
memory7550504000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2017-07-12.13:13:41
[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/M_musculus.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/M_musculus.fasta-enc.2.ngm
[SEQPROV] Reading 4838 Mbp from disk took 39.88s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/M_musculus.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 77.76s
[PREPROCESS] Max. k-mer frequency set so 1007!
[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: 132.892319s)
systime146.424
time134.558995962
usrtime4.34
working_directory/project/teaser/genometeaser/tests_generated/a720be7d510db33409606fe9fee515f7

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

command/project/teaser/genometeaser/software/ngm/ngm --output out_ngm.sam --ref /project/teaser/genometeaser/references/M_musculus.fasta -t 4 --qry reads.fastq --no-progress
memory7945452000
return0
status1
stderr
None
stdout
[MAIN] NextGenMap 0.5.3
[MAIN] Startup : x64 (build Jan 31 2017 16:20:45)
[MAIN] Starting time: 2017-07-12.13:15:57
[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/M_musculus.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/M_musculus.fasta-enc.2.ngm
[SEQPROV] Reading 4838 Mbp from disk took 36.72s
[PREPROCESS] Reading RefTable from /project/teaser/genometeaser/references/M_musculus.fasta-ht-13-2.3.ngm
[PREPROCESS] Reading from disk took 88.98s
[PREPROCESS] Max. k-mer frequency set so 1007!
[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: 806404
[INPUT] Average read length: 60 (min: 60, max: 62)
[INPUT] Corridor width: 14
[INPUT] Average kmer hits pro read: 13.105459
[INPUT] Max possible kmer hit: 16
[INPUT] Estimated sensitivity: 0.819091
[INPUT] Estimating parameter took 6.899s
[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: 792113
[MAIN] Done (791137 reads mapped (98.11%), 15267 reads not mapped, 806404 lines written)(elapsed: 281.974945s)
systime227.708
time287.737012148
usrtime420.228
working_directory/project/teaser/genometeaser/tests_generated/a720be7d510db33409606fe9fee515f7