@@ -180,7 +180,7 @@ def test_bam_compare_ZoverZ():
180180 _foo = open (outfile , 'r' )
181181 resp = _foo .readlines ()
182182 _foo .close ()
183- expected = ['3R\t 50\t 100\t -1\n ' , '3R\t 100\t 150\t 0\n ' , '3R\t 150\t 200\t -0.5849624 \n ' ]
183+ expected = ['3R\t 50\t 100\t -1\n ' , '3R\t 100\t 150\t 0\n ' , '3R\t 150\t 200\t -0.58 \n ' ]
184184 assert f"{ resp } " == f"{ expected } " , f"{ resp } != { expected } "
185185 unlink (outfile )
186186
@@ -293,7 +293,7 @@ def test_bam_compare_scale_factors_ratio():
293293 (scale factors [1,0.5]) (1+1)/(1+1*0.5)=1.33
294294 """
295295
296- expected = ['3R\t 0\t 50\t 1\n ' , '3R\t 50\t 100\t 0.6666667 \n ' , '3R\t 100\t 150\t 1.3333334 \n ' , '3R\t 150\t 200\t 1\n ' ]
296+ expected = ['3R\t 0\t 50\t 1\n ' , '3R\t 50\t 100\t 0.67 \n ' , '3R\t 100\t 150\t 1.33 \n ' , '3R\t 150\t 200\t 1\n ' ]
297297 assert f"{ resp } " == f"{ expected } " , f"{ resp } != { expected } "
298298 unlink (outfile )
299299
@@ -360,22 +360,22 @@ def test_bam_coverage_filter_blacklist():
360360
361361 expected = [
362362 "3R\t 0\t 100\t 0\n " ,
363- "3R\t 100\t 150\t 1.7521904 \n " ,
364- "3R\t 150\t 200\t 5.8406343 \n " ,
365- "3R\t 200\t 250\t 6.035322 \n " ,
366- "3R\t 250\t 300\t 3.6990685 \n " ,
367- "3R\t 300\t 400\t 2.5309415 \n " ,
368- "3R\t 400\t 450\t 4.283132 \n " ,
369- "3R\t 450\t 500\t 4.47782 \n " ,
370- "3R\t 500\t 550\t 2.141566 \n " ,
371- "3R\t 550\t 600\t 2.9203172 \n " ,
372- "3R\t 600\t 650\t 5.4512587 \n " ,
373- "3R\t 650\t 750\t 3.8937564 \n " ,
374- "3R\t 750\t 800\t 2.7256293 \n " ,
363+ "3R\t 100\t 150\t 1.75 \n " ,
364+ "3R\t 150\t 200\t 5.84 \n " ,
365+ "3R\t 200\t 250\t 6.04 \n " ,
366+ "3R\t 250\t 300\t 3.7 \n " ,
367+ "3R\t 300\t 400\t 2.53 \n " ,
368+ "3R\t 400\t 450\t 4.28 \n " ,
369+ "3R\t 450\t 500\t 4.48 \n " ,
370+ "3R\t 500\t 550\t 2.14 \n " ,
371+ "3R\t 550\t 600\t 2.92 \n " ,
372+ "3R\t 600\t 650\t 5.45 \n " ,
373+ "3R\t 650\t 750\t 3.89 \n " ,
374+ "3R\t 750\t 800\t 2.73 \n " ,
375375 "3R\t 800\t 900\t 0\n " ,
376- "3R\t 900\t 950\t 0.5840634 \n " ,
377- "3R\t 950\t 1000\t 1.3628147 \n " ,
378- "3R\t 1000\t 1050\t 0.77875125 \n " ,
376+ "3R\t 900\t 950\t 0.58 \n " ,
377+ "3R\t 950\t 1000\t 1.36 \n " ,
378+ "3R\t 1000\t 1050\t 0.78 \n " ,
379379 "3R\t 1050\t 1500\t 0\n "
380380 ]
381381
@@ -517,24 +517,23 @@ def test_bam_compare_filter_blacklist():
517517 _foo .close ()
518518 expected = [
519519 "3R\t 0\t 100\t 0\n " ,
520- "3R\t 100\t 150\t -0.13392672\n " ,
521- "3R\t 150\t 200\t -0.097087175\n " ,
522- "3R\t 200\t 250\t 0.034750275\n " ,
523- "3R\t 250\t 300\t 0.21310511\n " ,
524- "3R\t 300\t 350\t 0.09064066\n " ,
525- "3R\t 350\t 400\t -0.10200438\n " ,
526- "3R\t 400\t 450\t -0.009409866\n " ,
527- "3R\t 450\t 500\t 0.04944959\n " ,
528- "3R\t 500\t 550\t 0.106793426\n " ,
529- "3R\t 550\t 600\t 0.0784064\n " ,
530- "3R\t 600\t 650\t -0.009581514\n " ,
531- "3R\t 650\t 700\t -0.009331797\n " ,
532- "3R\t 700\t 750\t -0.13486266\n " ,
533- "3R\t 750\t 800\t -0.17890581\n " ,
520+ "3R\t 100\t 150\t -0.13\n " ,
521+ "3R\t 150\t 200\t -0.1\n " ,
522+ "3R\t 200\t 250\t 0.03\n " ,
523+ "3R\t 250\t 300\t 0.21\n " ,
524+ "3R\t 300\t 350\t 0.09\n " ,
525+ "3R\t 350\t 400\t -0.1\n " ,
526+ "3R\t 400\t 450\t -0.01\n " ,
527+ "3R\t 450\t 500\t 0.05\n " ,
528+ "3R\t 500\t 550\t 0.11\n " ,
529+ "3R\t 550\t 600\t 0.08\n " ,
530+ "3R\t 600\t 700\t -0.01\n " ,
531+ "3R\t 700\t 750\t -0.13\n " ,
532+ "3R\t 750\t 800\t -0.18\n " ,
534533 "3R\t 800\t 900\t 0\n " ,
535- "3R\t 900\t 950\t 0.3157759 \n " ,
536- "3R\t 950\t 1000\t 0.16194484 \n " ,
537- "3R\t 1000\t 1050\t -0.006837439 \n " ,
534+ "3R\t 900\t 950\t 0.32 \n " ,
535+ "3R\t 950\t 1000\t 0.16 \n " ,
536+ "3R\t 1000\t 1050\t -0.01 \n " ,
538537 "3R\t 1050\t 1500\t 0\n "
539538 ]
540539 assert f"{ resp } " == f"{ expected } " , f"{ resp } != { expected } "
0 commit comments