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 Mapped29950295.47%
Wrongly Mapped138924.428%
Not Mapped3180.101%
Total313712100.0%
Read Failure Statistics
Not mapped3180.101%
Missing in mapper output00.0%
Mapped to wrong chromosome107023.411%
Mapped to wrong position31901.017%
Mapped to wrong strand00.0%
Advanced Statistics
Missing in comparison alignment4
Secondary Alignments0
F-Measure0.976800
Precision0.955700
Recall0.998900
Timing
Raw Mapping Time22.279s
Effective Mapping Time19.775s
Effective Init Time2.504s
Effective Time MeasureWall clock
Mapping Time (Wall)22.279s
Mapping Time (CPU)62.610s
Mapping Time (CPU User)58.927s
Mapping Time (CPU System)3.683s
Init Time (Wall)2.504s
Init Time (CPU)2.456s
Init Time (CPU User)1.143s
Init Time (CPU System)1.313s
Additional Information
Mapper Memory Usage7076 MB
Total Test Runtime (Wall)96.893s
Mapper Command Line: 
/software/ngm/dev/ngm-0.5.2/ngm --qry2 reads2.fastq --output out_ngm-0.5.2.sam --qry1 reads1.fastq --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --no-progress
Errors and Warnings
No problems were encountered.
Subprocess Log

/software/ngm/dev/ngm-0.5.2/ngm --output out_ngm-0.5.2.sam --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --qry reads_base.fastq --no-progress

command/software/ngm/dev/ngm-0.5.2/ngm --output out_ngm-0.5.2.sam --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --qry reads_base.fastq --no-progress
memory6004604000
return0
status1
stderr

None
stdout

[MAIN] NextGenMap 0.5.2

[MAIN] Startup : x64 (build Nov 24 2016 17:40:01)

[MAIN] Starting time: 2017-02-06.14:52:38

