Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,612,067 | 0 | A | G | 52.6% | 32.8 / 54.1 | 38 | E161G (GAG→GGG) | yjjJ | type II toxin‑antitoxin system HipA family toxin YjjJ |
Reads supporting (aligned to +/- strand): ref base A (3/15); new base G (20/0); total (23/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.27e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCA > NZ_CP009273/4611982‑4612152 | gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggagggg < 1:722470/90‑1 (MQ=255) ataCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGGCGGGGGGGGGGAt > 1:402503/1‑90 (MQ=255) gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATc < 1:830499/90‑1 (MQ=255) gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGAATTAtt > 2:390929/1‑89 (MQ=255) tctcTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCa > 2:625459/1‑90 (MQ=255) tctcTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGGATTATCa > 2:576218/1‑90 (MQ=255) ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAg > 1:821221/1‑90 (MQ=255) tctGGCAGGAAGGAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGc < 1:576703/90‑1 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGc < 1:621824/90‑1 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGc < 1:300818/90‑1 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag > 2:202166/1‑76 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag < 1:202166/76‑1 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAgag > 1:534149/1‑88 (MQ=255) gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGAtt < 1:390929/90‑1 (MQ=255) gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGa > 1:263119/1‑75 (MQ=255) gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa < 2:263119/75‑1 (MQ=255) tgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATTAGCGATGGATTACTGCacaac > 1:158556/1‑90 (MQ=255) gCTCTACGCCCTGACCGTATTTAACGGAGAATATACTGGCGGTTGGTTTGACAGGGGGGGGTAATAAGAGGGCAGAAGTACTGCacaaaa > 1:220099/1‑88 (MQ=255) tACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCt < 1:654671/90‑1 (MQ=255) tACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCt < 2:16635/90‑1 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTTGGGGAGGGGGATTTTCAGCGATGGATTACTTCCCAAAAccccg > 1:590383/1‑88 (MQ=255) acagtGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑g < 1:86920/86‑1 (MQ=255) cccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑gg < 1:466348/90‑1 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑ggga > 2:463610/1‑88 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGAATTATCAGGGATGGGTTACTGCACAACACCCTGC‑GGaa > 2:82584/1‑90 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCTGC‑GGgaa > 1:189261/1‑90 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAGCGATGGATTACTTCACAACACCCCGC‑GGgaa > 1:103006/1‑90 (MQ=255) gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAAt < 1:224021/90‑1 (MQ=255) tATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCctc < 1:788013/90‑1 (MQ=255) tATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCctc < 2:402503/90‑1 (MQ=255) tATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATTAGCGATGGGTTACTGCACAACACCCCGCGGGAAATTctc > 1:260791/1‑88 (MQ=255) tttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTg < 1:155477/90‑1 (MQ=255) tAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGAAATTCctcctgt > 2:86920/1‑86 (MQ=255) aCGGTGAATATACTGGCGGTTGGGTGGTGGGGGGGGGGGAATTTTAGCGAAGGATTACTGGACAACACCCTCC‑GGGAATTCCTTTGGGTc > 1:871659/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGGTCaaaaa > 2:585420/1‑90 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCa > 2:219356/1‑90 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAAATCACCCa > 2:7897/1‑90 (MQ=255) tGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGAAATTTCTCTGGATCAAAAACTCACCCa > 2:592496/1‑90 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACt > 2:48825/1‑80 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACt < 1:48825/80‑1 (MQ=255) tggttggtCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGGAATTTCTCTGGGTCAAAAAATCACCCATTACGaa > 1:39123/1‑90 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 1:576218/90‑1 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 1:576065/90‑1 (MQ=255) | GATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCA > NZ_CP009273/4611982‑4612152 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |