[go: nahoru, domu]

Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Core.aln all the same nucleotide? #252

Closed
foamy1881 opened this issue Dec 5, 2019 · 3 comments
Closed

Core.aln all the same nucleotide? #252

foamy1881 opened this issue Dec 5, 2019 · 3 comments
Assignees

Comments

@foamy1881
Copy link
foamy1881 commented Dec 5, 2019

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 tseemann self-assigned this Dec 5, 2019
@tseemann
Copy link
Owner
tseemann commented Dec 5, 2019

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
Copy link
Author

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

@tseemann
Copy link
Owner

@foamy1881 is uses core.aln

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants