@jbrowse/plugin-comparative-adapters 2.6.1

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (178) hide show
  1. package/LICENSE +201 -0
  2. package/dist/ChainAdapter/ChainAdapter.d.ts +19 -0
  3. package/dist/ChainAdapter/ChainAdapter.js +22 -0
  4. package/dist/ChainAdapter/ChainAdapter.js.map +1 -0
  5. package/dist/ChainAdapter/configSchema.d.ts +39 -0
  6. package/dist/ChainAdapter/configSchema.js +45 -0
  7. package/dist/ChainAdapter/configSchema.js.map +1 -0
  8. package/dist/ChainAdapter/index.d.ts +3 -0
  9. package/dist/ChainAdapter/index.js +42 -0
  10. package/dist/ChainAdapter/index.js.map +1 -0
  11. package/dist/ChainAdapter/util.d.ts +16 -0
  12. package/dist/ChainAdapter/util.js +134 -0
  13. package/dist/ChainAdapter/util.js.map +1 -0
  14. package/dist/DeltaAdapter/DeltaAdapter.d.ts +20 -0
  15. package/dist/DeltaAdapter/DeltaAdapter.js +22 -0
  16. package/dist/DeltaAdapter/DeltaAdapter.js.map +1 -0
  17. package/dist/DeltaAdapter/configSchema.d.ts +39 -0
  18. package/dist/DeltaAdapter/configSchema.js +45 -0
  19. package/dist/DeltaAdapter/configSchema.js.map +1 -0
  20. package/dist/DeltaAdapter/index.d.ts +3 -0
  21. package/dist/DeltaAdapter/index.js +42 -0
  22. package/dist/DeltaAdapter/index.js.map +1 -0
  23. package/dist/DeltaAdapter/util.d.ts +17 -0
  24. package/dist/DeltaAdapter/util.js +152 -0
  25. package/dist/DeltaAdapter/util.js.map +1 -0
  26. package/dist/MCScanAnchorsAdapter/MCScanAnchorsAdapter.d.ts +33 -0
  27. package/dist/MCScanAnchorsAdapter/MCScanAnchorsAdapter.js +96 -0
  28. package/dist/MCScanAnchorsAdapter/MCScanAnchorsAdapter.js.map +1 -0
  29. package/dist/MCScanAnchorsAdapter/configSchema.d.ts +40 -0
  30. package/dist/MCScanAnchorsAdapter/configSchema.js +48 -0
  31. package/dist/MCScanAnchorsAdapter/configSchema.js.map +1 -0
  32. package/dist/MCScanAnchorsAdapter/index.d.ts +3 -0
  33. package/dist/MCScanAnchorsAdapter/index.js +42 -0
  34. package/dist/MCScanAnchorsAdapter/index.js.map +1 -0
  35. package/dist/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.d.ts +41 -0
  36. package/dist/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.js +121 -0
  37. package/dist/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.js.map +1 -0
  38. package/dist/MCScanSimpleAnchorsAdapter/configSchema.d.ts +40 -0
  39. package/dist/MCScanSimpleAnchorsAdapter/configSchema.js +48 -0
  40. package/dist/MCScanSimpleAnchorsAdapter/configSchema.js.map +1 -0
  41. package/dist/MCScanSimpleAnchorsAdapter/index.d.ts +3 -0
  42. package/dist/MCScanSimpleAnchorsAdapter/index.js +42 -0
  43. package/dist/MCScanSimpleAnchorsAdapter/index.js.map +1 -0
  44. package/dist/MashMapAdapter/MashMapAdapter.d.ts +5 -0
  45. package/dist/MashMapAdapter/MashMapAdapter.js +42 -0
  46. package/dist/MashMapAdapter/MashMapAdapter.js.map +1 -0
  47. package/dist/MashMapAdapter/configSchema.d.ts +37 -0
  48. package/dist/MashMapAdapter/configSchema.js +45 -0
  49. package/dist/MashMapAdapter/configSchema.js.map +1 -0
  50. package/dist/MashMapAdapter/index.d.ts +3 -0
  51. package/dist/MashMapAdapter/index.js +42 -0
  52. package/dist/MashMapAdapter/index.js.map +1 -0
  53. package/dist/PAFAdapter/PAFAdapter.d.ts +20 -0
  54. package/dist/PAFAdapter/PAFAdapter.js +152 -0
  55. package/dist/PAFAdapter/PAFAdapter.js.map +1 -0
  56. package/dist/PAFAdapter/SyntenyFeature.d.ts +4 -0
  57. package/dist/PAFAdapter/SyntenyFeature.js +17 -0
  58. package/dist/PAFAdapter/SyntenyFeature.js.map +1 -0
  59. package/dist/PAFAdapter/configSchema.d.ts +37 -0
  60. package/dist/PAFAdapter/configSchema.js +45 -0
  61. package/dist/PAFAdapter/configSchema.js.map +1 -0
  62. package/dist/PAFAdapter/index.d.ts +3 -0
  63. package/dist/PAFAdapter/index.js +42 -0
  64. package/dist/PAFAdapter/index.js.map +1 -0
  65. package/dist/PAFAdapter/util.d.ts +20 -0
  66. package/dist/PAFAdapter/util.js +134 -0
  67. package/dist/PAFAdapter/util.js.map +1 -0
  68. package/dist/index.d.ts +6 -0
  69. package/dist/index.js +51 -0
  70. package/dist/index.js.map +1 -0
  71. package/dist/util.d.ts +16 -0
  72. package/dist/util.js +57 -0
  73. package/dist/util.js.map +1 -0
  74. package/esm/ChainAdapter/ChainAdapter.d.ts +19 -0
  75. package/esm/ChainAdapter/ChainAdapter.js +16 -0
  76. package/esm/ChainAdapter/ChainAdapter.js.map +1 -0
  77. package/esm/ChainAdapter/configSchema.d.ts +39 -0
  78. package/esm/ChainAdapter/configSchema.js +43 -0
  79. package/esm/ChainAdapter/configSchema.js.map +1 -0
  80. package/esm/ChainAdapter/index.d.ts +3 -0
  81. package/esm/ChainAdapter/index.js +14 -0
  82. package/esm/ChainAdapter/index.js.map +1 -0
  83. package/esm/ChainAdapter/util.d.ts +16 -0
  84. package/esm/ChainAdapter/util.js +130 -0
  85. package/esm/ChainAdapter/util.js.map +1 -0
  86. package/esm/DeltaAdapter/DeltaAdapter.d.ts +20 -0
  87. package/esm/DeltaAdapter/DeltaAdapter.js +16 -0
  88. package/esm/DeltaAdapter/DeltaAdapter.js.map +1 -0
  89. package/esm/DeltaAdapter/configSchema.d.ts +39 -0
  90. package/esm/DeltaAdapter/configSchema.js +43 -0
  91. package/esm/DeltaAdapter/configSchema.js.map +1 -0
  92. package/esm/DeltaAdapter/index.d.ts +3 -0
  93. package/esm/DeltaAdapter/index.js +14 -0
  94. package/esm/DeltaAdapter/index.js.map +1 -0
  95. package/esm/DeltaAdapter/util.d.ts +17 -0
  96. package/esm/DeltaAdapter/util.js +148 -0
  97. package/esm/DeltaAdapter/util.js.map +1 -0
  98. package/esm/MCScanAnchorsAdapter/MCScanAnchorsAdapter.d.ts +33 -0
  99. package/esm/MCScanAnchorsAdapter/MCScanAnchorsAdapter.js +94 -0
  100. package/esm/MCScanAnchorsAdapter/MCScanAnchorsAdapter.js.map +1 -0
  101. package/esm/MCScanAnchorsAdapter/configSchema.d.ts +40 -0
  102. package/esm/MCScanAnchorsAdapter/configSchema.js +46 -0
  103. package/esm/MCScanAnchorsAdapter/configSchema.js.map +1 -0
  104. package/esm/MCScanAnchorsAdapter/index.d.ts +3 -0
  105. package/esm/MCScanAnchorsAdapter/index.js +14 -0
  106. package/esm/MCScanAnchorsAdapter/index.js.map +1 -0
  107. package/esm/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.d.ts +41 -0
  108. package/esm/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.js +116 -0
  109. package/esm/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.js.map +1 -0
  110. package/esm/MCScanSimpleAnchorsAdapter/configSchema.d.ts +40 -0
  111. package/esm/MCScanSimpleAnchorsAdapter/configSchema.js +46 -0
  112. package/esm/MCScanSimpleAnchorsAdapter/configSchema.js.map +1 -0
  113. package/esm/MCScanSimpleAnchorsAdapter/index.d.ts +3 -0
  114. package/esm/MCScanSimpleAnchorsAdapter/index.js +14 -0
  115. package/esm/MCScanSimpleAnchorsAdapter/index.js.map +1 -0
  116. package/esm/MashMapAdapter/MashMapAdapter.d.ts +5 -0
  117. package/esm/MashMapAdapter/MashMapAdapter.js +36 -0
  118. package/esm/MashMapAdapter/MashMapAdapter.js.map +1 -0
  119. package/esm/MashMapAdapter/configSchema.d.ts +37 -0
  120. package/esm/MashMapAdapter/configSchema.js +43 -0
  121. package/esm/MashMapAdapter/configSchema.js.map +1 -0
  122. package/esm/MashMapAdapter/index.d.ts +3 -0
  123. package/esm/MashMapAdapter/index.js +14 -0
  124. package/esm/MashMapAdapter/index.js.map +1 -0
  125. package/esm/PAFAdapter/PAFAdapter.d.ts +20 -0
  126. package/esm/PAFAdapter/PAFAdapter.js +147 -0
  127. package/esm/PAFAdapter/PAFAdapter.js.map +1 -0
  128. package/esm/PAFAdapter/SyntenyFeature.d.ts +4 -0
  129. package/esm/PAFAdapter/SyntenyFeature.js +14 -0
  130. package/esm/PAFAdapter/SyntenyFeature.js.map +1 -0
  131. package/esm/PAFAdapter/configSchema.d.ts +37 -0
  132. package/esm/PAFAdapter/configSchema.js +43 -0
  133. package/esm/PAFAdapter/configSchema.js.map +1 -0
  134. package/esm/PAFAdapter/index.d.ts +3 -0
  135. package/esm/PAFAdapter/index.js +14 -0
  136. package/esm/PAFAdapter/index.js.map +1 -0
  137. package/esm/PAFAdapter/util.d.ts +20 -0
  138. package/esm/PAFAdapter/util.js +127 -0
  139. package/esm/PAFAdapter/util.js.map +1 -0
  140. package/esm/index.d.ts +6 -0
  141. package/esm/index.js +45 -0
  142. package/esm/index.js.map +1 -0
  143. package/esm/util.d.ts +16 -0
  144. package/esm/util.js +49 -0
  145. package/esm/util.js.map +1 -0
  146. package/package.json +63 -0
  147. package/src/ChainAdapter/ChainAdapter.ts +18 -0
  148. package/src/ChainAdapter/configSchema.ts +50 -0
  149. package/src/ChainAdapter/index.ts +18 -0
  150. package/src/ChainAdapter/util.ts +170 -0
  151. package/src/DeltaAdapter/DeltaAdapter.ts +18 -0
  152. package/src/DeltaAdapter/configSchema.ts +50 -0
  153. package/src/DeltaAdapter/index.ts +18 -0
  154. package/src/DeltaAdapter/util.ts +149 -0
  155. package/src/MCScanAnchorsAdapter/MCScanAnchorsAdapter.test.ts +45 -0
  156. package/src/MCScanAnchorsAdapter/MCScanAnchorsAdapter.ts +134 -0
  157. package/src/MCScanAnchorsAdapter/configSchema.ts +52 -0
  158. package/src/MCScanAnchorsAdapter/index.ts +20 -0
  159. package/src/MCScanAnchorsAdapter/test_data/grape.bed.gz +0 -0
  160. package/src/MCScanAnchorsAdapter/test_data/grape.peach.anchors.gz +0 -0
  161. package/src/MCScanAnchorsAdapter/test_data/peach.bed.gz +0 -0
  162. package/src/MCScanSimpleAnchorsAdapter/MCScanSimpleAnchorsAdapter.ts +157 -0
  163. package/src/MCScanSimpleAnchorsAdapter/configSchema.ts +52 -0
  164. package/src/MCScanSimpleAnchorsAdapter/index.ts +19 -0
  165. package/src/MashMapAdapter/MashMapAdapter.ts +40 -0
  166. package/src/MashMapAdapter/configSchema.ts +51 -0
  167. package/src/MashMapAdapter/index.ts +18 -0
  168. package/src/PAFAdapter/PAFAdapter.test.ts +37 -0
  169. package/src/PAFAdapter/PAFAdapter.ts +181 -0
  170. package/src/PAFAdapter/SyntenyFeature.ts +15 -0
  171. package/src/PAFAdapter/configSchema.ts +50 -0
  172. package/src/PAFAdapter/index.ts +18 -0
  173. package/src/PAFAdapter/test_data/grape.peach.anchors +14966 -0
  174. package/src/PAFAdapter/test_data/peach_grape.paf +30 -0
  175. package/src/PAFAdapter/util.test.ts +7 -0
  176. package/src/PAFAdapter/util.ts +176 -0
  177. package/src/index.ts +62 -0
  178. package/src/util.ts +66 -0
