Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,612,067 | 0 | A | G | 43.3% | 88.5 / 14.6 | 60 | E161G (GAG→GGG) | yjjJ | type II toxin‑antitoxin system HipA family toxin YjjJ |
Reads supporting (aligned to +/- strand): ref base A (9/25); new base G (26/0); total (35/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGT > NZ_CP009273/4611989‑4612155 | gTCTCTGGCAGGAAGAACATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGAAATATc > 1:1020887/1‑90 (MQ=255) tctcTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCa < 2:345587/90‑1 (MQ=255) ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAg > 1:178671/1‑90 (MQ=255) ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAg > 1:349631/1‑90 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGaa > 1:411445/1‑82 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGaa < 2:411445/82‑1 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGATGGTTGGTCGGGGAGGGGAATTATCAGc < 2:780459/90‑1 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCg > 1:803622/1‑90 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGAATATTAGCg > 2:441286/1‑90 (MQ=255) tGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa < 2:1047182/90‑1 (MQ=255) gaGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 2:1490802/89‑1 (MQ=255) cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 2:871724/90‑1 (MQ=255) gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGAtt < 1:289359/90‑1 (MQ=255) agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 1:677851/90‑1 (MQ=255) agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 2:178671/90‑1 (MQ=255) agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 2:909064/90‑1 (MQ=255) gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGATTTCt > 2:984387/1‑90 (MQ=255) aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGcac < 1:1061773/90‑1 (MQ=255) gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATTAGCGATGGATTAATGCacaaaa > 2:1210192/1‑88 (MQ=255) gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTTTCAGCGATGGGTTACTGCacaaca > 1:933736/1‑90 (MQ=255) gCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑g < 1:275948/90‑1 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGGTGGATTACTGCACAACACCCTGC‑ggga > 2:989637/1‑88 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAAAACCCTGC‑ggga > 1:1267134/1‑88 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCCGC‑ggga > 1:632870/1‑88 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCCGC‑ggga > 2:339685/1‑88 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAGCGGGGGATTACTGCACAAAACCCCGC‑GGaa > 1:644563/1‑90 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag > 1:174957/1‑47 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag < 2:174957/47‑1 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGgaa > 2:1082582/1‑90 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGaaa > 1:173716/1‑90 (MQ=255) gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAAt < 1:892408/90‑1 (MQ=255) cGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCCGCGGGAAATTc > 2:508698/1‑90 (MQ=255) cGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGAATTATTAGCGATGGGTTACTGCACAACACCCCGCGGGAAATTc > 2:117672/1‑90 (MQ=255) aTTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTTTAAGGGAGGGATTATTGCCAAACACCCTCC‑GGAAATTCctct > 2:113939/1‑90 (MQ=255) tttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTg < 2:803622/90‑1 (MQ=255) cGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCa < 1:1249960/90‑1 (MQ=255) gTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGGAATTTCTCTGGGTCaaa > 2:209493/1‑90 (MQ=255) gTGAATATACTGGCGGCTGGGTTGTCGGGGGGGGGAATTTTCAGCGATGGGTTAATGCAAAACACCCTCC‑GGGAAATACTCTGGGTTaaa > 2:1329633/1‑90 (MQ=255) tGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCCCAACACCCCGC‑GGGAAATTCTCTGGGTCaaaa > 2:317513/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGGAATTCCTCTGGGTCaaaaa > 1:796609/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAAATATCAGCGATGGATTACTGCACAACACCCCGC‑GGGAATTTCTCTGGGTCaaaaa > 2:148009/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 2:408194/53‑1 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc > 1:408194/1‑53 (MQ=255) gAATATACTGGCGGTTGGGTTGTCGGGGGGGGGGATTTTCAGCGATGGGTTTCTGCACAAAACCCTGC‑GGGAATTTCTCTGGATCaaaaa > 1:1087142/1‑90 (MQ=255) gAATATACTGGCGGTTGGGTGGTTGGGGGGGGGGATTTTTAGCGATGGGTTACTGCACAAAACCCCCC‑GGGAAATTCTCTGGGTTaaaaa > 2:1098614/1‑90 (MQ=255) tataCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTc < 2:883829/90‑1 (MQ=255) taCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCAc < 2:678858/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCAcc < 1:984387/90‑1 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGGGATCGATTACTGGACAAAACCCCGC‑GGTAATTCCTCTGGATCAAAAACTCAcccc > 2:562984/1‑89 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCa > 1:932975/1‑90 (MQ=255) tGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCa > 1:566996/1‑90 (MQ=255) tGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCa < 1:151128/90‑1 (MQ=255) tGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGGAAATCCTCTGGATCAAAAAATCACCCa > 2:658001/1‑90 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaa < 1:249789/48‑1 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaa > 2:249789/1‑48 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAAc > 2:537028/1‑79 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAAc < 1:537028/79‑1 (MQ=255) ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAc < 2:375822/90‑1 (MQ=255) ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAc < 2:1478295/90‑1 (MQ=255) tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGc > 2:57720/1‑90 (MQ=255) tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGc < 2:1273509/90‑1 (MQ=255) ggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCt > 1:1373335/1‑90 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 2:577206/90‑1 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 2:1020887/90‑1 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 2:463623/90‑1 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 1:1098841/90‑1 (MQ=255) ggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAg > 1:1127217/1‑90 (MQ=255) gaggggaATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGt > 1:49110/1‑90 (MQ=255) | GTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGT > NZ_CP009273/4611989‑4612155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |