Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,612,067 | 0 | A | G | 47.5% | 29.8 / 15.7 | 41 | E161G (GAG→GGG) | yjjJ | type II toxin‑antitoxin system HipA family toxin YjjJ |
Reads supporting (aligned to +/- strand): ref base A (6/15); new base G (19/0); total (25/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTG > NZ_CP009273/4611979‑4612156 | gatgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag < 1:525098/90‑1 (MQ=255) gatgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCggggag > 1:531649/1‑90 (MQ=255) tgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGGCgggggggg > 2:1464692/1‑90 (MQ=255) gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggagggg < 1:213410/90‑1 (MQ=255) aCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGGGGGGGAGGGGGATTa > 2:199169/1‑90 (MQ=255) tctcTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCgggga < 2:977408/78‑1 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGc < 2:1085935/90‑1 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTGGGGGGGGGGGATTATTAGCg > 2:144142/1‑90 (MQ=255) tGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa < 1:257796/90‑1 (MQ=255) gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACt < 1:926426/90‑1 (MQ=255) aagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTg < 1:980171/90‑1 (MQ=255) aagaTGTGCACTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTg < 1:1360189/90‑1 (MQ=255) gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCacaa > 1:993229/1‑90 (MQ=255) tgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCacaac > 2:731427/1‑90 (MQ=255) ctctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAacac < 2:993229/90‑1 (MQ=255) tctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTTGTTGGTCGGGGGGGGGGAATTTAAACGATGGATTACTGCACCAAAcc > 1:896598/1‑90 (MQ=255) ctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACAccc < 2:137515/90‑1 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTTCACAACACCCTg > 2:827701/1‑90 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATTAGGGATGGATTACTGCACAACAccccg > 1:900034/1‑88 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATTAGCGATGGGTTTCTGCACAACAccccg > 1:365994/1‑88 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGAGTATAAGCGGTGGATTACTGCACAACACCCCGCGGaa > 2:1275165/1‑90 (MQ=255) gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGTGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAAt < 2:900034/90‑1 (MQ=255) cGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCCGCGGGAAATcc > 1:371009/1‑90 (MQ=255) aTTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTTCACAACACCCTGCGGAAATTCctct > 2:1071244/1‑90 (MQ=255) cGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCa > 2:1142595/1‑90 (MQ=255) ggTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaa < 2:494902/90‑1 (MQ=255) gTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTTTTAGCGATGGATTACTGCACAACACCCTGCGGAAATTTCTCTGGATTaaa > 2:1200676/1‑90 (MQ=255) gTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATAAGCGATGGATTACTGCACAAAACCCTGCGGAAATTCCTCTTGATCaaa > 1:442638/1‑90 (MQ=255) gTGAATATACTGGCGGGTGGGTGGGCGGGGGGGGGGATTATTAGCGATGGATTAATGCACAACACCCCGCGGGAATTTCTCTGGATCaaa > 1:972737/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTTTCAGCGATGGGTTACTGCACAAAACCCTGCGGAAATTTCTCTGGATCaaaaa > 2:220141/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaaaaa > 2:491190/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTTTCAGCGATGGATTACTGCACAAAACCCTGCGGAAATTCCTCTTGATCaaaaa > 2:1206390/1‑90 (MQ=255) gAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATTAGCGATGGATTAATTGACAACACCCCTCGGGAAATCCTCTGGGTCaaaaa > 1:473449/1‑90 (MQ=255) tataCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTc < 2:1226988/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 2:365994/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCATGCGGAAATTCCTCTGGATCAAAAACTCAcc < 2:928794/90‑1 (MQ=255) tGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCa < 1:1200676/90‑1 (MQ=255) tGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGGTGGATTTCTGCACAACACCCCGCGGAAATTCCTCTGGATCAAAAAATCCCCCa > 2:619069/1‑90 (MQ=255) gCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCAtt > 2:793177/1‑90 (MQ=255) gCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc > 2:1331151/1‑43 (MQ=255) gCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 1:1331151/43‑1 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTAc > 2:862816/1‑90 (MQ=255) ttggttggTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGa < 1:731427/90‑1 (MQ=255) ttgggtggTCGGGGGGGGGGATTATTAGCGAGGGATTACTGCACAACACCCTGCGGGAAGTCCTCTGGGGTAAAAACTCACCCATTAcca > 1:504032/1‑88 (MQ=255) tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaac < 2:528107/42‑1 (MQ=255) tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaac > 1:528107/1‑42 (MQ=255) ggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAg < 1:320432/90‑1 (MQ=255) gaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGt > 2:474301/1‑90 (MQ=255) aggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTg < 1:1071244/90‑1 (MQ=255) | GATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTG > NZ_CP009273/4611979‑4612156 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |