Figures for Titus

kevlar-dev · Nov 15, 2018 · 151b4f1 · 151b4f1
1 parent 0e95536
commit 151b4f1
Show file tree

Hide file tree

Showing 3 changed files with 568 additions and 0 deletions.
diff --git a/notebook/bigsim/roc-data.json b/notebook/bigsim/roc-data.json
@@ -0,0 +1 @@
+{"kevlar": {"10": {"SNVs": [0, 5, 6, 6, 7, 9, 9, 10, 10, 11, 11, 12, 12, 12, 13, 13, 13, 13, 15, 17, 18, 19, 21, 21, 21, 21, 21, 22, 24, 27, 27, 27, 27, 28, 28, 28, 28, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32], "INDELs 1-10bp": [17, 19, 26, 27, 28, 31, 32, 33, 33, 33, 33, 33, 33, 34, 34, 34, 34, 34, 34, 34, 35], "INDELs 11-100bp": [10, 11, 11, 11, 14, 15, 16, 34, 53, 58, 62, 62, 62, 62, 62, 63, 63, 63, 63, 63, 63, 63], "INDELs 101-200bp": [10, 11, 12, 12, 15, 15, 18, 26, 29, 31, 33, 34, 34, 52], "INDELs 201-300bp": [7, 20, 21, 40], "INDELs 301-400bp": [38]}, "20": {"SNVs": [121, 130, 130, 131, 131, 139, 140, 141, 142, 143, 144, 147, 147, 147, 151, 151, 152, 152, 152, 152, 152, 153, 154, 154, 155, 155, 156, 157, 158, 160, 160, 161, 161, 162, 162, 163, 163, 163, 164, 165, 165, 167, 167, 167, 167, 167, 167, 168, 168, 170, 170, 170, 170, 170, 170, 170, 170, 171, 171, 171, 171, 171, 171, 172, 172, 172, 172, 172, 172, 172, 172, 172, 172, 172, 172, 172, 172, 173, 173, 173, 173, 173, 173, 173, 173, 173, 173, 173, 173, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 174, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 175, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 176, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177, 177], "INDELs 1-10bp": [20, 77, 175, 176, 176, 178, 178, 180, 180], "INDELs 11-100bp": [53, 61, 74, 119, 150, 168, 187, 204, 211, 211, 212, 212, 212, 212, 212, 212, 212], "INDELs 101-200bp": [30, 61, 102, 172, 200], "INDELs 201-300bp": [8, 29, 45, 45, 59, 90, 209], "INDELs 301-400bp": [85, 214]}, "30": {"SNVs": [193, 198, 210, 213, 216, 216, 217, 221, 221, 221, 222, 222, 222, 224, 224, 224, 226, 226, 226, 226, 226, 226, 226, 226, 227, 227, 229, 231, 232], "INDELs 1-10bp": [136, 221, 235, 240], "INDELs 11-100bp": [85, 244], "INDELs 101-200bp": [65, 252], "INDELs 201-300bp": [46, 105, 261], "INDELs 301-400bp": [79, 117, 132, 132, 253]}, "50": {"SNVs": [215, 216, 218, 223, 224, 225, 228, 228, 232, 234, 234, 234, 235, 235, 240, 240, 242], "INDELs 1-10bp": [93, 194, 220, 246, 246, 247, 247, 251], "INDELs 11-100bp": [68, 203, 215, 247], "INDELs 101-200bp": [73, 149, 228, 254], "INDELs 201-300bp": [52, 95, 262], "INDELs 301-400bp": [72, 74, 131, 131, 132, 199, 257]}}, "gatk": {"10": {"SNVs": [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1], "INDELs 1-10bp": [0, 0, 0, 0], "INDELs 11-100bp": [3], "INDELs 101-200bp": [1], "INDELs 201-300bp": [0], "INDELs 301-400bp": [0]}, "20": {"SNVs": [21, 29, 53, 54, 75, 96, 96, 103, 104, 108, 110, 123, 126, 132, 138, 138, 139, 139, 139, 141, 145, 145, 146, 147, 148, 149, 150, 151, 151, 152, 152, 152, 152, 152, 152, 154, 154, 154, 154, 154, 154, 154, 155, 155, 155, 155, 155, 156, 156, 157, 157, 157, 157, 158, 158, 158, 158, 158, 158, 158, 158, 158, 158, 158], "INDELs 1-10bp": [31, 39, 97, 98, 128, 128, 128, 128, 128, 129, 130, 132, 132, 135, 136, 138, 138, 141, 141, 142, 147, 151, 154, 155, 155, 155, 156, 156, 156, 157, 157, 158, 158, 160, 161, 161, 162, 164, 166, 166, 166, 172, 176, 177, 177, 177, 177, 178, 178, 178, 178, 178, 178, 178, 178, 178, 178, 178, 179], "INDELs 11-100bp": [52, 58, 59, 59, 73, 73, 80, 84, 84, 96, 99, 99, 101, 103, 104, 106, 111, 115, 118, 122, 126, 129, 129, 130, 137, 137, 141, 144, 153, 158, 165, 166, 167, 167, 167, 167, 167, 167, 167, 167, 168, 169, 169, 170, 171, 171, 172, 172, 172, 172, 172, 172, 172, 173, 173, 173, 173, 173, 174, 174, 174, 174, 175, 175, 175, 176, 177, 177, 177, 177, 177, 177, 178, 178, 178, 178, 179], "INDELs 101-200bp": [40, 89], "INDELs 201-300bp": [12], "INDELs 301-400bp": [0]}, "30": {"SNVs": [129, 139, 146, 156, 158, 172, 173, 176, 179, 182, 182, 188, 191, 200, 206, 208, 208, 210, 212, 216, 218, 221, 222, 231, 231, 233, 235, 235, 235, 235, 236, 236, 236, 236, 237, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 239, 239, 239, 239, 241, 241, 241, 241, 241, 242, 242, 242, 242, 242, 242, 242, 243, 244, 244, 245, 245, 245, 245, 245, 245, 245, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247], "INDELs 1-10bp": [2, 106, 107, 108, 118, 149, 149, 168, 178, 178, 195, 195, 201, 213, 218, 224, 230, 231, 232, 232, 237, 237, 238, 238, 238, 239, 239, 242, 244, 244, 244, 249, 250, 250, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 252, 252, 253, 253, 254, 255, 255, 255, 255, 255, 256, 256, 256, 256, 256, 256], "INDELs 11-100bp": [57, 72, 101, 104, 117, 118, 121, 121, 129, 129, 137, 138, 139, 139, 139, 141, 142, 150, 153, 154, 154, 154, 154, 157, 157, 157, 169, 174, 174, 174, 176, 176, 176, 181, 181, 182, 182, 182, 188, 189, 189, 190, 190, 190, 193, 193, 194, 196, 196, 196, 201, 201, 202, 202, 202, 204, 204, 205, 209, 211, 212, 212, 214, 214, 215, 215, 219, 219, 219, 220, 221, 221, 222, 223, 223, 224, 224, 224, 224, 224, 224, 225, 225, 225, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 227, 227, 227, 227, 227, 227, 228], "INDELs 101-200bp": [32, 35, 44, 69, 93, 93, 129, 129], "INDELs 201-300bp": [12], "INDELs 301-400bp": [0]}, "50": {"SNVs": [67, 222, 232, 233, 237, 239, 239, 240, 240, 240, 240, 243, 244, 244, 246, 246, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 250, 250, 250, 250, 250, 250, 250, 250, 250], "INDELs 1-10bp": [82, 174, 190, 194, 196, 202, 229, 230, 230, 238, 238, 251, 253, 253, 253, 253, 254, 254, 254, 255, 257, 257, 257, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 261], "INDELs 11-100bp": [55, 78, 112, 138, 142, 147, 147, 157, 157, 163, 163, 176, 176, 178, 185, 185, 185, 186, 201, 201, 205, 205, 206, 206, 206, 207, 211, 212, 212, 213, 215, 215, 215, 216, 216, 217, 217, 217, 217, 217, 218, 220, 222, 222, 222, 222, 222, 224, 224, 224, 224, 224, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 226, 227, 227, 227, 227, 227, 227, 228, 228, 230, 230, 230, 230, 230, 230, 230, 231, 231, 231, 232, 232, 232, 232, 232, 233, 233, 233, 234, 234, 234, 235, 235, 235, 235, 235, 235, 235, 235, 236, 236, 236, 238, 240], "INDELs 101-200bp": [35, 51, 103, 112, 120, 132], "INDELs 201-300bp": [21], "INDELs 301-400bp": [0]}}, "triodenovo": {"10": {"SNVs": [12, 14, 14, 22, 24, 25, 31, 33, 35, 35, 37, 42, 51, 51, 51, 52, 52, 52, 54, 54, 58, 61, 61, 61, 65, 68, 70, 70, 70, 73, 74, 76, 76, 77, 77, 77, 78, 78, 79, 80, 81, 81, 81, 81, 87, 88, 88, 88, 89, 90, 90, 90, 90, 90, 90, 90, 90, 91, 91, 91, 91, 91, 93, 93, 95, 98, 98, 98, 99, 101, 101, 102, 103, 103, 103, 104, 105, 105, 107, 107, 108, 110, 110, 110, 110, 110, 110, 110, 111, 111, 112, 114, 114, 114, 114, 118, 118, 118, 118, 118, 120, 121, 122], "INDELs 1-10bp": [17, 36, 40, 65, 87, 87, 98, 106, 108, 113, 114, 114, 114, 115, 115, 116, 116, 116, 121, 122, 122, 122, 124], "INDELs 11-100bp": [22, 26, 26, 26, 39, 44, 48, 48, 52, 54, 56, 56, 56, 61, 62, 63, 78, 78, 83, 104, 108, 108, 112, 116, 116, 116, 117], "INDELs 101-200bp": [44], "INDELs 201-300bp": [1], "INDELs 301-400bp": [0]}, "20": {"SNVs": [41, 41, 54, 77, 94, 96, 100, 104, 115, 124, 130, 142, 147, 154, 156, 178, 192, 192, 199, 199, 199, 208, 215, 218, 229, 229, 230, 234, 235, 235, 236, 236, 236, 236, 236, 236, 238, 238, 238, 238, 239, 239, 239, 239, 239, 239, 239], "INDELs 1-10bp": [3, 32, 33, 35, 43, 47, 48, 59, 66, 77, 98, 100, 100, 102, 102, 104, 104, 107, 107, 107, 107, 112, 112, 112, 116, 117, 119, 143, 144, 144, 144, 144, 145, 145, 149, 149, 149, 160, 163, 164, 167, 169, 169, 169, 174, 195, 231, 233, 248, 252, 252, 252, 254, 254, 256, 257], "INDELs 11-100bp": [7, 22, 24, 36, 36, 38, 39, 39, 39, 50, 50, 53, 55, 56, 57, 57, 61, 62, 63, 65, 67, 68, 70, 70, 71, 71, 79, 89, 108, 113, 116, 116, 117, 121, 122, 122, 122, 124, 126, 131, 136, 137, 139, 140, 140, 140, 142, 142, 142, 154, 154, 156, 162, 162, 162, 166, 167, 169, 169, 171, 172, 173, 175, 175, 177, 178, 181, 186, 187, 188, 196, 197, 197, 197, 198, 198, 202, 210, 211, 211, 211, 221, 224, 224, 224, 224, 225], "INDELs 101-200bp": [109, 117], "INDELs 201-300bp": [15], "INDELs 301-400bp": [0]}, "30": {"SNVs": [27, 34, 42, 46, 55, 55, 55, 59, 72, 72, 79, 84, 87, 96, 100, 104, 104, 113, 142, 147, 156, 156, 159, 174, 180, 188, 194, 201, 224, 224, 241, 242, 244, 246, 246, 247, 247, 247, 248, 248, 248, 248, 248, 248, 248, 248, 249, 250], "INDELs 1-10bp": [4, 4, 4, 9, 10, 11, 16, 20, 34, 34, 34, 38, 46, 46, 46, 60, 65, 69, 69, 69, 86, 86, 88, 90, 90, 105, 113, 113, 139, 144, 147, 151, 151, 157, 159, 168, 168, 183, 186, 188, 189, 199, 205, 206, 215, 217, 217, 217, 231, 242, 244, 246, 252, 252, 252, 255, 255, 255, 255, 255, 258], "INDELs 11-100bp": [24, 24, 25, 27, 29, 30, 30, 36, 36, 36, 42, 47, 48, 50, 50, 74, 74, 83, 85, 85, 86, 86, 87, 90, 92, 93, 94, 97, 103, 112, 115, 118, 119, 119, 121, 121, 125, 126, 126, 129, 129, 131, 131, 131, 134, 138, 141, 143, 143, 146, 146, 147, 147, 152, 153, 153, 157, 160, 160, 160, 160, 160, 160, 161, 162, 162, 162, 165, 171, 171, 171, 177, 180, 181, 183, 184, 186, 188, 202, 202, 207, 207, 207, 207, 210, 210, 211, 214, 215, 218, 218, 218, 221, 222, 222, 222, 223, 225, 226, 226, 227, 228, 228, 228, 228, 228, 228, 228, 228], "INDELs 101-200bp": [30, 61, 97, 102, 107, 118, 118, 120, 131], "INDELs 201-300bp": [13], "INDELs 301-400bp": [0]}, "50": {"SNVs": [48, 99, 109, 128, 128, 130, 138, 143, 151, 202, 202, 207, 207, 213, 220, 245, 245, 247, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 251, 251, 251], "INDELs 1-10bp": [0, 12, 12, 12, 13, 13, 17, 17, 17, 24, 26, 26, 26, 54, 55, 66, 67, 72, 78, 82, 91, 92, 92, 94, 94, 105, 110, 110, 118, 131, 144, 153, 158, 158, 159, 162, 181, 192, 192, 192, 223, 230, 235, 238, 238, 246, 246, 248, 251, 252, 255, 257, 258, 259, 259, 260, 260, 260, 260, 260, 260, 261], "INDELs 11-100bp": [17, 36, 51, 72, 74, 75, 77, 82, 84, 85, 90, 91, 92, 94, 94, 96, 99, 100, 100, 109, 110, 110, 110, 111, 115, 127, 128, 128, 128, 128, 130, 134, 134, 134, 135, 135, 135, 135, 137, 140, 144, 144, 146, 148, 148, 153, 153, 154, 154, 155, 155, 155, 158, 158, 158, 158, 160, 162, 164, 164, 164, 164, 164, 164, 169, 170, 171, 179, 183, 183, 183, 183, 183, 183, 183, 185, 187, 187, 188, 189, 191, 192, 194, 201, 202, 204, 207, 212, 212, 225, 226, 226, 226, 227, 236, 236, 236, 236, 236, 236, 237, 237, 239, 240, 240, 240, 240], "INDELs 101-200bp": [14, 41, 44, 125, 128, 132], "INDELs 201-300bp": [21], "INDELs 301-400bp": [0]}}}
diff --git a/notebook/bigsim/roc-json.ipynb b/notebook/bigsim/roc-json.ipynb
diff --git a/notebook/bigsim/roc-to-json.ipynb b/notebook/bigsim/roc-to-json.ipynb