python /illumina/software/Survirus/SurVirus-master/surveyor.py /illumina/Projects/TEST_VIRIUS/Rawdata/L526401A/L526401A_R1.fastq.gz,/illumina/Projects/TEST_VIRIUS/Rawdata/L526401A/L526401A_R2.fastq.gz /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa --fq --threads 10 --dust /illumina/software/sdust/sdust-0.1/sdust [M::bwa_idx_load_from_disk] read 0 ALT contigs [M::process] read 200000 sequences (9997717 bp)... [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (115, 70565, 15, 75) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (40, 78, 431) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 1213) [M::mem_pestat] mean and std.dev: (166.03, 233.64) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1604) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (150, 239, 1155) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3165) [M::mem_pestat] mean and std.dev: (569.52, 703.84) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4170) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (73, 122, 474) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 1276) [M::mem_pestat] mean and std.dev: (140.75, 135.71) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1677) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 86, 502) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 1430) [M::mem_pestat] mean and std.dev: (202.61, 319.11) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1894) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 200000 reads in 71.489 CPU sec, 7.320 real sec [main] Version: 0.7.15-r1140 [main] CMD: bwa mem -t 10 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/head_1.fq /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/head_2.fq [main] Real time: 40.740 sec; CPU: 88.767 sec Executing: /illumina/software/Survirus/SurVirus-master/isolate_relevant_pairs_fq /illumina/Projects/TEST_VIRIUS/Rawdata/L526401A/L526401A_R1.fastq.gz /illumina/Projects/TEST_VIRIUS/Rawdata/L526401A/L526401A_R2.fastq.gz /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ Executing: bwa mem -t 10 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs_1.fq /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs_2.fq | samtools view -b -F 2304 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs.remapped.bam [M::bwa_idx_load_from_disk] read 0 ALT contigs [M::process] read 2000278 sequences (100000089 bp)... [M::process] read 2000282 sequences (100000047 bp)... [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1066, 705906, 150, 862) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 77, 253) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 687) [M::mem_pestat] mean and std.dev: (103.40, 128.12) [M::mem_pestat] low and high boundaries for proper pairs: (1, 904) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (150, 237, 1133) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3099) [M::mem_pestat] mean and std.dev: (558.17, 688.70) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4082) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (15, 141, 689) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2037) [M::mem_pestat] mean and std.dev: (290.20, 411.80) [M::mem_pestat] low and high boundaries for proper pairs: (1, 2711) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 86, 258) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 698) [M::mem_pestat] mean and std.dev: (104.35, 119.47) [M::mem_pestat] low and high boundaries for proper pairs: (1, 918) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 2000278 reads in 658.228 CPU sec, 67.109 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1080, 705580, 200, 927) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (33, 71, 190) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 504) [M::mem_pestat] mean and std.dev: (83.22, 87.69) [M::mem_pestat] low and high boundaries for proper pairs: (1, 661) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (149, 235, 1120) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3062) [M::mem_pestat] mean and std.dev: (551.60, 680.15) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4033) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (40, 243, 1279) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3757) [M::mem_pestat] mean and std.dev: (589.58, 770.25) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4996) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 83, 267) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 723) [M::mem_pestat] mean and std.dev: (103.86, 123.73) [M::mem_pestat] low and high boundaries for proper pairs: (1, 951) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000272 sequences (100000015 bp)... [M::mem_process_seqs] Processed 2000282 reads in 671.599 CPU sec, 67.794 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1056, 705565, 191, 962) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 78, 211) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 561) [M::mem_pestat] mean and std.dev: (87.83, 92.00) [M::mem_pestat] low and high boundaries for proper pairs: (1, 736) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (148, 233, 1104) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3016) [M::mem_pestat] mean and std.dev: (544.70, 670.73) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3972) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (32, 232, 1446) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 4274) [M::mem_pestat] mean and std.dev: (631.07, 837.68) [M::mem_pestat] low and high boundaries for proper pairs: (1, 5688) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 82, 271) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 739) [M::mem_pestat] mean and std.dev: (108.75, 133.74) [M::mem_pestat] low and high boundaries for proper pairs: (1, 973) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000035 bp)... [M::mem_process_seqs] Processed 2000272 reads in 670.427 CPU sec, 67.304 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1099, 705819, 214, 984) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 71, 202) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 536) [M::mem_pestat] mean and std.dev: (87.59, 92.84) [M::mem_pestat] low and high boundaries for proper pairs: (1, 703) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (148, 232, 1101) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3007) [M::mem_pestat] mean and std.dev: (541.86, 668.44) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3960) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (26, 225, 1027) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3029) [M::mem_pestat] mean and std.dev: (487.43, 652.20) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4030) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (41, 83, 293) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 797) [M::mem_pestat] mean and std.dev: (118.84, 151.45) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1049) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000284 sequences (100000074 bp)... [M::mem_process_seqs] Processed 2000294 reads in 676.997 CPU sec, 67.860 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1126, 706040, 190, 1027) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 69, 196) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 518) [M::mem_pestat] mean and std.dev: (82.19, 85.84) [M::mem_pestat] low and high boundaries for proper pairs: (1, 679) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1091) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2979) [M::mem_pestat] mean and std.dev: (537.24, 663.21) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3923) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (46, 289, 948) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2752) [M::mem_pestat] mean and std.dev: (496.82, 653.63) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3654) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 85, 301) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 829) [M::mem_pestat] mean and std.dev: (115.81, 144.44) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1093) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000280 sequences (100000045 bp)... [M::mem_process_seqs] Processed 2000284 reads in 676.004 CPU sec, 67.733 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1160, 705404, 190, 974) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 68, 202) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 538) [M::mem_pestat] mean and std.dev: (84.02, 91.40) [M::mem_pestat] low and high boundaries for proper pairs: (1, 706) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 231, 1090) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2976) [M::mem_pestat] mean and std.dev: (536.88, 661.55) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3919) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (20, 119, 999) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2957) [M::mem_pestat] mean and std.dev: (446.82, 650.18) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3936) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 83, 247) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 669) [M::mem_pestat] mean and std.dev: (102.79, 117.32) [M::mem_pestat] low and high boundaries for proper pairs: (1, 880) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000278 sequences (100000033 bp)... [M::mem_process_seqs] Processed 2000280 reads in 672.914 CPU sec, 67.370 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1147, 706101, 209, 1013) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 76, 235) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 635) [M::mem_pestat] mean and std.dev: (95.46, 112.24) [M::mem_pestat] low and high boundaries for proper pairs: (1, 835) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1087) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2967) [M::mem_pestat] mean and std.dev: (536.20, 661.14) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3907) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (22, 188, 767) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2257) [M::mem_pestat] mean and std.dev: (339.09, 454.33) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3002) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 72, 206) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 546) [M::mem_pestat] mean and std.dev: (87.74, 95.26) [M::mem_pestat] low and high boundaries for proper pairs: (1, 716) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000298 sequences (100000078 bp)... [M::mem_process_seqs] Processed 2000278 reads in 688.977 CPU sec, 69.205 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1150, 706166, 188, 987) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 70, 180) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 472) [M::mem_pestat] mean and std.dev: (81.93, 81.19) [M::mem_pestat] low and high boundaries for proper pairs: (1, 618) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 229, 1085) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2961) [M::mem_pestat] mean and std.dev: (534.86, 659.72) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3899) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (22, 119, 951) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2809) [M::mem_pestat] mean and std.dev: (393.27, 637.97) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3738) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 79, 217) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 577) [M::mem_pestat] mean and std.dev: (95.47, 103.36) [M::mem_pestat] low and high boundaries for proper pairs: (1, 757) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000284 sequences (100000070 bp)... [M::mem_process_seqs] Processed 2000298 reads in 703.781 CPU sec, 70.657 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1131, 705622, 171, 1028) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 73, 240) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 652) [M::mem_pestat] mean and std.dev: (95.57, 115.99) [M::mem_pestat] low and high boundaries for proper pairs: (1, 858) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 228, 1071) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2921) [M::mem_pestat] mean and std.dev: (528.96, 652.40) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3846) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (25, 138, 1012) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2986) [M::mem_pestat] mean and std.dev: (420.39, 583.90) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3973) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 84, 302) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 828) [M::mem_pestat] mean and std.dev: (113.67, 141.07) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1091) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000290 sequences (100000041 bp)... [M::mem_process_seqs] Processed 2000284 reads in 695.134 CPU sec, 69.781 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1160, 705740, 218, 1010) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 72, 209) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 557) [M::mem_pestat] mean and std.dev: (87.76, 94.93) [M::mem_pestat] low and high boundaries for proper pairs: (1, 731) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 226, 1055) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2873) [M::mem_pestat] mean and std.dev: (521.29, 642.03) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3782) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (24, 138, 935) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2757) [M::mem_pestat] mean and std.dev: (397.11, 582.30) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3668) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (35, 77, 205) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 545) [M::mem_pestat] mean and std.dev: (83.69, 84.10) [M::mem_pestat] low and high boundaries for proper pairs: (1, 715) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000288 sequences (100000076 bp)... [M::mem_process_seqs] Processed 2000290 reads in 684.862 CPU sec, 68.921 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1120, 704098, 198, 1062) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 72, 220) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 592) [M::mem_pestat] mean and std.dev: (89.04, 99.37) [M::mem_pestat] low and high boundaries for proper pairs: (1, 778) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 227, 1067) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2909) [M::mem_pestat] mean and std.dev: (527.20, 650.62) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3830) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (24, 154, 1000) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2952) [M::mem_pestat] mean and std.dev: (443.58, 643.34) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3928) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (35, 78, 226) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 608) [M::mem_pestat] mean and std.dev: (89.84, 102.25) [M::mem_pestat] low and high boundaries for proper pairs: (1, 799) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000304 sequences (100000043 bp)... [M::mem_process_seqs] Processed 2000288 reads in 704.614 CPU sec, 70.786 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1160, 705185, 196, 984) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 69, 193) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 509) [M::mem_pestat] mean and std.dev: (84.91, 87.04) [M::mem_pestat] low and high boundaries for proper pairs: (1, 667) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (145, 224, 1055) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2875) [M::mem_pestat] mean and std.dev: (522.05, 644.55) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3785) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (20, 235, 1078) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3194) [M::mem_pestat] mean and std.dev: (535.45, 666.16) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4252) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (35, 75, 311) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 863) [M::mem_pestat] mean and std.dev: (116.34, 158.92) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1139) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000314 sequences (100000022 bp)... [M::mem_process_seqs] Processed 2000304 reads in 726.556 CPU sec, 73.319 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1105, 706235, 187, 987) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 71, 195) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 515) [M::mem_pestat] mean and std.dev: (79.36, 78.48) [M::mem_pestat] low and high boundaries for proper pairs: (1, 675) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 226, 1064) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2900) [M::mem_pestat] mean and std.dev: (525.27, 647.81) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3818) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (22, 208, 854) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2518) [M::mem_pestat] mean and std.dev: (429.40, 551.07) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3350) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 87, 272) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 740) [M::mem_pestat] mean and std.dev: (108.36, 129.95) [M::mem_pestat] low and high boundaries for proper pairs: (1, 974) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000274 sequences (100000087 bp)... [M::mem_process_seqs] Processed 2000314 reads in 748.234 CPU sec, 75.286 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1112, 706081, 170, 891) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 79, 242) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 654) [M::mem_pestat] mean and std.dev: (96.24, 107.81) [M::mem_pestat] low and high boundaries for proper pairs: (1, 860) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (150, 237, 1133) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3099) [M::mem_pestat] mean and std.dev: (557.98, 688.77) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4082) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (22, 166, 1237) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3667) [M::mem_pestat] mean and std.dev: (581.48, 766.54) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4882) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (40, 93, 338) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 934) [M::mem_pestat] mean and std.dev: (135.11, 171.16) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1232) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000099 bp)... [M::mem_process_seqs] Processed 2000274 reads in 752.549 CPU sec, 75.763 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1128, 705348, 173, 872) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 69, 188) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 492) [M::mem_pestat] mean and std.dev: (82.64, 84.93) [M::mem_pestat] low and high boundaries for proper pairs: (1, 644) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (149, 235, 1115) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3047) [M::mem_pestat] mean and std.dev: (549.68, 676.86) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4013) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (19, 138, 1035) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3067) [M::mem_pestat] mean and std.dev: (519.58, 747.76) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4083) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 86, 271) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 737) [M::mem_pestat] mean and std.dev: (109.47, 123.84) [M::mem_pestat] low and high boundaries for proper pairs: (1, 970) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000288 sequences (100000036 bp)... [M::mem_process_seqs] Processed 2000294 reads in 746.862 CPU sec, 75.273 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1131, 707434, 176, 989) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (33, 67, 202) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 540) [M::mem_pestat] mean and std.dev: (84.71, 94.29) [M::mem_pestat] low and high boundaries for proper pairs: (1, 709) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (148, 232, 1098) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2998) [M::mem_pestat] mean and std.dev: (541.09, 667.12) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3948) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (20, 170, 1024) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3032) [M::mem_pestat] mean and std.dev: (481.48, 696.57) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4036) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 84, 249) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 671) [M::mem_pestat] mean and std.dev: (107.08, 124.88) [M::mem_pestat] low and high boundaries for proper pairs: (1, 882) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000036 bp)... [M::mem_process_seqs] Processed 2000288 reads in 759.433 CPU sec, 76.333 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1153, 706398, 198, 973) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 71, 188) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 492) [M::mem_pestat] mean and std.dev: (85.10, 86.16) [M::mem_pestat] low and high boundaries for proper pairs: (1, 644) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1076) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2934) [M::mem_pestat] mean and std.dev: (531.97, 655.10) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3863) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (16, 144, 922) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2734) [M::mem_pestat] mean and std.dev: (462.18, 628.35) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3640) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 82, 321) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 887) [M::mem_pestat] mean and std.dev: (118.54, 160.10) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1170) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000278 sequences (100000078 bp)... [M::mem_process_seqs] Processed 2000294 reads in 740.582 CPU sec, 74.306 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1145, 705890, 169, 963) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 69, 183) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 481) [M::mem_pestat] mean and std.dev: (78.35, 77.43) [M::mem_pestat] low and high boundaries for proper pairs: (1, 630) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1072) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2922) [M::mem_pestat] mean and std.dev: (530.10, 651.45) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3847) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (37, 208, 1343) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3955) [M::mem_pestat] mean and std.dev: (633.42, 851.67) [M::mem_pestat] low and high boundaries for proper pairs: (1, 5261) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 80, 245) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 663) [M::mem_pestat] mean and std.dev: (102.35, 121.29) [M::mem_pestat] low and high boundaries for proper pairs: (1, 872) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000296 sequences (100000009 bp)... [M::mem_process_seqs] Processed 2000278 reads in 733.495 CPU sec, 73.664 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1146, 706349, 161, 1044) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 72, 230) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 620) [M::mem_pestat] mean and std.dev: (94.95, 110.18) [M::mem_pestat] low and high boundaries for proper pairs: (1, 815) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1072) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2922) [M::mem_pestat] mean and std.dev: (529.62, 651.03) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3847) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (36, 281, 1372) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 4044) [M::mem_pestat] mean and std.dev: (741.15, 987.31) [M::mem_pestat] low and high boundaries for proper pairs: (1, 5380) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (40, 84, 256) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 688) [M::mem_pestat] mean and std.dev: (101.42, 116.73) [M::mem_pestat] low and high boundaries for proper pairs: (1, 904) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000288 sequences (100000050 bp)... [M::mem_process_seqs] Processed 2000296 reads in 734.066 CPU sec, 73.896 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1107, 705331, 190, 967) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 71, 198) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 524) [M::mem_pestat] mean and std.dev: (84.61, 85.56) [M::mem_pestat] low and high boundaries for proper pairs: (1, 687) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1081) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2949) [M::mem_pestat] mean and std.dev: (533.81, 657.55) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3883) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (22, 180, 1147) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3397) [M::mem_pestat] mean and std.dev: (446.77, 607.01) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4522) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (40, 85, 264) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 712) [M::mem_pestat] mean and std.dev: (113.60, 134.74) [M::mem_pestat] low and high boundaries for proper pairs: (1, 936) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000286 sequences (100000092 bp)... [M::mem_process_seqs] Processed 2000288 reads in 704.547 CPU sec, 70.857 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1216, 706711, 171, 1022) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 70, 196) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 516) [M::mem_pestat] mean and std.dev: (82.55, 84.40) [M::mem_pestat] low and high boundaries for proper pairs: (1, 676) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 228, 1063) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2897) [M::mem_pestat] mean and std.dev: (524.76, 645.85) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3814) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (44, 206, 968) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2816) [M::mem_pestat] mean and std.dev: (520.80, 676.11) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3740) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 84, 270) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 732) [M::mem_pestat] mean and std.dev: (106.31, 128.21) [M::mem_pestat] low and high boundaries for proper pairs: (1, 963) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000015 bp)... [M::mem_process_seqs] Processed 2000286 reads in 689.794 CPU sec, 69.261 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1189, 705933, 189, 1058) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (33, 69, 199) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 531) [M::mem_pestat] mean and std.dev: (81.52, 86.15) [M::mem_pestat] low and high boundaries for proper pairs: (1, 697) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 227, 1062) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2894) [M::mem_pestat] mean and std.dev: (524.79, 646.60) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3810) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (22, 117, 686) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2014) [M::mem_pestat] mean and std.dev: (309.36, 452.59) [M::mem_pestat] low and high boundaries for proper pairs: (1, 2678) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (35, 76, 274) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 752) [M::mem_pestat] mean and std.dev: (105.67, 135.65) [M::mem_pestat] low and high boundaries for proper pairs: (1, 991) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000302 sequences (100000075 bp)... [M::mem_process_seqs] Processed 2000294 reads in 687.702 CPU sec, 69.101 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1139, 705406, 186, 1047) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 67, 193) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 507) [M::mem_pestat] mean and std.dev: (77.05, 76.51) [M::mem_pestat] low and high boundaries for proper pairs: (1, 664) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 226, 1052) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2864) [M::mem_pestat] mean and std.dev: (520.20, 641.01) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3770) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (20, 119, 905) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2675) [M::mem_pestat] mean and std.dev: (416.17, 624.69) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3560) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 81, 225) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 601) [M::mem_pestat] mean and std.dev: (95.12, 107.22) [M::mem_pestat] low and high boundaries for proper pairs: (1, 789) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000012 bp)... [M::mem_process_seqs] Processed 2000302 reads in 699.261 CPU sec, 70.388 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1093, 705757, 189, 1029) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 74, 214) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 572) [M::mem_pestat] mean and std.dev: (88.39, 97.32) [M::mem_pestat] low and high boundaries for proper pairs: (1, 751) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (145, 225, 1052) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2866) [M::mem_pestat] mean and std.dev: (520.26, 641.66) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3773) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (16, 100, 932) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2764) [M::mem_pestat] mean and std.dev: (396.57, 604.27) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3680) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (33, 76, 186) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 492) [M::mem_pestat] mean and std.dev: (80.67, 78.40) [M::mem_pestat] low and high boundaries for proper pairs: (1, 645) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000296 sequences (100000062 bp)... [M::mem_process_seqs] Processed 2000294 reads in 680.773 CPU sec, 68.461 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1135, 705840, 191, 1044) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 66, 170) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 438) [M::mem_pestat] mean and std.dev: (75.52, 68.14) [M::mem_pestat] low and high boundaries for proper pairs: (1, 572) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (145, 223, 1038) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2824) [M::mem_pestat] mean and std.dev: (514.01, 634.05) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3717) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (28, 129, 1101) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3247) [M::mem_pestat] mean and std.dev: (561.34, 818.94) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4320) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 79, 235) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 633) [M::mem_pestat] mean and std.dev: (99.21, 114.40) [M::mem_pestat] low and high boundaries for proper pairs: (1, 832) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000308 sequences (100000083 bp)... [M::mem_process_seqs] Processed 2000296 reads in 679.216 CPU sec, 68.252 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1154, 705987, 183, 938) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 70, 207) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 551) [M::mem_pestat] mean and std.dev: (86.72, 92.81) [M::mem_pestat] low and high boundaries for proper pairs: (1, 723) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 227, 1064) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2900) [M::mem_pestat] mean and std.dev: (525.93, 648.74) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3818) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (31, 207, 1020) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2998) [M::mem_pestat] mean and std.dev: (464.63, 658.12) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3987) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 88, 252) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 680) [M::mem_pestat] mean and std.dev: (100.87, 111.00) [M::mem_pestat] low and high boundaries for proper pairs: (1, 894) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000262 sequences (100000041 bp)... [M::mem_process_seqs] Processed 2000308 reads in 693.219 CPU sec, 69.578 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1081, 706244, 171, 820) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (39, 80, 302) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 828) [M::mem_pestat] mean and std.dev: (118.43, 152.74) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1091) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (149, 236, 1124) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3074) [M::mem_pestat] mean and std.dev: (554.39, 684.73) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4049) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (35, 313, 1230) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3620) [M::mem_pestat] mean and std.dev: (654.47, 856.52) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4815) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 87, 271) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 735) [M::mem_pestat] mean and std.dev: (104.26, 123.90) [M::mem_pestat] low and high boundaries for proper pairs: (1, 967) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000270 sequences (100000082 bp)... [M::mem_process_seqs] Processed 2000262 reads in 711.813 CPU sec, 71.565 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1113, 706388, 156, 843) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (39, 82, 238) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 636) [M::mem_pestat] mean and std.dev: (95.69, 100.80) [M::mem_pestat] low and high boundaries for proper pairs: (1, 835) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::process] read 2000282 sequences (100000042 bp)... [M::mem_pestat] (25, 50, 75) percentile: (149, 234, 1104) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3014) [M::mem_pestat] mean and std.dev: (544.91, 671.06) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3969) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (29, 199, 1180) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3482) [M::mem_pestat] mean and std.dev: (493.66, 707.45) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4633) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 89, 255) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 689) [M::mem_pestat] mean and std.dev: (107.28, 121.48) [M::mem_pestat] low and high boundaries for proper pairs: (1, 906) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 2000270 reads in 699.860 CPU sec, 70.170 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1086, 706230, 190, 847) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 69, 219) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 585) [M::mem_pestat] mean and std.dev: (91.12, 104.82) [M::mem_pestat] low and high boundaries for proper pairs: (1, 768) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (148, 233, 1100) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3004) [M::mem_pestat] mean and std.dev: (542.22, 668.05) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3956) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (32, 262, 1084) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3188) [M::mem_pestat] mean and std.dev: (575.52, 772.55) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4240) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (43, 91, 259) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 691) [M::mem_pestat] mean and std.dev: (111.72, 125.81) [M::mem_pestat] low and high boundaries for proper pairs: (1, 907) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000302 sequences (100000073 bp)... [M::mem_process_seqs] Processed 2000282 reads in 692.819 CPU sec, 69.578 real sec [M::process] read 2000278 sequences (100000050 bp)... [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1131, 706371, 179, 917) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 67, 190) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 500) [M::mem_pestat] mean and std.dev: (82.43, 84.80) [M::mem_pestat] low and high boundaries for proper pairs: (1, 655) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1076) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2934) [M::mem_pestat] mean and std.dev: (531.60, 653.49) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3863) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (28, 221, 928) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2728) [M::mem_pestat] mean and std.dev: (406.24, 508.37) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3628) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (35, 79, 232) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 626) [M::mem_pestat] mean and std.dev: (97.55, 110.48) [M::mem_pestat] low and high boundaries for proper pairs: (1, 823) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 2000302 reads in 691.797 CPU sec, 69.239 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1134, 707133, 192, 910) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 67, 174) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 450) [M::mem_pestat] mean and std.dev: (77.14, 72.25) [M::mem_pestat] low and high boundaries for proper pairs: (1, 588) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 231, 1085) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2961) [M::mem_pestat] mean and std.dev: (536.39, 660.62) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3899) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (46, 193, 995) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2893) [M::mem_pestat] mean and std.dev: (495.73, 664.42) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3842) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 86, 248) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 666) [M::mem_pestat] mean and std.dev: (100.51, 109.19) [M::mem_pestat] low and high boundaries for proper pairs: (1, 875) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000278 sequences (100000047 bp)... [M::mem_process_seqs] Processed 2000278 reads in 684.415 CPU sec, 68.467 real sec [M::process] read 2000272 sequences (100000089 bp)... [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1089, 705496, 185, 921) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (37, 73, 207) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 547) [M::mem_pestat] mean and std.dev: (87.99, 94.63) [M::mem_pestat] low and high boundaries for proper pairs: (1, 717) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1086) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2964) [M::mem_pestat] mean and std.dev: (536.49, 661.61) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3903) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (33, 226, 1045) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3069) [M::mem_pestat] mean and std.dev: (545.08, 712.08) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4081) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (41, 90, 298) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 812) [M::mem_pestat] mean and std.dev: (121.96, 151.82) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1069) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 2000278 reads in 685.841 CPU sec, 68.533 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1082, 706647, 189, 896) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (37, 72, 200) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 526) [M::mem_pestat] mean and std.dev: (84.03, 85.44) [M::mem_pestat] low and high boundaries for proper pairs: (1, 689) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 229, 1084) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2958) [M::mem_pestat] mean and std.dev: (534.01, 659.19) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3895) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (18, 124, 923) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2733) [M::mem_pestat] mean and std.dev: (444.53, 642.38) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3638) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 79, 259) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 701) [M::mem_pestat] mean and std.dev: (101.41, 121.91) [M::mem_pestat] low and high boundaries for proper pairs: (1, 922) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000082 bp)... [M::mem_process_seqs] Processed 2000272 reads in 682.708 CPU sec, 68.282 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1157, 706493, 183, 946) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 66, 176) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 460) [M::mem_pestat] mean and std.dev: (77.81, 76.68) [M::mem_pestat] low and high boundaries for proper pairs: (1, 602) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 228, 1082) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2952) [M::mem_pestat] mean and std.dev: (532.60, 657.92) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3887) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (25, 249, 1151) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3403) [M::mem_pestat] mean and std.dev: (554.77, 746.90) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4529) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (43, 92, 365) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 1009) [M::mem_pestat] mean and std.dev: (130.19, 180.36) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1331) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000296 sequences (100000038 bp)... [M::mem_process_seqs] Processed 2000294 reads in 687.201 CPU sec, 68.704 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1125, 707478, 199, 863) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 69, 175) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 457) [M::mem_pestat] mean and std.dev: (80.79, 77.99) [M::mem_pestat] low and high boundaries for proper pairs: (1, 598) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 226, 1051) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2861) [M::mem_pestat] mean and std.dev: (520.38, 640.13) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3766) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (18, 119, 1204) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3576) [M::mem_pestat] mean and std.dev: (590.75, 872.57) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4762) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 91, 377) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 1053) [M::mem_pestat] mean and std.dev: (142.44, 198.15) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1391) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000292 sequences (100000005 bp)... [M::mem_process_seqs] Processed 2000296 reads in 673.621 CPU sec, 67.769 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1103, 707871, 160, 878) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (37, 72, 209) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 553) [M::mem_pestat] mean and std.dev: (83.62, 88.14) [M::mem_pestat] low and high boundaries for proper pairs: (1, 725) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 226, 1063) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2897) [M::mem_pestat] mean and std.dev: (524.50, 647.83) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3814) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (28, 304, 1003) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2953) [M::mem_pestat] mean and std.dev: (528.61, 708.10) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3928) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 81, 284) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 778) [M::mem_pestat] mean and std.dev: (111.44, 138.56) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1025) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000294 sequences (100000005 bp)... [M::mem_process_seqs] Processed 2000292 reads in 677.984 CPU sec, 67.865 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1144, 706952, 169, 888) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 70, 197) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 519) [M::mem_pestat] mean and std.dev: (82.72, 83.06) [M::mem_pestat] low and high boundaries for proper pairs: (1, 680) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (145, 224, 1044) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2842) [M::mem_pestat] mean and std.dev: (515.94, 637.20) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3741) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (18, 181, 1015) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3009) [M::mem_pestat] mean and std.dev: (521.20, 735.82) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4006) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 82, 245) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 663) [M::mem_pestat] mean and std.dev: (101.04, 119.58) [M::mem_pestat] low and high boundaries for proper pairs: (1, 872) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000276 sequences (100000093 bp)... [M::mem_process_seqs] Processed 2000294 reads in 672.668 CPU sec, 67.276 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1071, 707148, 184, 880) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 71, 209) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 559) [M::mem_pestat] mean and std.dev: (91.49, 101.12) [M::mem_pestat] low and high boundaries for proper pairs: (1, 734) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 225, 1053) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2867) [M::mem_pestat] mean and std.dev: (520.65, 641.62) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3774) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (31, 190, 1136) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3346) [M::mem_pestat] mean and std.dev: (576.01, 802.77) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4451) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 86, 321) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 885) [M::mem_pestat] mean and std.dev: (124.64, 170.71) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1167) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000262 sequences (100000053 bp)... [M::mem_process_seqs] Processed 2000276 reads in 669.031 CPU sec, 67.040 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1085, 706221, 190, 811) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 76, 242) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 654) [M::mem_pestat] mean and std.dev: (99.99, 116.70) [M::mem_pestat] low and high boundaries for proper pairs: (1, 860) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (149, 235, 1117) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3053) [M::mem_pestat] mean and std.dev: (551.03, 679.54) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4021) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (25, 202, 1137) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3361) [M::mem_pestat] mean and std.dev: (543.82, 741.27) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4473) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (44, 96, 448) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 1256) [M::mem_pestat] mean and std.dev: (169.72, 235.13) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1660) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000288 sequences (100000062 bp)... [M::mem_process_seqs] Processed 2000262 reads in 695.257 CPU sec, 69.821 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1091, 705733, 175, 846) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 73, 227) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 611) [M::mem_pestat] mean and std.dev: (92.79, 107.93) [M::mem_pestat] low and high boundaries for proper pairs: (1, 803) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (148, 234, 1113) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3043) [M::mem_pestat] mean and std.dev: (548.68, 677.11) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4008) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (30, 132, 923) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2709) [M::mem_pestat] mean and std.dev: (419.15, 597.26) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3602) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (42, 94, 329) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 903) [M::mem_pestat] mean and std.dev: (134.65, 172.84) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1190) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000280 sequences (100000054 bp)... [M::mem_process_seqs] Processed 2000288 reads in 690.778 CPU sec, 69.225 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1195, 706664, 196, 854) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (37, 75, 204) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 538) [M::mem_pestat] mean and std.dev: (88.72, 93.86) [M::mem_pestat] low and high boundaries for proper pairs: (1, 705) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (148, 232, 1096) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2992) [M::mem_pestat] mean and std.dev: (539.87, 665.33) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3940) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (31, 120, 804) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2350) [M::mem_pestat] mean and std.dev: (412.79, 591.51) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3123) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 80, 258) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 696) [M::mem_pestat] mean and std.dev: (102.71, 122.48) [M::mem_pestat] low and high boundaries for proper pairs: (1, 915) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000290 sequences (100000090 bp)... [M::mem_process_seqs] Processed 2000280 reads in 686.197 CPU sec, 69.043 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1095, 705958, 187, 906) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 73, 194) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 514) [M::mem_pestat] mean and std.dev: (79.39, 78.22) [M::mem_pestat] low and high boundaries for proper pairs: (1, 674) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 229, 1067) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2907) [M::mem_pestat] mean and std.dev: (527.44, 649.99) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3827) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (25, 303, 1157) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3421) [M::mem_pestat] mean and std.dev: (560.46, 680.50) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4553) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 79, 247) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 663) [M::mem_pestat] mean and std.dev: (96.43, 112.62) [M::mem_pestat] low and high boundaries for proper pairs: (1, 871) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000288 sequences (100000055 bp)... [M::mem_process_seqs] Processed 2000290 reads in 672.927 CPU sec, 67.579 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1130, 706711, 200, 857) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (36, 78, 234) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 630) [M::mem_pestat] mean and std.dev: (88.35, 96.54) [M::mem_pestat] low and high boundaries for proper pairs: (1, 828) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 230, 1086) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2964) [M::mem_pestat] mean and std.dev: (534.88, 659.93) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3903) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (18, 236, 1243) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3693) [M::mem_pestat] mean and std.dev: (544.52, 803.76) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4918) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 79, 275) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 751) [M::mem_pestat] mean and std.dev: (110.70, 143.37) [M::mem_pestat] low and high boundaries for proper pairs: (1, 989) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000278 sequences (100000064 bp)... [M::mem_process_seqs] Processed 2000288 reads in 680.599 CPU sec, 68.199 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1122, 706618, 176, 906) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (39, 72, 196) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 510) [M::mem_pestat] mean and std.dev: (81.30, 78.33) [M::mem_pestat] low and high boundaries for proper pairs: (1, 667) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 229, 1076) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2934) [M::mem_pestat] mean and std.dev: (530.61, 653.45) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3863) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (26, 176, 953) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2807) [M::mem_pestat] mean and std.dev: (487.75, 673.19) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3734) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 84, 229) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 611) [M::mem_pestat] mean and std.dev: (94.53, 96.11) [M::mem_pestat] low and high boundaries for proper pairs: (1, 802) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000286 sequences (100000059 bp)... [M::mem_process_seqs] Processed 2000278 reads in 675.478 CPU sec, 67.442 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1145, 706211, 175, 889) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 70, 198) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 526) [M::mem_pestat] mean and std.dev: (83.96, 83.78) [M::mem_pestat] low and high boundaries for proper pairs: (1, 690) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (147, 228, 1070) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2916) [M::mem_pestat] mean and std.dev: (528.23, 650.54) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3839) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (16, 95, 1057) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3139) [M::mem_pestat] mean and std.dev: (461.09, 732.56) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4180) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 80, 305) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 841) [M::mem_pestat] mean and std.dev: (118.35, 160.60) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1109) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000284 sequences (100000071 bp)... [M::mem_process_seqs] Processed 2000286 reads in 677.724 CPU sec, 67.863 real sec [M::process] read 2000286 sequences (100000005 bp)... [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1087, 705790, 182, 863) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (38, 76, 226) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 602) [M::mem_pestat] mean and std.dev: (90.99, 99.62) [M::mem_pestat] low and high boundaries for proper pairs: (1, 790) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 228, 1068) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2912) [M::mem_pestat] mean and std.dev: (527.83, 650.59) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3834) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (24, 316, 1461) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 4335) [M::mem_pestat] mean and std.dev: (853.20, 1154.94) [M::mem_pestat] low and high boundaries for proper pairs: (1, 5772) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (38, 81, 280) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 764) [M::mem_pestat] mean and std.dev: (114.67, 146.79) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1006) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 2000284 reads in 678.809 CPU sec, 68.007 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1077, 706864, 161, 961) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (34, 73, 208) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 556) [M::mem_pestat] mean and std.dev: (86.76, 90.42) [M::mem_pestat] low and high boundaries for proper pairs: (1, 730) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 226, 1060) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2888) [M::mem_pestat] mean and std.dev: (524.18, 646.76) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3802) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (24, 232, 1201) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3555) [M::mem_pestat] mean and std.dev: (606.57, 802.06) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4732) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 82, 249) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 675) [M::mem_pestat] mean and std.dev: (96.16, 107.70) [M::mem_pestat] low and high boundaries for proper pairs: (1, 888) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000302 sequences (100000091 bp)... [M::mem_process_seqs] Processed 2000286 reads in 670.494 CPU sec, 67.328 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1141, 706306, 180, 874) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (38, 79, 255) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 689) [M::mem_pestat] mean and std.dev: (103.76, 122.92) [M::mem_pestat] low and high boundaries for proper pairs: (1, 906) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (146, 225, 1059) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2885) [M::mem_pestat] mean and std.dev: (522.17, 645.50) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3798) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (17, 145, 1204) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3578) [M::mem_pestat] mean and std.dev: (558.20, 834.92) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4765) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (39, 85, 288) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 786) [M::mem_pestat] mean and std.dev: (111.89, 145.85) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1035) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 2000296 sequences (100000096 bp)... [M::mem_process_seqs] Processed 2000302 reads in 671.457 CPU sec, 67.429 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (1125, 707177, 168, 930) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (35, 74, 198) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 524) [M::mem_pestat] mean and std.dev: (83.34, 82.34) [M::mem_pestat] low and high boundaries for proper pairs: (1, 687) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (145, 224, 1050) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2860) [M::mem_pestat] mean and std.dev: (519.47, 642.01) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3765) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (17, 163, 1084) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3218) [M::mem_pestat] mean and std.dev: (541.78, 729.47) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4285) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (36, 81, 206) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 546) [M::mem_pestat] mean and std.dev: (88.43, 89.28) [M::mem_pestat] low and high boundaries for proper pairs: (1, 716) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::process] read 1496866 sequences (74832153 bp)... [M::mem_process_seqs] Processed 2000296 reads in 676.911 CPU sec, 67.989 real sec [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (795, 529531, 127, 647) [M::mem_pestat] analyzing insert size distribution for orientation FF... [M::mem_pestat] (25, 50, 75) percentile: (37, 75, 197) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 517) [M::mem_pestat] mean and std.dev: (83.22, 86.30) [M::mem_pestat] low and high boundaries for proper pairs: (1, 677) [M::mem_pestat] analyzing insert size distribution for orientation FR... [M::mem_pestat] (25, 50, 75) percentile: (145, 222, 1040) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 2830) [M::mem_pestat] mean and std.dev: (513.73, 634.32) [M::mem_pestat] low and high boundaries for proper pairs: (1, 3725) [M::mem_pestat] analyzing insert size distribution for orientation RF... [M::mem_pestat] (25, 50, 75) percentile: (20, 387, 1057) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 3131) [M::mem_pestat] mean and std.dev: (594.73, 716.63) [M::mem_pestat] low and high boundaries for proper pairs: (1, 4168) [M::mem_pestat] analyzing insert size distribution for orientation RR... [M::mem_pestat] (25, 50, 75) percentile: (37, 86, 285) [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 781) [M::mem_pestat] mean and std.dev: (116.77, 144.91) [M::mem_pestat] low and high boundaries for proper pairs: (1, 1029) [M::mem_pestat] skip orientation FF [M::mem_pestat] skip orientation RF [M::mem_pestat] skip orientation RR [M::mem_process_seqs] Processed 1496866 reads in 481.664 CPU sec, 48.280 real sec [main] Version: 0.7.15-r1140 [main] CMD: bwa mem -t 10 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs_1.fq /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs_2.fq [main] Real time: 3496.795 sec; CPU: 34486.047 sec Executing: samtools sort -@ 10 /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs.remapped.bam -o /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//retained-pairs.remapped.cs.bam [bam_sort_core] merging from 20 files and 10 in-memory blocks... Executing: /illumina/software/Survirus/SurVirus-master/extract_clips /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ Executing: bwa aln -t 10 /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.fa -f /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.sai [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 [main] Version: 0.7.15-r1140 [main] CMD: bwa aln -t 10 -f /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.sai /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.fa [main] Real time: 12.409 sec; CPU: 2.155 sec Executing: bwa samse /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.sai /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.fa | samtools view -b -F 2304 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.full.bam [main] Version: 0.7.15-r1140 [main] CMD: bwa samse /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.sai /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.fa [main] Real time: 0.031 sec; CPU: 0.003 sec Executing: samtools view -b -F 4 /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.full.bam > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.aln.bam Executing: samtools fasta -f 4 /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.full.bam > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.unmapped.fa [M::bam2fq_mainloop] discarded 0 singletons [M::bam2fq_mainloop] processed 0 reads Executing: bwa mem -t 10 /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.unmapped.fa | samtools view -b -F 2308 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.mem.bam [M::bwa_idx_load_from_disk] read 0 ALT contigs [main] Version: 0.7.15-r1140 [main] CMD: bwa mem -t 10 /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.unmapped.fa [main] Real time: 12.001 sec; CPU: 4.321 sec Executing: samtools cat /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.aln.bam /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.mem.bam -o /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-clips.bam Executing: bwa aln -t 10 /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.fa -f /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.sai [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... 56.28 sec [bwa_aln_core] write to the disk... 0.02 sec [bwa_aln_core] 262144 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 47.94 sec [bwa_aln_core] write to the disk... 0.02 sec [bwa_aln_core] 524288 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 56.48 sec [bwa_aln_core] write to the disk... 0.01 sec [bwa_aln_core] 786432 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 48.88 sec [bwa_aln_core] write to the disk... 0.03 sec [bwa_aln_core] 1048576 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 51.53 sec [bwa_aln_core] write to the disk... 0.03 sec [bwa_aln_core] 1310720 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 50.40 sec [bwa_aln_core] write to the disk... 0.02 sec [bwa_aln_core] 1572864 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 45.68 sec [bwa_aln_core] write to the disk... 0.03 sec [bwa_aln_core] 1835008 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 49.25 sec [bwa_aln_core] write to the disk... 0.03 sec [bwa_aln_core] 2097152 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 48.59 sec [bwa_aln_core] write to the disk... 0.01 sec [bwa_aln_core] 2359296 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 45.01 sec [bwa_aln_core] write to the disk... 0.02 sec [bwa_aln_core] 2621440 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 48.72 sec [bwa_aln_core] write to the disk... 0.01 sec [bwa_aln_core] 2883584 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 51.61 sec [bwa_aln_core] write to the disk... 0.03 sec [bwa_aln_core] 3145728 sequences have been processed. [bwa_aln_core] calculate SA coordinate... 10.90 sec [bwa_aln_core] write to the disk... 0.00 sec [bwa_aln_core] 3206365 sequences have been processed. [main] Version: 0.7.15-r1140 [main] CMD: bwa aln -t 10 -f /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.sai /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.fa [main] Real time: 75.153 sec; CPU: 614.150 sec Executing: bwa samse /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.sai /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.fa | samtools view -b -F 2304 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.full.bam [bwa_aln_core] convert to sequence coordinate... 0.05 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.16 sec [bwa_aln_core] 262144 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.04 sec [bwa_aln_core] refine gapped alignments... 0.05 sec [bwa_aln_core] print alignments... 0.16 sec [bwa_aln_core] 524288 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.04 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.16 sec [bwa_aln_core] 786432 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.05 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.16 sec [bwa_aln_core] 1048576 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.06 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.18 sec [bwa_aln_core] 1310720 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.06 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.16 sec [bwa_aln_core] 1572864 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.06 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.16 sec [bwa_aln_core] 1835008 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.04 sec [bwa_aln_core] refine gapped alignments... 0.03 sec [bwa_aln_core] print alignments... 0.14 sec [bwa_aln_core] 2097152 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.08 sec [bwa_aln_core] refine gapped alignments... 0.06 sec [bwa_aln_core] print alignments... 0.15 sec [bwa_aln_core] 2359296 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.04 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.17 sec [bwa_aln_core] 2621440 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.04 sec [bwa_aln_core] refine gapped alignments... 0.04 sec [bwa_aln_core] print alignments... 0.17 sec [bwa_aln_core] 2883584 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.05 sec [bwa_aln_core] refine gapped alignments... 0.03 sec [bwa_aln_core] print alignments... 0.15 sec [bwa_aln_core] 3145728 sequences have been processed. [bwa_aln_core] convert to sequence coordinate... 0.01 sec [bwa_aln_core] refine gapped alignments... 0.01 sec [bwa_aln_core] print alignments... 0.03 sec [bwa_aln_core] 3206365 sequences have been processed. [main] Version: 0.7.15-r1140 [main] CMD: bwa samse /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.sai /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.fa [main] Real time: 15.165 sec; CPU: 5.299 sec Executing: samtools view -b -F 4 /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.full.bam > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.aln.bam Executing: samtools fasta -f 4 /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.full.bam > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.unmapped.fa [M::bam2fq_mainloop] discarded 0 singletons [M::bam2fq_mainloop] processed 3127162 reads Executing: bwa mem -t 10 /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.unmapped.fa | samtools view -b -F 2308 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.mem.bam [M::bwa_idx_load_from_disk] read 0 ALT contigs [M::process] read 3127162 sequences (71794164 bp)... [M::mem_process_seqs] Processed 3127162 reads in 71.249 CPU sec, 7.762 real sec [main] Version: 0.7.15-r1140 [main] CMD: bwa mem -t 10 /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.unmapped.fa [main] Real time: 10.794 sec; CPU: 73.812 sec Executing: samtools cat /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.aln.bam /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.mem.bam -o /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-clips.bam Executing: /illumina/software/Survirus/SurVirus-master/reads_categorizer /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ Executing: bwa mem -t 10 -h 762 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-side.fq | samtools view -b -F 2308 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-side.bam [M::bwa_idx_load_from_disk] read 0 ALT contigs [main] Version: 0.7.15-r1140 [main] CMD: bwa mem -t 10 -h 762 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//virus-side.fq [main] Real time: 3.990 sec; CPU: 3.985 sec Executing: bwa mem -t 10 -h 100 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-side.fq | samtools view -b -F 2308 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-side.bam [M::bwa_idx_load_from_disk] read 0 ALT contigs [M::process] read 78720 sequences (3935282 bp)... [M::mem_process_seqs] Processed 78720 reads in 8.115 CPU sec, 0.818 real sec [main] Version: 0.7.15-r1140 [main] CMD: bwa mem -t 10 -h 100 /illumina/genomes/Virus_Human/BWAindex/hg38_viruses_757.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0//host-side.fq [main] Real time: 5.418 sec; CPU: 11.862 sec Executing: /illumina/software/Survirus/SurVirus-master/merge_retained_reads /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ Executing: /illumina/software/Survirus/SurVirus-master/build_region-reads_associations /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ Executing: /illumina/software/Survirus/SurVirus-master/remapper /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.txt 2> /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/log.txt sh: line 1: 6222 Segmentation fault /illumina/software/Survirus/SurVirus-master/remapper /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.txt 2> /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/log.txt Executing: /illumina/software/Survirus/SurVirus-master/bp_region_consensus_builder /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/genomes/Virus/BWAindex/viruses_757.fasta /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/bam_0/ Executing: /illumina/software/sdust/sdust-0.1/sdust /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/host_bp_seqs.fa > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/host_bp_seqs.masked.bed Executing: /illumina/software/sdust/sdust-0.1/sdust /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/virus_bp_seqs.fa > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/virus_bp_seqs.masked.bed Executing: /illumina/software/Survirus/SurVirus-master/filter /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.t1.txt sh: line 1: 6235 Segmentation fault /illumina/software/Survirus/SurVirus-master/filter /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.t1.txt Executing: /illumina/software/Survirus/SurVirus-master/filter /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus --remapped > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.remapped.t1.txt sh: line 1: 6237 Segmentation fault /illumina/software/Survirus/SurVirus-master/filter /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus --remapped > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.remapped.t1.txt Executing: /illumina/software/Survirus/SurVirus-master/filter /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus --print-rejected > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.discarded.txt sh: line 1: 6239 Segmentation fault /illumina/software/Survirus/SurVirus-master/filter /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus --print-rejected > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/results.discarded.txt Finding alternative locations... Executing: bwa mem -t 10 -h 1000 /illumina/genomes/Human_hg38/BWAindex/hg38.fa /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/host_bp_seqs.fa | samtools view -b -F 2308 > /illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/host_bp_seqs.bam [M::bwa_idx_load_from_disk] read 0 ALT contigs [E::main_mem] fail to open file `/illumina/Projects/TEST_VIRIUS/Workdir/L526401A/Survirus/host_bp_seqs.fa'. Traceback (most recent call last): File "/illumina/software/Survirus/SurVirus-master/surveyor.py", line 230, in with pysam.AlignmentFile("%s/host_bp_seqs.bam" % cmd_args.workdir) as bp_seqs_bam, \ File "pysam/libcalignmentfile.pyx", line 734, in pysam.libcalignmentfile.AlignmentFile.__cinit__ File "pysam/libcalignmentfile.pyx", line 983, in pysam.libcalignmentfile.AlignmentFile._open ValueError: file has no sequences defined (mode='r') - is it SAM/BAM format? Consider opening with check_sq=False