| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 700,996 | 0 | C | A | 49.2% | 59.3 / 38.6 | 65 | P533T (CCG→ACG) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
| Reads supporting (aligned to +/- strand): ref base C (9/24); new base A (32/0); total (41/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.24e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGA > NZ_CP009273/700907‑701081 | cGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAAc > 1:2169937/1‑90 (MQ=255) ccGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACg > 2:163095/1‑90 (MQ=255) gATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGAc < 2:2149915/90‑1 (MQ=255) aTTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACa > 2:401760/1‑90 (MQ=255) aTTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGGAAACGACa > 1:1553703/1‑90 (MQ=255) ttACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACCACAg > 2:926190/1‑90 (MQ=255) ttACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAACCCAcca > 2:2068903/1‑88 (MQ=255) ggTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGCCAGATaaa > 2:1483017/1‑90 (MQ=255) tGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATaaaaa > 2:151074/1‑89 (MQ=255) gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGGAAACGACAGATaaaaaa > 1:2762950/1‑88 (MQ=255) gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAtc < 1:330080/90‑1 (MQ=255) aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCTACAGATAAAAtcg < 2:1893414/90‑1 (MQ=255) aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAtcg < 1:2592921/90‑1 (MQ=255) tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAtcgt < 2:2169937/90‑1 (MQ=255) tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGGAAACCACAGATaaaaaaga > 1:1860237/1‑86 (MQ=255) tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAacgtcg > 2:745779/1‑90 (MQ=255) gTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTa < 1:1207096/90‑1 (MQ=255) tGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTAt < 2:1867270/90‑1 (MQ=255) ggCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGa > 2:611438/1‑71 (MQ=255) ggCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGa < 1:611438/71‑1 (MQ=255) cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCa > 2:2803741/1‑90 (MQ=255) cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGTATCACCa > 2:1538755/1‑90 (MQ=255) tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTTGAAACCG‑CATACAAAAATGTCGGATCACCac > 1:2249144/1‑89 (MQ=255) cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa < 1:401760/90‑1 (MQ=255) ggTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGGATCACCCGCCGCCgg > 2:474680/1‑90 (MQ=255) ttCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGGAAACGACAGATAAAATTGTCGTATTACCAGCCGCAgggg > 2:2801293/1‑89 (MQ=255) tCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGAc < 1:622075/90‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGGAAACGACAGATAAAAACGTCGTATTACCAACCGCCGGGGCa > 1:1004861/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGAc < 2:1161835/56‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAg < 1:2361799/79‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCAg > 2:2361799/1‑79 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGTATCACCAACCCCAGGGGCa > 2:1790609/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGAAAAAAAAGTCGTATCACCAGCCCCCggggga > 1:1227900/1‑87 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACAACAGGTAAAAAAGGCGGATTACCCGCCGCCggggga > 1:2127695/1‑87 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAAACGAc > 1:1161835/1‑56 (MQ=255) cTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGGAAACCACAAAAAAAAAAATCGGATCACCCGCCGCAGGGACaa > 1:1955342/1‑90 (MQ=255) tGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAAt > 2:1008273/1‑90 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGGATCACCAGCCGCAGGGGCAATCGTg > 2:2614456/1‑90 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGGATCACCAGCCGCCGGGGCAAACggg > 2:890162/1‑88 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGGGa > 2:726194/1‑90 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACCACAGATAAAAACGTGGTTTCACCAACCCCAGGGGGAATCgaga > 1:1221656/1‑90 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:2564332/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:819936/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:2019625/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:79516/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTCGGTGACGGTGTGGCGG‑GGAAAACGACAGATAAAAAAGTCGTATCACCAACCGCAGGGGGAAACGTGgaaa > 1:590388/1‑90 (MQ=255) acagcGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:2320275/86‑1 (MQ=255) ttCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAAACGACAGATAAAAAAGTCGTATCACCAGCCGCAGGGGCAAACGTGaaaaac > 1:1016937/1‑88 (MQ=255) tCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACCACAGAAAAAAAAGTCGTATCACCAGCCGCAGGGGCAAAAGTGAAAAACt > 1:462663/1‑90 (MQ=255) cGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAAACGACAGATAAAAACGTCGGATCACCCGCCGCAGGGGGAAACGTGGAAAACtt > 1:525822/1‑90 (MQ=255) gCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTc < 1:726194/90‑1 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa < 1:2562412/90‑1 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa < 2:2512950/90‑1 (MQ=255) aGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa < 2:2606353/90‑1 (MQ=255) aGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa < 2:2237764/90‑1 (MQ=255) gCGGTGGGTGACGGTGTGGCGG‑GGAAAACGGCAAAAAAAAAAGTCGTATTACCAGCCGCCGGGGCAAACGTGGAAAAATTCAAAacccgg > 1:2320275/1‑87 (MQ=255) gtgggtgACGGTGTGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGGGAAAAACTTCAACACCAACCAc > 2:2612965/1‑90 (MQ=255) tgACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGtt < 2:65601/90‑1 (MQ=255) aCGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAAACTTCAACACCAACCACGCGTTCt > 2:481698/1‑90 (MQ=255) gtgtGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAATCGTGAAAATCTTCAACACCAACCACGCGGTCTGcc > 2:2312593/1‑90 (MQ=255) gtgGCGG‑TGAAAACCACAAATAAAAAAGTCGTATCACCAGCCGCAGGGGGAAACGTGAAAAAATTTAAAAACCACCACCCGCTCTGCCTg > 1:1569865/1‑90 (MQ=255) ggCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa < 2:1291368/90‑1 (MQ=255) ggCGG‑TGAAAACGACAGATAAAAAAGTCGTATCACCAGCCCCAGGGGCAAACGGGAAAAAATTCAAAACCAACCCCCCGGTCTGCCTggg > 1:2216166/1‑89 (MQ=255) gAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa > 2:1891294/1‑90 (MQ=255) gAAACCGACAGATAAAATCGGCGTGTCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa > 2:644936/1‑90 (MQ=255) gAAAACGACAGAAAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAAAATTTAACACCCACCACCCGTTCTGCCTGGGAACCGa > 1:1688698/1‑90 (MQ=255) | CGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGA > NZ_CP009273/700907‑701081 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |