@@ -180,7 +180,7 @@ def test_bam_compare_ZoverZ():
180
180
_foo = open (outfile , 'r' )
181
181
resp = _foo .readlines ()
182
182
_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 ' ]
184
184
assert f"{ resp } " == f"{ expected } " , f"{ resp } != { expected } "
185
185
unlink (outfile )
186
186
@@ -293,7 +293,7 @@ def test_bam_compare_scale_factors_ratio():
293
293
(scale factors [1,0.5]) (1+1)/(1+1*0.5)=1.33
294
294
"""
295
295
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 ' ]
297
297
assert f"{ resp } " == f"{ expected } " , f"{ resp } != { expected } "
298
298
unlink (outfile )
299
299
@@ -360,22 +360,22 @@ def test_bam_coverage_filter_blacklist():
360
360
361
361
expected = [
362
362
"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 " ,
375
375
"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 " ,
379
379
"3R\t 1050\t 1500\t 0\n "
380
380
]
381
381
@@ -517,24 +517,23 @@ def test_bam_compare_filter_blacklist():
517
517
_foo .close ()
518
518
expected = [
519
519
"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 " ,
534
533
"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 " ,
538
537
"3R\t 1050\t 1500\t 0\n "
539
538
]
540
539
assert f"{ resp } " == f"{ expected } " , f"{ resp } != { expected } "
0 commit comments