Run the bin status calling on each sample.
- -Calculating statuses' bp per each bin for BS_007JTNB8 ...Calculating statuses' bp per each bin for BS_00TRPEQX ...Calculating statuses' bp per each bin for BS_01DQH017 ...Calculating statuses' bp per each bin for BS_01Y5F4PN ...Calculating statuses' bp per each bin for BS_02W5H7K5 ...Calculating statuses' bp per each bin for BS_02YBZSBY ...Calculating statuses' bp per each bin for BS_02ZSVZCB ...Calculating statuses' bp per each bin for BS_03SSN1N2 ...Calculating statuses' bp per each bin for BS_042DVDQM ...Calculating statuses' bp per each bin for BS_04B8X54S ...Calculating statuses' bp per each bin for BS_05S9WJW6 ...Calculating statuses' bp per each bin for BS_060QB7XB ...Calculating statuses' bp per each bin for BS_06GEF6DB ...Calculating statuses' bp per each bin for BS_06GPPDHZ ...Calculating statuses' bp per each bin for BS_090GV72Q ...Calculating statuses' bp per each bin for BS_0ATJ22QA ...Calculating statuses' bp per each bin for BS_0DKPGQWD ...Calculating statuses' bp per each bin for BS_0EN5R1AP ...Calculating statuses' bp per each bin for BS_0FQKT8EY ...Calculating statuses' bp per each bin for BS_0J5GC803 ...Calculating statuses' bp per each bin for BS_0J5MCBZV ...Calculating statuses' bp per each bin for BS_0Q8G4N8H ...Calculating statuses' bp per each bin for BS_0QEBZX3M ...Calculating statuses' bp per each bin for BS_0RF0G9SQ ...Calculating statuses' bp per each bin for BS_0S617SMM ...Calculating statuses' bp per each bin for BS_0SFXXBQ5 ...Calculating statuses' bp per each bin for BS_0W8AWY10 ...Calculating statuses' bp per each bin for BS_0XD58YD9 ...Calculating statuses' bp per each bin for BS_0XR7GYAV ...Calculating statuses' bp per each bin for BS_0YVR8Q4E ...Calculating statuses' bp per each bin for BS_0ZR4XA69 ...Calculating statuses' bp per each bin for BS_10AD3EFX ...Calculating statuses' bp per each bin for BS_10HD18B8 ...Calculating statuses' bp per each bin for BS_11322HD1 ...Calculating statuses' bp per each bin for BS_1135HC0V ...Calculating statuses' bp per each bin for BS_12ZB7R6A ...Calculating statuses' bp per each bin for BS_154FP4ZQ ...Calculating statuses' bp per each bin for BS_164KRQJ0 ...Calculating statuses' bp per each bin for BS_16TBX6G5 ...Calculating statuses' bp per each bin for BS_17AXPP1Y ...Calculating statuses' bp per each bin for BS_17FVWMNV ...Calculating statuses' bp per each bin for BS_185Q218S ...Calculating statuses' bp per each bin for BS_18NCV5QZ ...Calculating statuses' bp per each bin for BS_19EJ85F8 ...Calculating statuses' bp per each bin for BS_1AZ8YJSH ...Calculating statuses' bp per each bin for BS_1B4Q7CNG ...Calculating statuses' bp per each bin for BS_1C4WCNQF ...Calculating statuses' bp per each bin for BS_1CZPVCXR ...Calculating statuses' bp per each bin for BS_1EQGGHHQ ...Calculating statuses' bp per each bin for BS_1EYKDS0Z ...Calculating statuses' bp per each bin for BS_1GK2EMFY ...Calculating statuses' bp per each bin for BS_1JGQPJH3 ...Calculating statuses' bp per each bin for BS_1JMTTKMK ...Calculating statuses' bp per each bin for BS_1M63B97V ...Calculating statuses' bp per each bin for BS_1MME7FBS ...Calculating statuses' bp per each bin for BS_1Q524P3B ...Calculating statuses' bp per each bin for BS_1RF75MK2 ...Calculating statuses' bp per each bin for BS_1RFBH1SP ...Calculating statuses' bp per each bin for BS_1S2BHJ8K ...Calculating statuses' bp per each bin for BS_1S5V7WWC ...Calculating statuses' bp per each bin for BS_1VMBA3BN ...Calculating statuses' bp per each bin for BS_1YTHM07J ...Calculating statuses' bp per each bin for BS_20TBZG09 ...Calculating statuses' bp per each bin for BS_20X874J5 ...Calculating statuses' bp per each bin for BS_22VCR7DF ...Calculating statuses' bp per each bin for BS_230S192D ...Calculating statuses' bp per each bin for BS_235QEYMQ ...Calculating statuses' bp per each bin for BS_23ZE42Q3 ...Calculating statuses' bp per each bin for BS_245CTSH9 ...Calculating statuses' bp per each bin for BS_25VZM81R ...Calculating statuses' bp per each bin for BS_2AWCKNB3 ...Calculating statuses' bp per each bin for BS_2CD1XWVP ...Calculating statuses' bp per each bin for BS_2CQ6K2JW ...Calculating statuses' bp per each bin for BS_2E81A3FT ...Calculating statuses' bp per each bin for BS_2F70R10Z ...Calculating statuses' bp per each bin for BS_2HNHCSNB ...Calculating statuses' bp per each bin for BS_2J4FG4HV ...Calculating statuses' bp per each bin for BS_2JDVY4F7 ...Calculating statuses' bp per each bin for BS_2NN6ME5A ...Calculating statuses' bp per each bin for BS_2P7MJGM7 ...Calculating statuses' bp per each bin for BS_2RJN5Z4W ...Calculating statuses' bp per each bin for BS_2RX2NF1Q ...Calculating statuses' bp per each bin for BS_2SG8JNMB ...Calculating statuses' bp per each bin for BS_2VB8649B ...Calculating statuses' bp per each bin for BS_2VHFP7M8 ...Calculating statuses' bp per each bin for BS_2VSN5YC5 ...Calculating statuses' bp per each bin for BS_2W200YK5 ...Calculating statuses' bp per each bin for BS_2Z9JQDPV ...Calculating statuses' bp per each bin for BS_30W38X77 ...Calculating statuses' bp per each bin for BS_3108RS7G ...Calculating statuses' bp per each bin for BS_32HAVF6S ...Calculating statuses' bp per each bin for BS_32JF8TPP ...Calculating statuses' bp per each bin for BS_32VQRFDS ...Calculating statuses' bp per each bin for BS_354DTKHE ...Calculating statuses' bp per each bin for BS_36HCZ3KW ...Calculating statuses' bp per each bin for BS_376Q9QKG ...Calculating statuses' bp per each bin for BS_39Q80NT3 ...Calculating statuses' bp per each bin for BS_3A0CJMRH ...Calculating statuses' bp per each bin for BS_3AERRZQT ...Calculating statuses' bp per each bin for BS_3BMGF6GN ...Calculating statuses' bp per each bin for BS_3C4YAMKJ ...Calculating statuses' bp per each bin for BS_3CG2N4PD ...Calculating statuses' bp per each bin for BS_3DYRAZE7 ...Calculating statuses' bp per each bin for BS_3E5C1PN1 ...Calculating statuses' bp per each bin for BS_3J4T2YYW ...Calculating statuses' bp per each bin for BS_3JE2T7AR ...Calculating statuses' bp per each bin for BS_3MYQP7P5 ...Calculating statuses' bp per each bin for BS_3NX3RBCX ...Calculating statuses' bp per each bin for BS_3P96RT4E ...Calculating statuses' bp per each bin for BS_3QDRWJCF ...Calculating statuses' bp per each bin for BS_3R0KN6Z4 ...Calculating statuses' bp per each bin for BS_3VKW5988 ...Calculating statuses' bp per each bin for BS_3Y44RJXS ...Calculating statuses' bp per each bin for BS_3Z40EZHD ...Calculating statuses' bp per each bin for BS_3ZASRA3A ...Calculating statuses' bp per each bin for BS_3ZMH8578 ...Calculating statuses' bp per each bin for BS_40FX713C ...Calculating statuses' bp per each bin for BS_40SJZGT2 ...Calculating statuses' bp per each bin for BS_41SWNWAG ...Calculating statuses' bp per each bin for BS_41TD9XET ...Calculating statuses' bp per each bin for BS_42A66SH4 ...Calculating statuses' bp per each bin for BS_45913DTA ...Calculating statuses' bp per each bin for BS_47454JJ6 ...Calculating statuses' bp per each bin for BS_47X3VAM0 ...Calculating statuses' bp per each bin for BS_49JD1V3S ...Calculating statuses' bp per each bin for BS_4DTX295Z ...Calculating statuses' bp per each bin for BS_4DY83R02 ...Calculating statuses' bp per each bin for BS_4G9Y2VC6 ...Calculating statuses' bp per each bin for BS_4GKH983E ...Calculating statuses' bp per each bin for BS_4M0ZMCDC ...Calculating statuses' bp per each bin for BS_4NPEMNAD ...Calculating statuses' bp per each bin for BS_4P25NNZD ...Calculating statuses' bp per each bin for BS_4PAE19ZC ...Calculating statuses' bp per each bin for BS_4QFSH7C4 ...Calculating statuses' bp per each bin for BS_4RS1SC48 ...Calculating statuses' bp per each bin for BS_4T8J7JEP ...Calculating statuses' bp per each bin for BS_4XPPZTGG ...Calculating statuses' bp per each bin for BS_4Z6F1HJZ ...Calculating statuses' bp per each bin for BS_4ZKN0WGS ...Calculating statuses' bp per each bin for BS_50XV7CTQ ...Calculating statuses' bp per each bin for BS_53TV75NN ...Calculating statuses' bp per each bin for BS_55N7KFHW ...Calculating statuses' bp per each bin for BS_58X7SKJS ...Calculating statuses' bp per each bin for BS_59FR1NC2 ...Calculating statuses' bp per each bin for BS_5B6XZ7YP ...Calculating statuses' bp per each bin for BS_5BMNK8SY ...Calculating statuses' bp per each bin for BS_5BNV7MCF ...Calculating statuses' bp per each bin for BS_5D24XV4T ...Calculating statuses' bp per each bin for BS_5J5VH3X0 ...Calculating statuses' bp per each bin for BS_5JC116NM ...Calculating statuses' bp per each bin for BS_5KFAEMQJ ...Calculating statuses' bp per each bin for BS_5KK9P7XD ...Calculating statuses' bp per each bin for BS_5NAKYS1S ...Calculating statuses' bp per each bin for BS_5R0HHQ1Y ...Calculating statuses' bp per each bin for BS_5RD8KEED ...Calculating statuses' bp per each bin for BS_5S8VXASX ...Calculating statuses' bp per each bin for BS_5TT6TT4K ...Calculating statuses' bp per each bin for BS_5V2FBYQG ...Calculating statuses' bp per each bin for BS_5V2XCW17 ...Calculating statuses' bp per each bin for BS_5WK2KXYW ...Calculating statuses' bp per each bin for BS_5XZP7F4Q ...Calculating statuses' bp per each bin for BS_5Y2Q63Y4 ...Calculating statuses' bp per each bin for BS_5Z4XQC9X ...Calculating statuses' bp per each bin for BS_5ZRZC3ZM ...Calculating statuses' bp per each bin for BS_62XRDTM6 ...Calculating statuses' bp per each bin for BS_682Z7WH6 ...Calculating statuses' bp per each bin for BS_68TZMZH1 ...Calculating statuses' bp per each bin for BS_69VS8PS1 ...Calculating statuses' bp per each bin for BS_6F49F7WH ...Calculating statuses' bp per each bin for BS_6FV7BDE3 ...Calculating statuses' bp per each bin for BS_6FYTRKJJ ...Calculating statuses' bp per each bin for BS_6GN1FGTB ...Calculating statuses' bp per each bin for BS_6GTKRGNK ...Calculating statuses' bp per each bin for BS_6GV08HTE ...Calculating statuses' bp per each bin for BS_6HE8H9HY ...Calculating statuses' bp per each bin for BS_6JBE0947 ...Calculating statuses' bp per each bin for BS_6KEYD1S4 ...Calculating statuses' bp per each bin for BS_6NPSZJ4C ...Calculating statuses' bp per each bin for BS_6Q9PMNWY ...Calculating statuses' bp per each bin for BS_6SA5B35V ...Calculating statuses' bp per each bin for BS_6T87HSPX ...Calculating statuses' bp per each bin for BS_6TMPZKSZ ...Calculating statuses' bp per each bin for BS_6VKKMM6M ...Calculating statuses' bp per each bin for BS_6X2BQQGE ...Calculating statuses' bp per each bin for BS_6X3Z0AVR ...Calculating statuses' bp per each bin for BS_6XHDR1WS ...Calculating statuses' bp per each bin for BS_6XP91JM2 ...Calculating statuses' bp per each bin for BS_6Y7FX7QJ ...Calculating statuses' bp per each bin for BS_6Z213H2V ...Calculating statuses' bp per each bin for BS_6Z4WHYDG ...Calculating statuses' bp per each bin for BS_6ZP2KRPW ...Calculating statuses' bp per each bin for BS_70FKSKM9 ...Calculating statuses' bp per each bin for BS_71HPXFE7 ...Calculating statuses' bp per each bin for BS_742CDKAS ...Calculating statuses' bp per each bin for BS_74A1TB03 ...Calculating statuses' bp per each bin for BS_75VHXMZN ...Calculating statuses' bp per each bin for BS_78TQS4Y3 ...Calculating statuses' bp per each bin for BS_797V83YM ...Calculating statuses' bp per each bin for BS_79NQJZ09 ...Calculating statuses' bp per each bin for BS_79SYEHY3 ...Calculating statuses' bp per each bin for BS_7BNQBB26 ...Calculating statuses' bp per each bin for BS_7F07M7JG ...Calculating statuses' bp per each bin for BS_7GKF6M85 ...Calculating statuses' bp per each bin for BS_7HC4W4SA ...Calculating statuses' bp per each bin for BS_7M5272WN ...Calculating statuses' bp per each bin for BS_7M6F9GQD ...Calculating statuses' bp per each bin for BS_7M7JNG00 ...Calculating statuses' bp per each bin for BS_7PF3C1P7 ...Calculating statuses' bp per each bin for BS_7PKPNFFK ...Calculating statuses' bp per each bin for BS_7RQCH5Y7 ...Calculating statuses' bp per each bin for BS_7XZNDFV4 ...Calculating statuses' bp per each bin for BS_7Y5E5KQW ...Calculating statuses' bp per each bin for BS_80078QDG ...Calculating statuses' bp per each bin for BS_80X7AVCP ...Calculating statuses' bp per each bin for BS_8226HD4X ...Calculating statuses' bp per each bin for BS_823V5X6Z ...Calculating statuses' bp per each bin for BS_83MSWR59 ...Calculating statuses' bp per each bin for BS_84AQCZE4 ...Calculating statuses' bp per each bin for BS_853PNV7P ...Calculating statuses' bp per each bin for BS_859AV1DB ...Calculating statuses' bp per each bin for BS_85Q5P8GF ...Calculating statuses' bp per each bin for BS_85Y70BHE ...Calculating statuses' bp per each bin for BS_86QSEFS7 ...Calculating statuses' bp per each bin for BS_872Q3FK2 ...Calculating statuses' bp per each bin for BS_886M7JMG ...Calculating statuses' bp per each bin for BS_895FJ34T ...Calculating statuses' bp per each bin for BS_89RCN9PF ...Calculating statuses' bp per each bin for BS_8AY2GM4G ...Calculating statuses' bp per each bin for BS_8B1YF57H ...Calculating statuses' bp per each bin for BS_8BD8WD38 ...Calculating statuses' bp per each bin for BS_8BM1WTGR ...Calculating statuses' bp per each bin for BS_8BS10XE2 ...Calculating statuses' bp per each bin for BS_8CA214YR ...Calculating statuses' bp per each bin for BS_8D52JK1Q ...Calculating statuses' bp per each bin for BS_8EB1S8MK ...Calculating statuses' bp per each bin for BS_8FSJDG2T ...Calculating statuses' bp per each bin for BS_8H2XX24V ...Calculating statuses' bp per each bin for BS_8JFMP1T1 ...Calculating statuses' bp per each bin for BS_8Q8CAY84 ...Calculating statuses' bp per each bin for BS_8Q8N6AHW ...Calculating statuses' bp per each bin for BS_8S2CZBA6 ...Calculating statuses' bp per each bin for BS_8SETK1J4 ...Calculating statuses' bp per each bin for BS_8SYN7GXG ...Calculating statuses' bp per each bin for BS_8TFBNHT0 ...Calculating statuses' bp per each bin for BS_8W0FNBWF ...Calculating statuses' bp per each bin for BS_8WJR4AA3 ...Calculating statuses' bp per each bin for BS_8YQGSMDK ...Calculating statuses' bp per each bin for BS_922YMFYK ...Calculating statuses' bp per each bin for BS_928HWFFM ...Calculating statuses' bp per each bin for BS_93B38HPS ...Calculating statuses' bp per each bin for BS_93BV8AY9 ...Calculating statuses' bp per each bin for BS_947CK40E ...Calculating statuses' bp per each bin for BS_94AP61Q5 ...Calculating statuses' bp per each bin for BS_95RRF4NE ...Calculating statuses' bp per each bin for BS_969K7ZM1 ...Calculating statuses' bp per each bin for BS_96S0VQBN ...Calculating statuses' bp per each bin for BS_99PPRCW4 ...Calculating statuses' bp per each bin for BS_9BF58P80 ...Calculating statuses' bp per each bin for BS_9BGCH82S ...Calculating statuses' bp per each bin for BS_9BN45DFK ...Calculating statuses' bp per each bin for BS_9CMJ3QF6 ...Calculating statuses' bp per each bin for BS_9D39MFHT ...Calculating statuses' bp per each bin for BS_9E632VPJ ...Calculating statuses' bp per each bin for BS_9E6Q11G2 ...Calculating statuses' bp per each bin for BS_9F2J6VS1 ...Calculating statuses' bp per each bin for BS_9G7ETYHR ...Calculating statuses' bp per each bin for BS_9GJHMA3J ...Calculating statuses' bp per each bin for BS_9GN1QA3Q ...Calculating statuses' bp per each bin for BS_9GSWPZ44 ...Calculating statuses' bp per each bin for BS_9HDPJT21 ...Calculating statuses' bp per each bin for BS_9HEMB7RK ...Calculating statuses' bp per each bin for BS_9JVGXA2W ...Calculating statuses' bp per each bin for BS_9KDA3HG6 ...Calculating statuses' bp per each bin for BS_9M49SGAQ ...Calculating statuses' bp per each bin for BS_9N3B3HZB ...Calculating statuses' bp per each bin for BS_9NF39QAW ...Calculating statuses' bp per each bin for BS_9PQPZPGM ...Calculating statuses' bp per each bin for BS_9QMVZ4C0 ...Calculating statuses' bp per each bin for BS_9R82A3VT ...Calculating statuses' bp per each bin for BS_9SVMVJ45 ...Calculating statuses' bp per each bin for BS_9W6FK6X4 ...Calculating statuses' bp per each bin for BS_9X9E86DY ...Calculating statuses' bp per each bin for BS_9YPJANGX ...Calculating statuses' bp per each bin for BS_9ZFXXJPK ...Calculating statuses' bp per each bin for BS_A3EV6FXK ...Calculating statuses' bp per each bin for BS_A4KYP5H0 ...Calculating statuses' bp per each bin for BS_A70G7S2W ...Calculating statuses' bp per each bin for BS_A94A7EJ4 ...Calculating statuses' bp per each bin for BS_ABENB1H2 ...Calculating statuses' bp per each bin for BS_ABZ3BK38 ...Calculating statuses' bp per each bin for BS_ACE6QKVF ...Calculating statuses' bp per each bin for BS_ACTR13BQ ...Calculating statuses' bp per each bin for BS_AEMRNSQX ...Calculating statuses' bp per each bin for BS_AF5D41PD ...Calculating statuses' bp per each bin for BS_AFBPM6CN ...Calculating statuses' bp per each bin for BS_AG0BZN6F ...Calculating statuses' bp per each bin for BS_AG6M7TJE ...Calculating statuses' bp per each bin for BS_AGAW7M7W ...Calculating statuses' bp per each bin for BS_AGD2ATY1 ...Calculating statuses' bp per each bin for BS_AGGH3QGS ...Calculating statuses' bp per each bin for BS_AHAXPFG3 ...Calculating statuses' bp per each bin for BS_AJ0T7DJ7 ...Calculating statuses' bp per each bin for BS_AJAH8KG5 ...Calculating statuses' bp per each bin for BS_AJFG2QDD ...Calculating statuses' bp per each bin for BS_AK9BV52G ...Calculating statuses' bp per each bin for BS_APAQ2X2S ...Calculating statuses' bp per each bin for BS_AQECMDPP ...Calculating statuses' bp per each bin for BS_AQMKA8NC ...Calculating statuses' bp per each bin for BS_AQQY8HMR ...Calculating statuses' bp per each bin for BS_AR2BH0Z4 ...Calculating statuses' bp per each bin for BS_ARE4EM81 ...Calculating statuses' bp per each bin for BS_ASXY97FK ...Calculating statuses' bp per each bin for BS_AVS4DSZW ...Calculating statuses' bp per each bin for BS_AWDQEKK1 ...Calculating statuses' bp per each bin for BS_AWV1EDYB ...Calculating statuses' bp per each bin for BS_AXM96K20 ...Calculating statuses' bp per each bin for BS_AXW0XBH8 ...Calculating statuses' bp per each bin for BS_AYPBAHJ6 ...Calculating statuses' bp per each bin for BS_B1X2SB4P ...Calculating statuses' bp per each bin for BS_B3RH9V17 ...Calculating statuses' bp per each bin for BS_B4DY7ET3 ...Calculating statuses' bp per each bin for BS_B4P1NP87 ...Calculating statuses' bp per each bin for BS_B558FJXC ...Calculating statuses' bp per each bin for BS_B5Z2XV03 ...Calculating statuses' bp per each bin for BS_B6K0N6HS ...Calculating statuses' bp per each bin for BS_B6KCJF2B ...Calculating statuses' bp per each bin for BS_B73JK7T1 ...Calculating statuses' bp per each bin for BS_B91XGSA5 ...Calculating statuses' bp per each bin for BS_B94G701Q ...Calculating statuses' bp per each bin for BS_B9QP40ER ...Calculating statuses' bp per each bin for BS_BA0M71QZ ...Calculating statuses' bp per each bin for BS_BAVWP5JG ...Calculating statuses' bp per each bin for BS_BBCM52QE ...Calculating statuses' bp per each bin for BS_BBHFKBE7 ...Calculating statuses' bp per each bin for BS_BBTT7DH4 ...Calculating statuses' bp per each bin for BS_BD4RQ1G0 ...Calculating statuses' bp per each bin for BS_BFDEZK1C ...Calculating statuses' bp per each bin for BS_BM95DGCQ ...Calculating statuses' bp per each bin for BS_BMHYK93N ...Calculating statuses' bp per each bin for BS_BNFCXVYZ ...Calculating statuses' bp per each bin for BS_BQ81D2BP ...Calculating statuses' bp per each bin for BS_BQWMQM2Y ...Calculating statuses' bp per each bin for BS_BR1X3R68 ...Calculating statuses' bp per each bin for BS_BS5X4H0Y ...Calculating statuses' bp per each bin for BS_BSM3ZHW4 ...Calculating statuses' bp per each bin for BS_BVD7PWP5 ...Calculating statuses' bp per each bin for BS_BWZTMWTM ...Calculating statuses' bp per each bin for BS_BX8AR74E ...Calculating statuses' bp per each bin for BS_BYJDRXMK ...Calculating statuses' bp per each bin for BS_BZ4GKDP4 ...Calculating statuses' bp per each bin for BS_BZWM455T ...Calculating statuses' bp per each bin for BS_C0YDDQKB ...Calculating statuses' bp per each bin for BS_C19DFVR5 ...Calculating statuses' bp per each bin for BS_C1RRDNF5 ...Calculating statuses' bp per each bin for BS_C2NH5FDT ...Calculating statuses' bp per each bin for BS_C3KDA5ZC ...Calculating statuses' bp per each bin for BS_C3X4AEK5 ...Calculating statuses' bp per each bin for BS_C5RDMCCS ...Calculating statuses' bp per each bin for BS_C6GBYR17 ...Calculating statuses' bp per each bin for BS_C6NB914B ...Calculating statuses' bp per each bin for BS_C9TNGEA0 ...Calculating statuses' bp per each bin for BS_CBMAWSAR ...Calculating statuses' bp per each bin for BS_CDKC1E2B ...Calculating statuses' bp per each bin for BS_CFEWVDG1 ...Calculating statuses' bp per each bin for BS_CFPWSZB5 ...Calculating statuses' bp per each bin for BS_CGCNC41M ...Calculating statuses' bp per each bin for BS_CGXTFM67 ...Calculating statuses' bp per each bin for BS_CHAFNK64 ...Calculating statuses' bp per each bin for BS_CK7CPM73 ...Calculating statuses' bp per each bin for BS_CKY2TTXG ...Calculating statuses' bp per each bin for BS_CN7PB0DN ...Calculating statuses' bp per each bin for BS_CP435MAC ...Calculating statuses' bp per each bin for BS_CPYV7AHD ...Calculating statuses' bp per each bin for BS_CPZZR2R7 ...Calculating statuses' bp per each bin for BS_CQ8C6X4X ...Calculating statuses' bp per each bin for BS_CRKBDAYZ ...Calculating statuses' bp per each bin for BS_CSY11X7E ...Calculating statuses' bp per each bin for BS_CWQ35N2K ...Calculating statuses' bp per each bin for BS_CXCKHP3Z ...Calculating statuses' bp per each bin for BS_CXDKK6KS ...Calculating statuses' bp per each bin for BS_CXN0498B ...Calculating statuses' bp per each bin for BS_CZRA594T ...Calculating statuses' bp per each bin for BS_D09VWF8H ...Calculating statuses' bp per each bin for BS_D2BZ25FW ...Calculating statuses' bp per each bin for BS_D368BNRD ...Calculating statuses' bp per each bin for BS_D3KVGMRJ ...Calculating statuses' bp per each bin for BS_D3XVXF05 ...Calculating statuses' bp per each bin for BS_D3ZV6MMY ...Calculating statuses' bp per each bin for BS_D59VEERR ...Calculating statuses' bp per each bin for BS_D6STCMQS ...Calculating statuses' bp per each bin for BS_D6XHKZDZ ...Calculating statuses' bp per each bin for BS_DAEM3WRX ...Calculating statuses' bp per each bin for BS_DC7XYYG8 ...Calculating statuses' bp per each bin for BS_DE26D072 ...Calculating statuses' bp per each bin for BS_DGG7YE8R ...Calculating statuses' bp per each bin for BS_DH86PF5A ...Calculating statuses' bp per each bin for BS_DJEVCH34 ...Calculating statuses' bp per each bin for BS_DMBNB0EP ...Calculating statuses' bp per each bin for BS_DT92H1MJ ...Calculating statuses' bp per each bin for BS_DVDT4VXQ ...Calculating statuses' bp per each bin for BS_DVQYQSVW ...Calculating statuses' bp per each bin for BS_DW1CYEXP ...Calculating statuses' bp per each bin for BS_E0AAK0V0 ...Calculating statuses' bp per each bin for BS_E0S2Y0TS ...Calculating statuses' bp per each bin for BS_E1VJCEWB ...Calculating statuses' bp per each bin for BS_E2RXQ3JW ...Calculating statuses' bp per each bin for BS_E42SSQGJ ...Calculating statuses' bp per each bin for BS_E5B2YZJH ...Calculating statuses' bp per each bin for BS_E5H6CFYT ...Calculating statuses' bp per each bin for BS_E7HEQZ1K ...Calculating statuses' bp per each bin for BS_E8MC6R2K ...Calculating statuses' bp per each bin for BS_E8SJ36A6 ...Calculating statuses' bp per each bin for BS_E9ATGWME ...Calculating statuses' bp per each bin for BS_E9M7TDB6 ...Calculating statuses' bp per each bin for BS_EDCQ4M0R ...Calculating statuses' bp per each bin for BS_EDZCAHR2 ...Calculating statuses' bp per each bin for BS_EE73VE7V ...Calculating statuses' bp per each bin for BS_EFBBP9H3 ...Calculating statuses' bp per each bin for BS_EHBH0NYX ...Calculating statuses' bp per each bin for BS_EJ1H9PZY ...Calculating statuses' bp per each bin for BS_EJV0N3BX ...Calculating statuses' bp per each bin for BS_EKR4NSNW ...Calculating statuses' bp per each bin for BS_EP25TTAG ...Calculating statuses' bp per each bin for BS_EPM623G4 ...Calculating statuses' bp per each bin for BS_EQE19CB0 ...Calculating statuses' bp per each bin for BS_ERAWW3H7 ...Calculating statuses' bp per each bin for BS_ERFMPQN3 ...Calculating statuses' bp per each bin for BS_ETTZM63Y ...Calculating statuses' bp per each bin for BS_EXTEGB51 ...Calculating statuses' bp per each bin for BS_F064HMTP ...Calculating statuses' bp per each bin for BS_F0GNWEJJ ...Calculating statuses' bp per each bin for BS_F1B71HTZ ...Calculating statuses' bp per each bin for BS_F48RA7QF ...Calculating statuses' bp per each bin for BS_F4D169Q7 ...Calculating statuses' bp per each bin for BS_F5641KTR ...Calculating statuses' bp per each bin for BS_F5KRHGKC ...Calculating statuses' bp per each bin for BS_F5V8CWXT ...Calculating statuses' bp per each bin for BS_F7KYPE79 ...Calculating statuses' bp per each bin for BS_F8K4VQMF ...Calculating statuses' bp per each bin for BS_F8KR1QWE ...Calculating statuses' bp per each bin for BS_F969DNHD ...Calculating statuses' bp per each bin for BS_F9ZTZ4YS ...Calculating statuses' bp per each bin for BS_FBJ516WW ...Calculating statuses' bp per each bin for BS_FBP5NZNR ...Calculating statuses' bp per each bin for BS_FCXEAPJ1 ...Calculating statuses' bp per each bin for BS_FF73TT6D ...Calculating statuses' bp per each bin for BS_FFD25TTK ...Calculating statuses' bp per each bin for BS_FG649AWK ...Calculating statuses' bp per each bin for BS_FHE6Z73G ...Calculating statuses' bp per each bin for BS_FHFVCQ41 ...Calculating statuses' bp per each bin for BS_FJEBFQHG ...Calculating statuses' bp per each bin for BS_FJEZ3ASV ...Calculating statuses' bp per each bin for BS_FK3B5SDH ...Calculating statuses' bp per each bin for BS_FKQ7F6D1 ...Calculating statuses' bp per each bin for BS_FNF1TCSW ...Calculating statuses' bp per each bin for BS_FNM4BB83 ...Calculating statuses' bp per each bin for BS_FQ2ZBMPA ...Calculating statuses' bp per each bin for BS_FSM39DXA ...Calculating statuses' bp per each bin for BS_FT95T4YP ...Calculating statuses' bp per each bin for BS_FVPMPMRJ ...Calculating statuses' bp per each bin for BS_FVYBGMG1 ...Calculating statuses' bp per each bin for BS_FWCKP4X1 ...Calculating statuses' bp per each bin for BS_FXDY881S ...Calculating statuses' bp per each bin for BS_FYJMNRGK ...Calculating statuses' bp per each bin for BS_G12HNXYX ...Calculating statuses' bp per each bin for BS_G16NM5WJ ...Calculating statuses' bp per each bin for BS_G2016BE5 ...Calculating statuses' bp per each bin for BS_G4EX0H98 ...Calculating statuses' bp per each bin for BS_G61CN2V0 ...Calculating statuses' bp per each bin for BS_G9MQM1KK ...Calculating statuses' bp per each bin for BS_GAA327FP ...Calculating statuses' bp per each bin for BS_GCQEB05M ...Calculating statuses' bp per each bin for BS_GDPZN3SD ...Calculating statuses' bp per each bin for BS_GFS9JEE6 ...Calculating statuses' bp per each bin for BS_GKN0B47F ...Calculating statuses' bp per each bin for BS_GKPVZXYP ...Calculating statuses' bp per each bin for BS_GMQZNQH8 ...Calculating statuses' bp per each bin for BS_GRECE8Q9 ...Calculating statuses' bp per each bin for BS_GRS3T427 ...Calculating statuses' bp per each bin for BS_GS9JCK55 ...Calculating statuses' bp per each bin for BS_GTAJFTTR ...Calculating statuses' bp per each bin for BS_GWPJKM7G ...Calculating statuses' bp per each bin for BS_GXACEYKB ...Calculating statuses' bp per each bin for BS_H0638GZT ...Calculating statuses' bp per each bin for BS_H174820R ...Calculating statuses' bp per each bin for BS_H1K33JVK ...Calculating statuses' bp per each bin for BS_H1MNXNHB ...Calculating statuses' bp per each bin for BS_H2K7XX4X ...Calculating statuses' bp per each bin for BS_H50HR85Y ...Calculating statuses' bp per each bin for BS_H83DTMT2 ...Calculating statuses' bp per each bin for BS_H8E4K1DT ...Calculating statuses' bp per each bin for BS_H8MAKDRB ...Calculating statuses' bp per each bin for BS_H8NWA41N ...Calculating statuses' bp per each bin for BS_HBAKYKQM ...Calculating statuses' bp per each bin for BS_HC4T1EEP ...Calculating statuses' bp per each bin for BS_HEJ72V3F ...Calculating statuses' bp per each bin for BS_HGP0PV0D ...Calculating statuses' bp per each bin for BS_HJ95JNFN ...Calculating statuses' bp per each bin for BS_HM5GFJN8 ...Calculating statuses' bp per each bin for BS_HM5N78GV ...Calculating statuses' bp per each bin for BS_HN8DE43A ...Calculating statuses' bp per each bin for BS_HNE6BPMH ...Calculating statuses' bp per each bin for BS_HPMX6CT4 ...Calculating statuses' bp per each bin for BS_HQFNQHVW ...Calculating statuses' bp per each bin for BS_HSYF7G2D ...Calculating statuses' bp per each bin for BS_HTGDSPQA ...Calculating statuses' bp per each bin for BS_HW4J7JM8 ...Calculating statuses' bp per each bin for BS_HYGK88B0 ...Calculating statuses' bp per each bin for BS_HYKV2TH9 ...Calculating statuses' bp per each bin for BS_HZHEF797 ...Calculating statuses' bp per each bin for BS_HZNKSQ17 ...Calculating statuses' bp per each bin for BS_HZV4WDTB ...Calculating statuses' bp per each bin for BS_J037K4BM ...Calculating statuses' bp per each bin for BS_J33YVP27 ...Calculating statuses' bp per each bin for BS_J4E9SW51 ...Calculating statuses' bp per each bin for BS_J52D4XQY ...Calculating statuses' bp per each bin for BS_J5J0EA7Y ...Calculating statuses' bp per each bin for BS_J6VRW67C ...Calculating statuses' bp per each bin for BS_J8EK6RNF ...Calculating statuses' bp per each bin for BS_JA3XXP1Y ...Calculating statuses' bp per each bin for BS_JB3J82ZK ...Calculating statuses' bp per each bin for BS_JBHQ6G3J ...Calculating statuses' bp per each bin for BS_JBY9S5D4 ...Calculating statuses' bp per each bin for BS_JC2187NA ...Calculating statuses' bp per each bin for BS_JDMM2XG4 ...Calculating statuses' bp per each bin for BS_JDS0SYST ...Calculating statuses' bp per each bin for BS_JE1DF581 ...Calculating statuses' bp per each bin for BS_JGKRN7NA ...Calculating statuses' bp per each bin for BS_JHGFJR3D ...Calculating statuses' bp per each bin for BS_JJNSP29S ...Calculating statuses' bp per each bin for BS_JJWC13GG ...Calculating statuses' bp per each bin for BS_JJZPY3BG ...Calculating statuses' bp per each bin for BS_JK5GTHVG ...Calculating statuses' bp per each bin for BS_JM979K1M ...Calculating statuses' bp per each bin for BS_JNS7Q60A ...Calculating statuses' bp per each bin for BS_JQNF0FQP ...Calculating statuses' bp per each bin for BS_JR56RQDC ...Calculating statuses' bp per each bin for BS_JRFVST47 ...Calculating statuses' bp per each bin for BS_JSMAT49K ...Calculating statuses' bp per each bin for BS_JSNJZERZ ...Calculating statuses' bp per each bin for BS_JTBM5TSE ...Calculating statuses' bp per each bin for BS_JVTQD8HH ...Calculating statuses' bp per each bin for BS_JX1PK0VW ...Calculating statuses' bp per each bin for BS_K049HVGR ...Calculating statuses' bp per each bin for BS_K07KNTFY ...Calculating statuses' bp per each bin for BS_K3PPH522 ...Calculating statuses' bp per each bin for BS_K3Y06EKR ...Calculating statuses' bp per each bin for BS_K3YT05M7 ...Calculating statuses' bp per each bin for BS_K6A9Z04J ...Calculating statuses' bp per each bin for BS_K6BMXTW5 ...Calculating statuses' bp per each bin for BS_K89V7JNJ ...Calculating statuses' bp per each bin for BS_KAD49R68 ...Calculating statuses' bp per each bin for BS_KAZENYZE ...Calculating statuses' bp per each bin for BS_KB9GJDCS ...Calculating statuses' bp per each bin for BS_KBPBCWA0 ...Calculating statuses' bp per each bin for BS_KCN7NG92 ...Calculating statuses' bp per each bin for BS_KFRBXXV7 ...Calculating statuses' bp per each bin for BS_KG1YMTA4 ...Calculating statuses' bp per each bin for BS_KH0NFCWF ...Calculating statuses' bp per each bin for BS_KH3859M5 ...Calculating statuses' bp per each bin for BS_KHSYAB3J ...Calculating statuses' bp per each bin for BS_KJQXZ8GK ...Calculating statuses' bp per each bin for BS_KJTW1RBN ...Calculating statuses' bp per each bin for BS_KKWQW01N ...Calculating statuses' bp per each bin for BS_KM878Y43 ...Calculating statuses' bp per each bin for BS_KP0YNTAK ...Calculating statuses' bp per each bin for BS_KPWRDRFP ...Calculating statuses' bp per each bin for BS_KRYYDC2X ...Calculating statuses' bp per each bin for BS_KSHETTQC ...Calculating statuses' bp per each bin for BS_KSRYKCDE ...Calculating statuses' bp per each bin for BS_KT9H4B1N ...Calculating statuses' bp per each bin for BS_KV3XF9WV ...Calculating statuses' bp per each bin for BS_KVV3G6NP ...Calculating statuses' bp per each bin for BS_KVYCH9PB ...Calculating statuses' bp per each bin for BS_KXGJ1HMG ...Calculating statuses' bp per each bin for BS_KY1CTGJ3 ...Calculating statuses' bp per each bin for BS_KYG3H6C6 ...Calculating statuses' bp per each bin for BS_KYPMZZSY ...Calculating statuses' bp per each bin for BS_KYPZRS4S ...Calculating statuses' bp per each bin for BS_M0B42FPR ...Calculating statuses' bp per each bin for BS_M16CDR44 ...Calculating statuses' bp per each bin for BS_M18CDDGK ...Calculating statuses' bp per each bin for BS_M1TDS9CT ...Calculating statuses' bp per each bin for BS_M1YHV93N ...Calculating statuses' bp per each bin for BS_M28HG0NT ...Calculating statuses' bp per each bin for BS_M2MSWKBC ...Calculating statuses' bp per each bin for BS_M4E4H6NG ...Calculating statuses' bp per each bin for BS_M5FM63EB ...Calculating statuses' bp per each bin for BS_M659G06J ...Calculating statuses' bp per each bin for BS_M8EBKZDG ...Calculating statuses' bp per each bin for BS_M9YZK3SP ...Calculating statuses' bp per each bin for BS_MAVSJH2Y ...Calculating statuses' bp per each bin for BS_MBYEK6HV ...Calculating statuses' bp per each bin for BS_MCM78YPC ...Calculating statuses' bp per each bin for BS_MDVMT75X ...Calculating statuses' bp per each bin for BS_ME97J1GA ...Calculating statuses' bp per each bin for BS_MEPW36WF ...Calculating statuses' bp per each bin for BS_MFTAVQ5C ...Calculating statuses' bp per each bin for BS_MGBEMNYR ...Calculating statuses' bp per each bin for BS_MGY2V5N4 ...Calculating statuses' bp per each bin for BS_MHT5C2AP ...Calculating statuses' bp per each bin for BS_MJJZJMTK ...Calculating statuses' bp per each bin for BS_MRX1SPFB ...Calculating statuses' bp per each bin for BS_MS87PMR7 ...Calculating statuses' bp per each bin for BS_MSB1RDAJ ...Calculating statuses' bp per each bin for BS_MSGYZ69T ...Calculating statuses' bp per each bin for BS_MT5ZSRG4 ...Calculating statuses' bp per each bin for BS_MVC1NCZX ...Calculating statuses' bp per each bin for BS_MW16Y7ME ...Calculating statuses' bp per each bin for BS_MW63FCZC ...Calculating statuses' bp per each bin for BS_MW6YMRBJ ...Calculating statuses' bp per each bin for BS_MW74GF38 ...Calculating statuses' bp per each bin for BS_MWKJF6XH ...Calculating statuses' bp per each bin for BS_MWXDJFWW ...Calculating statuses' bp per each bin for BS_MWZCP1XW ...Calculating statuses' bp per each bin for BS_MZJJRDZF ...Calculating statuses' bp per each bin for BS_N0Z6V2YG ...Calculating statuses' bp per each bin for BS_N2JD3VX6 ...Calculating statuses' bp per each bin for BS_N4B3WB8W ...Calculating statuses' bp per each bin for BS_N5VYY66W ...Calculating statuses' bp per each bin for BS_N8D9FH6M ...Calculating statuses' bp per each bin for BS_N8WXTFN4 ...Calculating statuses' bp per each bin for BS_N9BF56FP ...Calculating statuses' bp per each bin for BS_N9N1X78X ...Calculating statuses' bp per each bin for BS_N9SMBR24 ...Calculating statuses' bp per each bin for BS_NASADC3P ...Calculating statuses' bp per each bin for BS_NB509KST ...Calculating statuses' bp per each bin for BS_NCGGFG6H ...Calculating statuses' bp per each bin for BS_NDSWS9TS ...Calculating statuses' bp per each bin for BS_NEX92KQH ...Calculating statuses' bp per each bin for BS_NFT9VHMJ ...Calculating statuses' bp per each bin for BS_NG0CR0EP ...Calculating statuses' bp per each bin for BS_NHMPQ8Q1 ...Calculating statuses' bp per each bin for BS_NHRNZRQV ...Calculating statuses' bp per each bin for BS_NJFK43N3 ...Calculating statuses' bp per each bin for BS_NJV6A585 ...Calculating statuses' bp per each bin for BS_NMJPX8VQ ...Calculating statuses' bp per each bin for BS_NNFDFAFM ...Calculating statuses' bp per each bin for BS_NS8G0465 ...Calculating statuses' bp per each bin for BS_NTXHP7QR ...Calculating statuses' bp per each bin for BS_NWYBD9CA ...Calculating statuses' bp per each bin for BS_NXZSD2KK ...Calculating statuses' bp per each bin for BS_NYVVN7SD ...Calculating statuses' bp per each bin for BS_NZRJE7QA ...Calculating statuses' bp per each bin for BS_P0KQFXNC ...Calculating statuses' bp per each bin for BS_P0QJ1QAH ...Calculating statuses' bp per each bin for BS_P36C0ABE ...Calculating statuses' bp per each bin for BS_P3EVK3C2 ...Calculating statuses' bp per each bin for BS_P3KD5EQG ...Calculating statuses' bp per each bin for BS_P3P96WCD ...Calculating statuses' bp per each bin for BS_P3PF53V8 ...Calculating statuses' bp per each bin for BS_P42V737Y ...Calculating statuses' bp per each bin for BS_P4K6WK9Y ...Calculating statuses' bp per each bin for BS_P6NWC51S ...Calculating statuses' bp per each bin for BS_P7CY04E7 ...Calculating statuses' bp per each bin for BS_P8WN8XEQ ...Calculating statuses' bp per each bin for BS_P99S85CY ...Calculating statuses' bp per each bin for BS_PC2YQ4F6 ...Calculating statuses' bp per each bin for BS_PDRNE606 ...Calculating statuses' bp per each bin for BS_PEDGRK8S ...Calculating statuses' bp per each bin for BS_PEFRDKDZ ...Calculating statuses' bp per each bin for BS_PFXDWV9R ...Calculating statuses' bp per each bin for BS_PJ08ZW9A ...Calculating statuses' bp per each bin for BS_PKZ1HWNB ...Calculating statuses' bp per each bin for BS_PMFR8NPN ...Calculating statuses' bp per each bin for BS_PMSZKZP7 ...Calculating statuses' bp per each bin for BS_PNXDHDF3 ...Calculating statuses' bp per each bin for BS_PP4WJBNR ...Calculating statuses' bp per each bin for BS_PPY0EWPG ...Calculating statuses' bp per each bin for BS_PV130A29 ...Calculating statuses' bp per each bin for BS_PVYMDM18 ...Calculating statuses' bp per each bin for BS_PVZJCXS1 ...Calculating statuses' bp per each bin for BS_PYVPKH11 ...Calculating statuses' bp per each bin for BS_PZDTK87A ...Calculating statuses' bp per each bin for BS_Q0XBZP88 ...Calculating statuses' bp per each bin for BS_Q3ANJDVG ...Calculating statuses' bp per each bin for BS_Q50P97Y0 ...Calculating statuses' bp per each bin for BS_Q5XKTB1V ...Calculating statuses' bp per each bin for BS_Q6GVRAWK ...Calculating statuses' bp per each bin for BS_Q807ENGY ...Calculating statuses' bp per each bin for BS_Q88059DJ ...Calculating statuses' bp per each bin for BS_Q9NHTR7G ...Calculating statuses' bp per each bin for BS_QA0W86WP ...Calculating statuses' bp per each bin for BS_QBZDQX7A ...Calculating statuses' bp per each bin for BS_QC5FNXGH ...Calculating statuses' bp per each bin for BS_QD22MJJ5 ...Calculating statuses' bp per each bin for BS_QDGHHS4S ...Calculating statuses' bp per each bin for BS_QE5B0T7R ...Calculating statuses' bp per each bin for BS_QE7WAAKC ...Calculating statuses' bp per each bin for BS_QG6V29H7 ...Calculating statuses' bp per each bin for BS_QGX93WPF ...Calculating statuses' bp per each bin for BS_QHTDJG77 ...Calculating statuses' bp per each bin for BS_QK61ST4T ...Calculating statuses' bp per each bin for BS_QM1Z8R38 ...Calculating statuses' bp per each bin for BS_QMG7XNEW ...Calculating statuses' bp per each bin for BS_QMHND3QK ...Calculating statuses' bp per each bin for BS_QMY84KF4 ...Calculating statuses' bp per each bin for BS_QQP17PDQ ...Calculating statuses' bp per each bin for BS_QRH793AV ...Calculating statuses' bp per each bin for BS_QSDKJM8T ...Calculating statuses' bp per each bin for BS_QSMFVHSB ...Calculating statuses' bp per each bin for BS_QWM9BPDY ...Calculating statuses' bp per each bin for BS_QWS88QXE ...Calculating statuses' bp per each bin for BS_QX754ADQ ...Calculating statuses' bp per each bin for BS_QXZTRMWM ...Calculating statuses' bp per each bin for BS_QYWGB52C ...Calculating statuses' bp per each bin for BS_QZRP3NSG ...Calculating statuses' bp per each bin for BS_QZT3SPY7 ...Calculating statuses' bp per each bin for BS_QZTPJ8YW ...Calculating statuses' bp per each bin for BS_R1CV479M ...Calculating statuses' bp per each bin for BS_R1G6SG0K ...Calculating statuses' bp per each bin for BS_R1RMKH1B ...Calculating statuses' bp per each bin for BS_R1ZBTEAG ...Calculating statuses' bp per each bin for BS_R3WB0PP7 ...Calculating statuses' bp per each bin for BS_R5G3QP11 ...Calculating statuses' bp per each bin for BS_R6CKWZW6 ...Calculating statuses' bp per each bin for BS_R726FKHP ...Calculating statuses' bp per each bin for BS_R7YMKNJN ...Calculating statuses' bp per each bin for BS_R82MYPKB ...Calculating statuses' bp per each bin for BS_R8ZKDFWH ...Calculating statuses' bp per each bin for BS_R913E3VB ...Calculating statuses' bp per each bin for BS_R922P3XJ ...Calculating statuses' bp per each bin for BS_R94EDP5H ...Calculating statuses' bp per each bin for BS_R9V82BC1 ...Calculating statuses' bp per each bin for BS_RA0VCMWC ...Calculating statuses' bp per each bin for BS_RB8VSEVG ...Calculating statuses' bp per each bin for BS_RCAZ62J7 ...Calculating statuses' bp per each bin for BS_RG6JK7ZN ...Calculating statuses' bp per each bin for BS_RGSJMSRE ...Calculating statuses' bp per each bin for BS_RH9PMH3G ...Calculating statuses' bp per each bin for BS_RHEBGZF8 ...Calculating statuses' bp per each bin for BS_RQ5Z7KNW ...Calculating statuses' bp per each bin for BS_RR3DKGKY ...Calculating statuses' bp per each bin for BS_RXP2ZRQT ...Calculating statuses' bp per each bin for BS_RXSBJ929 ...Calculating statuses' bp per each bin for BS_RZNNEX1Y ...Calculating statuses' bp per each bin for BS_S0T3CQ97 ...Calculating statuses' bp per each bin for BS_S34SCYW4 ...Calculating statuses' bp per each bin for BS_S791VC80 ...Calculating statuses' bp per each bin for BS_S88G7BXH ...Calculating statuses' bp per each bin for BS_SCQNN2A2 ...Calculating statuses' bp per each bin for BS_SD28E0SC ...Calculating statuses' bp per each bin for BS_SDT1J3A7 ...Calculating statuses' bp per each bin for BS_SE82EQF8 ...Calculating statuses' bp per each bin for BS_SFH8TGKT ...Calculating statuses' bp per each bin for BS_SFZ3A07S ...Calculating statuses' bp per each bin for BS_SK4H5MJQ ...Calculating statuses' bp per each bin for BS_SK949N33 ...Calculating statuses' bp per each bin for BS_SKRY0BJ4 ...Calculating statuses' bp per each bin for BS_SMCPASF0 ...Calculating statuses' bp per each bin for BS_SNC5WMVN ...Calculating statuses' bp per each bin for BS_SQSH5G7A ...Calculating statuses' bp per each bin for BS_SSYNTXDC ...Calculating statuses' bp per each bin for BS_ST3Z2B9B ...Calculating statuses' bp per each bin for BS_STNH7YSX ...Calculating statuses' bp per each bin for BS_SYGKQD8M ...Calculating statuses' bp per each bin for BS_SYWRCE3Z ...Calculating statuses' bp per each bin for BS_T1SNKF50 ...Calculating statuses' bp per each bin for BS_T54TAE3R ...Calculating statuses' bp per each bin for BS_T5HPW4Z0 ...Calculating statuses' bp per each bin for BS_T5T21G7N ...Calculating statuses' bp per each bin for BS_T7WMJ08W ...Calculating statuses' bp per each bin for BS_T8C13KNH ...Calculating statuses' bp per each bin for BS_T8TCEBB1 ...Calculating statuses' bp per each bin for BS_T97EP8H0 ...Calculating statuses' bp per each bin for BS_T99M7842 ...Calculating statuses' bp per each bin for BS_T9T3R2J0 ...Calculating statuses' bp per each bin for BS_TE8QFF7T ...Calculating statuses' bp per each bin for BS_TF5TTEXH ...Calculating statuses' bp per each bin for BS_TF6H437Q ...Calculating statuses' bp per each bin for BS_THH0CE5A ...Calculating statuses' bp per each bin for BS_TNQKJZCP ...Calculating statuses' bp per each bin for BS_TPX7YY57 ...Calculating statuses' bp per each bin for BS_TQ0J7WJQ ...Calculating statuses' bp per each bin for BS_TRGECFNP ...Calculating statuses' bp per each bin for BS_TX2WGF8K ...Calculating statuses' bp per each bin for BS_TX8C5VAJ ...Calculating statuses' bp per each bin for BS_TZPJHKQR ...Calculating statuses' bp per each bin for BS_V2MDX7HG ...Calculating statuses' bp per each bin for BS_V3SC6NWH ...Calculating statuses' bp per each bin for BS_V3ZNXRHK ...Calculating statuses' bp per each bin for BS_V6HKXX1Z ...Calculating statuses' bp per each bin for BS_V79WFJ6N ...Calculating statuses' bp per each bin for BS_V8G6BZRR ...Calculating statuses' bp per each bin for BS_VA07453Z ...Calculating statuses' bp per each bin for BS_VDSR4X6N ...Calculating statuses' bp per each bin for BS_VF099E8S ...Calculating statuses' bp per each bin for BS_VHD4ARSG ...Calculating statuses' bp per each bin for BS_VK46N40F ...Calculating statuses' bp per each bin for BS_VKH9KYDB ...Calculating statuses' bp per each bin for BS_VMZMGB7B ...Calculating statuses' bp per each bin for BS_VPJ9CARN ...Calculating statuses' bp per each bin for BS_VQC2PE2B ...Calculating statuses' bp per each bin for BS_VQEPFE1P ...Calculating statuses' bp per each bin for BS_VTTTQYQA ...Calculating statuses' bp per each bin for BS_VVD8TJNE ...Calculating statuses' bp per each bin for BS_VW4XN9Y7 ...Calculating statuses' bp per each bin for BS_VW96PFKQ ...Calculating statuses' bp per each bin for BS_VXDGXQKZ ...Calculating statuses' bp per each bin for BS_VXNT1TDB ...Calculating statuses' bp per each bin for BS_W23ERYW3 ...Calculating statuses' bp per each bin for BS_W2QCHQ7E ...Calculating statuses' bp per each bin for BS_W3R2DQ86 ...Calculating statuses' bp per each bin for BS_W50WEJE7 ...Calculating statuses' bp per each bin for BS_W5P7SPDH ...Calculating statuses' bp per each bin for BS_W6B48DD5 ...Calculating statuses' bp per each bin for BS_W74D279Q ...Calculating statuses' bp per each bin for BS_W8DA9AXT ...Calculating statuses' bp per each bin for BS_W9FW1PAE ...Calculating statuses' bp per each bin for BS_W9P0PFDH ...Calculating statuses' bp per each bin for BS_WB8P21ED ...Calculating statuses' bp per each bin for BS_WCNPE0AC ...Calculating statuses' bp per each bin for BS_WD37M8SD ...Calculating statuses' bp per each bin for BS_WD7AQV83 ...Calculating statuses' bp per each bin for BS_WDTT7PG2 ...Calculating statuses' bp per each bin for BS_WEGGXVYA ...Calculating statuses' bp per each bin for BS_WJ7GMQTF ...Calculating statuses' bp per each bin for BS_WJB33V17 ...Calculating statuses' bp per each bin for BS_WJF6ZXWN ...Calculating statuses' bp per each bin for BS_WMEYJGDC ...Calculating statuses' bp per each bin for BS_WP2D11SK ...Calculating statuses' bp per each bin for BS_WP9J88EB ...Calculating statuses' bp per each bin for BS_WPCBK1EG ...Calculating statuses' bp per each bin for BS_WQ19Z6MD ...Calculating statuses' bp per each bin for BS_WRJNMEX2 ...Calculating statuses' bp per each bin for BS_WSFB4SA2 ...Calculating statuses' bp per each bin for BS_WXKCR17Q ...Calculating statuses' bp per each bin for BS_WY5KYSHJ ...Calculating statuses' bp per each bin for BS_WYF8NGSD ...Calculating statuses' bp per each bin for BS_X1DSETVE ...Calculating statuses' bp per each bin for BS_X1R0N8QS ...Calculating statuses' bp per each bin for BS_X5VN0FW0 ...Calculating statuses' bp per each bin for BS_X6361T19 ...Calculating statuses' bp per each bin for BS_X7QJCVJB ...Calculating statuses' bp per each bin for BS_X9PJ715R ...Calculating statuses' bp per each bin for BS_XBBBE8V2 ...Calculating statuses' bp per each bin for BS_XBE002WV ...Calculating statuses' bp per each bin for BS_XCFNEZ8E ...Calculating statuses' bp per each bin for BS_XCNGQ25H ...Calculating statuses' bp per each bin for BS_XDY9BYRK ...Calculating statuses' bp per each bin for BS_XEVMEYFS ...Calculating statuses' bp per each bin for BS_XGMXJKFX ...Calculating statuses' bp per each bin for BS_XGVSX8BX ...Calculating statuses' bp per each bin for BS_XH5JP8DB ...Calculating statuses' bp per each bin for BS_XHAQXDR0 ...Calculating statuses' bp per each bin for BS_XJZQA3AP ...Calculating statuses' bp per each bin for BS_XKE77TR4 ...Calculating statuses' bp per each bin for BS_XKM2EAT4 ...Calculating statuses' bp per each bin for BS_XKP747AS ...Calculating statuses' bp per each bin for BS_XKXDH0YJ ...Calculating statuses' bp per each bin for BS_XM302FF4 ...Calculating statuses' bp per each bin for BS_XMNBSBGN ...Calculating statuses' bp per each bin for BS_XMP9XNR9 ...Calculating statuses' bp per each bin for BS_XNTQBZAJ ...Calculating statuses' bp per each bin for BS_XPSD2NRJ ...Calculating statuses' bp per each bin for BS_XVBX7ZK1 ...Calculating statuses' bp per each bin for BS_XVNNWN2Z ...Calculating statuses' bp per each bin for BS_XVYRXN8J ...Calculating statuses' bp per each bin for BS_XW26N96W ...Calculating statuses' bp per each bin for BS_XWQTX5HH ...Calculating statuses' bp per each bin for BS_Y2C78JC4 ...Calculating statuses' bp per each bin for BS_Y3472DYJ ...Calculating statuses' bp per each bin for BS_Y3ETG5AE ...Calculating statuses' bp per each bin for BS_Y43VC4VG ...Calculating statuses' bp per each bin for BS_Y4ZX0YX7 ...Calculating statuses' bp per each bin for BS_Y5DA59M3 ...Calculating statuses' bp per each bin for BS_Y74XAFJX ...Calculating statuses' bp per each bin for BS_Y7CA8B5R ...Calculating statuses' bp per each bin for BS_Y9GHDZKS ...Calculating statuses' bp per each bin for BS_YB99Q6TN ...Calculating statuses' bp per each bin for BS_YBCQBCQC ...Calculating statuses' bp per each bin for BS_YETTZ1NC ...Calculating statuses' bp per each bin for BS_YEWX69X4 ...Calculating statuses' bp per each bin for BS_YGQS4660 ...Calculating statuses' bp per each bin for BS_YH4DYNQ2 ...Calculating statuses' bp per each bin for BS_YH60H2FS ...Calculating statuses' bp per each bin for BS_YHXMYDBN ...Calculating statuses' bp per each bin for BS_YKQP4SBX ...Calculating statuses' bp per each bin for BS_YMDPFK0V ...Calculating statuses' bp per each bin for BS_YPCTK927 ...Calculating statuses' bp per each bin for BS_YRC54VPP ...Calculating statuses' bp per each bin for BS_YVFFHJ1Q ...Calculating statuses' bp per each bin for BS_YXCVR90M ...Calculating statuses' bp per each bin for BS_YYAPSA5P ...Calculating statuses' bp per each bin for BS_YZ0P5E8H ...Calculating statuses' bp per each bin for BS_YZ2YNEJJ ...Calculating statuses' bp per each bin for BS_YZD4SSMA ...Calculating statuses' bp per each bin for BS_Z14Z2JMW ...Calculating statuses' bp per each bin for BS_Z1BBZ01F ...Calculating statuses' bp per each bin for BS_Z34PQDYT ...Calculating statuses' bp per each bin for BS_Z4Q8RH5M ...Calculating statuses' bp per each bin for BS_Z4S81HG1 ...Calculating statuses' bp per each bin for BS_Z5VJCFNP ...Calculating statuses' bp per each bin for BS_Z7R92XX6 ...Calculating statuses' bp per each bin for BS_Z9NCH2XA ...Calculating statuses' bp per each bin for BS_Z9PKZ4RT ...Calculating statuses' bp per each bin for BS_ZAH4FHV2 ...Calculating statuses' bp per each bin for BS_ZAT0XJE5 ...Calculating statuses' bp per each bin for BS_ZBWQG7VR ...Calculating statuses' bp per each bin for BS_ZCY4ECNV ...Calculating statuses' bp per each bin for BS_ZECKZTFC ...Calculating statuses' bp per each bin for BS_ZH6PC8KM ...Calculating statuses' bp per each bin for BS_ZHJFX6C9 ...Calculating statuses' bp per each bin for BS_ZJSMC22A ...Calculating statuses' bp per each bin for BS_ZJT1REVB ...Calculating statuses' bp per each bin for BS_ZJZFT0AJ ...Calculating statuses' bp per each bin for BS_ZKFEMJZ8 ...Calculating statuses' bp per each bin for BS_ZKZACPVV ...Calculating statuses' bp per each bin for BS_ZP89F1JY ...Calculating statuses' bp per each bin for BS_ZPF89720 ...Calculating statuses' bp per each bin for BS_ZR75EKKX ...Calculating statuses' bp per each bin for BS_ZS1QRMXS ...Calculating statuses' bp per each bin for BS_ZSH09N84 ...Calculating statuses' bp per each bin for BS_ZV21J6YW ...Parsed with column specification:
+cols(
+ .default = col_character()
+)Calculating statuses' bp per each bin for BS_ZZJF26C4 ...See spec(...) for full column specifications.sessionInfo()R version 3.6.0 (2019-04-26)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Debian GNU/Linux 9 (stretch)
@@ -6316,13 +3580,13 @@ Session Info
[37] S4Vectors_0.24.4 jsonlite_1.6 XVector_0.26.0
[40] rjson_0.2.20 ggplot2_3.2.0 hms_0.4.2
[43] png_0.1-7 digest_0.6.20 stringi_1.4.3
-[46] dplyr_0.8.3 clue_0.3-57 GenomicRanges_1.38.0
-[49] grid_3.6.0 tools_3.6.0 bitops_1.0-6
+[46] dplyr_0.8.3 GenomicRanges_1.38.0 grid_3.6.0
+[49] clue_0.3-57 tools_3.6.0 bitops_1.0-6
[52] magrittr_1.5 RCurl_1.95-4.12 lazyeval_0.2.2
[55] tibble_2.1.3 cluster_2.1.0 crayon_1.3.4
-[58] tidyr_0.8.3 pkgconfig_2.0.2 ellipsis_0.2.0.1
-[61] data.table_1.12.2 assertthat_0.2.1 rmarkdown_1.13
-[64] R6_2.4.0 compiler_3.6.0 # Magrittr pipe
`%>%` <- dplyr::`%>%`Registered S3 method overwritten by 'dplyr':
- method from
- print.rowwise_df # Path to input directory
input_dir <- file.path("..", "..", "data")
scratch_dir <- file.path("..", "..", "scratch")
@@ -2983,7 +2978,7 @@ Directories and Files
Unzip and set up GISTIC results folder.
-
+
gistic_zip <- file.path(input_dir, "pbta-cnv-consensus-gistic.zip")
gistic_dir <- file.path(input_dir, "pbta-cnv-consensus-gistic")
@@ -3000,39 +2995,31 @@ Directories and Files
Read in data
-
+
# Read in metadata
metadata <-
- readr::read_tsv(file.path(input_dir, "pbta-histologies.tsv"))
+ readr::read_tsv(file.path(input_dir, "pbta-histologies.tsv"), guess_max = 10000)
-
+
Parsed with column specification:
cols(
.default = col_character(),
- OS_days = [32mcol_double()[39m,
- age_last_update_days = [32mcol_double()[39m,
- normal_fraction = [32mcol_double()[39m,
- tumor_fraction = [32mcol_double()[39m,
- tumor_ploidy = [32mcol_double()[39m,
- molecular_subtype = [33mcol_logical()[39m
-)
-See spec(...) for full column specifications.
-493 parsing failures.
- row col expected actual file
-2334 molecular_subtype 1/0/T/F/TRUE/FALSE Group3 '../../data/pbta-histologies.tsv'
-2335 molecular_subtype 1/0/T/F/TRUE/FALSE Group4 '../../data/pbta-histologies.tsv'
-2336 molecular_subtype 1/0/T/F/TRUE/FALSE Group3 '../../data/pbta-histologies.tsv'
-2337 molecular_subtype 1/0/T/F/TRUE/FALSE Group3 '../../data/pbta-histologies.tsv'
-2338 molecular_subtype 1/0/T/F/TRUE/FALSE Group3 '../../data/pbta-histologies.tsv'
-.... ................. .................. ...... .................................
-See problems(...) for more details.
-
+ OS_days = col_double(),
+ age_last_update_days = col_double(),
+ normal_fraction = col_double(),
+ tumor_fraction = col_double(),
+ tumor_ploidy = col_double()
+)See spec(...) for full column specifications.Read in the GISTIC scores file and format the data.
- +gistic_scores <- data.table::fread(file.path(gistic_dir,
"scores.gistic"),
data.table = FALSE) %>%
@@ -3052,40 +3039,55 @@ Read in data
)
)Registered S3 method overwritten by 'data.table':
- method from
- print.data.table
-Error in data.table::fread(file.path(gistic_dir, "scores.gistic"), data.table = FALSE) :
- File '../../data/pbta-cnv-consensus-gistic/scores.gistic' does not exist or is non-readable. getwd()=='/home/rstudio/kitematic/OpenPBTA-analysis/analyses/cnv-chrom-plot'Read in the seg file so we can plot data by histology.
rr gistic_plot <- ggbio::autoplot(gistic_ranges, ggplot2::aes(y = score, fill = mcols.Type), geom = , scales = _x, space = _x
-) + ggplot2::theme_classic() + ggplot2::theme(axis.text.x = ggplot2::element_text(size = 3, angle = 45, hjust = 1)) + colorblindr::scale_fill_OkabeIto(name = ) + ggplot2::ylab(scores)
# TODO: update file path when the consensus file has been added to the data release.
+# Read in the segment copy number data
+seg_data <- data.table::fread(file.path(
+ "..",
+ "copy_number_consensus_call",
+ "results",
+ "pbta-cnv-consensus.seg.gz"),
+ data.table = FALSE
+) %>%
+ # Join the histology column to this data
+ dplyr::inner_join(dplyr::select(
+ metadata,
+ "Kids_First_Biospecimen_ID",
+ "short_histology"
+ ),
+ by = c("ID" = "Kids_First_Biospecimen_ID")
+ ) %>%
+ # Reformat the chromosome variable to drop the "chr"
+ dplyr::mutate(chrom = factor(gsub("chr", "", chrom),
+ levels = c(1:22, "X", "Y")))
Registered S3 method overwritten by 'R.oo':
+ method from
+ throw.default R.methodsS3rr ggplot2::ggsave( plot = gistic_plot@ggplot, filename = file.path(plots_dir, _plot.png), width = 10, height = 3, units =
-)
gistic_ranges <- GenomicRanges::GRanges(
+ seqnames = gistic_scores$Chromosome,
+ ranges = IRanges::IRanges(
+ start = gistic_scores$Start,
+ end = gistic_scores$End
+ ),
+ score = gistic_scores$gscore,
+ mcols = gistic_scores
+)rr # Make a plot for each histology group lapply(unique(seg_data$short_histology), function(histology_group) {
# How many samples are of this histology label? sample_num <- seg_data %>% dplyr::filter(short_histology == histology_group) %>% dplyr::pull(ID) %>% unique() %>% length()
-# Isolate the data to the histology group hist_df <- seg_data %>% dplyr::filter(short_histology == histology_group) %>% dplyr::group_by(short_histology, chrom, loc.start, loc.end) %>% dplyr::summarize(hist.mean = mean(seg.mean)) %>% # Make Del/Amp variable dplyr::mutate(Type = dplyr::case_when( hist.mean < 0 ~ , hist.mean > 0 ~ , hist.mean == 0 ~
-))
# Turn into a GRanges for easier mapping hist_ranges <- GenomicRanges::GRanges( seqnames = hist_df\(chrom, ranges = IRanges::IRanges( start = hist_df\)loc.start, end = hist_df\(loc.end ), score = hist_df\)hist.mean, mcols = hist_df$Type ) # Map this on a plot hist_plot <- ggbio::autoplot(hist_ranges, ggplot2::aes(y = score, fill = mcols), geom = , scales = _x, space = _x
-) + ggplot2::theme_classic() + ggplot2::theme(axis.text.x = ggplot2::element_text(size = 3, angle = 45, hjust = 1)) + colorblindr::scale_fill_OkabeIto(name = ) + ggplot2::labs( title = histology_group, y = segment mean, subtitle = paste(=, sample_num) )
# Save the plot to a png # Have to use @ggplot to extract the portion of the ggbio object with the plot ggplot2::ggsave( plot = hist_plot@ggplot, filename = file.path(plots_dir, paste0(histology_group, _plot.png)), width = 10, height = 3, units =
-) })
gistic_plot <- ggbio::autoplot(gistic_ranges, ggplot2::aes(y = score, fill = mcols.Type),
+ geom = "bar", scales = "free_x", space = "free_x"
+) +
+ ggplot2::theme_classic() +
+ ggplot2::theme(axis.text.x = ggplot2::element_text(size = 3, angle = 45, hjust = 1)) +
+ colorblindr::scale_fill_OkabeIto(name = "Type") +
+ ggplot2::ylab("G-scores")Registered S3 method overwritten by 'GGally':
+ method from
+ +.gg ggplot2# Print out plot here
+gistic_plot@ggplot
rr sessionInfo()
ggplot2::ggsave(
+ plot = gistic_plot@ggplot,
+ filename = file.path(plots_dir, "gistic_plot.png"),
+ width = 10,
+ height = 3,
+ units = "in"
+)R version 3.6.0 (2019-04-26)
-Platform: x86_64-pc-linux-gnu (64-bit)
-Running under: Debian GNU/Linux 9 (stretch)
-
-Matrix products: default
-BLAS/LAPACK: /usr/lib/libopenblasp-r0.2.19.so
-
-locale:
- [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8
- [4] LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8 LC_MESSAGES=C
- [7] LC_PAPER=en_US.UTF-8 LC_NAME=C LC_ADDRESS=C
-[10] LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
-
-attached base packages:
-[1] stats graphics grDevices utils datasets methods base
-
-loaded via a namespace (and not attached):
- [1] colorspace_1.4-1 rjson_0.2.20 rprojroot_1.3-2
- [4] biovizBase_1.34.1 circlize_0.4.8 htmlTable_1.13.1
- [7] XVector_0.26.0 GenomicRanges_1.38.0 GlobalOptions_0.1.1
- [10] base64enc_0.1-3 dichromat_2.0-0 clue_0.3-57
- [13] rstudioapi_0.10 bit64_0.9-7 AnnotationDbi_1.48.0
- [16] R.methodsS3_1.7.1 splines_3.6.0 ggbio_1.34.0
- [19] knitr_1.23 Formula_1.2-3 jsonlite_1.6
- [22] colorblindr_0.1.0 Rsamtools_2.2.1 cluster_2.1.0
- [25] dbplyr_1.4.2 R.oo_1.22.0 png_0.1-7
- [28] graph_1.62.0 BiocManager_1.30.4 readr_1.3.1
- [31] compiler_3.6.0 httr_1.4.0 backports_1.1.4
- [34] assertthat_0.2.1 Matrix_1.2-17 lazyeval_0.2.2
- [37] cli_1.1.0 acepack_1.4.1 htmltools_0.3.6
- [40] prettyunits_1.0.2 tools_3.6.0 gtable_0.3.0
- [43] glue_1.3.1 GenomeInfoDbData_1.2.2 reshape2_1.4.3
- [46] dplyr_0.8.3 rappdirs_0.3.1 Rcpp_1.0.1
- [49] Biobase_2.46.0 styler_1.1.1 Biostrings_2.54.0
- [52] rtracklayer_1.46.0 xfun_0.8 stringr_1.4.0
- [55] ensembldb_2.10.2 XML_3.98-1.20 zlibbioc_1.32.0
- [58] scales_1.0.0 BSgenome_1.54.0 VariantAnnotation_1.32.0
- [61] hms_0.4.2 ProtGenerics_1.18.0 RBGL_1.62.1
- [64] parallel_3.6.0 SummarizedExperiment_1.16.1 AnnotationFilter_1.10.0
- [67] rematch2_2.0.1 RColorBrewer_1.1-2 ComplexHeatmap_2.2.0
- [70] yaml_2.2.0 curl_3.3 memoise_1.1.0
- [73] gridExtra_2.3 ggplot2_3.2.0 biomaRt_2.42.0
- [76] rpart_4.1-15 reshape_0.8.8 latticeExtra_0.6-28
- [79] stringi_1.4.3 RSQLite_2.1.1 S4Vectors_0.24.1
- [82] checkmate_1.9.4 GenomicFeatures_1.38.0 BiocGenerics_0.32.0
- [85] BiocParallel_1.20.1 shape_1.4.4 GenomeInfoDb_1.22.0
- [88] rlang_0.4.0 pkgconfig_2.0.2 matrixStats_0.55.0
- [91] bitops_1.0-6 evaluate_0.14 lattice_0.20-38
- [94] purrr_0.3.2 GenomicAlignments_1.22.1 htmlwidgets_1.3
- [97] labeling_0.3 cowplot_0.9.4 bit_1.1-14
-[100] tidyselect_0.2.5 GGally_1.4.0 plyr_1.8.4
-[103] magrittr_1.5 R6_2.4.0 IRanges_2.20.1
-[106] Hmisc_4.2-0 DelayedArray_0.12.2 DBI_1.0.0
-[109] pillar_1.4.2 foreign_0.8-71 withr_2.1.2
-[112] survival_2.44-1.1 RCurl_1.95-4.12 nnet_7.3-12
-[115] tibble_2.1.3 crayon_1.3.4 OrganismDbi_1.28.0
-[118] BiocFileCache_1.10.2 rmarkdown_1.13 GetoptLong_0.1.7
-[121] progress_1.2.2 grid_3.6.0 data.table_1.12.2
-[124] blob_1.1.1 digest_0.6.20 R.utils_2.9.0
-[127] openssl_1.4 stats4_3.6.0 munsell_0.5.0
-[130] askpass_1.1 # Make a plot for each histology group
+lapply(unique(seg_data$short_histology), function(histology_group) {
+
+ # How many samples are of this histology label?
+ sample_num <- seg_data %>%
+ dplyr::filter(short_histology == histology_group) %>%
+ dplyr::pull(ID) %>%
+ unique() %>%
+ length()
+
+ # Isolate the data to the histology group
+ hist_df <- seg_data %>%
+ dplyr::filter(short_histology == histology_group) %>%
+ dplyr::group_by(short_histology, chrom, loc.start, loc.end) %>%
+ dplyr::summarize(hist.mean = mean(seg.mean)) %>%
+ # Make Del/Amp variable
+ dplyr::mutate(Type = dplyr::case_when(
+ hist.mean < 0 ~ "Del",
+ hist.mean > 0 ~ "Amp",
+ hist.mean == 0 ~ "Neutral"
+ ))
+
+ # Turn into a GRanges for easier mapping
+ hist_ranges <- GenomicRanges::GRanges(
+ seqnames = hist_df$chrom,
+ ranges = IRanges::IRanges(
+ start = hist_df$loc.start,
+ end = hist_df$loc.end
+ ),
+ score = hist_df$hist.mean,
+ mcols = hist_df$Type
+ )
+ # Map this on a plot
+ hist_plot <- ggbio::autoplot(hist_ranges,
+ ggplot2::aes(y = score, fill = mcols),
+ geom = "bar", scales = "free_x", space = "free_x"
+ ) +
+ ggplot2::theme_classic() +
+ ggplot2::theme(axis.text.x = ggplot2::element_text(size = 3, angle = 45, hjust = 1)) +
+ colorblindr::scale_fill_OkabeIto(name = "Type") +
+ ggplot2::labs(
+ title = histology_group,
+ y = "Average segment mean",
+ subtitle = paste("n =", sample_num)
+ )
+
+ # Save the plot to a png
+ # Have to use @ggplot to extract the portion of the ggbio object with the plot
+ ggplot2::ggsave(
+ plot = hist_plot@ggplot,
+ filename = file.path(plots_dir, paste0(histology_group, "_plot.png")),
+ width = 10,
+ height = 3,
+ units = "in"
+ )
+})Warning: Removed 1053 rows containing missing values (geom_rect).Warning: Removed 3572 rows containing missing values (geom_rect).Warning: Removed 1166 rows containing missing values (geom_rect).Warning: Removed 183 rows containing missing values (geom_rect).Warning: Removed 398 rows containing missing values (geom_rect).Warning: Removed 1732 rows containing missing values (geom_rect).Warning: Removed 234 rows containing missing values (geom_rect).Warning: Removed 458 rows containing missing values (geom_rect).Warning: Removed 350 rows containing missing values (geom_rect).Warning: Removed 230 rows containing missing values (geom_rect).Warning: Removed 383 rows containing missing values (geom_rect).Warning in max(xlim): no non-missing arguments to max; returning -InfWarning: Removed 181 rows containing missing values (geom_rect).Warning: Removed 238 rows containing missing values (geom_rect).Warning: Removed 275 rows containing missing values (geom_rect).Warning: Removed 219 rows containing missing values (geom_rect).Warning: Removed 185 rows containing missing values (geom_rect).Warning: Removed 582 rows containing missing values (geom_rect).Warning: Removed 182 rows containing missing values (geom_rect).Warning: Removed 233 rows containing missing values (geom_rect).Warning: Removed 426 rows containing missing values (geom_rect).Warning: Removed 209 rows containing missing values (geom_rect).Warning: Removed 183 rows containing missing values (geom_rect).Warning: Removed 187 rows containing missing values (geom_rect).Warning: Removed 183 rows containing missing values (geom_rect).Warning: Removed 311 rows containing missing values (geom_rect).Warning: Removed 275 rows containing missing values (geom_rect).Warning in max(xlim): no non-missing arguments to max; returning -InfWarning: Removed 181 rows containing missing values (geom_rect).Warning: Removed 189 rows containing missing values (geom_rect).Warning: Removed 266 rows containing missing values (geom_rect).Warning: Removed 245 rows containing missing values (geom_rect).Warning: Removed 189 rows containing missing values (geom_rect).Warning: Removed 186 rows containing missing values (geom_rect).Warning: Removed 175 rows containing missing values (geom_rect).Warning: Removed 200 rows containing missing values (geom_rect).Warning: Removed 227 rows containing missing values (geom_rect).Warning: Removed 242 rows containing missing values (geom_rect).Warning: Removed 243 rows containing missing values (geom_rect).Warning: Removed 190 rows containing missing values (geom_rect).Warning: Removed 196 rows containing missing values (geom_rect).Warning: Removed 185 rows containing missing values (geom_rect).Warning: Removed 200 rows containing missing values (geom_rect).Warning: Removed 189 rows containing missing values (geom_rect).Warning in max(xlim): no non-missing arguments to max; returning -InfWarning: Removed 181 rows containing missing values (geom_rect).sessionInfo()R version 3.6.0 (2019-04-26)
+Platform: x86_64-pc-linux-gnu (64-bit)
+Running under: Debian GNU/Linux 9 (stretch)
+
+Matrix products: default
+BLAS/LAPACK: /usr/lib/libopenblasp-r0.2.19.so
+
+locale:
+ [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
+ [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
+ [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=C
+ [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
+ [9] LC_ADDRESS=C LC_TELEPHONE=C
+[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
+
+attached base packages:
+[1] stats graphics grDevices utils datasets methods base
+
+loaded via a namespace (and not attached):
+ [1] ProtGenerics_1.18.0 bitops_1.0-6
+ [3] matrixStats_0.54.0 bit64_0.9-7
+ [5] RColorBrewer_1.1-2 progress_1.2.2
+ [7] httr_1.4.0 GenomeInfoDb_1.22.1
+ [9] tools_3.6.0 backports_1.1.4
+ [11] R6_2.4.0 rpart_4.1-15
+ [13] Hmisc_4.2-0 DBI_1.0.0
+ [15] lazyeval_0.2.2 BiocGenerics_0.32.0
+ [17] colorspace_1.4-1 nnet_7.3-12
+ [19] tidyselect_0.2.5 gridExtra_2.3
+ [21] prettyunits_1.0.2 GGally_1.4.0
+ [23] curl_3.3 bit_1.1-14
+ [25] compiler_3.6.0 graph_1.64.0
+ [27] Biobase_2.46.0 htmlTable_1.13.1
+ [29] DelayedArray_0.12.3 labeling_0.3
+ [31] rtracklayer_1.46.0 ggbio_1.34.0
+ [33] scales_1.0.0 checkmate_1.9.4
+ [35] readr_1.3.1 RBGL_1.62.1
+ [37] askpass_1.1 rappdirs_0.3.1
+ [39] stringr_1.4.0 digest_0.6.20
+ [41] Rsamtools_2.2.3 foreign_0.8-71
+ [43] rmarkdown_1.13 R.utils_2.9.0
+ [45] XVector_0.26.0 base64enc_0.1-3
+ [47] dichromat_2.0-0 pkgconfig_2.0.2
+ [49] htmltools_0.3.6 ensembldb_2.10.2
+ [51] BSgenome_1.54.0 dbplyr_1.4.2
+ [53] htmlwidgets_1.3 rlang_0.4.0
+ [55] rstudioapi_0.10 RSQLite_2.1.1
+ [57] jsonlite_1.6 BiocParallel_1.20.1
+ [59] acepack_1.4.1 dplyr_0.8.3
+ [61] R.oo_1.22.0 VariantAnnotation_1.32.0
+ [63] RCurl_1.95-4.12 magrittr_1.5
+ [65] GenomeInfoDbData_1.2.2 Formula_1.2-3
+ [67] Matrix_1.2-17 Rcpp_1.0.1
+ [69] munsell_0.5.0 S4Vectors_0.24.4
+ [71] R.methodsS3_1.7.1 stringi_1.4.3
+ [73] yaml_2.2.0 SummarizedExperiment_1.16.1
+ [75] zlibbioc_1.32.0 BiocFileCache_1.10.2
+ [77] plyr_1.8.4 grid_3.6.0
+ [79] blob_1.1.1 parallel_3.6.0
+ [81] crayon_1.3.4 lattice_0.20-38
+ [83] Biostrings_2.54.0 splines_3.6.0
+ [85] GenomicFeatures_1.38.2 hms_0.4.2
+ [87] knitr_1.23 pillar_1.4.2
+ [89] GenomicRanges_1.38.0 reshape2_1.4.3
+ [91] biomaRt_2.42.1 stats4_3.6.0
+ [93] XML_3.98-1.20 glue_1.3.1
+ [95] evaluate_0.14 biovizBase_1.34.1
+ [97] latticeExtra_0.6-28 BiocManager_1.30.4
+ [99] data.table_1.12.2 colorblindr_0.1.0
+[101] gtable_0.3.0 openssl_1.4
+[103] purrr_0.3.2 reshape_0.8.8
+[105] assertthat_0.2.1 ggplot2_3.2.0
+[107] xfun_0.8 AnnotationFilter_1.10.0
+[109] survival_2.44-1.1 OrganismDbi_1.28.0
+[111] tibble_2.1.3 GenomicAlignments_1.22.1
+[113] AnnotationDbi_1.48.0 memoise_1.1.0
+[115] IRanges_2.20.2 cluster_2.1.0 seg.mean agree with status?
seg.mean agree with status?# Here, we arbitrarily define recurrence to be greater than 3 based on
# the plot above and on a similar decision made in
# `analyses/fusion_filtering/06-recurrent-fusions-per-histology.R`
@@ -2971,9 +2971,23 @@ Cutoffs
Libraries and functions
-
+
library(tidyverse)
+
+── Attaching packages ────────────────────────────────── tidyverse 1.2.1 ──
+
+
+✔ ggplot2 3.2.0 ✔ purrr 0.3.2
+✔ tibble 2.1.3 ✔ dplyr 0.8.3
+✔ tidyr 0.8.3 ✔ stringr 1.4.0
+✔ readr 1.3.1 ✔ forcats 0.4.0
+
+
+── Conflicts ───────────────────────────────────── tidyverse_conflicts() ──
+✖ dplyr::filter() masks stats::filter()
+✖ dplyr::lag() masks stats::lag()
+
@@ -2981,7 +2995,7 @@ Libraries and functions
Files and directories
-
+
results_dir <- "results"
plots_dir <- "plots"
@@ -2994,36 +3008,38 @@ Read in files
Read in focal CN dominant status calls prepared in 05-define-most-focal-cn-units.Rmd.
-
+
focal_cn_df <-
read_tsv(file.path(results_dir, "consensus_seg_most_focal_cn_status.tsv.gz")) %>%
# Filter out the `uncallable` statuses
filter(status != "uncallable")
-
+
Parsed with column specification:
cols(
- Kids_First_Biospecimen_ID = [31mcol_character()[39m,
- status = [31mcol_character()[39m,
- region = [31mcol_character()[39m,
- region_type = [31mcol_character()[39m
+ Kids_First_Biospecimen_ID = col_character(),
+ status = col_character(),
+ region = col_character(),
+ region_type = col_character()
)
-
-
+
+
metadata <- read_tsv(file.path("..", "..", "data", "pbta-histologies.tsv"), guess_max = 10000)
-
+
Parsed with column specification:
cols(
.default = col_character(),
- OS_days = [32mcol_double()[39m,
- age_last_update_days = [32mcol_double()[39m,
- normal_fraction = [32mcol_double()[39m,
- tumor_fraction = [32mcol_double()[39m,
- tumor_ploidy = [32mcol_double()[39m
-)
-See spec(...) for full column specifications.
-
+ OS_days = col_double(),
+ age_last_update_days = col_double(),
+ normal_fraction = col_double(),
+ tumor_fraction = col_double(),
+ tumor_ploidy = col_double()
+)See spec(...) for full column specifications.# Tally the number of samples that have the same status call by region
final_df_tally <- focal_cn_df %>%
group_by(status, region, region_type) %>%
@@ -3046,13 +3062,11 @@ Tally the sample counts of each unique status call by region.
final_df_tally %>%
arrange(desc(count))# Tally the number of samples that have the same status call by region
final_df_tally_histology <- focal_cn_df %>%
# Join the metadata by biospecimen IDs to get short histology data
@@ -3077,13 +3091,11 @@ Tally the sample counts of each unique status call by region and short histo
final_df_tally_histology %>%
arrange(desc(count))
# Plot the distribution of counts for all histologies
ggplot(final_df_tally, aes(x = count)) +
geom_histogram(bins = 100) +
@@ -3109,22 +3121,23 @@ Plot the distribution of counts
angle = 90,
vjust = 1.2,
size = 4
- )
-
-# Save plot
-ggsave(file.path(plots_dir, "all_histology_cn_counts.png"), plot = last_plot())Saving 6.67 x 4.11 in image
# Save plot
+ggsave(file.path(plots_dir, "all_histology_cn_counts.png"), plot = last_plot())Saving 7 x 5 in image# Plot the distribution of counts for histologies separated
ggplot(final_df_tally_histology, aes(x = count)) +
geom_histogram(bins = 100) +
@@ -3142,17 +3155,18 @@ Plot the distribution of counts
angle = 90,
vjust = 1.2,
size = 4
- )
-
-# Save plot
-ggsave(file.path(plots_dir, "separate_histology_cn_counts.png"), plot = last_plot())Saving 6.67 x 4.11 in image
# Save plot
+ggsave(file.path(plots_dir, "separate_histology_cn_counts.png"), plot = last_plot())Saving 7 x 5 in imageThe recurrence_cutoff defined as 3 appears to be reasonable as seen in the distribution of counts plot above.
# Create a data frame of only observations that have counts greater than the `recurrence_cutoff`
recurrent_df <- final_df_tally %>%
filter(count > recurrence_cutoff)
@@ -3171,13 +3185,11 @@ All histologies represented
# Display table
recurrent_df# Create a data frame of only observations that have counts greater than the `recurrence_cutoff`
recurrent_df_histology <- final_df_tally_histology %>%
filter(count > recurrence_cutoff)
@@ -3193,13 +3205,11 @@ Separated by histology
# Display table
recurrent_df_histologyfinal_df_recurrent <- recurrent_df %>%
bind_rows(recurrent_df_histology) %>%
arrange(region) %>%
@@ -3216,13 +3226,11 @@ Combine all histology recurrent data with the data separated by histology
recurrent calls
-
+
# Save list of genes that have `recurrent` status calls for use in downstream
# analyses
top_recurrent_genes <- final_df_recurrent %>%
@@ -3259,7 +3267,7 @@ Save list of genes that have recurrent calls
Write final data frame to file
-
+
# Write the recurrent focal CN calls by region and histology data frame to
# file
write_tsv(final_df_recurrent, file.path(results_dir, "consensus_seg_recurrent_focal_cn_units.tsv"))
@@ -3272,10 +3280,10 @@ Write final data frame to file
Session Info
-
+
sessionInfo()
-
+
R version 3.6.0 (2019-04-26)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Debian GNU/Linux 9 (stretch)
@@ -3284,17 +3292,20 @@ Session Info
BLAS/LAPACK: /usr/lib/libopenblasp-r0.2.19.so
locale:
- [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8
- [4] LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8 LC_MESSAGES=C
- [7] LC_PAPER=en_US.UTF-8 LC_NAME=C LC_ADDRESS=C
-[10] LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
+ [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
+ [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
+ [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=C
+ [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
+ [9] LC_ADDRESS=C LC_TELEPHONE=C
+[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
-[1] forcats_0.4.0 stringr_1.4.0 dplyr_0.8.3 purrr_0.3.2 readr_1.3.1
-[6] tidyr_0.8.3 tibble_2.1.3 ggplot2_3.2.0 tidyverse_1.2.1
+[1] forcats_0.4.0 stringr_1.4.0 dplyr_0.8.3 purrr_0.3.2
+[5] readr_1.3.1 tidyr_0.8.3 tibble_2.1.3 ggplot2_3.2.0
+[9] tidyverse_1.2.1
loaded via a namespace (and not attached):
[1] Rcpp_1.0.1 cellranger_1.1.0 pillar_1.4.2 compiler_3.6.0
@@ -3306,9 +3317,9 @@ Session Info
[25] hms_0.4.2 generics_0.0.2 grid_3.6.0 tidyselect_0.2.5
[29] glue_1.3.1 R6_2.4.0 readxl_1.3.1 rmarkdown_1.13
[33] modelr_0.1.4 magrittr_1.5 backports_1.1.4 scales_1.0.0
-[37] htmltools_0.3.6 rsconnect_0.8.13 rvest_0.3.4 assertthat_0.2.1
-[41] colorspace_1.4-1 labeling_0.3 stringi_1.4.3 lazyeval_0.2.2
-[45] munsell_0.5.0 broom_0.5.2 crayon_1.3.4
+[37] htmltools_0.3.6 rvest_0.3.4 assertthat_0.2.1 colorspace_1.4-1
+[41] labeling_0.3 stringi_1.4.3 lazyeval_0.2.2 munsell_0.5.0
+[45] broom_0.5.2 crayon_1.3.4