[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-0.5.2.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 /project2/NextGenMapTest/Teaser/references/GRCh37.fa --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-0.5.2.sam

[SEQPROV] Reading encoded reference from /project2/NextGenMapTest/Teaser/references/GRCh37.fa-enc.2.ngm

[SEQPROV] Reading 3137 Mbp from disk took 13.84s

[PREPROCESS] Reading RefTable from /project2/NextGenMapTest/Teaser/references/GRCh37.fa-ht-13-2.3.ngm

[PREPROCESS] Reading from disk took 43.34s

[PREPROCESS] Max. k-mer frequency set so 879!

[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.002s

[INPUT] Input is Fastq

Setting of real/effective user Id to 0/0 failed

FATAL: Module fglrx not found.

Error! Fail to load fglrx kernel module! Maybe you can switch to root user to load kernel module directly

Setting of real/effective user Id to 0/0 failed

FATAL: Module fglrx not found.

Error! Fail to load fglrx kernel module! Maybe you can switch to root user to load kernel module directly

[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: Intel(R) Core(TM) i5-2500K CPU @ 3.30GHz (Driver: 1214.3 (sse2,avx))

[OPENCL] 4 CPU cores available.

[MAIN] Alignments computed: 1

[MAIN] Done (1 reads mapped (100.00%), 0 reads not mapped, 1 lines written)(elapsed: 59.464256s)
systime2.544613
time59.9092600346
usrtime1.218814
working_directory/project2/NextGenMapTest/Teaser/tests_generated/H1_n

/software/ngm/dev/ngm-0.5.2/ngm --output out_ngm-0.5.2.sam --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --qry reads_base.fastq --no-progress

command/software/ngm/dev/ngm-0.5.2/ngm --output out_ngm-0.5.2.sam --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --qry reads_base.fastq --no-progress
memory6002048000
return0
status1
stderr

None
stdout

[MAIN] NextGenMap 0.5.2

[MAIN] Startup : x64 (build Nov 24 2016 17:40:01)

[MAIN] Starting time: 2017-02-06.14:53:38

[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-0.5.2.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 /project2/NextGenMapTest/Teaser/references/GRCh37.fa --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-0.5.2.sam

[SEQPROV] Reading encoded reference from /project2/NextGenMapTest/Teaser/references/GRCh37.fa-enc.2.ngm

[SEQPROV] Reading 3137 Mbp from disk took 0.28s

[PREPROCESS] Reading RefTable from /project2/NextGenMapTest/Teaser/references/GRCh37.fa-ht-13-2.3.ngm

[PREPROCESS] Reading from disk took 0.90s

[PREPROCESS] Max. k-mer frequency set so 879!

[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.000s

[INPUT] Input is Fastq

Setting of real/effective user Id to 0/0 failed

FATAL: Module fglrx not found.

Error! Fail to load fglrx kernel module! Maybe you can switch to root user to load kernel module directly

Setting of real/effective user Id to 0/0 failed

FATAL: Module fglrx not found.

Error! Fail to load fglrx kernel module! Maybe you can switch to root user to load kernel module directly

[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: Intel(R) Core(TM) i5-2500K CPU @ 3.30GHz (Driver: 1214.3 (sse2,avx))

[OPENCL] 4 CPU cores available.

[MAIN] Alignments computed: 1

[MAIN] Done (1 reads mapped (100.00%), 0 reads not mapped, 1 lines written)(elapsed: 2.328681s)
systime1.3128
time2.50399899483
usrtime1.142826
working_directory/project2/NextGenMapTest/Teaser/tests_generated/H1_n

/software/ngm/dev/ngm-0.5.2/ngm --qry2 reads2.fastq --output out_ngm-0.5.2.sam --qry1 reads1.fastq --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --no-progress

command/software/ngm/dev/ngm-0.5.2/ngm --qry2 reads2.fastq --output out_ngm-0.5.2.sam --qry1 reads1.fastq --ref /project2/NextGenMapTest/Teaser/references/GRCh37.fa -t 4 --no-progress
memory7076064000
return0
status1
stderr

None
stdout

[MAIN] NextGenMap 0.5.2

[MAIN] Startup : x64 (build Nov 24 2016 17:40:01)

[MAIN] Starting time: 2017-02-06.14:53: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-0.5.2.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 /project2/NextGenMapTest/Teaser/references/GRCh37.fa --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-0.5.2.sam

[SEQPROV] Reading encoded reference from /project2/NextGenMapTest/Teaser/references/GRCh37.fa-enc.2.ngm

[SEQPROV] Reading 3137 Mbp from disk took 0.26s

[PREPROCESS] Reading RefTable from /project2/NextGenMapTest/Teaser/references/GRCh37.fa-ht-13-2.3.ngm

[PREPROCESS] Reading from disk took 0.84s

[PREPROCESS] Max. k-mer frequency set so 879!

[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: 156858

[INPUT] Average read length: 100 (min: 100, max: 102)

[INPUT] Corridor width: 20

[INPUT] Average kmer hits pro read: 19.633987

[INPUT] Max possible kmer hit: 29

[INPUT] Estimated sensitivity: 0.677034

[INPUT] Estimating parameter took 0.962s

[INPUT] Input is Fastq

[INPUT] Input is Fastq

Setting of real/effective user Id to 0/0 failed

FATAL: Module fglrx not found.

Error! Fail to load fglrx kernel module! Maybe you can switch to root user to load kernel module directly

Setting of real/effective user Id to 0/0 failed

FATAL: Module fglrx not found.

Error! Fail to load fglrx kernel module! Maybe you can switch to root user to load kernel module directly

[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: Intel(R) Core(TM) i5-2500K CPU @ 3.30GHz (Driver: 1214.3 (sse2,avx))

[OPENCL] 4 CPU cores available.

[MAIN] Valid pairs found: 95.55%

[MAIN] Estimated insert size: 500 bp

[MAIN] Alignments computed: 313501

[MAIN] Done (313395 reads mapped (99.90%), 321 reads not mapped, 313716 lines written)(elapsed: 20.883675s)
systime3.68344
time22.2793629169
usrtime58.927041
working_directory/project2/NextGenMapTest/Teaser/tests_generated/H1_n