@@ -0,0 +1,147 @@
1
+ import { BaseFeatureDataAdapter, } from '@jbrowse/core/data_adapters/BaseAdapter';
2
+ import { doesIntersect2 } from '@jbrowse/core/util/range';
3
+ import { openLocation } from '@jbrowse/core/util/io';
4
+ import { ObservableCreate } from '@jbrowse/core/util/rxjs';
5
+ import { readConfObject, } from '@jbrowse/core/configuration';
6
+ import { unzip } from '@gmod/bgzf-filehandle';
7
+ import { MismatchParser } from '@jbrowse/plugin-alignments';
8
+ // locals
9
+ import SyntenyFeature from './SyntenyFeature';
10
+ import { isGzip, parseLineByLine } from '../util';
11
+ import { getWeightedMeans, flipCigar, swapIndelCigar, parsePAFLine, } from './util';
12
+ const { parseCigar } = MismatchParser;
13
+ class PAFAdapter extends BaseFeatureDataAdapter {
14
+ async setup(opts) {
15
+ if (!this.setupP) {
16
+ this.setupP = this.setupPre(opts).catch(e => {
17
+ this.setupP = undefined;
18
+ throw e;
19
+ });
20
+ }
21
+ return this.setupP;
22
+ }
23
+ async setupPre(opts) {
24
+ const pm = this.pluginManager;
25
+ const pafLocation = openLocation(this.getConf('pafLocation'), pm);
26
+ const buffer = (await pafLocation.readFile(opts));
27
+ const buf = isGzip(buffer) ? await unzip(buffer) : buffer;
28
+ return parseLineByLine(buf, parsePAFLine);
29
+ }
30
+ async hasDataForRefName() {
31
+ // determining this properly is basically a call to getFeatures
32
+ // so is not really that important, and has to be true or else
33
+ // getFeatures is never called (BaseAdapter filters it out)
34
+ return true;
35
+ }
36
+ getAssemblyNames() {
37
+ const assemblyNames = this.getConf('assemblyNames');
38
+ if (assemblyNames.length === 0) {
39
+ const query = this.getConf('queryAssembly');
40
+ const target = this.getConf('targetAssembly');
41
+ return [query, target];
42
+ }
43
+ return assemblyNames;
44
+ }
45
+ async getRefNames(opts = {}) {
46
+ var _a;
47
+ // @ts-expect-error
48
+ const r1 = (_a = opts.regions) === null || _a === void 0 ? void 0 : _a[0].assemblyName;
49
+ const feats = await this.setup(opts);
50
+ const idx = this.getAssemblyNames().indexOf(r1);
51
+ if (idx !== -1) {
52
+ const set = new Set();
53
+ for (const feat of feats) {
54
+ set.add(idx === 0 ? feat.qname : feat.tname);
55
+ }
56
+ return [...set];
57
+ }
58
+ console.warn('Unable to do ref renaming on adapter');
59
+ return [];
60
+ }
61
+ getFeatures(query, opts = {}) {
62
+ return ObservableCreate(async (observer) => {
63
+ let pafRecords = await this.setup(opts);
64
+ const { config } = opts;
65
+ // note: this is not the adapter config, it is responding to a display
66
+ // setting passed in via the opts parameter
67
+ if (config && readConfObject(config, 'colorBy') === 'meanQueryIdentity') {
68
+ pafRecords = getWeightedMeans(pafRecords);
69
+ }
70
+ const assemblyNames = this.getAssemblyNames();
71
+ // The index of the assembly name in the query list corresponds to the
72
+ // adapter in the subadapters list
73
+ const index = assemblyNames.indexOf(query.assemblyName);
74
+ const { start: qstart, end: qend, refName: qref, assemblyName } = query;
75
+ if (index === -1) {
76
+ console.warn(`${assemblyName} not found in this adapter`);
77
+ observer.complete();
78
+ }
79
+ for (let i = 0; i < pafRecords.length; i++) {
80
+ const r = pafRecords[i];
81
+ let start = 0;
82
+ let end = 0;
83
+ let refName = '';
84
+ let mateName = '';
85
+ let mateStart = 0;
86
+ let mateEnd = 0;
87
+ const flip = index === 0;
88
+ const assemblyName = assemblyNames[+!flip];
89
+ if (index === 0) {
90
+ start = r.qstart;
91
+ end = r.qend;
92
+ refName = r.qname;
93
+ mateName = r.tname;
94
+ mateStart = r.tstart;
95
+ mateEnd = r.tend;
96
+ }
97
+ else {
98
+ start = r.tstart;
99
+ end = r.tend;
100
+ refName = r.tname;
101
+ mateName = r.qname;
102
+ mateStart = r.qstart;
103
+ mateEnd = r.qend;
104
+ }
105
+ const { extra, strand } = r;
106
+ if (refName === qref && doesIntersect2(qstart, qend, start, end)) {
107
+ const { numMatches = 0, blockLen = 1, cg, ...rest } = extra;
108
+ let CIGAR = extra.cg;
109
+ if (extra.cg) {
110
+ if (flip && strand === -1) {
111
+ CIGAR = flipCigar(parseCigar(extra.cg)).join('');
112
+ }
113
+ else if (flip) {
114
+ CIGAR = swapIndelCigar(extra.cg);
115
+ }
116
+ }
117
+ observer.next(new SyntenyFeature({
118
+ uniqueId: i + assemblyName,
119
+ assemblyName,
120
+ start,
121
+ end,
122
+ type: 'match',
123
+ refName,
124
+ strand,
125
+ ...rest,
126
+ CIGAR,
127
+ syntenyId: i,
128
+ identity: numMatches / blockLen,
129
+ numMatches,
130
+ blockLen,
131
+ mate: {
132
+ start: mateStart,
133
+ end: mateEnd,
134
+ refName: mateName,
135
+ assemblyName: assemblyNames[+flip],
136
+ },
137
+ }));
138
+ }
139
+ }
140
+ observer.complete();
141
+ });
142
+ }
143
+ freeResources( /* { query } */) { }
144
+ }
145
+ PAFAdapter.capabilities = ['getFeatures', 'getRefNames'];
146
+ export default PAFAdapter;
147
+ //# sourceMappingURL=PAFAdapter.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"PAFAdapter.js","sourceRoot":"","sources":["../../src/PAFAdapter/PAFAdapter.ts"],"names":[],"mappings":"AAAA,OAAO,EACL,sBAAsB,GAEvB,MAAM,yCAAyC,CAAA;AAEhD,OAAO,EAAE,cAAc,EAAE,MAAM,0BAA0B,CAAA;AACzD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAA;AACpD,OAAO,EAAE,gBAAgB,EAAE,MAAM,yBAAyB,CAAA;AAE1D,OAAO,EAEL,cAAc,GACf,MAAM,6BAA6B,CAAA;AACpC,OAAO,EAAE,KAAK,EAAE,MAAM,uBAAuB,CAAA;AAC7C,OAAO,EAAE,cAAc,EAAE,MAAM,4BAA4B,CAAA;AAE3D,SAAS;AACT,OAAO,cAAc,MAAM,kBAAkB,CAAA;AAC7C,OAAO,EAAE,MAAM,EAAE,eAAe,EAAE,MAAM,SAAS,CAAA;AACjD,OAAO,EACL,gBAAgB,EAChB,SAAS,EACT,cAAc,EACd,YAAY,GAEb,MAAM,QAAQ,CAAA;AAEf,MAAM,EAAE,UAAU,EAAE,GAAG,cAAc,CAAA;AAMrC,MAAqB,UAAW,SAAQ,sBAAsB;IAK5D,KAAK,CAAC,KAAK,CAAC,IAAkB;QAC5B,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE;YAChB,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE;gBAC1C,IAAI,CAAC,MAAM,GAAG,SAAS,CAAA;gBACvB,MAAM,CAAC,CAAA;YACT,CAAC,CAAC,CAAA;SACH;QACD,OAAO,IAAI,CAAC,MAAM,CAAA;IACpB,CAAC;IAED,KAAK,CAAC,QAAQ,CAAC,IAAkB;QAC/B,MAAM,EAAE,GAAG,IAAI,CAAC,aAAa,CAAA;QAC7B,MAAM,WAAW,GAAG,YAAY,CAAC,IAAI,CAAC,OAAO,CAAC,aAAa,CAAC,EAAE,EAAE,CAAC,CAAA;QACjE,MAAM,MAAM,GAAG,CAAC,MAAM,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAW,CAAA;QAC3D,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAA;QACzD,OAAO,eAAe,CAAC,GAAG,EAAE,YAAY,CAAC,CAAA;IAC3C,CAAC;IAED,KAAK,CAAC,iBAAiB;QACrB,+DAA+D;QAC/D,8DAA8D;QAC9D,2DAA2D;QAC3D,OAAO,IAAI,CAAA;IACb,CAAC;IAED,gBAAgB;QACd,MAAM,aAAa,GAAG,IAAI,CAAC,OAAO,CAAC,eAAe,CAAa,CAAA;QAC/D,IAAI,aAAa,CAAC,MAAM,KAAK,CAAC,EAAE;YAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,eAAe,CAAW,CAAA;YACrD,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,gBAAgB,CAAW,CAAA;YACvD,OAAO,CAAC,KAAK,EAAE,MAAM,CAAC,CAAA;SACvB;QACD,OAAO,aAAa,CAAA;IACtB,CAAC;IAED,KAAK,CAAC,WAAW,CAAC,OAAoB,EAAE;;QACtC,mBAAmB;QACnB,MAAM,EAAE,GAAG,MAAA,IAAI,CAAC,OAAO,0CAAG,CAAC,EAAE,YAAY,CAAA;QACzC,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAA;QAEpC,MAAM,GAAG,GAAG,IAAI,CAAC,gBAAgB,EAAE,CAAC,OAAO,CAAC,EAAE,CAAC,CAAA;QAC/C,IAAI,GAAG,KAAK,CAAC,CAAC,EAAE;YACd,MAAM,GAAG,GAAG,IAAI,GAAG,EAAU,CAAA;YAC7B,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE;gBACxB,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;aAC7C;YACD,OAAO,CAAC,GAAG,GAAG,CAAC,CAAA;SAChB;QACD,OAAO,CAAC,IAAI,CAAC,sCAAsC,CAAC,CAAA;QACpD,OAAO,EAAE,CAAA;IACX,CAAC;IAED,WAAW,CAAC,KAAa,EAAE,OAAmB,EAAE;QAC9C,OAAO,gBAAgB,CAAU,KAAK,EAAC,QAAQ,EAAC,EAAE;YAChD,IAAI,UAAU,GAAG,MAAM,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAA;YACvC,MAAM,EAAE,MAAM,EAAE,GAAG,IAAI,CAAA;YAEvB,sEAAsE;YACtE,2CAA2C;YAC3C,IAAI,MAAM,IAAI,cAAc,CAAC,MAAM,EAAE,SAAS,CAAC,KAAK,mBAAmB,EAAE;gBACvE,UAAU,GAAG,gBAAgB,CAAC,UAAU,CAAC,CAAA;aAC1C;YACD,MAAM,aAAa,GAAG,IAAI,CAAC,gBAAgB,EAAE,CAAA;YAE7C,sEAAsE;YACtE,kCAAkC;YAClC,MAAM,KAAK,GAAG,aAAa,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAA;YACvD,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,YAAY,EAAE,GAAG,KAAK,CAAA;YACvE,IAAI,KAAK,KAAK,CAAC,CAAC,EAAE;gBAChB,OAAO,CAAC,IAAI,CAAC,GAAG,YAAY,4BAA4B,CAAC,CAAA;gBACzD,QAAQ,CAAC,QAAQ,EAAE,CAAA;aACpB;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBAC1C,MAAM,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAA;gBACvB,IAAI,KAAK,GAAG,CAAC,CAAA;gBACb,IAAI,GAAG,GAAG,CAAC,CAAA;gBACX,IAAI,OAAO,GAAG,EAAE,CAAA;gBAChB,IAAI,QAAQ,GAAG,EAAE,CAAA;gBACjB,IAAI,SAAS,GAAG,CAAC,CAAA;gBACjB,IAAI,OAAO,GAAG,CAAC,CAAA;gBACf,MAAM,IAAI,GAAG,KAAK,KAAK,CAAC,CAAA;gBACxB,MAAM,YAAY,GAAG,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,CAAA;gBAC1C,IAAI,KAAK,KAAK,CAAC,EAAE;oBACf,KAAK,GAAG,CAAC,CAAC,MAAM,CAAA;oBAChB,GAAG,GAAG,CAAC,CAAC,IAAI,CAAA;oBACZ,OAAO,GAAG,CAAC,CAAC,KAAK,CAAA;oBACjB,QAAQ,GAAG,CAAC,CAAC,KAAK,CAAA;oBAClB,SAAS,GAAG,CAAC,CAAC,MAAM,CAAA;oBACpB,OAAO,GAAG,CAAC,CAAC,IAAI,CAAA;iBACjB;qBAAM;oBACL,KAAK,GAAG,CAAC,CAAC,MAAM,CAAA;oBAChB,GAAG,GAAG,CAAC,CAAC,IAAI,CAAA;oBACZ,OAAO,GAAG,CAAC,CAAC,KAAK,CAAA;oBACjB,QAAQ,GAAG,CAAC,CAAC,KAAK,CAAA;oBAClB,SAAS,GAAG,CAAC,CAAC,MAAM,CAAA;oBACpB,OAAO,GAAG,CAAC,CAAC,IAAI,CAAA;iBACjB;gBACD,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,CAAC,CAAA;gBAC3B,IAAI,OAAO,KAAK,IAAI,IAAI,cAAc,CAAC,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,CAAC,EAAE;oBAChE,MAAM,EAAE,UAAU,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,IAAI,EAAE,GAAG,KAAK,CAAA;oBAE3D,IAAI,KAAK,GAAG,KAAK,CAAC,EAAE,CAAA;oBACpB,IAAI,KAAK,CAAC,EAAE,EAAE;wBACZ,IAAI,IAAI,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE;4BACzB,KAAK,GAAG,SAAS,CAAC,UAAU,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAA;yBACjD;6BAAM,IAAI,IAAI,EAAE;4BACf,KAAK,GAAG,cAAc,CAAC,KAAK,CAAC,EAAE,CAAC,CAAA;yBACjC;qBACF;oBAED,QAAQ,CAAC,IAAI,CACX,IAAI,cAAc,CAAC;wBACjB,QAAQ,EAAE,CAAC,GAAG,YAAY;wBAC1B,YAAY;wBACZ,KAAK;wBACL,GAAG;wBACH,IAAI,EAAE,OAAO;wBACb,OAAO;wBACP,MAAM;wBACN,GAAG,IAAI;wBACP,KAAK;wBACL,SAAS,EAAE,CAAC;wBACZ,QAAQ,EAAE,UAAU,GAAG,QAAQ;wBAC/B,UAAU;wBACV,QAAQ;wBACR,IAAI,EAAE;4BACJ,KAAK,EAAE,SAAS;4BAChB,GAAG,EAAE,OAAO;4BACZ,OAAO,EAAE,QAAQ;4BACjB,YAAY,EAAE,aAAa,CAAC,CAAC,IAAI,CAAC;yBACnC;qBACF,CAAC,CACH,CAAA;iBACF;aACF;YAED,QAAQ,CAAC,QAAQ,EAAE,CAAA;QACrB,CAAC,CAAC,CAAA;IACJ,CAAC;IAED,aAAa,EAAC,eAAe,IAAS,CAAC;;AA/IzB,uBAAY,GAAG,CAAC,aAAa,EAAE,aAAa,CAAC,CAAA;eAHxC,UAAU"}
@@ -0,0 +1,4 @@
1
+ import { SimpleFeature } from '@jbrowse/core/util';
2
+ export default class SyntenyFeature extends SimpleFeature {
3
+ get(arg: string): any;
4
+ }
@@ -0,0 +1,14 @@
1
+ import { SimpleFeature } from '@jbrowse/core/util';
2
+ import { MismatchParser } from '@jbrowse/plugin-alignments';
3
+ // locals
4
+ const { getMismatches } = MismatchParser;
5
+ export default class SyntenyFeature extends SimpleFeature {
6
+ // eslint-disable-next-line @typescript-eslint/no-explicit-any
7
+ get(arg) {
8
+ if (arg === 'mismatches') {
9
+ return getMismatches(this.get('CIGAR'));
10
+ }
11
+ return super.get(arg);
12
+ }
13
+ }
14
+ //# sourceMappingURL=SyntenyFeature.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"SyntenyFeature.js","sourceRoot":"","sources":["../../src/PAFAdapter/SyntenyFeature.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAA;AAClD,OAAO,EAAE,cAAc,EAAE,MAAM,4BAA4B,CAAA;AAE3D,SAAS;AACT,MAAM,EAAE,aAAa,EAAE,GAAG,cAAc,CAAA;AAExC,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,aAAa;IACvD,8DAA8D;IAC9D,GAAG,CAAC,GAAW;QACb,IAAI,GAAG,KAAK,YAAY,EAAE;YACxB,OAAO,aAAa,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAA;SACxC;QACD,OAAO,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;IACvB,CAAC;CACF"}
@@ -0,0 +1,37 @@
1
+ declare const PAFAdapter: import("@jbrowse/core/configuration/configurationSchema").ConfigurationSchemaType<{
2
+ /**
3
+ * #slot
4
+ */
5
+ assemblyNames: {
6
+ type: string;
7
+ defaultValue: never[];
8
+ description: string;
9
+ };
10
+ /**
11
+ * #slot
12
+ */
13
+ targetAssembly: {
14
+ type: string;
15
+ defaultValue: string;
16
+ description: string;
17
+ };
18
+ /**
19
+ * #slot
20
+ */
21
+ queryAssembly: {
22
+ type: string;
23
+ defaultValue: string;
24
+ description: string;
25
+ };
26
+ /**
27
+ * #slot
28
+ */
29
+ pafLocation: {
30
+ type: string;
31
+ defaultValue: {
32
+ uri: string;
33
+ locationType: string;
34
+ };
35
+ };
36
+ }, import("@jbrowse/core/configuration/configurationSchema").ConfigurationSchemaOptions<undefined, undefined>>;
37
+ export default PAFAdapter;
@@ -0,0 +1,43 @@
1
+ import { ConfigurationSchema } from '@jbrowse/core/configuration';
2
+ /**
3
+ * #config PAFAdapter
4
+ */
5
+ function x() { } // eslint-disable-line @typescript-eslint/no-unused-vars
6
+ const PAFAdapter = ConfigurationSchema('PAFAdapter', {
7
+ /**
8
+ * #slot
9
+ */
10
+ assemblyNames: {
11
+ type: 'stringArray',
12
+ defaultValue: [],
13
+ description: 'Array of assembly names to use for this file. The target assembly name is the first value in the array, query assembly name is the second',
14
+ },
15
+ /**
16
+ * #slot
17
+ */
18
+ targetAssembly: {
19
+ type: 'string',
20
+ defaultValue: '',
21
+ description: 'Alternative to assemblyNames: the target assembly name',
22
+ },
23
+ /**
24
+ * #slot
25
+ */
26
+ queryAssembly: {
27
+ type: 'string',
28
+ defaultValue: '',
29
+ description: 'Alternative to assemblyNames: the query assembly name',
30
+ },
31
+ /**
32
+ * #slot
33
+ */
34
+ pafLocation: {
35
+ type: 'fileLocation',
36
+ defaultValue: {
37
+ uri: '/path/to/file.paf',
38
+ locationType: 'UriLocation',
39
+ },
40
+ },
41
+ }, { explicitlyTyped: true });
42
+ export default PAFAdapter;
43
+ //# sourceMappingURL=configSchema.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"configSchema.js","sourceRoot":"","sources":["../../src/PAFAdapter/configSchema.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,mBAAmB,EAAE,MAAM,6BAA6B,CAAA;AAEjE;;GAEG;AACH,SAAS,CAAC,KAAI,CAAC,CAAC,wDAAwD;AAExE,MAAM,UAAU,GAAG,mBAAmB,CACpC,YAAY,EACZ;IACE;;OAEG;IACH,aAAa,EAAE;QACb,IAAI,EAAE,aAAa;QACnB,YAAY,EAAE,EAAE;QAChB,WAAW,EACT,2IAA2I;KAC9I;IACD;;OAEG;IACH,cAAc,EAAE;QACd,IAAI,EAAE,QAAQ;QACd,YAAY,EAAE,EAAE;QAChB,WAAW,EAAE,wDAAwD;KACtE;IACD;;OAEG;IACH,aAAa,EAAE;QACb,IAAI,EAAE,QAAQ;QACd,YAAY,EAAE,EAAE;QAChB,WAAW,EAAE,uDAAuD;KACrE;IACD;;OAEG;IACH,WAAW,EAAE;QACX,IAAI,EAAE,cAAc;QACpB,YAAY,EAAE;YACZ,GAAG,EAAE,mBAAmB;YACxB,YAAY,EAAE,aAAa;SAC5B;KACF;CACF,EACD,EAAE,eAAe,EAAE,IAAI,EAAE,CAC1B,CAAA;AAED,eAAe,UAAU,CAAA"}
@@ -0,0 +1,3 @@
1
+ import PluginManager from '@jbrowse/core/PluginManager';
2
+ declare const _default: (pluginManager: PluginManager) => void;
3
+ export default _default;
@@ -0,0 +1,14 @@
1
+ import AdapterType from '@jbrowse/core/pluggableElementTypes/AdapterType';
2
+ import configSchema from './configSchema';
3
+ export default (pluginManager) => {
4
+ pluginManager.addAdapterType(() => new AdapterType({
5
+ name: 'PAFAdapter',
6
+ displayName: 'PAF adapter',
7
+ configSchema,
8
+ adapterMetadata: {
9
+ hiddenFromGUI: true,
10
+ },
11
+ getAdapterClass: () => import('./PAFAdapter').then(r => r.default),
12
+ }));
13
+ };
14
+ //# sourceMappingURL=index.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/PAFAdapter/index.ts"],"names":[],"mappings":"AACA,OAAO,WAAW,MAAM,iDAAiD,CAAA;AACzE,OAAO,YAAY,MAAM,gBAAgB,CAAA;AAEzC,eAAe,CAAC,aAA4B,EAAE,EAAE;IAC9C,aAAa,CAAC,cAAc,CAC1B,GAAG,EAAE,CACH,IAAI,WAAW,CAAC;QACd,IAAI,EAAE,YAAY;QAClB,WAAW,EAAE,aAAa;QAC1B,YAAY;QACZ,eAAe,EAAE;YACf,aAAa,EAAE,IAAI;SACpB;QACD,eAAe,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC;KACnE,CAAC,CACL,CAAA;AACH,CAAC,CAAA"}
@@ -0,0 +1,20 @@
1
+ export interface PAFRecord {
2
+ qname: string;
3
+ qstart: number;
4
+ qend: number;
5
+ tname: string;
6
+ tstart: number;
7
+ tend: number;
8
+ strand: number;
9
+ extra: {
10
+ cg?: string;
11
+ blockLen?: number;
12
+ mappingQual: number;
13
+ numMatches?: number;
14
+ meanScore?: number;
15
+ };
16
+ }
17
+ export declare function getWeightedMeans(ret: PAFRecord[]): PAFRecord[];
18
+ export declare function parsePAFLine(line: string): PAFRecord;
19
+ export declare function flipCigar(cigar: string[]): string[];
20
+ export declare function swapIndelCigar(cigar: string): string;
@@ -0,0 +1,127 @@
1
+ import { zip } from '../util';
2
+ // based on "weighted mean" method from https://github.com/tpoorten/dotPlotly
3
+ // License reproduced here
4
+ //
5
+ // MIT License
6
+ // Copyright (c) 2017 Tom Poorten
7
+ // Permission is hereby granted, free of charge, to any person obtaining a copy
8
+ // of this software and associated documentation files (the "Software"), to deal
9
+ // in the Software without restriction, including without limitation the rights
10
+ // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
11
+ // copies of the Software, and to permit persons to whom the Software is
12
+ // furnished to do so, subject to the following conditions:
13
+ // The above copyright notice and this permission notice shall be included in all
14
+ // copies or substantial portions of the Software.
15
+ // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
16
+ // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
17
+ // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
18
+ // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
19
+ // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
20
+ // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
21
+ // SOFTWARE.
22
+ //
23
+ // Notes: in the weighted mean longer alignments factor in more heavily of all
24
+ // the fragments of a query vs the reference that it mapped to
25
+ //
26
+ // this uses a combined key query+'-'+ref to iteratively map all the alignments
27
+ // that match a particular ref from a particular query (so 1d array of what
28
+ // could be a 2d map)
29
+ //
30
+ // the result is a single number that says e.g. chr5 from human mapped to chr5
31
+ // on mouse with 0.8 quality, and that0.8 is then attached to all the pieces of
32
+ // chr5 on human that mapped to chr5 on mouse. if chr5 on human also more
33
+ // weakly mapped to chr6 on mouse, then it would have another value e.g. 0.6.
34
+ // this can show strong and weak levels of synteny, especially in polyploidy
35
+ // situations
36
+ export function getWeightedMeans(ret) {
37
+ const scoreMap = {};
38
+ for (let i = 0; i < ret.length; i++) {
39
+ const entry = ret[i];
40
+ const query = entry.qname;
41
+ const target = entry.tname;
42
+ const key = query + '-' + target;
43
+ if (!scoreMap[key]) {
44
+ scoreMap[key] = { quals: [], len: [] };
45
+ }
46
+ scoreMap[key].quals.push(entry.extra.mappingQual);
47
+ scoreMap[key].len.push(entry.extra.blockLen || 1);
48
+ }
49
+ const meanScoreMap = Object.fromEntries(Object.entries(scoreMap).map(([key, val]) => {
50
+ const vals = zip(val.quals, val.len);
51
+ return [key, weightedMean(vals)];
52
+ }));
53
+ for (let i = 0; i < ret.length; i++) {
54
+ const entry = ret[i];
55
+ const query = entry.qname;
56
+ const target = entry.tname;
57
+ const key = query + '-' + target;
58
+ entry.extra.meanScore = meanScoreMap[key];
59
+ }
60
+ let min = 10000;
61
+ let max = 0;
62
+ for (let i = 0; i < ret.length; i++) {
63
+ const entry = ret[i];
64
+ min = Math.min(entry.extra.meanScore || 0, min);
65
+ max = Math.max(entry.extra.meanScore || 0, max);
66
+ }
67
+ for (let i = 0; i < ret.length; i++) {
68
+ const entry = ret[i];
69
+ const b = entry.extra.meanScore || 0;
70
+ entry.extra.meanScore = (b - min) / (max - min);
71
+ }
72
+ return ret;
73
+ }
74
+ // https://gist.github.com/stekhn/a12ed417e91f90ecec14bcfa4c2ae16a
75
+ function weightedMean(tuples) {
76
+ // eslint-disable-next-line unicorn/no-array-reduce
77
+ const [valueSum, weightSum] = tuples.reduce(([valueSum, weightSum], [value, weight]) => [
78
+ valueSum + value * weight,
79
+ weightSum + weight,
80
+ ], [0, 0]);
81
+ return valueSum / weightSum;
82
+ }
83
+ export function parsePAFLine(line) {
84
+ const [qname, , qstart, qend, strand, tname, , tstart, tend, numMatches, blockLen, mappingQual, ...fields] = line.split('\t');
85
+ const rest = Object.fromEntries(fields.map(field => {
86
+ const r = field.indexOf(':');
87
+ const fieldName = field.slice(0, r);
88
+ const fieldValue = field.slice(r + 3);
89
+ return [fieldName, fieldValue];
90
+ }));
91
+ return {
92
+ tname,
93
+ tstart: +tstart,
94
+ tend: +tend,
95
+ qname,
96
+ qstart: +qstart,
97
+ qend: +qend,
98
+ strand: strand === '-' ? -1 : 1,
99
+ extra: {
100
+ numMatches: +numMatches,
101
+ blockLen: +blockLen,
102
+ mappingQual: +mappingQual,
103
+ ...rest,
104
+ },
105
+ };
106
+ }
107
+ export function flipCigar(cigar) {
108
+ const arr = [];
109
+ for (let i = cigar.length - 2; i >= 0; i -= 2) {
110
+ arr.push(cigar[i]);
111
+ const op = cigar[i + 1];
112
+ if (op === 'D') {
113
+ arr.push('I');
114
+ }
115
+ else if (op === 'I') {
116
+ arr.push('D');
117
+ }
118
+ else {
119
+ arr.push(op);
120
+ }
121
+ }
122
+ return arr;
123
+ }
124
+ export function swapIndelCigar(cigar) {
125
+ return cigar.replaceAll('D', 'K').replaceAll('I', 'D').replaceAll('K', 'I');
126
+ }
127
+ //# sourceMappingURL=util.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"util.js","sourceRoot":"","sources":["../../src/PAFAdapter/util.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,GAAG,EAAE,MAAM,SAAS,CAAA;AAkB7B,6EAA6E;AAC7E,0BAA0B;AAC1B,EAAE;AACF,cAAc;AAEd,iCAAiC;AAEjC,+EAA+E;AAC/E,gFAAgF;AAChF,+EAA+E;AAC/E,4EAA4E;AAC5E,wEAAwE;AACxE,2DAA2D;AAE3D,iFAAiF;AACjF,kDAAkD;AAElD,6EAA6E;AAC7E,2EAA2E;AAC3E,8EAA8E;AAC9E,yEAAyE;AACzE,gFAAgF;AAChF,gFAAgF;AAChF,YAAY;AACZ,EAAE;AACF,8EAA8E;AAC9E,8DAA8D;AAC9D,EAAE;AACF,+EAA+E;AAC/E,2EAA2E;AAC3E,qBAAqB;AACrB,EAAE;AACF,8EAA8E;AAC9E,+EAA+E;AAC/E,yEAAyE;AACzE,6EAA6E;AAC7E,4EAA4E;AAC5E,aAAa;AAEb,MAAM,UAAU,gBAAgB,CAAC,GAAgB;IAC/C,MAAM,QAAQ,GAA0D,EAAE,CAAA;IAC1E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACnC,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAA;QACpB,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAA;QACzB,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAA;QAC1B,MAAM,GAAG,GAAG,KAAK,GAAG,GAAG,GAAG,MAAM,CAAA;QAChC,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE;YAClB,QAAQ,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,CAAA;SACvC;QACD,QAAQ,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,WAAW,CAAC,CAAA;QACjD,QAAQ,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,QAAQ,IAAI,CAAC,CAAC,CAAA;KAClD;IAED,MAAM,YAAY,GAAG,MAAM,CAAC,WAAW,CACrC,MAAM,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE;QAC1C,MAAM,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,GAAG,CAAC,CAAA;QACpC,OAAO,CAAC,GAAG,EAAE,YAAY,CAAC,IAAI,CAAC,CAAC,CAAA;IAClC,CAAC,CAAC,CACH,CAAA;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACnC,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAA;QACpB,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAA;QACzB,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAA;QAC1B,MAAM,GAAG,GAAG,KAAK,GAAG,GAAG,GAAG,MAAM,CAAA;QAChC,KAAK,CAAC,KAAK,CAAC,SAAS,GAAG,YAAY,CAAC,GAAG,CAAC,CAAA;KAC1C;IAED,IAAI,GAAG,GAAG,KAAK,CAAA;IACf,IAAI,GAAG,GAAG,CAAC,CAAA;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACnC,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAA;QACpB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,SAAS,IAAI,CAAC,EAAE,GAAG,CAAC,CAAA;QAC/C,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,SAAS,IAAI,CAAC,EAAE,GAAG,CAAC,CAAA;KAChD;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACnC,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAA;QACpB,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,SAAS,IAAI,CAAC,CAAA;QACpC,KAAK,CAAC,KAAK,CAAC,SAAS,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAA;KAChD;IAED,OAAO,GAAG,CAAA;AACZ,CAAC;AAED,kEAAkE;AAClE,SAAS,YAAY,CAAC,MAA0B;IAC9C,mDAAmD;IACnD,MAAM,CAAC,QAAQ,EAAE,SAAS,CAAC,GAAG,MAAM,CAAC,MAAM,CACzC,CAAC,CAAC,QAAQ,EAAE,SAAS,CAAC,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC,EAAE,EAAE,CAAC;QAC1C,QAAQ,GAAG,KAAK,GAAG,MAAM;QACzB,SAAS,GAAG,MAAM;KACnB,EACD,CAAC,CAAC,EAAE,CAAC,CAAC,CACP,CAAA;IACD,OAAO,QAAQ,GAAG,SAAS,CAAA;AAC7B,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,IAAY;IACvC,MAAM,CACJ,KAAK,EACL,AADM,EAEN,MAAM,EACN,IAAI,EACJ,MAAM,EACN,KAAK,EACL,AADM,EAEN,MAAM,EACN,IAAI,EACJ,UAAU,EACV,QAAQ,EACR,WAAW,EACX,GAAG,MAAM,CACV,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAA;IAEpB,MAAM,IAAI,GAAG,MAAM,CAAC,WAAW,CAC7B,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE;QACjB,MAAM,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAA;QAC5B,MAAM,SAAS,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;QACnC,MAAM,UAAU,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACrC,OAAO,CAAC,SAAS,EAAE,UAAU,CAAC,CAAA;IAChC,CAAC,CAAC,CACH,CAAA;IAED,OAAO;QACL,KAAK;QACL,MAAM,EAAE,CAAC,MAAM;QACf,IAAI,EAAE,CAAC,IAAI;QACX,KAAK;QACL,MAAM,EAAE,CAAC,MAAM;QACf,IAAI,EAAE,CAAC,IAAI;QACX,MAAM,EAAE,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/B,KAAK,EAAE;YACL,UAAU,EAAE,CAAC,UAAU;YACvB,QAAQ,EAAE,CAAC,QAAQ;YACnB,WAAW,EAAE,CAAC,WAAW;YACzB,GAAG,IAAI;SACR;KACW,CAAA;AAChB,CAAC;AAED,MAAM,UAAU,SAAS,CAAC,KAAe;IACvC,MAAM,GAAG,GAAG,EAAE,CAAA;IACd,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE;QAC7C,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAA;QAClB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACvB,IAAI,EAAE,KAAK,GAAG,EAAE;YACd,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAA;SACd;aAAM,IAAI,EAAE,KAAK,GAAG,EAAE;YACrB,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAA;SACd;aAAM;YACL,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAA;SACb;KACF;IACD,OAAO,GAAG,CAAA;AACZ,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,KAAa;IAC1C,OAAO,KAAK,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,CAAA;AAC7E,CAAC"}
package/esm/index.d.ts ADDED
@@ -0,0 +1,6 @@
1
+ import Plugin from '@jbrowse/core/Plugin';
2
+ import PluginManager from '@jbrowse/core/PluginManager';
3
+ export default class ComparativeAdaptersPlugin extends Plugin {
4
+ name: string;
5
+ install(pluginManager: PluginManager): void;
6
+ }
package/esm/index.js ADDED
@@ -0,0 +1,45 @@
1
+ import Plugin from '@jbrowse/core/Plugin';
2
+ import PAFAdapterF from './PAFAdapter';
3
+ import MCScanAnchorsAdapterF from './MCScanAnchorsAdapter';
4
+ import MCScanSimpleAnchorsAdapterF from './MCScanSimpleAnchorsAdapter';
5
+ import MashMapAdapterF from './MashMapAdapter';
6
+ import DeltaAdapterF from './DeltaAdapter';
7
+ import ChainAdapterF from './ChainAdapter';
8
+ import { getFileName, } from '@jbrowse/core/util/tracks';
9
+ export default class ComparativeAdaptersPlugin extends Plugin {
10
+ constructor() {
11
+ super(...arguments);
12
+ this.name = 'ComparativeAdaptersPlugin';
13
+ }
14
+ install(pluginManager) {
15
+ PAFAdapterF(pluginManager);
16
+ DeltaAdapterF(pluginManager);
17
+ ChainAdapterF(pluginManager);
18
+ MCScanAnchorsAdapterF(pluginManager);
19
+ MCScanSimpleAnchorsAdapterF(pluginManager);
20
+ MashMapAdapterF(pluginManager);
21
+ pluginManager.addToExtensionPoint('Core-guessAdapterForLocation', (adapterGuesser) => {
22
+ return (file, index, adapterHint) => {
23
+ const regexGuess = /\.paf/i;
24
+ const adapterName = 'PAFAdapter';
25
+ const fileName = getFileName(file);
26
+ if (regexGuess.test(fileName) || adapterHint === adapterName) {
27
+ return {
28
+ type: adapterName,
29
+ pafLocation: file,
30
+ };
31
+ }
32
+ return adapterGuesser(file, index, adapterHint);
33
+ };
34
+ });
35
+ pluginManager.addToExtensionPoint('Core-guessTrackTypeForLocation', (trackTypeGuesser) => {
36
+ return (adapterName) => {
37
+ if (adapterName === 'PAFAdapter') {
38
+ return 'SyntenyTrack';
39
+ }
40
+ return trackTypeGuesser(adapterName);
41
+ };
42
+ });
43
+ }
44
+ }
45
+ //# sourceMappingURL=index.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,sBAAsB,CAAA;AAIzC,OAAO,WAAW,MAAM,cAAc,CAAA;AACtC,OAAO,qBAAqB,MAAM,wBAAwB,CAAA;AAC1D,OAAO,2BAA2B,MAAM,8BAA8B,CAAA;AACtE,OAAO,eAAe,MAAM,kBAAkB,CAAA;AAC9C,OAAO,aAAa,MAAM,gBAAgB,CAAA;AAC1C,OAAO,aAAa,MAAM,gBAAgB,CAAA;AAE1C,OAAO,EACL,WAAW,GAGZ,MAAM,2BAA2B,CAAA;AAElC,MAAM,CAAC,OAAO,OAAO,yBAA0B,SAAQ,MAAM;IAA7D;;QACE,SAAI,GAAG,2BAA2B,CAAA;IA2CpC,CAAC;IAzCC,OAAO,CAAC,aAA4B;QAClC,WAAW,CAAC,aAAa,CAAC,CAAA;QAC1B,aAAa,CAAC,aAAa,CAAC,CAAA;QAC5B,aAAa,CAAC,aAAa,CAAC,CAAA;QAC5B,qBAAqB,CAAC,aAAa,CAAC,CAAA;QACpC,2BAA2B,CAAC,aAAa,CAAC,CAAA;QAC1C,eAAe,CAAC,aAAa,CAAC,CAAA;QAE9B,aAAa,CAAC,mBAAmB,CAC/B,8BAA8B,EAC9B,CAAC,cAA8B,EAAE,EAAE;YACjC,OAAO,CACL,IAAkB,EAClB,KAAoB,EACpB,WAAoB,EACpB,EAAE;gBACF,MAAM,UAAU,GAAG,QAAQ,CAAA;gBAC3B,MAAM,WAAW,GAAG,YAAY,CAAA;gBAChC,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,CAAC,CAAA;gBAClC,IAAI,UAAU,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,WAAW,KAAK,WAAW,EAAE;oBAC5D,OAAO;wBACL,IAAI,EAAE,WAAW;wBACjB,WAAW,EAAE,IAAI;qBAClB,CAAA;iBACF;gBACD,OAAO,cAAc,CAAC,IAAI,EAAE,KAAK,EAAE,WAAW,CAAC,CAAA;YACjD,CAAC,CAAA;QACH,CAAC,CACF,CAAA;QACD,aAAa,CAAC,mBAAmB,CAC/B,gCAAgC,EAChC,CAAC,gBAAkC,EAAE,EAAE;YACrC,OAAO,CAAC,WAAmB,EAAE,EAAE;gBAC7B,IAAI,WAAW,KAAK,YAAY,EAAE;oBAChC,OAAO,cAAc,CAAA;iBACtB;gBACD,OAAO,gBAAgB,CAAC,WAAW,CAAC,CAAA;YACtC,CAAC,CAAA;QACH,CAAC,CACF,CAAA;IACH,CAAC;CACF"}
package/esm/util.d.ts ADDED
@@ -0,0 +1,16 @@
1
+ /// <reference types="node" />
2
+ import { BaseOptions } from '@jbrowse/core/data_adapters/BaseAdapter';
3
+ import { GenericFilehandle } from 'generic-filehandle';
4
+ import { PAFRecord } from './PAFAdapter/util';
5
+ export declare function isGzip(buf: Buffer): boolean;
6
+ export declare function parseBed(text: string): Map<string, {
7
+ refName: string;
8
+ start: number;
9
+ end: number;
10
+ score: number;
11
+ name: string;
12
+ strand: number;
13
+ }>;
14
+ export declare function readFile(file: GenericFilehandle, opts?: BaseOptions): Promise<string>;
15
+ export declare function zip(a: number[], b: number[]): [number, number][];
16
+ export declare function parseLineByLine(buffer: Buffer, cb: (line: string) => PAFRecord): PAFRecord[];
package/esm/util.js ADDED
@@ -0,0 +1,49 @@
1
+ import { unzip } from '@gmod/bgzf-filehandle';
2
+ export function isGzip(buf) {
3
+ return buf[0] === 31 && buf[1] === 139 && buf[2] === 8;
4
+ }
5
+ export function parseBed(text) {
6
+ return new Map(text
7
+ .split(/\n|\r\n|\r/)
8
+ .filter(f => !!f || f.startsWith('#'))
9
+ .map(line => {
10
+ const [refName, start, end, name, score, strand] = line.split('\t');
11
+ return [
12
+ name,
13
+ {
14
+ refName,
15
+ start: +start,
16
+ end: +end,
17
+ score: +score,
18
+ name,
19
+ strand: strand === '-' ? -1 : 1,
20
+ },
21
+ ];
22
+ }));
23
+ }
24
+ export async function readFile(file, opts) {
25
+ const buffer = (await file.readFile(opts));
26
+ return new TextDecoder('utf8', { fatal: true }).decode(isGzip(buffer) ? await unzip(buffer) : buffer);
27
+ }
28
+ export function zip(a, b) {
29
+ return a.map((e, i) => [e, b[i]]);
30
+ }
31
+ const decoder = typeof TextDecoder !== 'undefined' ? new TextDecoder('utf8') : undefined;
32
+ export function parseLineByLine(buffer, cb) {
33
+ let blockStart = 0;
34
+ const entries = [];
35
+ while (blockStart < buffer.length) {
36
+ const n = buffer.indexOf('\n', blockStart);
37
+ if (n === -1) {
38
+ break;
39
+ }
40
+ const b = buffer.slice(blockStart, n);
41
+ const line = ((decoder === null || decoder === void 0 ? void 0 : decoder.decode(b)) || b.toString()).trim();
42
+ if (line) {
43
+ entries.push(cb(line));
44
+ }
45
+ blockStart = n + 1;
46
+ }
47
+ return entries;
48
+ }
49
+ //# sourceMappingURL=util.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"util.js","sourceRoot":"","sources":["../src/util.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,KAAK,EAAE,MAAM,uBAAuB,CAAA;AAG7C,MAAM,UAAU,MAAM,CAAC,GAAW;IAChC,OAAO,GAAG,CAAC,CAAC,CAAC,KAAK,EAAE,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAA;AACxD,CAAC;AAED,MAAM,UAAU,QAAQ,CAAC,IAAY;IACnC,OAAO,IAAI,GAAG,CACZ,IAAI;SACD,KAAK,CAAC,YAAY,CAAC;SACnB,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;SACrC,GAAG,CAAC,IAAI,CAAC,EAAE;QACV,MAAM,CAAC,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAA;QACnE,OAAO;YACL,IAAI;YACJ;gBACE,OAAO;gBACP,KAAK,EAAE,CAAC,KAAK;gBACb,GAAG,EAAE,CAAC,GAAG;gBACT,KAAK,EAAE,CAAC,KAAK;gBACb,IAAI;gBACJ,MAAM,EAAE,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aAChC;SACF,CAAA;IACH,CAAC,CAAC,CACL,CAAA;AACH,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,QAAQ,CAAC,IAAuB,EAAE,IAAkB;IACxE,MAAM,MAAM,GAAG,CAAC,MAAM,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAW,CAAA;IACpD,OAAO,IAAI,WAAW,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,CAAC,CAAC,MAAM,CACpD,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAC9C,CAAA;AACH,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAW,EAAE,CAAW;IAC1C,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAqB,CAAC,CAAA;AACvD,CAAC;AAED,MAAM,OAAO,GACX,OAAO,WAAW,KAAK,WAAW,CAAC,CAAC,CAAC,IAAI,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,SAAS,CAAA;AAE1E,MAAM,UAAU,eAAe,CAC7B,MAAc,EACd,EAA+B;IAE/B,IAAI,UAAU,GAAG,CAAC,CAAA;IAClB,MAAM,OAAO,GAAG,EAAE,CAAA;IAClB,OAAO,UAAU,GAAG,MAAM,CAAC,MAAM,EAAE;QACjC,MAAM,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC,IAAI,EAAE,UAAU,CAAC,CAAA;QAC1C,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE;YACZ,MAAK;SACN;QACD,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,CAAA;QACrC,MAAM,IAAI,GAAG,CAAC,CAAA,OAAO,aAAP,OAAO,uBAAP,OAAO,CAAE,MAAM,CAAC,CAAC,CAAC,KAAI,CAAC,CAAC,QAAQ,EAAE,CAAC,CAAC,IAAI,EAAE,CAAA;QACxD,IAAI,IAAI,EAAE;YACR,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAA;SACvB;QAED,UAAU,GAAG,CAAC,GAAG,CAAC,CAAA;KACnB;IACD,OAAO,OAAO,CAAA;AAChB,CAAC"}
package/package.json ADDED
@@ -0,0 +1,63 @@
1
+ {
2
+ "name": "@jbrowse/plugin-comparative-adapters",
3
+ "version": "2.6.1",
4
+ "description": "JBrowse 2 comparative adapters",
5
+ "keywords": [
6
+ "jbrowse",
7
+ "jbrowse2"
8
+ ],
9
+ "license": "Apache-2.0",
10
+ "homepage": "https://jbrowse.org",
11
+ "bugs": "https://github.com/GMOD/jbrowse-components/issues",
12
+ "repository": {
13
+ "type": "git",
14
+ "url": "https://github.com/GMOD/jbrowse-components.git",
15
+ "directory": "plugins/comparative-adapters"
16
+ },
17
+ "author": "JBrowse Team",
18
+ "distMain": "dist/index.js",
19
+ "srcMain": "src/index.ts",
20
+ "main": "dist/index.js",
21
+ "files": [
22
+ "dist",
23
+ "src",
24
+ "esm"
25
+ ],
26
+ "scripts": {
27
+ "build": "npm-run-all build:*",
28
+ "test": "cd ../..; jest plugins/comparative-adapters",
29
+ "prepublishOnly": "yarn test",
30
+ "prepack": "yarn build && yarn useDist",
31
+ "postpack": "yarn useSrc",
32
+ "useDist": "node ../../scripts/useDist.js",
33
+ "useSrc": "node ../../scripts/useSrc.js",
34
+ "prebuild": "npm run clean",
35
+ "build:esm": "tsc --build tsconfig.build.esm.json",
36
+ "build:es5": "tsc --build tsconfig.build.es5.json",
37
+ "clean": "rimraf dist esm *.tsbuildinfo"
38
+ },
39
+ "dependencies": {
40
+ "@gmod/bgzf-filehandle": "^1.4.3",
41
+ "generic-filehandle": "^3.0.0"
42
+ },
43
+ "peerDependencies": {
44
+ "@jbrowse/core": "^2.0.0",
45
+ "@jbrowse/plugin-alignments": "^2.0.0",
46
+ "@jbrowse/plugin-linear-genome-view": "^2.0.0",
47
+ "@mui/material": "^5.0.0",
48
+ "mobx": "^6.0.0",
49
+ "mobx-react": "^7.0.0",
50
+ "mobx-state-tree": "^5.0.0",
51
+ "prop-types": "^15.0.0",
52
+ "react": ">=16.8.0",
53
+ "react-dom": ">=16.8.0",
54
+ "rxjs": "^7.0.0"
55
+ },
56
+ "distModule": "esm/index.js",
57
+ "srcModule": "src/index.ts",
58
+ "module": "esm/index.js",
59
+ "publishConfig": {
60
+ "access": "public"
61
+ },
62
+ "gitHead": "1cbe7ba097fb2d2763c776e5e429e4670cdd583c"
63
+ }