Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,612,067 | 0 | A | G | 51.6% | 28.4 / 46.1 | 31 | E161G (GAG→GGG) | yjjJ | type II toxin‑antitoxin system HipA family toxin YjjJ |
Reads supporting (aligned to +/- strand): ref base A (3/12); new base G (16/0); total (19/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCA > NZ_CP009273/4611991‑4612152 | ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAg > 1:795670/1‑90 (MQ=255) tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTTGGCGGGGGGGGGAATTTTCAGc > 2:743957/1‑90 (MQ=255) ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGAt < 1:200729/90‑1 (MQ=255) cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 2:948711/90‑1 (MQ=255) cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 1:743957/90‑1 (MQ=255) cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 1:250045/90‑1 (MQ=255) gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGCTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGAtt < 2:376918/90‑1 (MQ=255) agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 2:360753/90‑1 (MQ=255) gaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGca < 2:834287/90‑1 (MQ=255) aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGcac > 1:853483/1‑90 (MQ=255) gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGGCGGGGGGGGGGATTATCAGCGATTGGTTTATGCacaa > 1:463646/1‑90 (MQ=255) gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCacaa > 1:403942/1‑90 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTTTTAGCGGTGGGTTACTGCACAACACCCCGCGGaa > 2:197879/1‑90 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGGTGGGTGGTCGGGGGGGGGGAATATCAGGGATGGGTTTCTGAACAACACCCCGCGGgaa > 2:942683/1‑90 (MQ=255) gACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCTGCGGAAAt > 1:413865/1‑90 (MQ=255) tttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTg < 2:164103/90‑1 (MQ=255) tttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTg < 2:172351/90‑1 (MQ=255) gTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTTCTCTGGATCaaa > 2:55408/1‑90 (MQ=255) gTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAAAACCCTGAGGGAATTCCTCTGGATCaaa > 1:943244/1‑90 (MQ=255) tGAATATACTGGCGGTTTGTTGGGCCGGGGGGGGGAATATCAGCGGTGGATTACTGGACAACACCCCGCGGAAAATCCCCTTGATCaaaa > 2:723173/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaaaaa > 1:844132/1‑90 (MQ=255) gAATATACTGGCGGCTGGTTGGCCGGGGGGGGGGATTATCAACGATGGATTACCGCCCAACACCCCCTGGAAAATCCTCCGGATCaaaaa > 2:216021/1‑90 (MQ=255) atatACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACt < 2:209929/90‑1 (MQ=255) taCTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGGAAATCCTCTGGGTCAAAAACTCAc > 1:456323/1‑90 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 2:341261/90‑1 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCTGCGGGAATTTCTCTGGATCAAAAAATCAcccc > 2:97115/1‑89 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGGAATTTTAGCGATGGGTTACTGCACAAAACCCCGCGGGAATTTCTCTGGGTTAAAAAATCAcccc > 2:699950/1‑89 (MQ=255) gCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGGTTACTGCACAACACCCCGCGGGAATTTCTCTGGATCAAAAACTCACCCAtt > 2:690188/1‑90 (MQ=255) ggttggttGGGCGGGGGGGGGGATTATCAGCGATGGATTACTGGACAACACCCCGGGGGAATTTCTCTGGATCAAAAAATCACCCATTAc > 2:577773/1‑90 (MQ=255) ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAc < 2:844132/90‑1 (MQ=255) gtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTg > 1:480488/1‑90 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 1:46543/90‑1 (MQ=255) gggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaa > 1:787701/1‑62 (MQ=255) gggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaa < 2:787701/62‑1 (MQ=255) | CTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCA > NZ_CP009273/4611991‑4612152 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |