Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 56,379 | 0 | T | G | 57.6% | 25.4 / 42.0 | 33 | Y244S (TAC→TCC) | lptD | LPS assembly protein LptD |
Reads supporting (aligned to +/- strand): ref base T (1/13); new base G (19/0); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.44e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG > NZ_CP009273/56297‑56468 | caTGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTc < 2:801786/90‑1 (MQ=255) atgatgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCa < 2:359458/90‑1 (MQ=255) tgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt < 1:857206/90‑1 (MQ=255) tgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGACGCaaaat > 2:274854/1‑88 (MQ=255) gTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAagta > 1:610388/1‑90 (MQ=255) tGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAagtagt > 1:68940/1‑90 (MQ=255) gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt < 2:185841/90‑1 (MQ=255) gacgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggtgg > 1:477099/1‑90 (MQ=255) acgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggtggt < 1:908003/90‑1 (MQ=255) aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTggtggtggt < 1:256564/90‑1 (MQ=255) aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTggtgggggt > 1:581072/1‑90 (MQ=255) aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAAACAAAGTAGTtggggggggg > 2:374972/1‑89 (MQ=255) aTGCATATAATGCGGCGTGATGGT‑GCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAAGAGTtgggggggggg > 1:542377/1‑88 (MQ=255) aaTGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAGAGTTGGGGGGGGTGTACTTgg > 2:202274/1‑90 (MQ=255) ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTT‑GTGGTGGTGTACTTGGCGTTCg < 1:508441/90‑1 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTTGGGGTGGTGGACTTTGCGTTCg > 1:29986/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTGGGGGTGGTGTACTTGGCGTTCg > 1:326195/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGTGGGGTACTTGGCGGTCg > 1:831301/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGGGGTGTACTTGGCGTTCg > 2:319314/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGTTGGGGGGGGGGTACTTGGGGTTCg > 2:242873/1‑90 (MQ=255) tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa < 1:202274/90‑1 (MQ=255) cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAAc < 2:2669/90‑1 (MQ=255) tCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCa < 2:899555/90‑1 (MQ=255) ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTGGTGGGGGGGTACTTGGCGTTCGGGGTCAAGAAACCag > 1:514457/1‑90 (MQ=255) ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGTTGGGGGGGGTGTACTTGGCGGTCGGGGTTAAGAAACCag > 2:177666/1‑90 (MQ=255) tatTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc < 2:386898/90‑1 (MQ=255) tatTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGGGGGGTACTTTGCGTTCGGGATCAAGAAACCAGAGc > 2:789283/1‑90 (MQ=255) gcgcGATGTTCCAGTAATATGGCAGGGAGAACTCAAAAAAGTTGGTGGGGGGGTACTTGGCGGTCGGGGTCAAGAAACCAGAGCGACGtt > 2:87723/1‑90 (MQ=255) gcgATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGGTCAAGAAACCAGAGCGACGTTTg > 1:910777/1‑90 (MQ=255) aTGTTCCAGTAATATGGCAGGGAGAACTCAAAGGAGTTGGTGGGGGGGTACTTGGGGTTCGGGGTCAAGAAAACAGAGCGGCGGTTGTca > 2:27700/1‑90 (MQ=255) atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCaa < 1:637166/90‑1 (MQ=255) ggCAGGTAGAACTCAAAGTAGTTGGTTGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGc < 2:406063/90‑1 (MQ=255) ggCAGGGAGAACTCAAAATAGTTGGTGGGGGGGGACTTGGCGTTCGGGGTCAAGAAAACAAAGCGACGTTTGTCACCCACCGGCCACTGc > 2:342676/1‑90 (MQ=255) tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg < 1:719443/90‑1 (MQ=255) tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg < 1:431803/90‑1 (MQ=255) | CATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG > NZ_CP009273/56297‑56468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |