tseemann / nullarbor

:floppy_disk: :page_with_curl: "Reads to report" for public health and clinical microbiology
GNU General Public License v2.0
134 stars 37 forks source link

Core.aln all the same nucleotide? #252

Closed foamy1881 closed 4 years ago

foamy1881 commented 4 years ago

Hi finished running nullarbor2 and when I looked at the core SNP and distance matrix output most of them are 0 and 1. I read numerous post saying please removed unaligned >200,000 bp and snippy must run on isolates of more than 10x coverage. I check those and it seem fine and I manage to run snippy by itself but it still produces the same outcome. Can someone point my mistake?

This is snippy-core 4.4.5
Obtained from http://github.com/tseemann/snippy
Enabling bundled tools.
Found any2fasta - /home/yfong/miniconda3/envs/snippy/binaries/noarch/any2fasta
Found samtools - /home/yfong/miniconda3/envs/snippy/bin/samtools
Found minimap2 - /home/yfong/miniconda3/envs/snippy/bin/minimap2
Found bedtools - /home/yfong/miniconda3/envs/snippy/bin/bedtools
Found snp-sites - /home/yfong/miniconda3/envs/snippy/bin/snp-sites
Saving reference FASTA: core.ref.fa
This is any2fasta 0.4.2
Opening 'M001/ref.fa'
Detected FASTA format
Read 73527 lines from 'M001/ref.fa'
Wrote 1 sequences from FASTA file.
Processed 1 files.
Done.
Loaded 1 sequences totalling 4411532 bp.
Will mask 0 regions totalling 0 bp ~ 0.00%
0   M001    snp=622 del=76  ins=50  het=108 unaligned=133260
1   M004    snp=1420    del=231 ins=205 het=395 unaligned=99174
2   M005    snp=1343    del=168 ins=149 het=267 unaligned=110619
3   M006    snp=1448    del=421 ins=362 het=617 unaligned=86679
4   M007    snp=1391    del=198 ins=178 het=352 unaligned=105004
5   M008    snp=1304    del=154 ins=136 het=446 unaligned=118573
6   M009    snp=1332    del=195 ins=185 het=345 unaligned=114121
7   M011    snp=2062    del=668 ins=337 het=821 unaligned=86549
8   M012    snp=1353    del=195 ins=165 het=456 unaligned=105708
9   M013    snp=1455    del=279 ins=256 het=979 unaligned=96660
10  M014    snp=1015    del=90  ins=58  het=161 unaligned=118938
11  M016    snp=2009    del=617 ins=452 het=750 unaligned=91697
12  M017    snp=1046    del=175 ins=137 het=1893    unaligned=94427
13  M018    snp=1214    del=136 ins=130 het=613 unaligned=102553
14  M019    snp=1344    del=165 ins=139 het=529 unaligned=104828
15  M020    snp=1331    del=294 ins=359 het=1183    unaligned=86500
16  M023    snp=1309    del=149 ins=139 het=482 unaligned=117782
17  M024    snp=880 del=111 ins=74  het=724 unaligned=121227
18  M025    snp=1359    del=285 ins=414 het=988 unaligned=62649
19  M026    snp=1991    del=358 ins=243 het=388 unaligned=97976
20  M027    snp=686 del=71  ins=125 het=3782    unaligned=70180
21  M028    snp=844 del=82  ins=45  het=51  unaligned=117821
22  M032    snp=924 del=55  ins=60  het=42  unaligned=117118
23  M036    snp=1494    del=624 ins=451 het=771 unaligned=82055
24  M045    snp=277 del=4   ins=5   het=143 unaligned=112022
25  M047    snp=301 del=18  ins=13  het=170 unaligned=78764
26  M050    snp=1414    del=264 ins=227 het=305 unaligned=103399
27  M051    snp=2029    del=331 ins=307 het=540 unaligned=94108
28  M052    snp=1390    del=203 ins=198 het=741 unaligned=107805
29  M053    snp=1431    del=420 ins=418 het=806 unaligned=87381
30  M056    snp=1396    del=250 ins=275 het=636 unaligned=97441
31  M057    snp=1467    del=539 ins=374 het=689 unaligned=87344
32  M058    snp=936 del=177 ins=140 het=553 unaligned=84978
33  M059    snp=1988    del=242 ins=237 het=461 unaligned=107036
34  M060    snp=1448    del=382 ins=344 het=552 unaligned=91235
35  M061    snp=1362    del=175 ins=157 het=689 unaligned=98914
36  M063    snp=2069    del=386 ins=406 het=860 unaligned=93972
37  M064    snp=1405    del=274 ins=341 het=544 unaligned=98250
38  M065    snp=1970    del=263 ins=280 het=734 unaligned=96985
39  M066    snp=1471    del=296 ins=366 het=1215    unaligned=101900
40  M067    snp=1406    del=263 ins=280 het=434 unaligned=99535
41  M069    snp=788 del=159 ins=114 het=540 unaligned=73572
42  M070    snp=1293    del=143 ins=149 het=343 unaligned=91754
43  M072    snp=1457    del=334 ins=400 het=1193    unaligned=89888
44  M077    snp=937 del=94  ins=104 het=576 unaligned=102341
45  M078    snp=756 del=78  ins=88  het=250 unaligned=76888
46  M081    snp=1436    del=374 ins=331 het=322 unaligned=91983
47  M082    snp=1322    del=222 ins=172 het=405 unaligned=90869
48  M083    snp=1434    del=390 ins=316 het=1236    unaligned=96484
49  M084    snp=1983    del=280 ins=305 het=1123    unaligned=76589
50  M085    snp=1444    del=251 ins=230 het=1945    unaligned=93179
51  M086    snp=1948    del=261 ins=276 het=605 unaligned=92812
52  M087    snp=1355    del=154 ins=149 het=223 unaligned=101555
53  M088    snp=734 del=75  ins=72  het=182 unaligned=82023
54  M089    snp=676 del=87  ins=132 het=2644    unaligned=72026
55  M090    snp=1877    del=258 ins=256 het=333 unaligned=96738
56  M094    snp=1428    del=265 ins=346 het=836 unaligned=100277
57  M097    snp=1928    del=314 ins=245 het=392 unaligned=90026
58  M098    snp=1379    del=144 ins=114 het=389 unaligned=118076
59  M099    snp=1198    del=130 ins=139 het=179 unaligned=112555
60  M100    snp=1185    del=178 ins=110 het=279 unaligned=114223
61  M101    snp=1047    del=108 ins=83  het=636 unaligned=106960
62  M103    snp=1386    del=207 ins=138 het=366 unaligned=110462
63  M104    snp=1447    del=460 ins=427 het=921 unaligned=89684
64  M105    snp=1356    del=219 ins=122 het=208 unaligned=109469
65  M106    snp=1991    del=183 ins=314 het=325 unaligned=100987
66  M107    snp=1347    del=187 ins=139 het=339 unaligned=119335
67  M108    snp=1327    del=170 ins=123 het=246 unaligned=111021
68  M109    snp=1372    del=183 ins=192 het=284 unaligned=104246
69  M110    snp=775 del=70  ins=96  het=222 unaligned=79096
70  M111    snp=1322    del=136 ins=131 het=317 unaligned=105293
71  M113    snp=1368    del=193 ins=164 het=599 unaligned=102652
72  M114    snp=1276    del=152 ins=121 het=357 unaligned=98546
73  M116    snp=1296    del=164 ins=133 het=436 unaligned=120980
74  M117    snp=1398    del=248 ins=215 het=426 unaligned=103640
75  M120    snp=1390    del=210 ins=179 het=434 unaligned=100949
76  M122    snp=1358    del=190 ins=160 het=386 unaligned=111393
77  M123    snp=1362    del=166 ins=178 het=360 unaligned=107263
78  M124    snp=1324    del=166 ins=139 het=523 unaligned=121729
79  M125    snp=1349    del=185 ins=119 het=238 unaligned=121240
80  M126    snp=1367    del=188 ins=192 het=467 unaligned=111955
81  M127    snp=1263    del=205 ins=159 het=899 unaligned=99601
82  M128    snp=1332    del=233 ins=137 het=406 unaligned=127194
83  M130    snp=1389    del=188 ins=142 het=399 unaligned=113170
84  M131    snp=1248    del=157 ins=118 het=665 unaligned=102704
85  M132    snp=1318    del=190 ins=127 het=328 unaligned=112309
86  M133    snp=1384    del=220 ins=170 het=352 unaligned=122480
87  M134    snp=1978    del=307 ins=259 het=621 unaligned=88984
88  M135    snp=728 del=129 ins=159 het=1009    unaligned=69485
89  M136    snp=297 del=22  ins=27  het=1801    unaligned=68391
90  M137    snp=2018    del=324 ins=356 het=1458    unaligned=85552
91  M138    snp=1949    del=276 ins=249 het=773 unaligned=88086
92  M139    snp=1988    del=309 ins=272 het=575 unaligned=94038
93  M141    snp=1393    del=229 ins=178 het=283 unaligned=112047
94  M142    snp=1910    del=240 ins=237 het=715 unaligned=90220
95  M143    snp=1997    del=264 ins=261 het=414 unaligned=86079
96  M144    snp=1278    del=164 ins=122 het=254 unaligned=129216
97  M146    snp=1380    del=171 ins=201 het=466 unaligned=102088
98  M147    snp=1279    del=150 ins=117 het=190 unaligned=112690
99  M148    snp=1366    del=232 ins=165 het=312 unaligned=106273
100 M149    snp=1399    del=296 ins=230 het=382 unaligned=101541
101 M150    snp=1331    del=165 ins=143 het=477 unaligned=105291
102 M152    snp=1327    del=240 ins=203 het=588 unaligned=107936
103 M153    snp=2053    del=385 ins=423 het=688 unaligned=79050
104 M154    snp=1367    del=203 ins=215 het=249 unaligned=105392
105 M155    snp=1354    del=176 ins=117 het=728 unaligned=107582
106 M157    snp=1347    del=209 ins=149 het=361 unaligned=107663
107 M158    snp=1182    del=212 ins=107 het=256 unaligned=111772
108 M159    snp=1377    del=246 ins=146 het=379 unaligned=106873
109 M160    snp=1287    del=150 ins=107 het=208 unaligned=120343
110 M161    snp=1384    del=192 ins=190 het=151 unaligned=105985
111 M162    snp=1162    del=88  ins=63  het=94  unaligned=115239
112 M163    snp=1375    del=152 ins=117 het=357 unaligned=114957
113 M164    snp=1365    del=153 ins=118 het=251 unaligned=114578
114 M165    snp=1394    del=174 ins=137 het=332 unaligned=118497
115 M166    snp=1956    del=262 ins=225 het=467 unaligned=85211
116 M167    snp=1386    del=151 ins=128 het=370 unaligned=108039
117 M168    snp=1314    del=95  ins=100 het=258 unaligned=115305
118 M169    snp=1368    del=145 ins=122 het=333 unaligned=108055
119 M170    snp=1322    del=146 ins=101 het=198 unaligned=84492
120 M171    snp=1332    del=143 ins=89  het=180 unaligned=110865
121 M172    snp=1199    del=94  ins=67  het=200 unaligned=120312
122 M174    snp=1371    del=120 ins=113 het=202 unaligned=112145
123 M175    snp=1252    del=118 ins=53  het=174 unaligned=111925
124 M176    snp=1340    del=128 ins=106 het=199 unaligned=111180
125 M177    snp=912 del=123 ins=83  het=179 unaligned=97489
126 M178    snp=1295    del=141 ins=108 het=148 unaligned=114522
127 M179    snp=2019    del=184 ins=189 het=464 unaligned=107556
128 M180    snp=1274    del=133 ins=100 het=349 unaligned=106539
129 M181    snp=1393    del=145 ins=130 het=325 unaligned=112241
130 M182    snp=1355    del=186 ins=134 het=802 unaligned=108332
131 M184    snp=1410    del=173 ins=160 het=414 unaligned=109908
132 M185    snp=1410    del=153 ins=127 het=470 unaligned=106679
133 M186    snp=1331    del=110 ins=76  het=221 unaligned=109636
134 M187    snp=2007    del=180 ins=178 het=395 unaligned=106899
135 M188    snp=1839    del=127 ins=114 het=283 unaligned=104542
136 M190    snp=1188    del=102 ins=103 het=243 unaligned=120401
137 M191    snp=904 del=54  ins=50  het=263 unaligned=112981
138 M192    snp=1364    del=184 ins=161 het=318 unaligned=113868
139 M193    snp=1999    del=256 ins=227 het=422 unaligned=103993
140 M195    snp=1373    del=155 ins=126 het=246 unaligned=111584
141 M196    snp=717 del=30  ins=114 het=4018    unaligned=76032
142 M197    snp=1140    del=80  ins=62  het=170 unaligned=125148
143 M198    snp=485 del=21  ins=25  het=790 unaligned=83400
144 M199    snp=785 del=73  ins=64  het=131 unaligned=85827
145 M200    snp=1290    del=118 ins=58  het=144 unaligned=119710
146 M201    snp=1758    del=119 ins=96  het=269 unaligned=104199
147 M202    snp=1769    del=90  ins=107 het=131 unaligned=105284
148 M203    snp=1346    del=143 ins=99  het=260 unaligned=118190
149 M204    snp=756 del=55  ins=32  het=128 unaligned=86818
150 M205    snp=1372    del=139 ins=158 het=248 unaligned=112983
151 M206    snp=1296    del=102 ins=117 het=242 unaligned=110550
152 M207    snp=1397    del=125 ins=150 het=406 unaligned=110273
153 M208    snp=1073    del=62  ins=78  het=147 unaligned=115979
154 M209    snp=690 del=38  ins=54  het=925 unaligned=86578
155 M210    snp=530 del=17  ins=41  het=535 unaligned=90229
156 M211    snp=1117    del=122 ins=56  het=259 unaligned=113822
157 M220    snp=1351    del=224 ins=211 het=277 unaligned=117030
158 M246    snp=1469    del=263 ins=202 het=524 unaligned=98123
159 M247    snp=1448    del=323 ins=253 het=268 unaligned=93023
160 M248    snp=1325    del=229 ins=249 het=1088    unaligned=89036
161 M249    snp=1476    del=257 ins=298 het=2720    unaligned=100588
162 M250    snp=1450    del=236 ins=247 het=352 unaligned=90797
163 M251    snp=1475    del=340 ins=248 het=540 unaligned=93882
164 M252    snp=1436    del=220 ins=215 het=264 unaligned=92322
165 M253    snp=1390    del=310 ins=385 het=552 unaligned=93391
166 M254    snp=1425    del=337 ins=250 het=540 unaligned=98770
167 M262    snp=1167    del=220 ins=192 het=1557    unaligned=72729
168 M266    snp=884 del=92  ins=122 het=286 unaligned=111270
169 M268    snp=1307    del=248 ins=191 het=179 unaligned=113602
170 M269    snp=1349    del=163 ins=236 het=502 unaligned=106949
171 M271    snp=1325    del=175 ins=180 het=164 unaligned=111226
172 M272    snp=1429    del=380 ins=432 het=351 unaligned=95957
173 M273    snp=1367    del=286 ins=305 het=217 unaligned=102389
174 M275    snp=1487    del=254 ins=376 het=702 unaligned=103994
175 M276    snp=1442    del=189 ins=276 het=297 unaligned=112154
176 M277    snp=1384    del=229 ins=206 het=173 unaligned=117257
177 M278    snp=1369    del=322 ins=239 het=677 unaligned=92782
178 M279    snp=1404    del=303 ins=209 het=451 unaligned=95434
179 M281    snp=1416    del=218 ins=289 het=231 unaligned=108482
180 M282    snp=1458    del=250 ins=299 het=516 unaligned=106271
181 M283    snp=1301    del=206 ins=282 het=258 unaligned=99717
182 M284    snp=1413    del=204 ins=235 het=153 unaligned=111235
183 M285    snp=1448    del=338 ins=471 het=452 unaligned=91711
184 M288    snp=1405    del=196 ins=221 het=313 unaligned=105677
185 M290    snp=1212    del=198 ins=164 het=257 unaligned=116121
186 M291    snp=1968    del=311 ins=295 het=651 unaligned=92568
187 M292    snp=1474    del=294 ins=394 het=356 unaligned=96377
188 M293    snp=1426    del=274 ins=227 het=327 unaligned=90885
189 M294    snp=1466    del=312 ins=376 het=432 unaligned=93511
190 M295    snp=1448    del=232 ins=288 het=282 unaligned=96298
191 M296    snp=2163    del=476 ins=546 het=728 unaligned=82100
192 M297    snp=1484    del=249 ins=299 het=324 unaligned=100420
193 M298    snp=2078    del=323 ins=457 het=606 unaligned=77867
194 M299    snp=1419    del=246 ins=317 het=314 unaligned=91170
195 M300    snp=1439    del=214 ins=274 het=320 unaligned=101838
196 M301    snp=1448    del=369 ins=286 het=609 unaligned=60818
197 M312    snp=1219    del=137 ins=125 het=272 unaligned=119871
198 M313    snp=959 del=179 ins=37  het=183 unaligned=139459
199 M314    snp=1311    del=142 ins=118 het=140 unaligned=102772
200 M315    snp=1106    del=161 ins=94  het=345 unaligned=101919
201 M316    snp=517 del=23  ins=15  het=135 unaligned=175825
202 M318    snp=1061    del=113 ins=65  het=277 unaligned=118584
203 M319    snp=1113    del=125 ins=118 het=482 unaligned=116494
204 M320    snp=1760    del=168 ins=194 het=493 unaligned=69201
205 M321    snp=904 del=107 ins=52  het=883 unaligned=101537
206 M322    snp=743 del=41  ins=48  het=925 unaligned=88172
207 M325    snp=1694    del=141 ins=121 het=701 unaligned=87386
208 M327    snp=866 del=87  ins=64  het=1138    unaligned=79505
209 M328    snp=1395    del=251 ins=199 het=231 unaligned=90591
210 M330    snp=2063    del=349 ins=319 het=529 unaligned=66670
211 M332    snp=931 del=70  ins=52  het=378 unaligned=130058
212 M333    snp=1478    del=479 ins=286 het=277 unaligned=84092
213 M334    snp=790 del=84  ins=34  het=125 unaligned=151486
214 M336    snp=1221    del=120 ins=166 het=128 unaligned=100734
215 M337    snp=1420    del=213 ins=286 het=197 unaligned=87815
216 M338    snp=1222    del=123 ins=85  het=485 unaligned=104272
217 M340    snp=1816    del=215 ins=268 het=495 unaligned=83224
218 M341    snp=720 del=85  ins=33  het=266 unaligned=122954
219 M343    snp=880 del=103 ins=53  het=416 unaligned=105425
220 M344    snp=862 del=72  ins=29  het=427 unaligned=98288
221 M346    snp=717 del=11  ins=16  het=134 unaligned=127435
222 M347    snp=1410    del=346 ins=331 het=403 unaligned=83688
223 M349    snp=1287    del=178 ins=183 het=199 unaligned=94069
224 M350    snp=535 del=51  ins=22  het=89  unaligned=137350
225 M356    snp=1044    del=96  ins=118 het=91  unaligned=111180
226 M357    snp=839 del=91  ins=155 het=102 unaligned=111477
227 M358    snp=777 del=84  ins=30  het=47  unaligned=132965
228 M359    snp=1962    del=248 ins=208 het=223 unaligned=94107
229 M361    snp=1246    del=249 ins=96  het=215 unaligned=93057
230 M362    snp=1432    del=297 ins=255 het=261 unaligned=86174
231 M364    snp=224 del=10  ins=15  het=110 unaligned=102702
232 M365    snp=1200    del=186 ins=98  het=328 unaligned=106300
233 M366    snp=1803    del=310 ins=210 het=609 unaligned=70770
234 M367    snp=654 del=126 ins=47  het=330 unaligned=110876
235 M368    snp=839 del=50  ins=71  het=180 unaligned=123685
236 M369    snp=879 del=216 ins=70  het=86  unaligned=94907
237 M370    snp=1406    del=202 ins=157 het=216 unaligned=101385
238 M371    snp=416 del=36  ins=33  het=373 unaligned=97100
239 M373    snp=1423    del=164 ins=280 het=209 unaligned=85214
240 M374    snp=1914    del=219 ins=257 het=174 unaligned=99659
241 M386    snp=1401    del=313 ins=242 het=212 unaligned=91981
242 M402    snp=884 del=162 ins=62  het=126 unaligned=86438
243 M403    snp=1884    del=285 ins=237 het=150 unaligned=79801
244 M404    snp=1079    del=162 ins=89  het=188 unaligned=86387
245 M406    snp=1523    del=202 ins=100 het=172 unaligned=81783
246 M407    snp=943 del=102 ins=97  het=297 unaligned=92702
247 M408    snp=1291    del=228 ins=159 het=83  unaligned=85956
248 M409    snp=1044    del=99  ins=112 het=122 unaligned=86084
249 M410    snp=1250    del=218 ins=90  het=254 unaligned=82976
250 M411    snp=891 del=206 ins=259 het=150 unaligned=87970
251 M412    snp=1380    del=267 ins=155 het=243 unaligned=83015
252 M413    snp=853 del=136 ins=59  het=191 unaligned=87816
253 M416    snp=1162    del=434 ins=285 het=342 unaligned=83346
254 M417    snp=1479    del=570 ins=453 het=276 unaligned=80895
255 M418    snp=1240    del=211 ins=169 het=181 unaligned=85885
256 M420    snp=824 del=139 ins=65  het=121 unaligned=88915
257 M422    snp=1453    del=349 ins=324 het=586 unaligned=76843
258 M423    snp=1104    del=159 ins=132 het=151 unaligned=85150
259 M426    snp=1404    del=211 ins=299 het=170 unaligned=94974
260 M430    snp=2048    del=547 ins=398 het=446 unaligned=65365
261 M433    snp=2150    del=721 ins=574 het=530 unaligned=65703
262 M462    snp=742 del=33  ins=43  het=185 unaligned=130258
263 M487    snp=803 del=142 ins=131 het=468 unaligned=78263
264 M490    snp=785 del=95  ins=42  het=654 unaligned=123300
265 M491    snp=521 del=86  ins=51  het=638 unaligned=128814
266 M492    snp=742 del=173 ins=105 het=513 unaligned=131446
267 M493    snp=779 del=84  ins=37  het=339 unaligned=95018
268 M494    snp=1587    del=205 ins=183 het=820 unaligned=94271
269 M495    snp=1073    del=104 ins=97  het=1153    unaligned=89866
270 M496    snp=744 del=160 ins=132 het=685 unaligned=118411
271 M497    snp=866 del=204 ins=106 het=727 unaligned=109548
272 M498    snp=1000    del=116 ins=39  het=620 unaligned=106143
273 M499    snp=619 del=117 ins=95  het=620 unaligned=80526
274 M500    snp=1115    del=145 ins=86  het=696 unaligned=116229
275 M501    snp=785 del=165 ins=145 het=383 unaligned=116032
276 M502    snp=912 del=98  ins=66  het=259 unaligned=117165
277 M503    snp=819 del=145 ins=89  het=333 unaligned=117423
278 M504    snp=1735    del=212 ins=144 het=239 unaligned=95726
279 M505    snp=1352    del=162 ins=142 het=133 unaligned=116052
280 M506    snp=1289    del=150 ins=120 het=537 unaligned=100102
281 M507    snp=1375    del=202 ins=243 het=609 unaligned=80246
282 M508    snp=996 del=271 ins=140 het=692 unaligned=100548
283 M509    snp=1379    del=235 ins=310 het=793 unaligned=72125
284 M510    snp=1928    del=338 ins=303 het=587 unaligned=84790
285 M511    snp=1718    del=357 ins=361 het=970 unaligned=83848
286 M512    snp=1776    del=313 ins=200 het=565 unaligned=97898
287 M513    snp=1899    del=443 ins=250 het=740 unaligned=85226
288 M515    snp=1298    del=209 ins=103 het=598 unaligned=108834
289 M516    snp=1378    del=176 ins=213 het=384 unaligned=100807
290 M518    snp=1002    del=119 ins=91  het=507 unaligned=107193
291 M519    snp=1471    del=181 ins=138 het=747 unaligned=91392
292 M520    snp=1828    del=280 ins=263 het=1057    unaligned=92027
293 M525    snp=1393    del=346 ins=372 het=284 unaligned=90448
294 M526    snp=1419    del=333 ins=293 het=295 unaligned=89866
295 M528    snp=1396    del=441 ins=428 het=252 unaligned=92203
296 M529    snp=1442    del=542 ins=454 het=581 unaligned=89755
297 M530    snp=1447    del=423 ins=485 het=467 unaligned=92268
298 M531    snp=1428    del=312 ins=412 het=194 unaligned=95146
299 M534    snp=1467    del=430 ins=481 het=553 unaligned=93297
300 M535    snp=627 del=220 ins=208 het=835 unaligned=77900
301 M536    snp=615 del=180 ins=190 het=340 unaligned=81101
302 M537    snp=1524    del=538 ins=454 het=1822    unaligned=88484
303 M538    snp=1333    del=411 ins=354 het=1163    unaligned=92429
304 M539    snp=1287    del=408 ins=353 het=367 unaligned=93206
305 M540    snp=1117    del=178 ins=344 het=120 unaligned=96620
306 M541    snp=1414    del=460 ins=425 het=527 unaligned=88988
307 M542    snp=1522    del=416 ins=207 het=275 unaligned=91707
308 M543    snp=1825    del=373 ins=249 het=794 unaligned=85411
309 M544    snp=1430    del=304 ins=392 het=1049    unaligned=98728
310 M547    snp=2009    del=656 ins=387 het=352 unaligned=80409
311 M548    snp=2080    del=605 ins=414 het=1204    unaligned=73706
312 M549    snp=1366    del=487 ins=406 het=264 unaligned=91765
313 M550    snp=1448    del=554 ins=451 het=646 unaligned=88454
314 M551    snp=1346    del=467 ins=403 het=338 unaligned=90434
315 M552    snp=1487    del=461 ins=433 het=728 unaligned=86708
316 M553    snp=1456    del=359 ins=368 het=169 unaligned=100108
317 M554    snp=1389    del=446 ins=464 het=345 unaligned=91451
318 M565    snp=593 del=81  ins=108 het=2438    unaligned=66262
319 M566    snp=1306    del=184 ins=175 het=581 unaligned=87280
320 M567    snp=1263    del=276 ins=205 het=642 unaligned=97084
321 M568    snp=677 del=161 ins=174 het=553 unaligned=135494
322 M569    snp=1396    del=266 ins=192 het=336 unaligned=101124
323 M570    snp=1981    del=307 ins=251 het=947 unaligned=83374
324 M572    snp=1355    del=218 ins=146 het=649 unaligned=104540
325 M576    snp=1357    del=164 ins=193 het=446 unaligned=96865
326 M577    snp=1301    del=228 ins=142 het=274 unaligned=107973
327 M579    snp=1364    del=190 ins=153 het=201 unaligned=113605
328 M580    snp=1133    del=129 ins=66  het=355 unaligned=129172
329 M581    snp=1957    del=244 ins=215 het=661 unaligned=96210
330 M582    snp=1353    del=205 ins=198 het=396 unaligned=104351
331 M583    snp=1342    del=204 ins=185 het=312 unaligned=112126
332 M584    snp=1400    del=229 ins=247 het=349 unaligned=98252
333 M585    snp=1357    del=223 ins=193 het=522 unaligned=97605
334 M586    snp=1382    del=207 ins=160 het=564 unaligned=100459
335 M587    snp=1941    del=258 ins=256 het=493 unaligned=95872
336 M588    snp=1314    del=191 ins=155 het=538 unaligned=78443
337 M590    snp=874 del=117 ins=101 het=294 unaligned=94391
338 M591    snp=1355    del=183 ins=169 het=404 unaligned=106011
339 M593    snp=1095    del=155 ins=132 het=620 unaligned=79839
340 M594    snp=1912    del=242 ins=194 het=322 unaligned=118435
341 M595    snp=1358    del=292 ins=145 het=303 unaligned=99081
342 M597    snp=1895    del=271 ins=205 het=238 unaligned=103243
343 M601    snp=1314    del=136 ins=105 het=374 unaligned=79737
344 M602    snp=1984    del=195 ins=253 het=794 unaligned=99224
345 M603    snp=1339    del=222 ins=164 het=196 unaligned=101581
346 M604    snp=1334    del=134 ins=151 het=447 unaligned=105988
347 M607    snp=1858    del=233 ins=187 het=213 unaligned=104419
Opening: core.tab
Opening: core.vcf
Processing contig: NC_000962
Generating core.full.aln
Creating TSV file: core.txt
Running: snp-sites -c -o core.aln core.full.aln
Have a suggestion? Tell me at http://github.com/tseemann/snippy/issues
Done.

I have add ".txt" at the end of the core.aln file so i can upload them. i'm puzzle why all the core alignment SNP are the same for all isolates? Maybe some stupid mistake by me? can someone help?

core.aln.txt

tseemann commented 4 years ago

why all the core alignment SNP are the same for all isolates

what do you mean?

SNP and distance matrix output most of them are 0 and 1

assuming you are using most of the reference genome in your core, this is a normal result if you have an outbreak or strongly related isolates.

foamy1881 commented 4 years ago

quick question when nullarbor2 generate the distance matrix which file does it use from snippy? core.aln or core.full.aln?

tseemann commented 4 years ago

@foamy1881 is uses core.aln

distances.tab : core.aln
  snp-dists -b $< > $@