Skip to content

Commit c9abafd

Browse files
committed
feat: add influenza-a diamond db
1 parent 22f7cec commit c9abafd

3 files changed

Lines changed: 65 additions & 7 deletions

File tree

README.md

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -12,18 +12,18 @@ Then create a config.yaml file (like the example) with the list of all INSDC ref
1212

1313
```yaml
1414
references:
15-
seg4_H1: U08903.1
16-
seg4_H2: CY005413.1
17-
seg3__H5N1: NC_007359.1
15+
seg4-H1: U08903.1
16+
seg4-H2: CY005413.1
17+
seg3--H5N1: NC_007359.1
1818
```
1919
20-
Note that the key e.g. `seg4_H1` will be the prefix of the `dataset` name that is returned in the diamond blastx results.tsv. (As DIAMOND matches proteins each CDS in the sequence receives its own identifier, in Loculus we map all sequences that match the protein `id|CDS{i}` to the sequence `id`. ). The keys should be the same as your sequenceName they are assigned to e.g. `{segment}-{reference}` or alternatively, you can add lists of accepted matches to the `config.accepted_dataset_matches` field.
20+
Note that the key e.g. `seg4-H1` will be the prefix of the `dataset` name that is returned in the diamond blastx results.tsv. (As DIAMOND matches proteins each CDS in the sequence receives its own identifier, in Loculus we map all sequences that match the protein `id|CDS{i}` to the sequence `id`. ). The keys should be the same as your sequenceName they are assigned to e.g. `{segment}-{reference}` or alternatively, you can add lists of accepted matches to the `config.accepted_dataset_matches` field.
2121

2222
```tsv
2323
seqName dataset pident ...
24-
MW874350.1 seg3__H5N1|CDS1 0.4829120176662018 ...
25-
MW874350.1 seg3__H1N1|CDS2 0.34937439846005774 ...
26-
MW874350.1 seg3__H3N2|CDS1 0.22552301255230126 ...
24+
MW874350.1 seg3-H5N1|CDS1 0.4829120176662018 ...
25+
MW874350.1 seg3-H1N1|CDS2 0.34937439846005774 ...
26+
MW874350.1 seg3-H3N2|CDS1 0.22552301255230126 ...
2727
```
2828

2929
```bash

results/influenza-a/config.yaml

Lines changed: 58 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,58 @@
1+
references:
2+
seg1-h1n1pdm: NC_002205.1
3+
seg1-h1n1: NC_002023.1
4+
seg1-h3n2: NC_007373.1
5+
seg1-h2n2: NC_007378.1
6+
seg2-h1n1pdm: NC_026435.1
7+
seg2-h1n1: NC_002021.1
8+
seg2-h3n2: NC_007372.1
9+
seg2-h2n2: NC_007375.1
10+
seg3-h1n1pdm: NC_026437.1
11+
seg3-h1n1: NC_002022.1
12+
seg3-h3n2: NC_007371.1
13+
seg3-h2n2: NC_007376.1
14+
seg4-h3_h3n8: CY028836
15+
seg4-h4_h4n6: CY181241
16+
seg4-h5_h5n1: NC_007362.1
17+
seg4-h5_h5n2: KU143256
18+
seg4-h6_h6n2: CY130030
19+
seg4-h7_h7n9: NC_026425.1
20+
seg4-h8_h8n4: CY136131
21+
seg4-h9_h9n2: NC_004908.1
22+
seg4-h10_h10n7: CY136094
23+
seg4-h11_h11n9: CY130070
24+
seg4-h12_h12n5: CY130078
25+
seg4-h13_h13n6: CY130086
26+
seg4-h14_h14n5: JN696314
27+
seg4-h15_h15n9: CY006010
28+
seg4-h16_h16n3: CY136630
29+
seg4-h17_h17n10: CY103876
30+
seg4-h18_h18n11: CY125945
31+
seg5-h1n1pdm: NC_026436.1
32+
seg5-h1n1: NC_002019.1
33+
seg5-h3n2: NC_007369.1
34+
seg5-h2n2: NC_007378.1
35+
seg6-n1_h5n1: NC_007361.1
36+
seg6-n2_h6n2: CY130032
37+
seg6-n2_h9n2: NC_004909.1
38+
seg6-n2_h5n2: KU143347
39+
seg6-n3_h16n3: CY136632
40+
seg6-n4_h8n4: CY136133
41+
seg6-n5_h12n5: CY130080
42+
seg6-n6_h4n6: CY181243
43+
seg6-n6_h13n6: CY130088
44+
seg6-n7_h10n7: CY136096
45+
seg6-n8_h3n8: CY028838
46+
seg6-n9_h7n9: NC_026429.1
47+
seg6-n9_h11n9: CY130072
48+
seg6-n9_h15n9: CY005407
49+
seg6-n10_h17n10: CY103878
50+
seg6-n11_h18n11: CY125947
51+
seg7-h1n1pdm: NC_026431.1
52+
seg7-h1n1: NC_002016.1
53+
seg7-h3n2: NC_007367.1
54+
seg7-h2n2: NC_007377.1
55+
seg8-h1n1pdm: NC_026432.1
56+
seg8-h1n1: NC_002020.1
57+
seg8-h3n2: NC_007370.1
58+
seg8-h2n2: NC_007380.1

results/influenza-a/diamond.dmnd

34.1 KB
Binary file not shown.

0 commit comments

Comments
 (0)