Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,970,032 | (G)6→7 | coding (1002/1047 nt) | wzzE → | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,970,026 | 1 | . | G | 97.6% | 142.2 / ‑3.4 | 42 | coding (996/1047 nt) | wzzE | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (20/21); total (21/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01 |
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑G‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA > NC_000913/3969888‑3970157 | ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc > 1:656260/1‑139 (MQ=255) ccATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑g < 2:196271/139‑1 (MQ=255) tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑gggg > 1:123480/1‑139 (MQ=255) tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTA‑TTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑ggggg > 2:841782/1‑139 (MQ=255) ttAAACCCCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑gggggg < 1:1316192/139‑2 (MQ=255) tGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTgg < 2:488430/139‑1 (MQ=255) gTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGgt < 1:1274932/139‑1 (MQ=255) ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTc > 1:1839495/1‑139 (MQ=255) ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTc > 1:525286/1‑139 (MQ=255) cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCAt > 2:403526/1‑139 (MQ=255) gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACccgc < 1:403526/139‑1 (MQ=255) cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc > 2:982578/1‑139 (MQ=255) cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc > 2:1792589/1‑139 (MQ=255) ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑g < 1:1840008/94‑1 (MQ=255) ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑g > 2:1840008/1‑94 (MQ=255) tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa < 2:123480/139‑1 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaacaac > 1:132874/1‑139 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaacaac > 1:967866/1‑139 (MQ=255) ttGCGTACGCCGGAAGAACCGGCAAAACGCAATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑GCGGGGGCTGATCGGAGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACtactgc > 1:621380/1‑139 (MQ=255) cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAt < 1:1609327/122‑1 (MQ=255) cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAt > 2:1609327/1‑122 (MQ=255) cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAgc > 2:652511/1‑139 (MQ=255) ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATT‑ACCCGCCGTTGCTCGAAATAACaa < 1:1563589/111‑1 (MQ=255) ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATT‑ACCCGCCGTTGCTCGAAATAACaa > 2:1563589/1‑111 (MQ=255) tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa < 1:521140/100‑1 (MQ=255) tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa > 2:521140/1‑100 (MQ=255) aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg < 2:132874/139‑1 (MQ=255) aCGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa > 2:1747900/1‑139 (MQ=255) cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt < 1:1747900/139‑1 (MQ=255) gtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATc > 2:1521452/1‑139 (MQ=255) tGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑gggggg < 1:549343/38‑2 (MQ=255) tGCCTTCCTGATGATTATGCGGGGCATTGTC‑g‑gggggg > 2:549343/1‑37 (MQ=255) cTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa < 2:1839495/139‑1 (MQ=255) ccTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg < 1:572536/133‑1 (MQ=255) ccTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg > 2:572536/1‑133 (MQ=255) ccTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAg > 2:832671/1‑139 (MQ=255) cTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt < 2:525286/139‑1 (MQ=255) tgatgaTTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa < 1:652511/139‑1 (MQ=255) gatTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga < 1:832671/139‑1 (MQ=255) gatTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga < 1:841782/139‑1 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca > 1:760557/1‑73 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca < 2:760557/73‑1 (MQ=255) aTGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGt < 1:1792589/139‑1 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa < 2:805275/127‑1 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa > 1:805275/4‑130 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa > 1:405511/4‑139 (MQ=255) gggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa < 2:761221/65‑1 (MQ=255) gggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa > 1:761221/1‑65 (MQ=255) ggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTgg > 2:1785386/1‑139 (MQ=255) gCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGt < 1:1521452/139‑1 (MQ=255) cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa < 2:967866/139‑1 (MQ=255) | GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑G‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA > NC_000913/3969888‑3970157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |