Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 56,379 | 0 | T | G | 59.5% | 8.2 / 35.1 | 37 | Y244S (TAC→TCC) | lptD | LPS assembly protein LptD |
Reads supporting (aligned to +/- strand): ref base T (4/11); new base G (22/0); total (26/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.60e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG > NZ_CP009273/56296‑56468 | acaTGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACt > 2:414418/1‑90 (MQ=255) atgatgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGACTCa > 2:46630/1‑90 (MQ=255) tgatgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaa < 2:723647/90‑1 (MQ=255) atgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAag < 2:588664/90‑1 (MQ=255) atgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAag > 1:689483/1‑90 (MQ=255) gTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAAATCAAagga > 1:324330/1‑88 (MQ=255) gTTGCCACGACGATGCACATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCCAagta > 1:78992/1‑90 (MQ=255) gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATA‑Gtt > 2:618153/1‑90 (MQ=255) gacgaTGCATATAATGCGGCGTGATGGGGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGTGGAATCAAAATATGTtggtg > 2:72946/1‑90 (MQ=255) aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTggtggtggt < 2:114609/90‑1 (MQ=255) tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATA‑GGTGGGGGTGGGGGACtttggg > 1:337602/1‑86 (MQ=255) tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTA‑GTTGGTGGGGGTGGACTTgggg > 2:310609/1‑88 (MQ=255) tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATA‑AGTGGTGGTGGGGGACTTgggg > 2:178387/1‑88 (MQ=255) tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAAA‑GTTGGTGGGGGGGTACtttggg > 2:365137/1‑86 (MQ=255) cggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGtt < 1:282453/90‑1 (MQ=255) cggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGAATCAAAATA‑GTTGGGGGGGGTGGACTTTGCggt > 2:648985/1‑88 (MQ=255) cggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGGA‑GTTGGGGGGGGGGTACTTgggggt > 2:438865/1‑86 (MQ=255) ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTc < 1:438865/90‑1 (MQ=255) ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAAA‑GGTGGTGGGGGGGTACTTTGCGTTc > 2:379874/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGAATCAAAATA‑GTTGGGGGTGGTGTTCTTGGGGTTCg > 2:262529/1‑90 (MQ=255) cgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAAATCAAAAAA‑GTTGGTGGTGGTGTACTTGGGGGTCgg > 1:341209/1‑90 (MQ=255) tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGa < 2:501254/90‑1 (MQ=255) ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGGGAACTCAAAAAA‑GTTGGGGGGGGTGTACTTGGGGCTCGGGGTCAaaaa > 2:488216/1‑90 (MQ=255) tatTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc < 2:41663/90‑1 (MQ=255) ttGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc < 1:414418/90‑1 (MQ=255) gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGGGTACTTGGCGGTTGGGAGCAAGAAAACAAAGCGACggt > 2:619587/1‑88 (MQ=255) gcgcGATGTTACAGTAATATGGCAGGTGGAAATCAAAAAA‑GTTGGTGGTGGGGGACTTGGGGGTCGGGGGCAAGAAACCAGAGCGACGtg > 2:640175/1‑89 (MQ=255) gATGTTCCAGTAATATGGCAGGGAGAACTCAAAATA‑GTTGGGGGTGGGGTACTTGGCGTTCGGGGTCAAGAAACCAGAGCGACGTTTGTc > 2:295383/1‑90 (MQ=255) gATGTTCCAGTAATATGGCAGGGAGAACTCAAAAGA‑GTTGGGGGGGGGGGACTTGGGGGTCGGGGTTAAGAAACCAGAGCGACGGTTTTc > 2:101372/1‑90 (MQ=255) aaTATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCa < 2:768699/90‑1 (MQ=255) tatGGCAGGGAGAAATCAAAGTA‑GTTGGGGGGGGGGTACTTGGGGTTCGGGATCAAAAAACCAAAGCGACGTTTGTCACCCACCGGCAAc > 2:135818/1‑90 (MQ=255) atGGCAGGGAGAACTCAAAGAA‑GGTGGGGGTGGGGTACTTGGCGGTCGGGGGCAAAAAAACAAAGCGGCGTTTGTCACCCACCGGCAAcc > 2:406156/1‑89 (MQ=255) ggCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCCCCCACCGGCAACTGc < 1:331717/90‑1 (MQ=255) aGGGAGAAGTCAAAGCA‑GTTGGGGGTGGGGTACTTGGGGTTCGGGGTCAAAAAAAAAGAGGGGCGTTTTTCACCCACCGGCAACGGCgaa > 2:116509/1‑90 (MQ=255) aGGGAGAACTCAAAGTA‑GTTGGTGGTGGGGTACTTGGCGTTCGGGGTCAAGAAAACAGAGCGACGTTTGTCACCCACCGGCAACTGGaaa > 2:678960/1‑90 (MQ=255) aGGGAGAACTCAAAGGA‑GTTGGGGGTGGGGGACTTGGCGTTCGGGGTCAAGAAACCAGAGCGGCGTTTGTCACCCACCGGCAACTGCaaa > 2:470890/1‑90 (MQ=255) gTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATa < 1:648985/90‑1 (MQ=255) tAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCCCCGGCAACTGCAAATAg < 2:341209/90‑1 (MQ=255) tAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg < 1:116994/90‑1 (MQ=255) | ACATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTA‑GTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG > NZ_CP009273/56296‑56468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |